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-talkReceived on Mon Mar 4 09:01:04 2013
This archive was generated by hypermail 2.1.8 : Thu Mar 07 2013 - 08:55:58 MST