[Thread Prev][Thread Next][Index]

[ferret_users] different results when calculate climatology



Dear ferreters,
     The following script provide two ways to calculate climatology. But I am confused with different results.
 
use atl.nc
set reg/y=43
show data
        V        Meridional Component of Velocit    1:61     1:31      1:23      1:156
        U        Zonal Component of Velocity         1:61     1:30      1:23      1:156
 climatological_axes
CANCEL DATA climatological_axes
 let v_clm = v[gt=month_reg@MOD]
 
! the first way to get Jan climatology
 fill v_clm[l=1]
 
! the second way to get Jan climatology (the most stupid way)
 let a=v[l=1]+v[l=13]+v[l=25]+v[l=37]+v[l=49];
 let b=v[l=61]+v[l=73]+v[l=85]+v[l=97]+v[l=109];
 let c=v[l=121]+v[l=133]+v[l=145]
 
let tt=(a+b+c)/13
 
fill tt
 
Of course, it has the similar results between the figure of "fill v_clm[l=1]"  and that of "fill tt", but still different.
 
Why? The calendar of the data is Gregorian, it might be the reason?
 
Dawn
 
 
 
 



Expecting? Get great news right away with email Auto-Check.
Try the Yahoo! Mail Beta.
[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement