Re: problem with climDayTLL

Hi Dennis,

Thank you for your reply.

Here is the input metadata as well as a value for one time slice (day 10, 10N, 10E):

ncl 9> printVarSummary(TA)

Variable: TA
Type: float
Total Size: 1132037280 bytes
            283009320 values
Number of Dimensions: 3
Dimensions and sizes: [time | 8030] x [lat | 198] x [lon | 178]
Coordinates:
            time: [701256..894096]
            lat: [-12.20433..37.77321]
            lon: [-20.74124..27.01102]
Number Of Attributes: 4
  long_name : Anemom Temp
  units : K
  missing_value : -1e+30
  _FillValue : -1e+30

; at day 10 at 10N 10E the value is....
ncl 14> print(TA(10,{10},{10}))

Variable: TA (subsection)
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
Number Of Attributes: 7
  lon : 10.01446
  lat : 10.04281
  time : 701496
  _FillValue : -1e+30
  missing_value : -1e+30
  units : K
  long_name : Anemom Temp
(0) 296.6686

Now,

Here is the command i execute and the output metadata (after using climDayTLL) as well as a value for the same point, same day climatology value:

;dataset has attributes of time lat lon and begins on year 1980 at day Jan 1 (001)
ncl 8> TA_daily=clmDayTLL(TA,19800001)
ncl 9> printVarSummary(TA_daily)

Variable: TA_daily
Type: float
Total Size: 51597216 bytes
            12899304 values
Number of Dimensions: 3
Dimensions and sizes: [year_day | 366] x [lat | 198] x [lon | 178]
Coordinates:
            year_day: [1..366]
            lat: [-12.20433..37.77321]
            lon: [-20.74124..27.01102]
Number Of Attributes: 5
  _FillValue : -1e+30
  long_name : Daily Climatology: Anemom Temp
  units : K
  information : Raw daily averages across all years
  smoothing : None

;looks to have executed with no problem, however plotting or looking at a variable value reveals....

ncl 11> print(TA_daily(10,{10},{10}))

Variable: TA_daily (subsection)
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
Number Of Attributes: 8
  lon : 10.01446
  lat : 10.04281
  year_day : 11
  smoothing : None
  information : Raw daily averages across all years
  units : K
  long_name : Daily Climatology: Anemom Temp
  _FillValue : -1e+30
(0) -1e+30

And this is for every timestep and for every point. Clearly, I am doing something wrong. If you have any ideas, it would be much appreciated.

-Marc

________________________________________
From: Dennis Shea [shea@ucar.edu]
Sent: Monday, March 04, 2013 10:40 AM
To: Marc Pace Marcella
Cc: ncl-talk@ucar.edu
Subject: Re: problem with climDayTLL

This does not tell us much other than you think 'clmDayTLL'
is not working.

Post a printVarSummary of a sample variable.

     f = addfile("RCM.nc", "r")
     x = f->Whatever
     printVarSummary(x)

--
If you are new to NCL, please read the mini-language manual at:
    http://www.ncl.ucar.edu/Document/Manuals/
On 3/4/13 8:25 AM, Marc Pace Marcella wrote:
> Hi NCLers,
>
> I'm new here to NCL and I am trying to use the climDayTLL function to take a netcdf of daily RCM variables for ~30 years and construct the daily climatology (eg the Jan 30th climate).  However, when I use the climDayTLL, the new variable I declare is filled with all -1E30's, the same value that is used for Missing Values.  My dataset begins on Jan 01, 1980 so I am using 19800001 as the YYYYDDD input and my input dataset is stamped as time and time is the left most axis in the form of Time, Lat, Lon.  So, Im not sure what my issue could be?? Any help anyone could give, would be greatly appreciated.  Thanks!
>
> Marc
> _______________________________________________
> 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