Re: vorticity plot problem: draw error

Sahidul,

There is an easy fix for this. Please go to:

http://www.ncl.ucar.edu/FAQ/#err_msgs

and search for "workspace reallocation"

--Mary

On Jul 28, 2010, at 6:39 AM, Sahidul wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 28 10:27:40 2010