Hi,
I have a wrfout file which contain 48 hour , I want to draw slp of every hour.
When run the script, I can only get pictures in 1-36 hour, the 37-48 hour picture can not be plot.
The error is as follows:
Troubles finding level 100 above ground.
Problems first occur at ( 1, 1)
Surface pressure = NaN hPa.
STOP Error_in_finding_100_hPa_up
I have transfer the data on the ftp. Can anyone give me some suggestions?
Thanks.
JinQ
-------------------------------------------------------------------------------------------------
ftp> put wrfout_d01_2010-09-18_00:00:00
-------------------------------------------------------------------------------------------------
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/WRF_contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
a = addfile("wrfout_d01_2010-09-18_00:00:00.nc","r")
wks=gsn_open_wks("ps","slp") ; Open PS file.
gsn_define_colormap(wks,"BlGrYeOrReVi200") ; Change color map.
times = wrf_user_list_times(a) ; get times in the file
ntimes = dimsizes(times)
res= True
pltres = True
mpres = True
do it=0,ntimes-1
slp = wrf_user_getvar(a,"slp",it)
opts=res
opts@cnFillOn = True
opts@TimeLabel = times(it)
contour_a = wrf_contour(a,wks,slp,opts)
plot = wrf_map_overlays(a,wks,(/contour_a/),pltres,mpres)
end do
end
-- ` `
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 22 21:51:36 2011
This archive was generated by hypermail 2.1.8 : Wed Nov 23 2011 - 16:46:19 MST