...
Note3: for TR FV to FV maps, the --mono options are ignored. monotone conservative maps are constructed using --in_np 1 (piecewise constant reconstruction). higher order maps (--in_np 2) are not monotone.
ncremap build built in FV to FV map file generator
oRRS18to6 → ne256pg2 status: GOOD
~3min with 4 threads on Cori frontend node
Weights are good. row sum and conservation error: 6e-13
oRRS18to6 → ne1024pg2
~40 min on Compy with 4 threads. 1.5 h on Cori-Haswell node
Weights are good. row sum error 2e-12, conservation error 1e-12
TR “mono” . (TR from TR from with new --correct_areas option)
oRRS18to6 → ne256pg2 status: GOOD
2.5h on Anvil (single processor)
Weights are good. Row sums error: 7e-16. Conservation error 7e-16
vortex test field l2/max error: 3.6145e-05 4.3881e-04
oRRS18to6 → ne512pg2: status GOOD
~9h to compute overlap map
weights: good. row sum error: 7e-16. Conservation error 7e-16
vortex: l2, linf: 1.7867e-04 2.2985e-03
oRRS18to6 → ne1024pg2: status: GOOD
15.5h on Anvil, using 20GB
weights: good. row sum error: 4e-16. conservation error: 1e15
MBTR “mono”. Status: NOT USABLE
oRRS18to6 → ne256pg2
Weights are good, but 316 MPAS grid points have weight=0. This is incorrect and corrupts other diagnostics.
ESMF “aave”
oRRS18to6 → ne256pg2. Status: GOOD
weights are good. row sum max 8.5e-12. Conservation error 1.5e-11
vortex test field l2/max error: 3.6144e-05 4.3881e-04
oRRS18to6 → ne512pg2. Status: GOOD
weights are good. row sum error 3e-11. conservation error 1e-11
vortex: l2, linf: 1.7866e-04 2.2985e-03
oRRS18to6 → ne1024pg2. Status: NOT USABLE
35min on Cori Haswell node
( /global/project/projectdirs/acme/taylorm/mapping/map_oRRS18to6v3_to_ne1024pg2_aave.nc.BAD )
weights are bad: max 2.6, row sum max 3.0. conservation error 4e-4
same issue with ESMF v8.
...
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 |
To check maps, use