...
./GenerateCSMesh --alt --res 30 --file NE30.g # the standard atmosphere SE grid, using the np4" GLL spectral element
./GenerateVolumetricMesh --in NE30.g --np 2 --out NE30pg2.g # An atmosphere "FV" grid (for the new pg2 option where physics is on a FV grid)
Note that it 1: It is important to use the "–alt" option, so the mesh ordering matches what is used internally by E3SM. This option was made the default in 2020/2.
Note 2: For SE atmosphere grids, the Exodus file is what is used by the model, and it contains the corners of all the spectral elements. The dual grid (SCRIP format) metadata files are not needed by TempestRemap.
...
Note: Some of the maps used for state variables in ACME v1 are labeled "_conserve.nc", which are conservative, high-order, non-monotone maps ("highorder" in the above proposed naming convention), which causes some confusion because "conserve" implies it's a replacement for an "aave" type map, when really that type of map is a replacement for "bilin". Thus I'm editing the naming convention (2017/11) to not use "conserve".
The Overlap Mesh
IMPORTANT!!!: MAPS grids often have holes where there are no mesh points. When computing the TempestRemap overlap grid (step 1 below), the grid with less coverage
must be the "a" grid, while the global grid should be the "b" grid. To generate the overlap grid:
...
TR options to generate the various maps:Paul Ullrich: which of these should we add "--correct_areas"?
- SE→FV mono. (conservative, monotone, 1st order)
- ./GenerateOfflineMap --in_mesh $atmgrid --out_mesh $ocngrid --ov_mesh overlap_mesh.g --in_type cgll --in_np 4 --out_type fv --mono --correct_areas --out_map map_atm2ocn_mono.nc
- SE→FV intbilin (nonconservative, monotone, higher order)
- ./GenerateOfflineMap --in_mesh $atmgrid --out_mesh $ocngrid --ov_mesh overlap_mesh.g --in_type cgll --in_np 4 --out_type fv --mono3 --noconserve --out_map map_atm2ocn_intbilin.nc
- Does these maps need "–correct_areas" ?
- SE→FV highorder. (conservative, nonmonotone)
- ./GenerateOfflineMap --in_mesh $atmgrid --out_mesh $ocngrid --ov_mesh overlap_mesh.g --in_type cgll --in_np 4 --out_type fv --correct_areas --out_map map_atm2ocn_highorder.nc
- FV→SE monotr. (conservative, monotone, 1st order)
- The transpose of the SE→FV mono map:
- ./GenerateTransposeMap --in map_atm2ocn_mono.nc --out map_ocn2atm_monotr.nc
- FV→SE highorder. (conservative, non-monotone)
- ./GenerateOfflineMap --in_mesh $ocngrid --out_mesh $atmgrid --ov_mesh overlap_mesh.g --in_type fv --in_np 2 --out_type cgll --out_np 4 --volumetric --out_map map_ocn2atm_highorder.nc
- Does these maps need "–correct_areas" ?
- FV→SE mono. (conservative, monotone, 1st order)
- ./GenerateOfflineMap --in_mesh $ocngrid --out_mesh $atmgrid --ov_mesh overlap_mesh.g --in_type fv --in_np 1 --out_type cgll --out_np 4 --volumetric --out_map map_ocn2atm_mono.nc
- Does these maps need "–correct_areas" ?
- FV→SE intbilintr (nonconservative, monotone, higher order)
- The transpose of the SE→FV intbilin map
- ./GenerateTransposeMap --in map_atm2ocn_intbilin.nc --out map_ocn2atm_intbilintr.nc
- FV→FV mono (conservative, monotone, 1st order)
- The classic cell integrated piecewise constant map. Equivalent to the ESMF "aave" map and NCO's built in map.
- ./GenerateOfflineMap --in_mesh $ocngrid --out_mesh $atmgrid --ov_mesh overlap_mesh.g --in_type fv --in_np 1 --out_type fv --out_np 1 -correct_areas --out_map map_ocn2atm_mono.nc
- FV→FV highorder. (conservative, non-monotone, high order)
- ./GenerateOfflineMap --in_mesh $atmgrid --out_mesh $ocngrid --ov_mesh overlap_mesh.g --in_np 2 --in_type fv --out_type fv --out_type fv -correct_areas --out_map map_atm2ocn_highorder.nc
...