Skip to content

b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274 #239

@gustavo-marques

Description

@gustavo-marques

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:


Sandbox:


Diagnostics:


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:

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions