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

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Mon Dec 06 2010 - 08:54:00 MST

The error is telling you exactly what the problem is. You are trying to reference the variable "rainfall_tot" which doesn't exist.

The error is on this line:

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

Perhaps you meant to use "rainfall_hour" here instead?

--Mary

On Dec 6, 2010, at 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","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:54:10 2010

This archive was generated by hypermail 2.1.8 : Mon Dec 06 2010 - 09:18:09 MST