Hi NCL help,
I am having WRF output files, saved in 1 hour interval. Total number of
files is 82.
I want to calculate vortcity at 850 hPa. following script is showing errors
like:
fatal:ContourPlotPreDraw: Workspace reallocation would exceed maximum size
33554432
fatal:ContourPlotDraw: draw error
warning:WorkstationDeactivate: workstation not active or not opened
In my .hluresfile: wsMaximumSize: 33554432
;========Script is here ===========
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
datadir = "/WRF/output/2010072800/"
FILES = systemfunc ("ls -1 " + datadir + "wrfout_d02* " )
numFILES = dimsizes(FILES)
res = True
res@InitTime = False
res@Footer = False
res@cnFillOn = True
res@FieldTitle = ""
res@cnLinesOn = False
res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = (/
-90.,-70.,-50.,-30.,0.,30.,50.,70.,90.,110.,130.,150./)
pltres = True
mpres = True
mpres@mpNationalLineThicknessF = 1.0
mpres@mpNationalLineColor = "black"
mpres@mpGeophysicalLineColor = "black"
mpres@mpGeophysicalLineThicknessF = 1.0
plevs = (/ 850./)
nlevels = dimsizes(plevs)
do ifil = 0,numFILES-1
a = addfile(FILES(ifil)+".nc","r")
times = wrf_user_list_times(a)
;;;===========================
day = systemfunc("echo "+times +"|cut -c9-10")
yr = systemfunc("echo "+times +"|cut -c1-4")
mon = systemfunc("date '+%b'|tr [a-z] [A-Z]")
hr = systemfunc("echo "+times +"|cut -c12-13")
;;;===========================
p = wrf_user_getvar(a, "pressure",0)
vo = wrf_user_getvar(a,"avo",0)
F = a->F
f1 = F(0,:,:)
pressure = plevs(0)
wks =
gsn_open_wks("ps","vorticity_"+pressure+"hPa"+"_"+day+mon+yr+"_"+hr+"Z")
gsn_define_colormap(wks,"gui_default")
vo_plane = wrf_user_intrp3d(vo,p,"h",850.,0.,False)
vortf = (vo_plane - f1)*10 ;;;
Vorticity
vort = wrf_contour(a,wks,vortf,res)
if(any(pressure.eq.plevs)) then
plot = wrf_map_overlays(a,wks,vort,pltres,mpres)
end if
end do
end
=============================
Thanks & regards,
Sahidul
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 28 06:39:43 2010
This archive was generated by hypermail 2.1.8 : Wed Jul 28 2010 - 10:42:06 MDT