print 'got grid_id:',grid_id
myaxes[2] = grid_id
mapnm = 'polar_stereographic'
params = [ "standard_parallel",
"latitude_of_projection_origin",
"false_easting",
"false_northing",
"straight_vertical_longitude_from_pole",
"scale_factor_at_projection_origin"]
punits = ["","","","","","" ]
pvalues = [71.,90.,0.,0.,-39., 1. ]
cmor.set_grid_mapping(grid_id=myaxes[2],
mapping_name = mapnm,
parameter_names = params,
parameter_values = pvalues,
parameter_units = punits)
cmor.set_table(tables[1])
pass_axes = [myaxes[2]]
data = f['orog'][0,:]
myvars[0] = cmor.variable( table_entry = te,
units = u,
axis_ids = pass_axes,
missing_value = data.missing,
history = '',
comment = ''
)
cmor.write(myvars[0], data)
def testCMIP6(self):
# Create x and y coordinates and bounds.
#
coordVals = numpy.arange(-533750.0, 533750.0 +
2500.0, 2500.0, numpy.float32)
coordBnds = numpy.zeros((coordVals.shape[0], 2), numpy.float32)
coordBnds[:, 0] = coordVals - 1250.0
coordBnds[:, 1] = coordVals + 1250.0
# Create longitude and latitude fields for a polar stereographic projection.
#
xgrid, ygrid = numpy.broadcast_arrays(
numpy.expand_dims(
coordVals, 0), numpy.expand_dims(
coordVals, 1))
rhogrid = numpy.sqrt(xgrid**2 + ygrid**2)
cgrid = 2.0 * numpy.arctan((0.5 / 6378137.0) * rhogrid)
latgrid = (180.0 / 3.141592654) * numpy.arcsin(numpy.cos(cgrid))
longrid = (180.0 / 3.141592654) * numpy.arctan2(xgrid, -ygrid)
# Set up CMOR with information from the CMOR config dictionary.
#
cmor.setup("Tables", netcdf_file_action=cmor.CMOR_REPLACE_4)
# Create the output CMOR dataset using the output configuration.
#
cmor.dataset_json("Test/common_user_input.json")
# Load the grid table.
#
cmor_table_obj = cmor.load_table("CMIP6_grids.json")
# Create ygre and xgre axes.
#
entry = {
'table_entry': 'y',
'units': 'm',
'coord_vals': coordVals,
'cell_bounds': coordBnds}
axis_ygre = cmor.axis(**entry)
entry = {
'table_entry': 'x',
'units': 'm',
'coord_vals': coordVals,
'cell_bounds': coordBnds}
axis_xgre = cmor.axis(**entry)
# Create the grid
#
grid_id = cmor.grid(
axis_ids=[
axis_ygre,
axis_xgre],
latitude=latgrid,
longitude=longrid)
# Set the CMOR grid mapping.
#
# mapnm = 'polar_stereographic'
param_dict = {
'latitude_of_projection_origin': [90.0, 'degrees_north'],
'straight_vertical_longitude_from_pole': [135.0, 'degrees_east'],
'standard_parallel': [70.0, 'degrees_north'],
'false_northing': [0.0, 'meters'],
'false_easting': [0.0, 'meters']
}
ierr = cmor.set_grid_mapping(
grid_id, 'polar_stereographic', param_dict)
self.assertEqual(ierr, 0)
print('got grid_id:', grid_id)
myaxes[2] = grid_id
mapnm = 'polar_stereographic'
params = ["standard_parallel",
"latitude_of_projection_origin",
"false_easting",
"false_northing",
"straight_vertical_longitude_from_pole",
"scale_factor_at_projection_origin"]
punits = ["", "", "", "", "", ""]
pvalues = [71., 90., 0., 0., -39., 1.]
cmor.set_grid_mapping(grid_id=myaxes[2],
mapping_name=mapnm,
parameter_names=params,
parameter_values=pvalues,
parameter_units=punits)
cmor.set_table(tables[1])
pass_axes = [myaxes[2]]
data = f['orog'][0, :]
myvars[0] = cmor.variable(table_entry=te,
units=u,
axis_ids=pass_axes,
missing_value=data.missing,
history='',
comment=''
)
cmor.write(myvars[0], data)