search

NCAR / Topo

NCAR Global Model Topography Generation Software for Unstructured Grids
26 stars 16 forks source link

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

Open PeterHjortLauritzen opened 2 years ago

PeterHjortLauritzen commented 2 years ago

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

Screen Shot 2022-05-20 at 2 20 26 PM

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 commented 2 years ago

For reference some plots:

PHIS, SGH and SGH30 for 60-10km mesh

Screen Shot 2022-05-20 at 4 14 59 PM

Screen Shot 2022-05-20 at 4 15 20 PM

Screen Shot 2022-05-20 at 4 15 36 PM

PeterHjortLauritzen commented 2 years ago

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