[Thread Prev][Thread Next][Index]

Re: [ferret_users] netcdf reading



Hi -
There are some useful NetCDF file-check utilities on the internet, including the CF-Convention checker,

http://puma.nerc.ac.uk/cgi-bin/cf-checker.pl



On 10/18/2011 5:03 PM, Ryo Furue wrote:
Namba-san,

First of all, your email message looks corrupted and is
hard to read.  See the attached text, which is a copy
of your message (with my quotation characters on top of
each line).

However, I did manage to read it, and I found your netCDF file
is in a "wrong" format.  To learn how to create a "correct"
netCDF file, you may want to look at examples.  For example,
this is the header of "coads_climatology.cdf", which you find
in your Ferret data directory:

[/usr/local/ferret/data]$ ncdump -h coads_climatology.cdf
netcdf coads_climatology {
dimensions:
    COADSX = 180 ;
    COADSY = 90 ;
    TIME = UNLIMITED ; // (12 currently)
variables:
    double COADSX(COADSX) ;
        COADSX:units = "degrees_east" ;
        COADSX:modulo = " " ;
        COADSX:point_spacing = "even" ;
    double COADSY(COADSY) ;
        COADSY:units = "degrees_north" ;
        COADSY:point_spacing = "even" ;
    double TIME(TIME) ;
        TIME:units = "hour since 0000-01-01 00:00:00" ;
        TIME:time_origin = "1-JAN-0000 00:00:00" ;
        TIME:modulo = " " ;
    float SST(TIME, COADSY, COADSX) ;
        SST:missing_value = -1.e+34f ;
        SST:_FillValue = -1.e+34f ;
  . . . . . . .

This is how the coordinates, dimensions, and variables have to be
related to each other.

Hope this helps,
Ryo

| Dear Sir
| I have the netcdf like
| ncdump -h 19930101.nc
| netcdf \19930101 {dimensions: longitude = 801 ; latitude = 251 ; day = 101 ;variables: float xlongitude(longitude, latitude, day) ; xlongitude:units = "degrees_north" ; float ylatitude(longitude, latitude, day) ; ylatitude:units = "degrees_east" ; float temperature(longitude, latitude, day) ; temperature:units = "celsius" ; float salinity(longitude, latitude, day) ; salinity:units = "psu" ; float mixlayer(longitude, latitude, day) ; mixlayer:units = "meter" ;}
| When I do in ferret
| yes? use 19930101.ncyes? show data currently SET data sets: 1> ./19930101.nc (default) name title I J K L XLONGITUDE 1:101 1:251 1:801 ... YLATITUDE 1:101 1:251 1:801 ... TEMPERATURE 1:101 1:251 1:801 ... SALINITY 1:101 1:251 1:801 ... MIXLAYER 1:101 1:251 1:801 ... | I expected and I would like to do that I is longitude 1:801 and J is latitude 1:251 and K is day 1:101 as the ncdump shows,
| Could you give me how to do it?
| Best Regards,
| Takaya Namba
|
|
|
|
|
|
|


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

Privacy Policy | Disclaimer | Accessibility Statement