[ncl-talk] Manage netCDF metadata

[Mark Chan]

Dear Dennis,
Thanks so much, and it works greatly.
Cheers,Mark 

     On Wednesday, December 10, 2014 4:18 PM, Dennis Shea <shea at ucar.edu> wrote:
   

 http://www.ncl.ucar.edu/Document/Manuals/language_man.pdf

See the section:    Writing scalars to netCDF

=========
height   = ...height!0 = "ncl_scalar"    ; <=== 'special' name

rotated_pole    = ...
rotated_pole!0 = "ncl_scalar"
 
On Wed, Dec 10, 2014 at 7:23 AM, Mark Chan <cym263 at yahoo.com> wrote:

Dear NCL group,
In the below is about the metadata of a netCDF file, in which each of the variables has the dimension info such as "bnds = 2". Is it any possible in NCL scripts that some dimension info could be omitted such as "skip" and "ncl5", so the variable "rotated_pole" would be as: ("ncl5" and "skip" won't show up)

    char rotated_pole ;
        rotated_pole:grid_mapping_name = "rotated_latitude_longitude" ;
        rotated_pole:grid_north_pole_latitude = 90.f ;
        rotated_pole:grid_north_pole_longitude = -180.f ;
Thanks very much in advance!

Best regards,Mark

===netCDF metadata generated by NCL===
dimensions:
    time = UNLIMITED ; // (12 currently)
    rlon = 213 ;
    rlat = 220 ;
    skip = 1 ;
    bnds = 2 ;
    ncl5 = 1 ;
variables:
    double time(time) ;
        time:long_name = "time" ;
        time:standard_name = "time" ;
        time:units = "days since 1949-12-01 00:00:00" ;
        time:calendar = "standard" ;
        time:axis = "T" ;
        time:bounds = "time_bnds" ;
    double rlon(rlon) ;
        rlon:long_name = "longitude in rotated pole grid" ;
        rlon:standard_name = "grid_longitude" ;
        rlon:units = "degrees" ;
        rlon:axis = "X" ;
    double rlat(rlat) ;
        rlat:long_name = "latitude in rotated pole grid" ;
        rlat:standard_name = "grid_latitude" ;
        rlat:units = "degrees" ;
        rlat:axis = "Y" ;
    char rotated_pole(skip) ;
        rotated_pole:grid_mapping_name = "rotated_latitude_longitude" ;
        rotated_pole:grid_north_pole_latitude = 90.f ;
        rotated_pole:grid_north_pole_longitude = -180.f ;
    float tas(time, rlat, rlon) ;
        tas:long_name = "Near-Surface Air Temperature" ;
        tas:standard_name = "air_temperature" ;
        tas:units = "K" ;
        tas:cell_methods = "time: mean" ;
        tas:_FillValue = 1.e+20f ;
        tas:missing_value = 1.e+20f ;
        tas:coordinates = "lon lat height" ;
    double time_bnds(time, bnds) ;
    double height(ncl5) ;
        height:long_name = "height" ;
        height:standard_name = "height" ;
        height:units = "m" ;
        height:positive = "up" ;
        height:axis = "Z" ;
    double lon(rlat, rlon) ;
        lon:long_name = "longitude" ;
        lon:standard_name = "longitude" ;
        lon:units = "degrees_east" ;
    double lat(rlat, rlon) ;
        lat:long_name = "latitude" ;
        lat:standard_name = "latitude" ;
        lat:units = "degrees_north" ;



_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk





   
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141210/cdd69a99/attachment.html