[Thread Prev][Thread Next][Index]

Re: [ferret_users] Ferret does not recognize lon/lat axis



Hi Olivier,
Are you saying that you are running pyferret 1.2 on a MAC?  We have not built or tested pyferret under Mac OSX.  So that may be the trouble when trying pyferret.

The OSX system I have available is 10.10.2, and running ferret v6.94, it's reading this file correctly:
 
yes? use hc300_WOA13_decav_Reg1_clim.nc
yes? sh dat
     currently SET data sets:
    1> hc300_WOA13_decav_Reg1_clim.nc  (default)
 name     title                             I         J         K         L         M         N
 HC300    Heat content 0-300m              1:360     1:180     ...       1:12      ...       ...
 
yes? sh grid HC300
    GRID GOL1
 name       axis              # pts   start                end
 LON       LONGITUDE          360mr   179.5W(-179.5)       179.5E
 LAT       LATITUDE           180 r   89.5S                89.5N
 normal    Z
 TIME_COUNTER T (month)        12 r   1                    12
 normal    E
 normal    F
The order in which dimensions or attributes appear in the file does not matter.  The file works fine under linux running both ferret v6.96 and pyferret 1.2

Ansley

On 2/17/2016 5:44 AM, Olivier Marti wrote:
Hi,

When plotting variable hc300 from hc300_WOA13_decav_Reg1_clim.nc (ncdump -h included), with ferret v6.94 and pyferret 1.2 (ferret 6.96) on my Mac (OS X 10.11.3), the longitude and latitude axis are not recognized. If I run ncks -v hc300_WOA13_decav_Reg1_clim.nc toto.nc, the axis in toto.nc (ncdump -h included) are correctly recognized and interpreted. The difference between the two files seems to be only the order of the variables and the order of the attributes.

With ferret v6.94 on a Linux system, the axis are correctly interpreted.

Is this a bug or a feature, a problem with my install, etc ... ?

The file is written by python/cdms2, and I have no real control of the order of the attributes.


Olivier


marti@Spip-14:24:06-[Unix]/Analyses 👉  ncdump -h hc300_WOA13_decav_Reg1_clim.nc
netcdf hc300_WOA13_decav_Reg1_clim {
dimensions:
time_counter = UNLIMITED ; // (12 currently)
lat = 180 ;
lon = 360 ;
variables:
float time_counter(time_counter) ;
time_counter:units = "month" ;
float lat(lat) ;
lat:units = "degrees_north" ;
lat:standard_name = "latitude" ;
lat:axis = "Y" ;
lat:long_name = "Latitude" ;
float lon(lon) ;
lon:modulo = 360. ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
lon:topology = "circular" ;
float hc300(time_counter, lat, lon) ;
hc300:name = "hc300" ;
hc300:coordinates = "time_counter lat lon" ;
hc300:long_name = "Heat content 0-300m" ;
hc300:standard_name = "integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" ;
hc300:units = "J/m2" ;
hc300:missing_value = 1.e+20f ;

// global attributes:
:Conventions = "CF-1.0" ;
}


marti@Spip-14:24:09-[Unix]/Analyses 👉  ncdump -h toto.nc
netcdf toto {
dimensions:
time_counter = UNLIMITED ; // (12 currently)
lat = 180 ;
lon = 360 ;
variables:
float hc300(time_counter, lat, lon) ;
hc300:name = "hc300" ;
hc300:coordinates = "time_counter lat lon" ;
hc300:long_name = "Heat content 0-300m" ;
hc300:standard_name = "integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" ;
hc300:units = "J/m2" ;
hc300:missing_value = 1.e+20f ;
float lat(lat) ;
lat:units = "degrees_north" ;
lat:standard_name = "latitude" ;
lat:axis = "Y" ;
lat:long_name = "Latitude" ;
float lon(lon) ;
lon:modulo = 360. ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
lon:topology = "circular" ;
float time_counter(time_counter) ;
time_counter:units = "month" ;

// global attributes:
:Conventions = "CF-1.0" ;
:history = "Wed Feb 17 14:23:27 2016: ncks -v hc300 hc300_WOA13_decav_Reg1_clim.nc toto.nc" ;
:NCO = "4.4.8" ;
}

-- 
Dr. Olivier Marti  - mailto:olivier.marti@xxxxxxxxxxxx
Laboratoire des Sciences du Climat et de l'Environnement
Institut Pierre Simon Laplace
Laboratoire CEA-CNRS-UVSQ  -  UMR 8212
Tel : +33 1 69 08 77 27 - Mob : +33 6 45 36 43 74
Fax : +33 1 69 08 30 73















[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement