vorticity plot problem: draw error

From: Sahidul <sahidul_at_nyahnyahspammersnyahnyah>
Date: Wed Jul 28 2010 - 06:39:38 MDT

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