g.e23_a16g.GJRAv4.TL319_t232_hycom1_N75.2024.046

[gustavo-marques]

Purpose:
Investigate the impact of limiting the layer thicknesses in the HYCOM1 grid.

---

Description:
Twin of 2024.018, but using KV = KD = 0, topo and channel edits from 025, and MAX_LAYER_THICKNESS_CONFIG that eliminates the initial layer thicknesses that are >= than the 90'th quantile of the initial layer thicknesses in 018.

---

Case directory:

• Locally (if still available):
  /glade/work/gmarques/cesm.cases/G/g.e23_a16g.GJRAv4.TL319_t232_hycom1_N75.2024.046

---

Sandbox:

• Locally (if still available): /glade/work/gmarques/cesm.sandboxes/cesm2_3_alpha16g

• On github: https://github.com/ESCOMP/cesm/tree/cesm2_3_alpha16g

MOM6 tag: dev/ncar_240311

MOM_interface tag: mi_240311

CICE6 tag: cice_20230707

CESM_CICE tag: cesm_cice6_4_1_10

---

Diagnostics:

• Diags (if available):
  TODO

---

Output:

• Raw output on derecho (if still available): /glade/derecho/scratch/gmarques/archive/g.e23_a16g.GJRAv4.TL319_t232_hycom1_N75.2024.046/

---

Contacts:
@gustavo-marques

---

Extra details:

--compset 2000_DATM%JRA-1p4-2018_SLND_CICE_MOM6_DROF%JRA-1p4-2018_SGLC_SWAV_SESP --res TL319_t232

The following included in user_nl_mom

MAXTRUNC = 100
INPUTDIR = "INPUT"
STOCH_EOS=True
STANLEY_COEFF=0.5
USE_STANLEY_PGF=True
USE_STANLEY_ML=True
USE_STANLEY_ISO=True
USE_STANLEY_GM=True
USE_NEUTRAL_DIFFUSION = True
NDIFF_TAPERING = True
KHTR_USE_EBT_STRUCT = True
MEKE_KHTR_FAC = 1.0
DO_GEOTHERMAL = True
GEOTHERMAL_SCALE = 1.0
GEOTHERMAL_FILE = "geothermal_davies2013_tx2_3_20231120_cdf5.nc"
GEOTHERMAL_VARNAME = "geothermal_hf"
!!!!! Franks's grid !!!!!!!!!!!
NK = 75
REGRIDDING_COORDINATE_MODE = "HYCOM1"
BOUNDARY_EXTRAPOLATION = True
ALE_COORDINATE_CONFIG = "HYBRID:hybrid_75layer_zstar2.50m-2020-11-23.nc,sigma2,dz"
REGRID_COMPRESSIBILITY_FRACTION = 0.01
MAXIMUM_INT_DEPTH_CONFIG = "FNC1:5,8000.0,1.0,.01"
MAX_LAYER_THICKNESS_CONFIG = "FILE:dz_max_90th_quantile.nc,dz"
COORD_CONFIG = "ALE"
!!! interp/remap scheme !!!!!!
INTERPOLATION_SCHEME = "PPM_CW"
REMAPPING_SCHEME = "PPM_CW"
HYCOM1_ONLY_IMPROVES = True
VELOCITY_REMAPPING_SCHEME = "PPM_CW"
!!! topo and channel edits !!!
TOPO_EDITS_FILE = "TopoEdit_tx2_3v2_025.nc"
CHANNEL_CONFIG = "list"
CHANNEL_LIST_FILE = "Channels_tx2_3v2_025"
FATAL_UNUSED_CHANNEL_WIDTHS = True
! background visc/diff
KD = 0.
KV = 0.
KD_MIN = 2.0E-07

File dz_max_90th_quantile.nc:

netcdf dz_max_90th_quantile {
  dimensions:
    z = 75 ;

  variables:
    double dz(z) ;

  // global attributes:
    :title = "Maximum thickness using h > 90th quantile of run 018" ;
    :author = "Gustavo Marques ([email protected])" ;
    :Created = "2024-04-17" ;


  data:
    dz = 2.5, 2.5, 2.5, 2.5, 4.65954471, 4.91196546, 5.19412109, 5.50946059, 5.86181952, 6.25546, 6.69511391, 7.18602944, 7.73402094, 8.34552182, 9.02764048, 9.78821859, 10.63589128, 11.58014822, 12.63139452, 13.80100975, 15.10140314, 16.54606232, 18.14959251, 19.92774221, 21.89741083, 24.07663286, 26.4845323, 29.14124041, 32.06776905, 35.28583156, 38.81760291, 42.68541118, 46.91135357, 51.51683185, 56.52200507, 61.94516131, 67.80201558, 74.1049473, 80.86219824, 88.07705989, 95.74708725, 103.86338297, 112.41000134, 121.36352345, 130.69285331, 140.35927768, 150.31682062, 160.51290687, 170.8893278, 181.38348107, 191.92983316, 202.46153511, 212.91210815, 223.21710955, 233.31569094, 243.15197092, 252.67615996, 261.84539616, 270.62427305, 278.98506237, 286.90765387, 294.37924899, 301.39385482, 307.95162933, 314.05812899, 319.72350586, 324.96169559, 329.78962955, 334.22649664, 338.29307219, 342.01112468, 345.40290465, 348.49071585, 351.29656501, 353.84188435 ;

}