g.e23_a16g.GJRAv4.TL319_t232_zstar_N75.2024.030

[gustavo-marques]

Purpose:

Test new 75 layer z-star grid (see #31)

---

Description:

This is a twin of 2024.024,, but with a new 75 layer grid.

---

Case directory:

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

---

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_zstar_N65.2024.030/

---

Contacts:
@gustavo-marques, @klindsay28, @fobryan3

---

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

! new zstar with 75 layers
NK = 75
ALE_COORDINATE_CONFIG = "FILE:zstar_75layer_2.5m_248.4m-2024-03-29.nc,dz"
! background visc/diff
KD = 0.
KV = 0.
KD_MIN = 2.0E-07
! interp/remap scheme
INTERPOLATION_SCHEME = "PPM_CW"
REMAPPING_SCHEME = "PPM_CW"
VELOCITY_REMAPPING_SCHEME = "PPM_CW"
!!! topo and channel edits !!!
TOPO_EDITS_FILE = "TopoEdit_tx2_3v2_024.nc"
CHANNEL_CONFIG = "list"
CHANNEL_LIST_FILE = "Channels_tx2_3v2_024"
FATAL_UNUSED_CHANNEL_WIDTHS = True
! others
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"

Topo edits (TopoEdit_tx2_3v2_024.nc):

netcdf TopoEdit_tx2_3v2_024 {
  dimensions:
    nEdits = UNLIMITED ; // (3 currently)

  variables:
    int iEdit(nEdits) ;
      iEdit:long_name = "i-index of edited depth" ;

    int jEdit(nEdits) ;
      jEdit:long_name = "j-index of edited depth" ;

    int ni ;

    int nj ;

    double zEdit(nEdits) ;
      zEdit:_FillValue = NaN ;
      zEdit:long_name = "New value of depth" ;
      zEdit:units = "meters" ;

  // global attributes:
    :title = "Topography Edits" ;
    :original\ grid = "topo.sub150.tx2_3v2.srtm.edit4.SmL1.0_C1.0.nc" ;
    :Edit\ History = "Date Created 03/26/2024" ;


  data:
    iEdit = 473, 469, 470 ;

    jEdit = 355, 353, 353 ;

    ni = 540 ;

    nj = 480 ;

    zEdit = 10, 97.6712265014648, 49.4505958557129 ;

} // group /

Channels_tx2_3v2_024:

U_width,    -5.75,    -5.50,    35.60,    36.30,    12000.0 ! St. of Gibralter
V_width,    28.75,    29.50,    41.10,    41.30,     5000.0 ! Bosphorus St.
U_width,    26.50,    26.80,    40.00,    40.70,     5000.0 ! Dardanelles