b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366

[cecilehannay]

Purpose:
Increase the taper depth to determine whether we can prevent Labrador Sea freezing while avoiding the reversal of the AMOC circulation within the Nordic Seas.

Same configuration as 360a, but with the taper depth increased to 2000 m.

Start from 360a at year 139.

---

Description:

Same as 360a but

• taper depth 1500. -> 2000.

MLE_TAPER_FN_DEPTH = 1500. -> 2000. 

• Start from 360a at year 139.

./xmlchange RUN_REFCASE=bg.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.360a
./xmlchange RUN_REFDATE=0139-01-01
./xmlchange RUN_TYPE=hybrid

**Compset
Compset: 1850C_CAM70%MT_CLM60%BGC-CROP_CICE_MOM6_MOSART_CISM2%GRIS-EVOLVE_WW3
res: ne30pg3_t233_wg37_gris4

---

Case directory:

• Locally (if still available):
  /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366

• On github:
  https://github.com/NCAR/cesm_dev/tree/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366

---

Sandbox:

• Locally (if still available): /glade/campaign/cesm/development/liwg/leguy/CESM3/code/cesm3_0_alpha09b/cime/scripts

• On github: https://github.com/ESCOMP/cesm/tree/cesm3_0_alpha09b

• add this block to ccs_config/modelgrid_aliases_nuopc.xml (next to the
  other ne30pg3_t233 aliases):

  <model_grid alias="ne30pg3_t233_wg37_gris4">
    <grid name="atm">ne30np4.pg3</grid>
    <grid name="lnd">ne30np4.pg3</grid>
    <grid name="ocnice">tx2_3v3</grid>
    <grid name="glc">gris4</grid>
    <grid name="wav">wgx3v7</grid>
    <mask>tx2_3v3</mask>
  </model_grid>

---

Diagnostics:

• Diags (if available):
  https://webext.cgd.ucar.edu/BMT1850/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366/

---

Output:

• Raw output on derecho (if still available): /glade/derecho/scratch/$USER/archive/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366/
• CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366

---

Contacts:
@cecilehannay @dlawrenncar @adamrher @PeterHjortLauritzen @JulioTBacmeister @swrneale @slevis-lmwg @wwieder @olyson @gustavo-marques @iangrooms @fobryan3 @klindsay28 @dabail10 @duvivier @tilmes @mbramberger @gunterl

---

Extra details:

• create_newcase

./create_newcase --compset 1850C_CAM70%MT_CLM60%BGC-CROP_CICE_MOM6_MOSART_CISM2%GRIS-EVOLVE_WW3 --res ne30pg3_t233_wg37_gris4 \
--case /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366  \
--run-unsupported  

• pe_layout

/bin/cp /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/bg.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.365/env_mach_pes.xml .

xmlchange

./xmlchange CAM_CONFIG_OPTS=' -pcols 9' --append
./xmlchange CASE_GIT_REPOSITORY=git@github.com:NCAR/cesm_dev.git
./xmlchange RUN_REFCASE=bg.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.360a
./xmlchange RUN_REFDATE=0139-01-01
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=true
./xmlchange RUN_REFDIR=cesm2_init
./xmlchange OCN_DIAG_MODE=spinup
./xmlchange JOB_PRIORITY=special --force
./xmlchange PROJECT=CESM0023,RESUBMIT=15,STOP_N=4,STOP_OPTION=nyears
./xmlchange PROJECT=CESM0023
./xmlchange JOB_WALLCLOCK_TIME=12:00:00  --subgroup case.run

• user_nl_cam

mfilt               =       0,      0,     20,      40,      12,      120,      1,   1
nhtfrq              =       0,      0,    -24,      -3,       0,       -2,      0,  -8760
ndens               =       2,       2,      2,       2,       2,        1,      2,   1
interpolate_output  =  .true.,  .false., .true.,  .true., .false.,   .true.,     .true.
interpolate_nlat    =     192,     192,    192,     192,      192,     192,      192
interpolate_nlon    =     288,     288,    288,     288,      288,     288,      288

empty_htapes = .true.

fincl1 = 'ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODDUSTdn','AODVIS', 'AODVISdn','BURDENBC',
'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT',
'BURDENSO4', 'BURDENSOA', 'CAPE', 'CDNUMC',  'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW',
'CLDMED', 'CLDTOT', 'CLOUD', 
'FCTI', 'FCTL', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT',
'FREQZM', 'FSDS', 'FSDSC', 'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LANDFRAC', 'LHFLX', 'LWCF',
'OCNFRAC', 'OMEGA', 'OMEGA500', 'PBLH', 'PHIS', 'PINT', 'PMID', 'PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL',  'Q',
'QFLX', 'RELHUM', 'SHFLX', 'SOLIN', 'SST',
'SWCF',
'T', 'TAUX', 'TAUY', 'TFIX', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'U', 'U10',
'UBOT', 
'V', 'VBOT', 'Z3', 
'AODVISdn', 'CDNUMC',  'NUMICE', 'NUMLIQ','OMEGA500',
 'AQSO4_H2O2','AQSO4_O3', 
'ncl_a1SFWET',
'ncl_a2SFWET', 'ncl_a3SFWET', 'pom_a1SFWET', 'pom_a4SFWET', 'so4_a1SFWET', 'so4_a2SFWET', 'so4_a3SFWET', 'soa_a1SFWET',
'soa_a2SFWET', 'bc_c1SFWET', 'bc_c4SFWET', 'dst_c1SFWET', 'dst_c2SFWET', 'dst_c3SFWET', 'ncl_c1SFWET', 'ncl_c2SFWET',
'ncl_c3SFWET', 'pom_c1SFWET', 'pom_c4SFWET', 'so4_c1SFWET', 'so4_c2SFWET', 'so4_c3SFWET', 'soa_c1SFWET', 'soa_c2SFWET',
'bc_a1DDF', 'bc_a4DDF', 'dst_a1DDF', 'dst_a2DDF', 'dst_a3DDF', 'ncl_a1DDF', 'ncl_a2DDF', 'ncl_a3DDF',
'pom_a1DDF', 'pom_a4DDF', 'so4_a1DDF', 'so4_a2DDF', 'so4_a3DDF', 'soa_a1DDF', 'soa_a2DDF',
'so4_a1_CLXF', 'so4_a2_CLXF', 'SFbc_a4', 'SFpom_a4', 'SFso4_a1', 'SFso4_a2',
'so4_a1_sfgaex1', 'so4_a2_sfgaex1', 'so4_a3_sfgaex1', 'soa_a1_sfgaex1', 'soa_a2_sfgaex1',
'SFdst_a1','SFdst_a2', 'SFdst_a3', 'SFncl_a1', 'SFncl_a2', 'SFncl_a3',
'num_a2_sfnnuc1', 'SFSO2', 'OCN_FLUX_DMS', 'SAD_SULFC', 'SAD_TROP', 'SAD_AERO',
'AODVISstdn', 'REFF_AERO', 'so4_a1', 'so4_a2', 'SAD_AERO', 'PS', 'TROP_P', 'FSNTOAC',
'AODSSdn', 'AODSSdn', 'BURDENSO4dn', 'AODVISdn', 'AODVISstdn',
'AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn', 'AODDUSTdn','AODVISdn', 'AODVISstdn','BURDENSO4dn', 'SFDMS','BURDENBCdn', 'BURDENDUSTdn', 
'BURDENPOMdn', 'BURDENSEASALTdn', 'BURDENSO4dn', 'BURDENSOAdn', 'AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn', 'AODDUSTdn','AODVISdn', 'AODVISstdn','BURDENSO4dn', 'SFDMS','BURDENBCdn', 'BURDENDUSTdn', 'BURDENPOMdn', 'BURDENSEASALTdn', 'BURDENSO4dn', 'BURDENSOAdn', 
'CCN3', 'CFC11', 'CFC12', 'CH4','CO2', 'EXTINCTdn', 'N2O', 
'EFIX'

fincl3 = 'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL', 'Z500'

fincl4 =  'PRECC','PRECL'


seasalt_emis_scale		= 1.6D0
tracer_cnst_datapath           = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/TracerCnstCyclical1850'

tracer_cnst_file               = 'tracer_cnst_MTt4s_halons_3Dmonthly_1850climo330_c20260423.nc'

prescribed_ozone_datapath      = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/ozone_strataero_day_5_zm'

prescribed_ozone_file          = 'ozone_strataero_b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix_1850_5_30_climo_zm_5day_c260424.nc'

prescribed_strataero_datapath  = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/ozone_strataero_day_5_zm'

 prescribed_strataero_file              = 'ozone_strataero_b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix_1850_5_30_climo_zm_5day_c260424.nc'

stream_ndep_data_filename              = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/NDepositionCyclical1850/fndep_clm_MTt4s_3Dmonthly_1850climo330_c20260424.nc'

clubb_gamma_coef                =  0.298     
clubb_gamma_coefb               =  0.298   

clubb_penta_solve_method = 2
clubb_tridiag_solve_method = 2

• user_nl_clm

use_init_interp = .true.
reseed_dead_plants = .true.
use_excess_ice_streams=.false. 

• user_nl_cpl

histaux_l2x1yrg = .true.
histaux_atm2med_file1_enabled = .true.
histaux_atm2med_file2_enabled = .true.
histaux_atm2med_file3_enabled = .true.
histaux_atm2med_file4_enabled = .true.
histaux_atm2med_file5_enabled = .true.
histaux_atm2med_file5_history_n = 1
histaux_atm2med_file5_history_option = 'ndays'
histaux_atm2med_file5_ntperfile = 1
histaux_rof2med_file1_enabled = .true.
histaux_rof2med_file1_ntperfile = 8

• user_nl_cice

r_snw = 0.75
restart_fsd = .true.
histfreq = "m", "x", "x", "x", "x"

• user_nl_cism

dew = 4000.
dns = 4000.
ewn = 421
nsn = 721
upn = 11

idiag = 136
jdiag = 294

dycore = 2
flow_law = 2
evolution = 3
temperature = 1
temp_init = 3
marine_margin = 1
calving_domain = 1
calving_init = 1
limit_marine_cliffs = .false.
basal_mass_balance = 1

gthf = 1
smb_input = 0
smb_input_function = 0
artm_input_function = 0
bmlt_float = 0
dm_dt_diag = 1

compute_blocks = 0
which_ho_powerlaw_c = 2
which_ho_babc = 9
which_ho_bwat = 0
which_ho_effecpress = 0
which_ho_efvs = 2
which_ho_thermal_timestep = 0
which_ho_resid = 4
which_ho_sparse = 3
which_ho_nonlinear = 1
which_ho_gradient = 0
which_ho_gradient_margin = 1
which_ho_assemble_beta = 1
which_ho_assemble_bfric = 1
which_ho_assemble_taud = 1
which_ho_assemble_lateral = 1
which_ho_approx = 4
diva_slope_correction = .true.
which_ho_precond = 3
which_ho_calving_front = 0
which_ho_ground = 2
which_ho_flotation_function = 2
glissade_maxiter = 100
block_inception = .false.
remove_ice_caps = .false.

flow_factor_ground = 1.
flow_factor_float = 1.
ice_limit = 1.
thck_gradient_ramp = 100.
marine_limit = -200
max_slope = 0.5
rhoi = 917.
rhoo = 1026.
grav = 9.80616
shci = 2117.27
lhci = 3.337e5
trpt = 273.16
pmp_offset = 10.
beta_grounded_min = 100.
beta_powerlaw_umax = 5000.
pseudo_plastic_q = 0.25
pseudo_plastic_u0 = 100.
p_ocean_penetration = 0.5
cliff_timescale = 0.0
taumax_cliff = 1.e6
powerlaw_m = 3.0
powerlaw_c_const = 20000.
powerlaw_c_max = 100000.
powerlaw_c_min = 100.

• user_nl_mom

DIAG_COORD_DEF_RHO2 = "FILE:ocean_rho2_190917.nc,interfaces=rho2"
MLE_TAPER_FN_DEPTH = 2000.
MLE_TAPER_FN_POWER = 4
BBL_EFFIC = 0.2 
EOF

-SourceMods

SourceMods/src.cam/ 
dimensions_mod.F90  
dyn_comp.F90  
fvm_mapping.F90  
global_norms_mod.F90  
prim_advance_mod.F90  
prim_driver_mod.F90  README  
viscosity_mod.F90

SourceMods/src.mom/ 
diag_table.unresolved  
MOM_mixed_layer_restrat.F90

[cecilehannay]

Spaghetti plot

MOC reversal in the Nordic Seas

See discussion: #370

[gustavo-marques]

I updated the discussion after the OMWG meeting. We think AMOC is okay in 366.