[Thread Prev][Thread Next][Index]

[ferret_users] How to find the number of seconds in a timestep



I'm sure I'm just missing something here .... 
 
I'm reading variables from a netcdf file that has a regular time axis in days (as I know by listing a few timesteps). 
I can also see by this listing that the timesteps are 2 hours, but neither the axis unit nor the timestep length is a number I can use (without listing and looking each time).
 
How can I reliably find the number of seconds in a timestep?
Alternately, how can I find the units of the time axis?
 
yes? sho axis time
name       axis              # pts   start                end
TIME      TIME              7149 r   17-JUL-2012 04:00    04-MAR-2014 20:00
T0 = 01-JAN-1970 00:00:00
 
let TT = t[gt=u]   ! listing a few TT values shows by eye that the time axis is days, and the timestep 2 hours, but these are not numbers I can use except manually
 
list TBOX[gt=u]  ! shows that TBOX=1/12 of the time axis unit, but I don't know that unit length except by eye-balling a listing
 
let TAXIS_UNIT = 1/TT[l=@ddf]   ! length of a time axis unit in seconds (because Ferret takes time derivatives in seconds)
 
let TSEC = TBOX[gt=u] * TAXIS_UNIT   ! this calculation finally tells me the length of a timestep, in seconds
 
Clearly Ferret knows this information but won't tell me without a calculation. Am I missing something obvious?
 
It seems that TAXIS_UNIT should be a keyword like TBOX, TBOXLO, TBOXHI. Those are all similar to index units in that they are only meaningful once you know what the geophysical time axis unit is. TAXIS_UNIT (or equivalently TSEC) would make these available for calculations. 
 
And note that my method above fails at the final timestep (because of @ddf). In the present case that doesn't matter because I know that the axis is regular, but if it wasn't I would have a hard time getting this information.
 
Thanks ... Billy K

[Thread Prev][Thread Next][Index]