...
MBTR FV->FV error - missing source points. It appears MBTR maps have a similar loss of precision as TR maps, but this could be corrupted du to the missing source points.
MBTR SE->SE: montone weights are bad, either NaN or >1
ESMF “aave” and “bilin” maps are good up to ne512pg2, but both produce too large weights with ne1024pg2 (both SE->SE and FV->SE)
AnalyzeMap
request: add min/max of map file frac_a, frac_b data to output
request: check for NaNs when printing weight min/max
for non global (MPAS) source grid, prints tens of thousands of useless consistency warnings
ncremap option is missing global attributes needed by CIME’s domain generation code. they can be added via:
ncatted -a grid_file_src,global,a,c,"$target" $mapname
ncatted -a grid_file_dst,global,a,c,"$atm_scrip" $mapname
Command line options used:
Produce pg2 SE template files:
Code Block |
---|
#atm grid
set NE=1024
set atmgridnp4 = TEMPEST_NE${NE}.g
set atmgrid = TEMPEST_NE${NE}pg2.g # used by TR, MBTR
set atm_scrip = TEMPEST_NE${NE}pg2.scrip.nc # used by ESMF, ncremap
GenerateCSMesh --alt --res $NE --out_format Netcdf4 --file $atmgridnp4
GenerateVolumetricMesh --in $atmgridnp4 --out $atmgrid --np 2 --uniform
ConvertExodusToSCRIP --in $atmgrid --out $atm_scrip
# ocean grid:
set target = ocean.oRRS18to6v3.scrip.181106.nc |
Produce conservative monotone mapping:
Code Block |
---|
# NCREMAP:
set mapname = map_oRRS18to6v3_to_ne${NE}pg2_nco.nc
ncremap -6 --thr_nbr=4 --grd_src=$target --grd_dst=$atm_scrip --map=$mapname
ncatted -a grid_file_src,global,a,c,"$target" $mapname
ncatted -a grid_file_dst,global,a,c,"$atm_scrip" $mapname
#ESMF:
set mapname = map_oRRS18to6v3_to_ne${NE}pg2_aave.nc
srun -N 1 -c 4 -C haswell ESMF_RegridWeightGen -s $target -d $atm_scrip -i --src_regional \
--method conserve --netcdf4 -w $mapname
#TR:
set mapname = map_oRRS18to6v3_to_ne${NE}pg2_mono.nc
GenerateOverlapMesh --b $atmgrid --a $target --out overlap.g
GenerateOfflineMap --out_mesh $atmgrid --in_mesh $target --ov_mesh overlap.g \
--in_type fv --in_np 1 --out_type fv --out_double --out_format Netcdf4 \
--correct_areas --out_map $mapname
|