Re: Error while trying run clmDayTLL

Thank you very much for the reply. It seems functions such as smthClmDayTLL
and band_pass_latlon_time don't accept missing values. May I know what
might be the best way to deal with it?

Best,
Abheera

On Sat, Aug 25, 2012 at 8:36 AM, Dennis Shea <shea@ucar.edu> wrote:

> The issue is that you did not name the dimensions.
> From the clmDayTLL documentation:
>
> x
>
> A three-dimensional array (time, lat, lon).
> The dimensions must be named. <====
> The Description section shows how to do this.
>
> Normally, variables from netCDF of HDF are named.
> Since, the data are from a binary file, the user
> must provide this information.
>
> A few comments:
>
> -------
> Rather than multiple do loops to create yyyyddd:
>
> remove:
> yyyyddd = new(1096,integer)
> and the 3 do loops. Use
>
> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/contributed.**ncl"
>
> date = yyyymmdd_time(2007,2009, "integer") ; yyyymmdd
> yyyyddd = yyyymmdd_to_yyyyddd( date )
> print("date="+date+" yyyyddd="+yyyyddd)
>
> http://www.ncl.ucar.edu/**Document/Functions/**Contributed/yyyymmdd_time.*
> *shtml<http://www.ncl.ucar.edu/Document/Functions/Contributed/yyyymmdd_time.shtml>
>
> http://www.ncl.ucar.edu/**Document/Functions/**Contributed/yyyymmdd_to_**
> yyyyddd.shtml<http://www.ncl.ucar.edu/Document/Functions/Contributed/yyyymmdd_to_yyyyddd.shtml>
> -------
>
> ntim = dimsizes(yyyyddd)
>
> and replace
>
>
> xin = new((/1096,27,121,240/),"**float")
>
> with
> xin = new((/ntim,27,121,240/),"**float")
>
> actually, I'd suggest not hard-wiring any constants
>
> klev = 27
> nlat = 121
> mlon = 240
> xin = new((/ntim,klev,nlat,mlon/),"**float")
>
> This will facilitate future reuse with different array sizes.
>
> ***********
> Further, name the dimensions. *Any* unique identifier can be used
> for the dimension names. Common, names are
>
> xin!0 = "time" ; "T", "temporal_dimension",......
> xin!1 = "lev"
> xin!2 = "lat"
> xin!3 = "lon"
>
> ***********
> The
>
> set_default_fillvalue("float", 9.999e20)
>
> should be placed before you read from the binary files
>
> =========
> Also, replace the levelel loop with
>
> nDAY = 366
> xclim = new((/nDAY,klev,nlat,mlon/),"**float")
>
> do i=0,klev-1
> print("i="+i) ; debug
> xout = clmDayTLL(xin(:,i,:,:),**yyyyddd)
> printVarSummary(xout) ; debug
>
> xclim(:,i,:,:) = smthClmDayTLL(xout, 3)
> end do
>
> ++++++++++++++++++++++++++++++**+++++++
> You can always look at the script being used
>
> %> less $NCARG_ROOT/lib/ncarg/**nclscripts/csm/contributed.ncl
>
> search for 'function clmDayTLL' or 'smthClmDayTLL'
>
> Any number of harmonics (up to 365/2)
>
>
>
>
>
>
> On 8/24/12 6:26 PM, abheera hazra wrote:
>
>> Hi,
>>
>> I am trying to use the function clmDayTLL, I am using binary data for
>> input, so I create the yyyyddd variable, since the function says it is
>> an integer array, that is how i declare it initially, here is my code;
>>
>> ****************
>>
>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/contributed.**ncl"
>>
>> pathi = new(10,string)
>> patho = new(10,string)
>>
>> yyyyddd = new(1096,integer)
>> xin = new((/1096,27,121,240/),"**float")
>> xcl = new((/1096,121,240/),"float")
>> xclim = new((/366,27,121,240/),"float"**)
>> ;xout = new((/366,121,240/),"float")
>> ;xout!0="year_day"
>> ;yyyyddd!0="year_day";xcl!0="**year_day"
>>
>> pathi(1)= "/project/dheat.mean.2007"
>> pathi(2)= "/project/dheat.mean.2008"
>> pathi(3)= "/project/monisi/dheat.mean.**2009"
>>
>> patho(1)= "/project/dheat.clim.leap"
>> patho(2)= "/project/dheat.clim"
>>
>> patho(3)= "/project/dheat.anom.2007/"
>> patho(4)= "/project/dheat.anom.2008/"
>> patho(5)= "/project/dheat.anom.2009/"
>>
>>
>> k=-1
>> do j=1,365
>> k=k+1
>> yyyyddd(k)= 2007*1000+j
>> ; print (k)
>> end do
>>
>> k=364
>> do j=1,366
>> k=k+1
>> yyyyddd(k)= 2008*1000+j
>> ; print (k)
>> end do
>>
>> k=730
>> do j=1,365
>> k=k+1
>> yyyyddd(k)= 2009*1000+j
>> ; print (k)
>> end do
>>
>> print ("yyyyddd")
>>
>> nrec = 0
>> dims = (/365,27,121,240/)
>> x1 = fbindirread(pathi(1), nrec, dims, "float")
>>
>> dims = (/366,27,121,240/)
>> x2 = fbindirread(pathi(2), nrec, dims, "float")
>>
>> dims = (/365,27,121,240/)
>> x3 = fbindirread(pathi(3), nrec, dims, "float")
>> set_default_fillvalue("float", 9.999e20)
>>
>> print ("dheat mean files read")
>>
>> xin(0:364,:,:,:) = x1(0:364,:,:,:)
>> xin(365:730,:,:,:) = x2(0:365,:,:,:)
>> xin(731:1095,:,:,:) = x3(0:364,:,:,:)
>>
>> print ("dbh rewritten")
>> set_default_fillvalue("**integer", 10000000000)
>>
>> ;
>>
>> do i=0,26
>> xcl(0:1095,0:120,0:239) = xin(0:1095,i,0:120,0:239)
>> print ("xcl <- xin")
>> xout = clmDayTLL(xcl,yyyyddd)
>> print("clm")
>> xoutsm = smthClmDayTLL(xout, 3)
>> print("smoothed")
>> xclim(:,i,:,:) = xoutsm(:,:,:)
>> end do
>>
>> print ("smoothed climatology")
>>
>> ****************
>>
>>
>> But I keep getting the following error
>>
>>
>> *****************
>>
>> Copyright (C) 1995-2011 - All Rights Reserved
>> University Corporation for Atmospheric Research
>> NCAR Command Language Version 6.0.0
>> The use of this software is governed by a License Agreement.
>> See http://www.ncl.ucar.edu/ for more details.
>> (0) yyyyddd
>> (0) dheat mean files read
>> (0) dbh rewritten
>> (0) xcl <- xin
>> warning:Attempt to reference attribute (year_day) which is undefined
>> (0) clm
>> fatal:(missing) is not a dimension name in variable (clmDay), could not
>> determine dimension number
>> fatal:Execute: Error occurred at or near line 3969 in file
>> $NCARG_ROOT/lib/ncarg/**nclscripts/csm/contributed.ncl
>>
>> fatal:Execute: Error occurred at or near line 84 in file clim.ncl
>>
>> *****************
>>
>> Would like to know where I'm going wrong,
>>
>> Also, for the function smthClmDayTLL, how many harmonics can one
>> include? The NCL site mentions generally 2 or 3, but is it possible to
>> have more harmonics included?
>>
>> Thanks,
>> Abheera
>>
>>
>> ______________________________**_________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<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 Aug 28 00:49:08 2012