Re: How to read this data correctly

This is not a particularly 'good' netCDF file.
The latitude & longitude units are non-standard.
The units "degrees" is not recognized.
'degrees' what ... temperature, angle, ?????

 From an NCL perspective, it can be modified upon input

    f = addfile("CPTdaily_minsmr1979.nc", "r")

    lat = f->vlat
    lat!0 = "lat"
    lat@units = "degrees_north"
    lat&lat = lat

    lon = f->vlon
    lon!0 = "lon"
    lon@units = "degrees_east"
    lon&lon = lon

    h = f->h
    h&lat = lat
    h&lon = lon
    printVarSummary(h)

On 10/22/2012 12:33 PM, Thomas Tobian wrote:
> Dear All,
>
> I am trying to plot the derived data that looks like as follow:
> When i tried to plot using NCL, I got wrong result and error.
> Anyone has suggestion how could i modify the format of this netcdf to
> the readable and correct one.
> Thanks.
>
> netcdf CPTdaily_minsmr1979 {
>
> dimensions:
> lat = 25 ;
> lon = 144 ;
> level = 1 ;
> time = UNLIMITED ; // (365 currently)
> variables:
> float h(time, level, lat, lon) ;
> h:units = "km" ;
> float T(time, level, lat, lon) ;
> T:units = "K" ;
> float p(time, level, lat, lon) ;
> p:units = "mbar" ;
> float SMR(time, level, lat, lon) ;
> SMR:units = "ppmv" ;
> *float vlon(lon) ;*
> *vlon:units = "degrees" ;*
> *vlon:no_button = 1 ;*
> *float vlat(lat) ;*
> *vlat:units = "degrees" ;*
> *vlat:no_button = 1 ;*
> *float time(time) ;*
> *time:units = "days" ;*
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Oct 22 14:21:45 2012