Hello NCL user,
Is it possible to include simulation time at the top of graphics file
generated by following script
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
a = addfile("/root/Desktop/wrfout_d03_2012-04-26_00:00:00.nc","r")
; We generate plots, but what kind do we prefer?
type = "pdf"
wks = gsn_open_wks(type,"pressure_wrfout")
res = True ; Set up some basic plot resources
res@MainTitle = "PRESSURE FILES"
res@Footer = False
pltres = True
mpres = True
opts = res
opts@cnFillOn = True
times = wrf_user_getvar(a,"times",-1) ; get times in the file
ntimes = dimsizes(times) ; number of times in the file
infoTimeStamp(wks,0.015,"infoTimeStamp test")
print(ntimes)
do it = 0,ntimes-1,4 ; TIME LOOP
res@TimeLabel=times(it)
print("Working on time: " +it+" "+ times(it) )
slp = wrf_user_getvar(a,"PSFC",it) ; Get sea level pressure from file
; Add time stamp.
infoTimeStamp(wks,0.015,"infoTimeStamp test")
contour = wrf_contour(a,wks,slp,opts)
plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
end do
end
Any help will be appreciated.
Thank you.
Sangeeta Maharjan.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Apr 6 00:45:45 2013
This archive was generated by hypermail 2.1.8 : Sun Apr 07 2013 - 21:13:30 MDT