...
Command line options used:
Produce pg2 SE template files: (need utilities from TR New versions required:
TR: code from or later
...
NCO 4.9.2-alpha05 or newer
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 # option "--alt" now the default and has to be removed from newer TR versions: GenerateCSMesh --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 |
...
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 “ncks --chk_map”:
weights: verify that all weights are between [0,1], none of them are NaNs.
row sums: The row sum consistency check measures the ability of the map to preserve the constant function.
ocean->atm maps: max value should be 1.0. min value will be between [0..1] at land and fractional land points. Error reported above is from the max value.
atm->ocean maps: min/max value should be 1.0 for all rows with nonzero weights (since atmosphere grid is global)
area weighted column sums: this measures if the map is conservative.
ocean->atm maps: area_b weighted column sum should equal area_a for all columns.
atm->ocean: area_b weighted column sum should equal area_a at all ocean points, and be between [0,1] at fractional (coastline) points.