[CF-metadata] How to include EPSG codes or WKT information in a CF file [SEC=UNCLASSIFIED]

From: Timothy Hume <T.Hume>
Date: Mon, 30 Jul 2012 11:36:19 +1000

Hi,

I am adding grid mapping information to my NetCDF files. My data are on a lon-lat grid on a spherical Earth. I can put in the standard CF attributes, grid_mapping_name etc, but also want to include the EPSG and WKT codes to help software applications that might want these. I have not been able to find out how this is done in NetCDF-CF, so am using this format:
=
int crs ;
                crs:grid_mapping_name = "latitude_longitude" ;
                crs:earth_radius = 6371007.0 ;
                crs:prime_meridian = 0.0 ;
                crs:epsg_code = 4047 ;
                crs:wkt = "GEOGCS[\"Unspecified datum based upon the GRS 1980 Authalic Sphere\",DATUM[\"Not_specified_based_on
_GRS_1980_Authalic_Sphere\",SPHEROID[\"GRS 1980 Authalic Sphere\",6371007,0,AUTHORITY[\"EPSG\",\"7048\"]],AUTHORITY[\"EPSG\",\
"6047\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4047\"]]" ;

What do other people do when including this extra grid information?

Cheers,

Tim Hume
Centre for Australian Weather and Climate Research
Bureau of Meteorology
Melbourne
Australia
Received on Sun Jul 29 2012 - 19:36:19 BST