Scientifically validate high uniform and high resolution variable resolution MPAS topo files

[PeterHjortLauritzen]

I have been generating several MPAS topo files that need to be scientifically validated:

1. currently doing rrfac manipulation. Should we? (subtract uniform resolution topo files with resolution of low res and high res areas; sanity check)
2. when MPAS reaches ~3km in refinement region; is the topo "well-behaved"? Would need to do higher resolution intermediate cubed-sphere resolution and check sensitivity to that)

Command used for 60-3km grid:

./cube_to_target --grid_descriptor_file='/glade/work/skamaroc/MPAS_meshes/CESM_MPAS/60-3km/West_Pacific/mpas_x20.835586.wpac.1_esmf.nc' --intermediate_cs_name='/glade/p/cgd/amp/pel/topo/cubedata/gmted2010_modis_bedmahcine-ncube3000-220518.nc' --output_grid='mpas_x20.835586.wpac' -u 'Peter Hjort Lauritzen, pel@ucar' -q 'output/' --smoothing_scale=50.0 -r -y 20 -v -m

-m = smoothing over ocean points
-v = rrfac manipulation
-y = refinement factor

Command used for 60-10km grid:

./cube_to_target --grid_descriptor_file='/glade/work/skamaroc/MPAS_meshes/CESM_MPAS/60-10km/mpas_x6.999426.wpac.1_esmf.nc' --intermediate_cs_name='/glade/p/cgd/amp/pel/topo/cubedata/gmted2010_modis_bedmahcine-ncube3000-220518.nc' --output_grid='mpas_x6.999426.wpac' -u 'Peter Hjort Lauritzen, pel@ucar' -q 'output/' --smoothing_scale=50.0 -r -y 6 -v -m

Topo files here:

/glade/p/cgd/amp/pel/topo/grids/mpas-60-3km-WestPacific/mpas_x20.835586.wpac_gmted2010_bedmachine_nc3000_Laplace0050_20220520.nc
/glade/p/cgd/amp/pel/topo/grids/mpas-60-10km-WestPacific/mpas_x6.999426.wpac_gmted2010_bedmachine_nc3000_Laplace0050_20220520.nc

[PeterHjortLauritzen]

For reference some plots:

PHIS, SGH and SGH30 for 60-10km mesh

[PeterHjortLauritzen]

Commands have changed (-r now turns ridge finding off) and we recommend not smoothing over ocean (turn off with -m option) and no manipulation of refinement factor for MPAS grids (turn on with -v option):

Recommended command for 60-3km grid:

./cube_to_target --grid_descriptor_file='/glade/work/skamaroc/MPAS_meshes/CESM_MPAS/60-3km/West_Pacific/mpas_x20.835586.wpac.1_esmf.nc' --intermediate_cs_name='/glade/p/cgd/amp/pel/topo/cubedata/gmted2010_modis_bedmachine-ncube3000-220518.nc' --output_grid='mpas_x20.835586.wpac' -u 'Peter Hjort Lauritzen, pel@ucar' -q 'output/' --smoothing_scale=50.0 -y 20

Command used for 60-10km grid:

./cube_to_target --grid_descriptor_file='/glade/work/skamaroc/MPAS_meshes/CESM_MPAS/60-10km/mpas_x6.999426.wpac.1_esmf.nc' --intermediate_cs_name='/glade/p/cgd/amp/pel/topo/cubedata/gmted2010_modis_bedmachine-ncube3000-220518.nc' --output_grid='mpas_x6.999426.wpac' -u 'Peter Hjort Lauritzen, pel@ucar' -q 'output/' --smoothing_scale=50.0 -y 6