Re: Error_in_finding_100_hPa_up

This is being forwarded to wrfhelp@ucar.edu

On 11/22/11 9:51 PM, JinQ Liu wrote:
> 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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 22 22:38:35 2011