Commit 88f7db70 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Added doc/Gridmapping

parent 02bbd0b4
=======================================
Example 1: rotated_latitude_longitude
=======================================
rlon = 106 ;
rlat = 111 ;
float lon(rlat, rlon) ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:standard_name = "grid_longitude" ;
float lat(rlat, rlon) ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:standard_name = "grid_latitude" ;
float rlat(rlat) ;
rlat:standard_name = "grid_latitude" ;
rlat:long_name = "rotated latitude" ;
rlat:units = "degrees" ;
float rlon(rlon) ;
rlon:standard_name = "grid_longitude" ;
rlon:long_name = "rotated longitude" ;
rlon:units = "degrees" ;
char rotated_pole ;
rotated_pole:long_name = "coordinates of the rotated North Pole" ;
rotated_pole:grid_mapping_name = "rotated_latitude_longitude" ;
rotated_pole:grid_north_pole_latitude = 39.25f ;
rotated_pole:grid_north_pole_longitude = -162.f ;
float T_2M_AV(time, height_2m, rlat, rlon) ;
T_2M_AV:standard_name = "air_temperature" ;
T_2M_AV:units = "K" ;
T_2M_AV:coordinates = "lon lat" ;
gridtype = GRIB_CURVILINEAR
xsize = 334
ysize = 334
xname = lon
xvals = ...
yname = lat
yvals = ...
xdimname = rlon
ydimname = rlat
projection = rotated_latitude_longitude
pname = rotated_pole
pxvals = ....
pyvals = ....
grid_north_pole_latitude = 39.25f
grid_north_pole_longitude = -162.f
long_name = "coordinates of the rotated North Pole"
=======================================
Example 2: lambert_azimuthal_equal_area
=======================================
x = 334 ;
y = 334 ;
double x(x) ;
x:standard_name = "projection_x_coordinate" ;
x:units = "m" ;
x:axis = "X" ;
double y(y) ;
y:standard_name = "projection_y_coordinate" ;
y:units = "m" ;
y:axis = "Y" ;
char laea ;
laea:grid_mapping_name = "lambert_azimuthal_equal_area" ;
laea:earth_radius = 6370997. ;
laea:longitude_of_projection_origin = 0. ;
laea:latitude_of_projection_origin = 90. ;
float CFC_FC_F(time, y, x) ;
CFC_FC_F:grid_mapping = "laea" ;
gridtype = GRIB_CURVILINEAR
xsize = 334
ysize = 334
projection = lambert_azimuthal_equal_area
pxvals = ....
pyvals = ....
earth_radius = 6370997.
longitude_of_projection_origin = 0.
latitude_of_projection_origin = 90.
cdo gen_geoloc ifile ofile
//cdo del_geoloc ifile ofile
gridtype = GRIB_CURVILINEAR
xsize = 334
ysize = 334
xvals = ....
yvals = ....
xbounds = ....
ybounds = ....
projection = lambert_azimuthal_equal_area
pxvals = ....
pyvals = ....
earth_radius = 6370997.
longitude_of_projection_origin = 0.
latitude_of_projection_origin = 90.
cdo del_proj ifile ofile
gridtype = GRIB_CURVILINEAR
xsize = 334
ysize = 334
xvals = ....
yvals = ....
xbounds = ....
ybounds = ....
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment