[Thread Prev][Thread Next][Index]

[ferret_users] @sum bug ?



Hi all,

I am just discovering that the @sum operator gives different results
from 6.61 and next release 6.62, 6.64.
I haven't any new about this serious bug in release notes.

My test are available from:
http://dods.ipsl.jussieu.fr/brocksce/ferret_bugs/err661_sum.jnl
http://dods.ipsl.jussieu.fr/brocksce/ferret_bugs/err661_sum_file.nc

And gives:
######################################"
$ ferret -script err661_sum.jnl
use err661_sum_file.nc
list/p=10 tas[i=@sum,j=@sum]
            VARIABLE : TAS
            FILENAME : err661_sum_file.nc
            X        : 0.5 to 96.5 (summed)
            Y        : 0.5 to 96.5 (summed)
         2526427.750
let a=xsequence(tas)
list/p=10 a[i=@sum]
            VARIABLE : XSEQUENCE(TAS)
            FILENAME : err661_sum_file.nc
            X        : 0.5 to 9216.5 (summed)
         2526427.750

And with 6.62 and 6.64
######################################"
use err661_sum_file.nc
list/p=10 tas[i=@sum,j=@sum]
            VARIABLE : TAS
            FILENAME : err661_sum_file.nc
            X        : 0.5 to 96.5 (summed)
            Y        : 0.5 to 96.5 (summed)
         2526427.250
let a=xsequence(tas)
list/p=10 a[i=@sum]
            VARIABLE : XSEQUENCE(TAS)
            FILENAME : err661_sum_file.nc
            X        : 0.5 to 9216.5 (summed)
         2526434.500

######################################"
From python/netCDF4/numpy
I know that the sum is really 2526434.500


Other question related: why tas[i=@sum,j=@sum]
and xsequence then i=@sum gives different results.
Both should give 2526434.500

Hope it a bad dream...
Let me know if you can confirm this.

Patrick


[Thread Prev][Thread Next][Index]

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

Privacy Policy | Disclaimer | Accessibility Statement