b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274

Open
#239 4 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Assessment

Difficulty
5/5
Estimated time
Over a week
Newbie friendliness
25/100
Issue type
Feature
Clarity
Mostly clear
Activity status
Stale
Domain
hpc

Research direction

Start with the referenced case directory and the cesm3_0_alpha07g sandbox; compare the setup with issue #232 and inspect the user_nl_cice and user_nl_mom entries. Run or locate diagnostics for case b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.274 and determine whether the Lab Sea freezes and later recovers from the G154 restart; completion is a documented result.

Written by the indexing model from the issue text.

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:


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:

Dominant language
No language data
Stars
11
Forks
3
PR merge metrics
No merged PRs in 30d

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

More from NCAR/cesm_dev

All issues in NCAR/cesm_dev

Similar issues

More HPC issues

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.