Memory Issue with WRF Functions

Dear WRF Users,

I have stumbled upon an apparent memory issue and am having difficulty
resolving it. What I want to do is produce a .PNG for each variable and
output time, but memory allocation goes up after each loop even though I
delete all variables. In fact, I delete every variable that the list_vars()
function returns and NCL still has memory allocated! It seems at least one
of the WRF functions (*wrf_contour* or *wrf_map_overlays*) is allocating
memory in a non-obvious way. Can someone provide some help or insight on
how to solve this issue?

Thanks,
Gus

Here is a snippet of code to give you a better idea of what I am doing:

;================================================================
; Load scripts
;================================================================

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/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

;================================================================
; Begin the code.
;================================================================

begin

print("Time Started: "+systemfunc("date"))
print("")

;================================================================
; User-defined.
;================================================================

mov_dir = "/Users/galaka/Documents/Research/mov/"

plot_type = "png"
plot_type@wkWidth = 1500
plot_type@wkHeight = 1500

system("rm "+mov_dir+"*"+plot_type+"*")

today = systemfunc("date +%Y%m%d")
today2 = systemfunc("date +%Y-%m-%d")

;================================================================
; Variable Plotting.
;================================================================

res = True
res@MainTitle = "Metgrid Files"
res@Footer = False
res@cnFillOn = True

pltres = True
mpres = True

opts = res

print("")
print("Domain #1")
f =
addfile("/maloney-scratch/datasets/WRF_output/"+today+"/wrfout_d01_"+today2+"_00:00:
00.nc","r")

;D01 Variables
lon = f->XLONG(0,0,:)
lat = f->XLAT(0,0,:)
times = chartostring(f->Times)
nt = dimsizes(times)

;Base Reflectivity
print("")
print(" Base Reflectivity")
do a = 1,nt-1
dbz = wrf_user_getvar(f,"dbz",a)

wks = gsn_open_wks(plot_type,mov_dir+"WRFout."+systemfunc("date
+%Y%m%d")+"00.d01.basedbz"+sprinti("%0.3i",a))
;wks = gsn_open_wks(plot_type,"testing")

res@TimeLabel = times(a)

contour = wrf_contour(f,wks,dbz(1,:,:),res)
plot = wrf_map_overlays(f,wks,(/contour/),pltres,mpres)

delete(wks)
delete(contour)
delete(plot)
delete(dbz)
end do
;delete(lat)
;delete(lon)
;delete(times)
;delete(opts)
;delete(pltres)
;delete(today2)
;delete(mov_dir)
;delete(nt)
;delete(mpres)
;delete(plot_type)
;delete(res)
;delete(a)
;delete(today)
;list_vars()
print("Sleeping")
sleep(30)
exit

-- 
Ghassan J. Alaka, Jr., Ph.D. Candidate
Department of Atmospheric Science
Colorado State University
1371 Campus Delivery
Fort Collins, CO 80523-1371
E: gusalaka@atmos.colostate.edu
P: (970) 491-1487

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 18 11:46:13 2012