One more time - making climatologies and anomalies...

[strandwg@xxxxxxxx (Gary Strand)]

I've read over the doc in the manual and on the web, and it still don't work!

To wit:

! Example data set: taux.nc (T42, 60 months of data, however, time_origin =
! "01-SEP-0000 00:00:00", a problem?)

SET DATA taux.nc
define grid/like=taux tx_grid

!
use climatological_axes
canc data climatological_axes
LET taux_climatology = taux[GT=month_reg@MOD]
!
save/file=taux_anom.nc taux_climatology
stat taux_climatology

             TAUX[GT=MONTH_REG@MOD]
             LONGITUDE: 1.4W(-1.4) to 1.4W(358.6)
             LATITUDE: 89.2S to 89.2N
             Z:  N/A
             TIME: 01-JAN      00:45 to 31-DEC      06:34
             DATA SET: ./taux.nc
 
 Total # of data points: 98304 (128*64*1*12)
 # flagged as bad  data: 0
 Minimum value: -2.1242
 Maximum value: 0.78542
 Mean    value: -0.011204    (unweighted average)
 Standard deviation: 0.11109
!
let taux_anomaly = taux - taux_climatology[gt=tx_grid]
save/append/file=taux_anom.nc taux_anomaly
stat taux_anomaly

             TAUX - TAUX_CLIMATOLOGY[GT=TX_GRID]
             LONGITUDE: 1.4W(-1.4) to 1.4W(358.6)
             LATITUDE: 89.2S to 89.2N
             Z:  N/A
             TIME: 15-JAN-0014 00:00 to 12-JAN-0019 12:00
             DATA SET: ./taux.nc
 
 Total # of data points: 491520 (128*64*1*60)
 # flagged as bad  data: 491520

-------------------------------------------------------------------------------

Not good! What happened?

--
Gary Strand
strandwg@ucar.edu                              http://www.cgd.ucar.edu/ccr/gary