NCAR / omwg_dev

Repo to centralize the OMWG development, including workflows and diagnostics
0 stars 0 forks source link

g.e23_a16g.GJRAv4.TL319_t232_hycom1_N75.2024.050 #43

Open gustavo-marques opened 7 months ago

gustavo-marques commented 7 months ago

Purpose:

Test the effect of bathymetric modifications near the Surigao Strait on the temperature bias in the Zulu Sea (see #41).


Description:

This is a twin of 2024.025, but with topo edits (make strait shallower) near the Surigao Strait (see below).


Case directory:


Sandbox:

MOM6 tag: dev/ncar_240311

MOM_interface tag: mi_240311

CICE6 tag: cice_20230707

CESM_CICE tag: cesm_cice6_4_1_10


Diagnostics:

--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

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 !!!!!!!!!!!
!!!!! 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 = "FNC1:400,31000.0,0.1,.01"
COORD_CONFIG = "ALE"
!!! interp/remap scheme !!!!!!
INTERPOLATION_SCHEME = "PPM_CW"
REMAPPING_SCHEME = "PPM_CW"
HYCOM1_ONLY_IMPROVES = True
VELOCITY_REMAPPING_SCHEME = "PPM_CW"
!!! regrid filter !!!!
REGRID_TIME_SCALE = 14400.0
REGRID_FILTER_SHALLOW_DEPTH = 100.0
REGRID_FILTER_DEEP_DEPTH = 125.0
!!! topo and channel edits !!!
TOPO_EDITS_FILE = "TopoEdit_tx2_3v2_050.nc"
CHANNEL_CONFIG = "list"
CHANNEL_LIST_FILE = "Channels_tx2_3v2_025"
FATAL_UNUSED_CHANNEL_WIDTHS = True

Topo edits (TopoEdit_tx2_3v2_049.nc):

netcdf TopoEdit_tx2_3v2_049 {
  dimensions:
    nEdits = UNLIMITED ; // (7 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 04/25/2024" ;

  data:
    iEdit = 78, 79, 79, 79, 473, 469, 470 ;

    jEdit = 294, 292, 293, 294, 355, 353, 353 ;

    ni = 540 ;

    nj = 480 ;

    zEdit = 269, 68, 82, 61, 10, 97.6712265014648, 49.4505958557129 ;

}

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