[Thread Prev][Thread Next][Index]
Re: file conflicts
Hi Matt,
Thanks for the trouble report.
It turns out that the problem was reported and fixed last October, so it V5.3 BETA already
has the problem fixed.
The problem only shows up if
* an axis name is identical in the two files
* both of the axes are the same length and have *irregularly* spaced points
Ferret version 5.22 mistakenly assumes that the axes are identical in the two files.
thanks for the report. Keep 'em coming! - steve
============
Matthew Harrison wrote:
> I've never encountered this using Ferret, so I'm wondering if anyone
> else has. Using v5.22 on IRIX6.5 and linux.
>
> Two files (attached) are identical with the exception of the time values
> and the filename. Ferret seems to be confusing the time and assigns the
> time values from the 1st file to the second file. If I open the files
> separately, ferret gets the time values correct. This problem does not
> duplicate in v4.91.
>
> ------------------------------------------------------------------------
>
> use file1.nc ! 1981 data
> use file2.nc ! 1982 data
>
> list t[d=1,g=temp]
>
> axis TIME
> DATA SET: ./file1.nc
> 16-JAN-1981 12 / 1: 747.
> 15-FEB-1981 00 / 2: 776.
> 16-MAR-1981 12 / 3: 806.
> 16-APR-1981 00 / 4: 836.
> 16-MAY-1981 12 / 5: 867.
> 16-JUN-1981 00 / 6: 897.
> 16-JUL-1981 12 / 7: 928.
> 16-AUG-1981 12 / 8: 959.
> 16-SEP-1981 00 / 9: 989.
> 16-OCT-1981 12 / 10: 1020.
> 16-NOV-1981 00 / 11: 1050.
> 16-DEC-1981 12 / 12: 1081.
>
> list t[d=2,g=temp]
>
> axis TIME
> DATA SET: ./file2.nc
> 16-JAN-1981 12 / 1: 747.
> 15-FEB-1981 00 / 2: 776.
> 16-MAR-1981 12 / 3: 806.
> 16-APR-1981 00 / 4: 836.
> 16-MAY-1981 12 / 5: 867.
> 16-JUN-1981 00 / 6: 897.
> 16-JUL-1981 12 / 7: 928.
> 16-AUG-1981 12 / 8: 959. !!! should be 1982 dates
> 16-SEP-1981 00 / 9: 989.
> 16-OCT-1981 12 / 10: 1020.
> 16-NOV-1981 00 / 11: 1050.
> 16-DEC-1981 12 / 12: 1081.
>
> --
> Matthew J. Harrison GFDL
> Physical Scientist P.O. Box 308
> NOAA Geophysical Fluid Dynamics Laboratory Princeton, NJ
> 08542
> ----------------------------------------------------------------------------
>
> e-mail: mjh@gfdl.gov
> Tel: (609)452-6579
> Fax: (609)987-5063
> http address: http://www.gfdl.gov/~mjh
>
> ------------------------------------------------------------------------
> netcdf file1 {
> dimensions:
> lat = 90 ;
> level = 17 ;
> lon = 144 ;
> time = UNLIMITED ; // (12 currently)
> variables:
> float lat(lat) ;
> lat:long_name = "latitude" ;
> lat:cartesian_axis = "Y" ;
> lat:units = "degrees_N" ;
> lat:edges = "latb" ;
> float level(level) ;
> level:long_name = "pressure" ;
> level:cartesian_axis = "Z" ;
> level:units = "hPa" ;
> float lon(lon) ;
> lon:long_name = "longitude" ;
> lon:cartesian_axis = "X" ;
> lon:units = "degrees_E" ;
> lon:edges = "lonb" ;
> short temp(time, level, lat, lon) ;
> temp:long_name = "temperature" ;
> temp:units = "deg_k" ;
> temp:valid_range = -32750s, 32750s ;
> temp:missing_value = -32751s ;
> temp:scale_factor = 0.004580152f ;
> temp:add_offset = 250.f ;
> temp:time_avg_info = "average_T1,average_T2,average_DT" ;
> double time(time) ;
> time:long_name = "time" ;
> time:cartesian_axis = "T" ;
> time:units = "days since 1979-01-01 00:00:00" ;
> time:calendar_type = "JULIAN" ;
>
> // global attributes:
> :filename = "1981feb01h00.land_month.nc" ;
> :MPP_IO_VERSION = "$Id: mpp_io.F90,v 5.5 2000/07/28 20:17:19 fms Exp $" ;
> :title = "Diagnostics from an FMS integration" ;
> :history = "Tue Jan 30 14:11:14 2001: ncks -v temp file1.nc out.nc" ;
> data:
>
> lat = -89, -87, -85, -83, -81, -79, -77, -75, -73, -71, -69, -67, -65, -63,
> -61, -59, -57, -55, -53, -51, -49, -47, -45, -43, -41, -39, -37, -35,
> -33, -31, -29, -27, -25, -23, -21, -19, -17, -15, -13, -11, -9, -7, -5,
> -3, -1, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33,
> 35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67, 69,
> 71, 73, 75, 77, 79, 81, 83, 85, 87, 89 ;
>
> level = 1000, 925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50,
> 30, 20, 10 ;
>
> lon = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25, 18.75, 21.25, 23.75,
> 26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75, 46.25, 48.75,
> 51.25, 53.75, 56.25, 58.75, 61.25, 63.75, 66.25, 68.75, 71.25, 73.75,
> 76.25, 78.75, 81.25, 83.75, 86.25, 88.75, 91.25, 93.75, 96.25, 98.75,
> 101.25, 103.75, 106.25, 108.75, 111.25, 113.75, 116.25, 118.75, 121.25,
> 123.75, 126.25, 128.75, 131.25, 133.75, 136.25, 138.75, 141.25, 143.75,
> 146.25, 148.75, 151.25, 153.75, 156.25, 158.75, 161.25, 163.75, 166.25,
> 168.75, 171.25, 173.75, 176.25, 178.75, 181.25, 183.75, 186.25, 188.75,
> 191.25, 193.75, 196.25, 198.75, 201.25, 203.75, 206.25, 208.75, 211.25,
> 213.75, 216.25, 218.75, 221.25, 223.75, 226.25, 228.75, 231.25, 233.75,
> 236.25, 238.75, 241.25, 243.75, 246.25, 248.75, 251.25, 253.75, 256.25,
> 258.75, 261.25, 263.75, 266.25, 268.75, 271.25, 273.75, 276.25, 278.75,
> 281.25, 283.75, 286.25, 288.75, 291.25, 293.75, 296.25, 298.75, 301.25,
> 303.75, 306.25, 308.75, 311.25, 313.75, 316.25, 318.75, 321.25, 323.75,
> 326.25, 328.75, 331.25, 333.75, 336.25, 338.75, 341.25, 343.75, 346.25,
> 348.75, 351.25, 353.75, 356.25, 358.75 ;
>
> time = 746.5, 776, 805.5, 836, 866.5, 897, 927.5, 958.5, 989, 1019.5, 1050,
> 1080.5 ;
> }
>
> ------------------------------------------------------------------------
> netcdf file2 {
> dimensions:
> lat = 90 ;
> level = 17 ;
> lon = 144 ;
> time = UNLIMITED ; // (12 currently)
> variables:
> float lat(lat) ;
> lat:long_name = "latitude" ;
> lat:cartesian_axis = "Y" ;
> lat:units = "degrees_N" ;
> lat:edges = "latb" ;
> float level(level) ;
> level:long_name = "pressure" ;
> level:cartesian_axis = "Z" ;
> level:units = "hPa" ;
> float lon(lon) ;
> lon:long_name = "longitude" ;
> lon:cartesian_axis = "X" ;
> lon:units = "degrees_E" ;
> lon:edges = "lonb" ;
> short temp(time, level, lat, lon) ;
> temp:long_name = "temperature" ;
> temp:units = "deg_k" ;
> temp:valid_range = -32750s, 32750s ;
> temp:missing_value = -32751s ;
> temp:scale_factor = 0.004580152f ;
> temp:add_offset = 250.f ;
> temp:time_avg_info = "average_T1,average_T2,average_DT" ;
> double time(time) ;
> time:long_name = "time" ;
> time:cartesian_axis = "T" ;
> time:units = "days since 1979-01-01 00:00:00" ;
> time:calendar_type = "JULIAN" ;
>
> // global attributes:
> :filename = "1982feb01h00.land_month.nc" ;
> :MPP_IO_VERSION = "$Id: mpp_io.F90,v 5.5 2000/07/28 20:17:19 fms Exp $" ;
> :title = "Diagnostics from an FMS integration" ;
> :history = "Tue Jan 30 14:11:34 2001: ncks -v temp file2.nc out.nc" ;
> data:
>
> lat = -89, -87, -85, -83, -81, -79, -77, -75, -73, -71, -69, -67, -65, -63,
> -61, -59, -57, -55, -53, -51, -49, -47, -45, -43, -41, -39, -37, -35,
> -33, -31, -29, -27, -25, -23, -21, -19, -17, -15, -13, -11, -9, -7, -5,
> -3, -1, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33,
> 35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67, 69,
> 71, 73, 75, 77, 79, 81, 83, 85, 87, 89 ;
>
> level = 1000, 925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50,
> 30, 20, 10 ;
>
> lon = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25, 18.75, 21.25, 23.75,
> 26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75, 46.25, 48.75,
> 51.25, 53.75, 56.25, 58.75, 61.25, 63.75, 66.25, 68.75, 71.25, 73.75,
> 76.25, 78.75, 81.25, 83.75, 86.25, 88.75, 91.25, 93.75, 96.25, 98.75,
> 101.25, 103.75, 106.25, 108.75, 111.25, 113.75, 116.25, 118.75, 121.25,
> 123.75, 126.25, 128.75, 131.25, 133.75, 136.25, 138.75, 141.25, 143.75,
> 146.25, 148.75, 151.25, 153.75, 156.25, 158.75, 161.25, 163.75, 166.25,
> 168.75, 171.25, 173.75, 176.25, 178.75, 181.25, 183.75, 186.25, 188.75,
> 191.25, 193.75, 196.25, 198.75, 201.25, 203.75, 206.25, 208.75, 211.25,
> 213.75, 216.25, 218.75, 221.25, 223.75, 226.25, 228.75, 231.25, 233.75,
> 236.25, 238.75, 241.25, 243.75, 246.25, 248.75, 251.25, 253.75, 256.25,
> 258.75, 261.25, 263.75, 266.25, 268.75, 271.25, 273.75, 276.25, 278.75,
> 281.25, 283.75, 286.25, 288.75, 291.25, 293.75, 296.25, 298.75, 301.25,
> 303.75, 306.25, 308.75, 311.25, 313.75, 316.25, 318.75, 321.25, 323.75,
> 326.25, 328.75, 331.25, 333.75, 336.25, 338.75, 341.25, 343.75, 346.25,
> 348.75, 351.25, 353.75, 356.25, 358.75 ;
>
> time = 1111.5, 1141, 1170.5, 1201, 1231.5, 1262, 1292.5, 1323.5, 1354,
> 1384.5, 1415, 1445.5 ;
> }
--
| NOAA/PMEL | ph. (206) 526-6080
Steve Hankin | 7600 Sand Point Way NE | FAX (206) 526-6744
| Seattle, WA 98115-0070 | hankin@pmel.noaa.gov
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement