Re: fatal:Undefined identifier: (rainfall_tot) is undefined, can't continue

The 1st place ranfall_tot is referenced is in line 31

  prcDay(nt/24,:,:) = dim_sum_n(rainfall_tot(ntStrt:ntLast,:,:), 0)

===
I think it should be

  prcDay(nt/24,:,:) = dim_sum_n(rainfall_hour(ntStrt:ntLast,:,:), 0)

===

Prior to that .... Why are you doing this???

> rainfall_hour = new ( (/ntim,nlat,mlon/), "float")
> do it = 0,ntimes-1
> rain_con = wrf_user_getvar(a,"RAINC",-1)
> rain_exp = wrf_user_getvar(a,"RAINNC",-1)
> rainfall_hour = rain_con + rain_exp
> end do

You should be able to dor

  rain_con = wrf_user_getvar(a,"RAINC",-1)
  rain_exp = wrf_user_getvar(a,"RAINNC",-1)
  rainfall_hour = rain_con + rain_exp

On 12/6/10 7:18 AM, wei huang wrote:
> Dear NCL users,
> I have 90 hours of wrf data and I want to convert it into number of
> days. During my conversion process the following error is appeared given
> below:
>
> fatal:Undefined identifier: (rainfall_tot) is undefined, can't continue
>
> fatal:Execute: Error occurred at or near line 31 in file rain10d1.ncl
>
> Any advice/suggestion will be appreciated
>
> My script is attached herewith
>
> Thanks
>
> wei
>
> *************************************
>
> 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/WRFUserARW.ncl"
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
>
> begin
>
> a = addfile("/wrfout.nc <http://wrfout.nc>","r")
>
> times= wrf_user_list_times(a); get times in the file
>
> ntimes = dimsizes(times)
>
> print(times)
>
> xlat = a->XLAT(0,:,:)
>
> xlon = a->XLONG(0,:,:)
>
> dimll= dimsizes( xlat )
>
> nlat = dimll(0)
>
> mlon = dimll(1)
>
> ntim = ntimes
>
> rainfall_hour = new ( (/ntim,nlat,mlon/), "float")
>
> do it = 0,ntimes-1
>
> rain_con = wrf_user_getvar(a,"RAINC",-1)
>
> rain_exp = wrf_user_getvar(a,"RAINNC",-1)
>
> rainfall_hour = rain_con + rain_exp
>
> end do
>
> ntJump = 24 ; 24 'samples' per day
>
> prcDay = rainfall_hour(::ntJump,:,:)
>
> ntStrt = 0
>
> ntLast = ntJump-1
>
> do nt=0,ntim-1,ntJump
>
> prcDay(nt/24,:,:) = dim_sum_n(rainfall_tot(ntStrt:ntLast,:,:), 0)
>
> ntStrt = ntStrt + ntJump
>
> ntLast = ntLast + ntJump
>
> end do
>
> prcDay@long_name = "Daily total precipitation"
>
> prcDay@units = rain_con@units
>
> copy_VarCoords(xlat,prcDay(0,:,:))
>
> printVarSummary(prcDay)
>
> 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 Mon Dec 6 08:56:09 2010