Hi Patrick, In Ferret, coordinate axes are stored in double precision internally, but all variables, like tt and x1 and x2 are stored and handled in single precision. So when you evaluate a time step and assign its value to a variable, the accuracy of the single-precision variables isn't enough to express the number. Single precision has only about 6 digits of accuracy, so even though the digits are written to the Ferret output, the last couple of digits are meaningless. If you write the data to a netCDF file, the double precision of the coordinates is preserved: def axis/T="01-JAN-2008":"31-DEC-2010":1/units="seconds"/t0="01-JAN-2008"/cal=GREGORIAN mytaxisEven in the file, the values of the variable tt are in single precision, so the values of tt do not match the double precision values of the coordinate variable mytaxis. We just don't have the ability in Ferret to handle or compute with double precision variables. Ansley Brockmann Patrick wrote: Hi all, |