-
Notifications
You must be signed in to change notification settings - Fork 3
Description
Purpose:
Check whether the new ocean baseline settings (268), starting from a corresponding G case, lead to Lab Sea freezing and potential recovery.
Description:
Same as 268 (#232) except:
Ocean and sea ice initial conditions from restarts from year 9 of a G case (G154). Here are the Lab Sea temperature and salinity drifts in this G case.
- add to user_nl_cice
ice_ic = "/glade/derecho/scratch/gmarques/archive/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154/rest/0009-01-01-00000/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154.cice.r.0009-01-01-00000.nc"
tr_fsd = .true.
restart_fsd = .true.
nfsd = 12
wave_spec_type = 'alt'
histfreq = "m", "x", "x", "x", "x"
histfreq_n = 1, 1, 0, 0, 0
f_fresh = 'mxxxx'
f_fresh_ai = 'mxxxx'
f_frzmlt = 'mdxxx'
f_afsd = "mdxxx"
f_afsdn = "mdxxx"
f_dafsd_latg = "mdxxx"
f_dafsd_latm = "mdxxx"
f_dafsd_newi = "mdxxx"
f_dafsd_wave = "mdxxx"
f_dafsd_weld = "mdxxx"
f_wave_sig_ht = "mdxxx"
f_diam_ww = "mdxxx"
f_hice_ww = "mdxxx"
f_fsdrad = "mdxxx"
f_fsdperim = "mdxxx"
f_uvel = "mdxxx"
f_vvel = "mdxxx"
r_snw = 1.0
- add to user_nl_mom
INPUTDIR = "/glade/derecho/scratch/gmarques/g.e30_a07g.GW1850MARBL_JRA.TL319_t232_wgx3_hycom1_N75.2025.152/run/INPUT"
MEKE_GEOMETRIC_ALPHA = 0.09
SMAGORINSKY_KH = False
MAXTRUNC = 1.0E6
! Leith+E Settings
TAPER_LEITHY = True
LEITH_BI_CONST = 85.0
LEITHY_CK = 1.0
! Background KH
USE_KH_BG_2D = True
KH_BG_2D_FILENAME = "KH_BG_2D.nc"
KH_BG_2D_VARNAME = "KH"
! Resolution scaling:
!! GM & Redi the same
!! Don't apply scaling to Laplacian viscosity, which is the 2D background
RESOLN_SCALED_KH = False
RESOLN_SCALED_KHTH = True ! Already using, here for clarity
RESOLN_SCALED_KHTR = True ! Already using, here for clarity
KH_RES_SCALE_COEF = 0.5 ! Tighter around the equator
KH_RES_FN_POWER = 6 ! Already using, here for clarity
VISC_RES_SCALE_COEF = 0.5 ! Used for KHTR
VISC_RES_FN_POWER = 6 ! Used for KHTR
INTERPOLATE_RES_FN = True ! Only computes resolution fn once, then interpolates
THICKNESS_CONFIG = thickness_file
VELOCITY_CONFIG = file
TS_CONFIG = file
INIT_LAYERS_FROM_Z_FILE = False
THICKNESS_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154/rest/0009-01-01-00000/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154.mom6.r.0009-01-01-00000.nc"
TS_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154/rest/0009-01-01-00000/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154.mom6.r.0009-01-01-00000.nc"
VELOCITY_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154/rest/0009-01-01-00000/g.e30_a07g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2025.154.mom6.r.0009-01-01-00000.nc"
TEMP_IC_VAR = "Temp"
SALT_IC_VAR = "Salt"
Case directory:
-
Locally (if still available):
/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274 -
On github:
https://github.com/NCAR/cesm_dev/tree/b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274
Sandbox:
-
Locally (if still available):
/glade/work/gmarques/cesm.sandboxes/cesm3_0_alpha07g -
On github: https://github.com/ESCOMP/cesm/tree/cesm3_0_alpha07g
Diagnostics:
- Diags (if available):
https://webext.cgd.ucar.edu/BLT1850/b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274/
Output:
- Raw output on derecho (if still available):
/glade/derecho/scratch/gmarques/archive/b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274/ - CAM Climos:
/glade/campaign/cgd/amp/amwg/climo/b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274
Contacts:
@cecilehannay @dlawrenncar @adamrher @PeterHjortLauritzen @JulioTBacmeister @swrneale @slevis-lmwg @wwieder @olyson @gustavo-marques @iangrooms @fobryan3 @klindsay28 @dabail10 @duvivier @tilmes @mbramberger @dmleung
Extra details: