[Thread Prev][Thread Next][Index]

bad flag gets reset



Hi,
I'm computing the following and it appears that the missing vlaue defined in my netcdf file is getting
reset when I compute the ration of @var/@ave. The missingvalue is 1.0000000E+20 but when I compute @var/@ave the bad flag is set to -9.9999998E+33 (= -0.1000E+35)
I would have expected the missing value to remain the same, ie 1.0000000E+20.
Any ideas?
Here's the Netcdf header:

 obs.globe.lan.1901-2000.pre.annual {
dimensions:
        lon = 720 ;
        lat = 360 ;
        time = UNLIMITED ; // (100 currently)
variables:
        float lon(lon) ;
                lon:units = "degrees_east" ;
                lon:valid_min = -180.f ;
                lon:valid_max = 180.f ;
                lon:title = "Longitude" ;
        float lat(lat) ;
                lat:units = "degrees_north" ;
                lat:valid_min = -90.f ;
                lat:valid_max = 90.f ;
                lat:title = "Latitude" ;
        float time(time) ;
                time:units = "years since 1900-01-01" ;
                time:title = "Time" ;
                time:long_name = "Time axis" ;
                time:time_origin = " 1900-01-15 00:00:00" ;
        float pre(time, lat, lon) ;
                pre:units = "mm/month" ;
                pre:long_name = "Precipitation" ;
                pre:missing_value = 1.e+20f

yes? set region/i=1:1/j=1:1      !look at a grid cell that I know is missing
yes? list/l=1:1/format=(e15.6)  pre[l=@ave]
             VARIABLE : Precipitation (mm/month)
             FILENAME : obs.globe.lan.1901-2000.pre.annual.nc
             BAD FLAG :  1.0000000E+20
             LONGITUDE: 179.8W(-179.8)
             LATITUDE : 89.8S
             TIME     : 02-JUL-1900 14:54 to 02-JUL-1901 20:43 (averaged)
   0.100000E+21  (OK)
yes? list/l=1:1/format=(e12.4)  pre[l=@var]
             VARIABLE : Variance of Precipitation ((mm/month)^2)
             FILENAME : obs.globe.lan.1901-2000.pre.annual.nc
             BAD FLAG :  1.0000000E+20
             LONGITUDE: 179.8W(-179.8)
             LATITUDE : 89.8S
             TIME     : 02-JUL-1900 14:54 to 02-JUL-1901 20:43 (variance)
  0.1000E+21 (OK)
yes? list/l=1:1/format=(e12.4)  pre[l=@var]/pre[l=@ave]
             VARIABLE : PRE[L=@VAR]/PRE[L=@AVE]
             FILENAME : obs.globe.lan.1901-2000.pre.annual.nc
             BAD FLAG : -9.9999998E+33 (aka -0.1000E+35)
             LONGITUDE: 179.8W(-179.8)
             LATITUDE : 89.8S
             TIME     : 02-JUL-1900 14:54 to 02-JUL-1901 20:43
 -0.1000E+35 (NOT OK)

-- 
_________________________
Steve Knox
NREL
Colorado State University
(970)491-1972

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement