[Thread Prev][Thread Next][Index]

Re: error in line definition



> 
> Dear Ferret users,
> 
> I use ferret for adapting inputs ECMWF to my grid.
> Here is the script I am using:
> 
> use grilles_gcm.nc
> use shum.94.nc
> define grid/like=shum[d=2]/x=grille_v[d=1]/y=grille_u[d=1] grille_h
> repeat/l=1181:1184 save/file=hur.nc/append shum[d=2,g=grille_h,i=1:97,j=1:73]
> 
> and here is the output:
> 
> !-> REPEAT: L=1181
>  LISTing to file hur.nc
>  **TMAP ERR: error in line definition
>              disordered output coordinate value: 1.74774E+07 Axis: TIME
> !-> REPEAT: L=1182
>  LISTing to file hur.nc
>  **TMAP ERR: error in line definition
>              disordered output coordinate value: 1.74774E+07 Axis: TIME
> !-> REPEAT: L=1183
>  LISTing to file hur.nc
>  **TMAP ERR: error in line definition
>              disordered output coordinate value: 1.74774E+07 Axis: TIME
> !-> REPEAT: L=1184
>  LISTing to file hur.nc
>  **TMAP ERR: error in line definition
>              disordered output coordinate value: 1.74774E+07 Axis: TIME
> 
> Thank you for your help.
> 

It's quite possible that you are running into precision problems. Consider
the following fortran90 program

write(*,*) nearest(1.74774e7,1.)
write(*,*) nearest(1.74774e7,-1.)
end

For 32 bit IEEE arithmetic I get the answer

17477402.
17477398.

So we see that the minimum increment is 4 in this region. It's possible that
ferret thinks that all your time points are identical. Try defining a new axis
with the time origin closer to the time of interest and map via subscript
association (@asn).

Cheers,
Russ


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement