Remapping Analyses and Initial Conditions#

Analysis remapping, supported by the available remapping types (i.e., ESMF, SLINT, etc.,), can be accomplished by using the remapper application. For optimal results using the remapper application, the user should first generate the tripolar grid projections using the gridspec for the respective grid-scale resolutions. and subsequently compute the respective interpolate algorithm (i.e., ESMF or SLINT) remapping coefficients via esmf_remap.

Application Configuration#

An example remapper application configuration for a \(0.25^{\circ}\) MOM6 analysis remapping to a \(1.0^{\circ}\) MOM6 tripolar grid projection follows while the respective variables are defined in the subsequent tables.

forecast_model: mom6
app_module: remapper.models.mom6.mom6
app_class: MOM6
output_netcdf: /run/MOM6.ic.1p0.nc
nlevs: 75
bathy_file: /work/gridspec/1p0/mom6/mom6.topo.1p0.nc
bathy_edits_file: /work/gridspec/1p0/mom6/mom6.topo-edits.1p0.nc
grids:
  source: !INC /ufs_tools/parm/remapper/grids/mom6.0p25.yaml
  destination: !INC /ufs_tools/parm/remapper/grids/mom6.1p0.yaml
remap:
  interp_scheme: esmf
  rotate_currents: true
  dstmass2dstuvel_bilinear: /work/esmf_remap/tripolar.1p0_1p0.bilinear.t2u_remap.nc
  dstmass2dstvvel_bilinear: /work/esmf_remap/tripolar.1p0_1p0.bilinear.t2v_remap.nc
  srcmass2dstmass_bilinear: /work/esmf_remap/tripolar.0p25_1p0.bilinear.t2t_remap.nc
  srcmass2dstmass_nrstnghbr: /work/esmf_remap/tripolar.0p25_1p0.nearest_s2d.t2t_remap.nc
  srcuvel2dstmass_bilinear: /work/esmf_remap/tripolar.0p25_1p0.bilinear.u2t_remap.nc
  srcvvel2dstmass_bilinear: /work/esmf_remap/tripolar.0p25_1p0.bilinear.v2t_remap.nc
  srcuvel2srcmass_bilinear: /work/esmf_remap/tripolar.0p25_0p25.bilinear.u2t_remap.nc
  srcvvel2srcmass_bilinear: /work/esmf_remap/tripolar.0p25_0p25.bilinear.v2t_remap.nc
variables: !INC /ufs_tools/parm/remapper/variables/mom62mom6.yaml

Configuration Variable

Description

bathy_edits_file

This is an optional variable specifying the file path to a netCDF-formatted file containing bathymetry value edits; if not being used this should be set to Null or commented out.

bathy_file

The file path to the netCDF formatted file containing the destination grid projection bathymetry.

forecast_model

The forecast model from which the analysis was produced; currently supported are MOM6, CICE, ORAS5, and GLORYS sea-ice.

grids

The respective source and destination grid informations; see the corresponding table below for further information.

interp_scheme

The interpolation scheme to be used for the remapping; currently supported schemes are ESMF and SLINT; Note that the files provided beneath the YAML remap key must be generated from the same interpolation scheme.

nlevs

The total number of vertical levels to be remapped to the destination grid projection.

output_netcdf

The netCDF-formatted file path to the remapped (MOM6) analysis variables.

remap

The remapping attributes; see the table below for further information.

variables

The analysis variable remapping attributes; see the table below for more information.

The configuration keys for the respective grids attributes are as follows.

Configuration Variable

Description

arakawa

The Arakawa type for the respective grid (e.g., destination or source).

grid_ncfile

The path to the netCDF-formatted file containing the respective grid geographical locations.

mass, uvel, vvel

The respective grid coordinate types for which variables may be defined on an Arakawa type grid; the nclat and nclon variables are the netCDF-formatted file variable names for the respective latitude and longitude geographical locations.

topo_ncfile

The path to the netCDF-formatted file containing the respective grid topography and corresponding landmask.

The configuration keys for the respective remap attributes are as follows.

Configuration Variable

Description

dstuvel2dstmass_bilinear

Bilinear remapping coefficients from the destination uvel grid locations to the destination grid mass grid locations.

dstuvel2dstmass_nrstnghbr

Nearest-neighbor remapping coefficients from the destination uvel grid locations to the destination grid mass grid locations.

dstvvel2dstmass_bilinear

Bilinear remapping coefficients from the destination vvel grid locations to the destination grid mass grid locations.

dstvvel2dstmass_nrstnghbr

Nearest-neighbor remapping coefficients from the destination vvel grid locations to the destination grid mass grid locations.

dstmass2dstuvel_bilinear

Bilinear remapping coefficients from the destination mass grid locations to the destination grid uvel grid locations.

dstmass2dstuvel_nrstnghbr

Nearest-neighbor remapping coefficients from the destination mass grid locations to the destination grid uvel grid locations.

dstmass2dstvvel_bilinear

Bilinear remapping coefficients from the destination mass grid locations to the destination grid vvel grid locations.

dstmass2dstvvel_nrstnghbr

Nearest-neighbor remapping coefficients from the destination mass grid locations to the destination grid vvel grid locations.

rotate_currents

A boolean variable specifying whether the current vector components must be rotated following remapping; this is True for all staggered Arakawa grids.

srcmass2dstmass_bilinear

Bilinear remapping coefficients from the source mass grid locations to the destination grid mass grid locations.

srcmass2dstmass_nrstnghbr

Nearest-neighbor remapping coefficients from the source mass grid locations to the destination grid mass grid locations.

srcmass2dstuvel_bilinear

Bilinear remapping coefficients from the source mass grid locations to the destination grid uvel grid locations.

srcmass2dstuvel_nrstnghbr

Nearest-neighbor remapping coefficients from the source mass grid locations to the destination grid uvel grid locations.

srcmass2dstvvel_bilinear

Bilinear remapping coefficients from the source mass grid locations to the destination grid vvel grid locations.

srcmass2dstvvel_nrstnghbr

Nearest-neighbor remapping coefficients from the source mass grid locations to the destination grid vvel grid locations.

srcuvel2dstmass_bilinear

Bilinear remapping coefficients from the source uvel grid locations to the destination grid mass grid locations.

srcuvel2dstmass_nrstnghbr

Nearest-neighbor remapping coefficients from the source uvel grid locations to the destination grid mass grid locations.

srcvvel2dstmass_bilinear

Bilinear remapping coefficients from the source vvel grid locations to the destination grid mass grid locations.

srcvvel2dstmass_nrstnghbr

Nearest-neighbor remapping coefficients from the source vvel grid locations to the destination grid mass grid locations.

The configuration keys for the respective variables attributes are as follows.

Configuration Variable

Description

dst_ncvarname

The destination grid netCDF variable name.

grid_stagger

The destination grid grid-staggered locations (e.g., mass, uvel, or vvel).

interp_type

The interpolation type to be used for remapping; may be either bilinear or nrstnghbr.

ncfilename

The netCDF-formatted file path containing the respective source grid variable values.

src_ncvarname

The source grid netCDF variable name.

xdim_name

The netCDF x-dimension coordinate variable name for the destination grid (i.e., remapped) varible.

ydim_name

The netCDF y-dimension coordinate variable name for the destination grid (i.e., remapped) varible.

zdim_name

The netCDF z-dimension coordinate variable name for the destination grid (i.e., remapped) varible; this is only required for 3-dimensional variables.

Example configuration files for both MOM6 and CICE remappings can be found here.

Launching the Application#

The remapper application may be launched within the UFS Tools package for the MOM6 example described in the previous section as follows.

user@host:$ cd scripts
user@host:$ ./compute_remapper.py --yaml /ufs_tools/parm/remapper/remapper.yaml

If successful the netCDF-formatted file path defined by output_netcdf will be created. The following are selected images illustrating the remapping of MOM6 and CICE \(0.25^{\circ}\) (left) resolution analyses to \(1.0^{\circ}\) (center) and \(5.0^{\circ}\) (right) tripolar grid projections.

MOM6 Sea-Surface Temperature#
_images/MOM6.sst.0p25.png _images/MOM6.sst.1p0.png _images/MOM6.sst.5p0.png
CICE Northern Hemisphere Ice Fraction#
_images/CICE.icefrac.0p25.npstere.png _images/CICE.icefrac.1p0.npstere.png _images/CICE.icefrac.5p0.npstere.png
CICE Southern Hemisphere Ice Fraction#
_images/CICE.icefrac.0p25.spstere.png _images/CICE.icefrac.1p0.spstere.png _images/CICE.icefrac.5p0.spstere.png