Error_in_finding_100_hPa_up

From: JinQ Liu <liujq_0912_at_nyahnyahspammersnyahnyah>
Date: Tue Nov 22 2011 - 21:51:25 MST

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