b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374

[cecilehannay]

Purpose:
Same as 373 but with tuning of backscatter

---

Description:
Same as 373 except

• add to user_nl_mom

LEITH_BI_CONST = 80 

• start from 371 at year 53

./xmlchange RUN_REFCASE=b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.371
./xmlchange RUN_REFDATE=0053-01-01
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=true
./xmlchange RUN_REFDIR=cesm2_init

---

Case directory:

• Locally (if still available):
  /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374

• On github:
  https://github.com/NCAR/cesm_dev/tree/b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374

---

Sandbox:

• Locally (if still available): /glade/work/hannay/cesm_tags/cesm3_0_alpha09e_MARBL

• On github:

# get the cesm tag
git clone -b cesm3_0_alpha09e https://github.com/ESCOMP/cesm.git cesm3_0_alpha09e_MARBL

# Make the following change in the [submodule "mom"] section of .gitmodules:
cd cesm3_0_alpha09e_MARBL
emacs -nw .gitmodules
-        fxtag = mi_260625
+        fxtag = mi_260715

# check out the tag
./bin/git-fleximod update

---

Diagnostics:

• Diags (if available):
  https://webext.cgd.ucar.edu/BMT1850/b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374/

---

Output:

• Raw output on derecho (if still available): /glade/derecho/scratch/hannay/archive/b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374/
• CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374

---

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

---

Extra details:

cd /glade/work/hannay/cesm_tags/cesm3_0_alpha09e_MARBL/cime/scripts

./create_newcase --compset  1850C_CAM70%MT_CLM60%BGC-CROP_CICE_MOM6%MARBL-BIO_MOSART_CISM2%GRIS-EVOLVE_WW3_SESP \
--res ne30pg3_t233_wg37_gris4 \
--case /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374  

• copy pe_layout

/bin/cp /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.366/env_mach_pes.xml .
./xmlchange NTASKS_OCN=2560 

• xmlchange

./xmlchange CAM_CONFIG_OPTS=' -pcols 9' --append
./xmlchange CASE_GIT_REPOSITORY=git@github.com:NCAR/cesm_dev.git

./xmlchange RUN_REFCASE=b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.371
./xmlchange RUN_REFDATE=0053-01-01
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=true
./xmlchange RUN_REFDIR=cesm2_init

• 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_datapath = '/glade/derecho/scratch/shawnh/Forcings_tmp/20260706/b.e30_alpha09b.B1850C_MTt4s.ne30_t233_wgx3.360/atm/proc/TracerCnstFileCyclic1850'

!tracer_cnst_file               = 'tracer_cnst_MTt4s_halons_3Dmonthly_1850climo330_c20260423.nc'
tracer_cnst_file = 'tracer_cnst_MTt4s_halons_3Dmonthly_1850climo360_c20260706.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_datapath = '/glade/derecho/scratch/shawnh/Forcings_tmp/20260706/b.e30_alpha09b.B1850C_MTt4s.ne30_t233_wgx3.360/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_ozone_file = 'ozone_strataero_b.e30_alpha09b.B1850C_MTt4s.ne30_t233_wgx3.360_1850_0005_0030_transient_zm_5day_c20260706.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_datapath = '/glade/derecho/scratch/shawnh/Forcings_tmp/20260706/b.e30_alpha09b.B1850C_MTt4s.ne30_t233_wgx3.360/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'
prescribed_strataero_file = 'ozone_strataero_b.e30_alpha09b.B1850C_MTt4s.ne30_t233_wgx3.360_1850_0005_0030_transient_zm_5day_c20260706.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'
stream_ndep_data_filename = '/glade/derecho/scratch/shawnh/Forcings_tmp/20260706/b.e30_alpha09b.B1850C_MTt4s.ne30_t233_wgx3.360/atm/proc/NDepositionCyclic1850/fndep_clm_MTt4s_3Dmonthly_1850climo330_c20260706.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

• 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.
which_ho_precond = 1

• user_nl_mom

DIAG_COORD_DEF_RHO2 = "FILE:ocean_rho2_190917.nc,interfaces=rho2"
BBL_EFFIC = 0.2 
MARBL_TRACERS_IC_FILE = "/glade/campaign/cesm/cesmdata/inputdata/ocn/mom/grid_indpt/ecosys_jan_IC_omip_latlon_1x1_180W_c250613.nc"
MARBL_TRACERS_IC_FILE_IS_Z = True
LEITH_BI_CONST = 80 

• user_nl_marbl

autotroph_settings(2)%kFe = 4e-05
grazing_relationship_settings(3,1)%z_grz = 2.3
autotroph_settings(4)%thetaN_max = 2.1
grazing_relationship_settings(1,2)%z_umax_0_per_day = 0.93
grazing_relationship_settings(1,2)%z_grz = 0.87
grazing_relationship_settings(1,2)%f_zoo_detr = 0.38
autotroph_settings(2)%mort_per_day = 0.14
grazing_relationship_settings(2,2)%graze_poc = 0.15
grazing_relationship_settings(2,1)%z_grz = 1.5
grazing_relationship_settings(3,2)%z_grz = 1.5

• SourceMods

/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09d_m.B1850C_MTso_Gris_Marbl.ne30_t233_wgx3.374/SourceMods/src.mom/
marbl_co2calc_mod.F90  

[cecilehannay]

We are adding extra diagnostics for mom-marbl at year 113:

cp /glade/u/home/klindsay/hannay/diag_table.unresolved.374 SourceMods/src.mom/diag_table.unresolved

From Keith:

We don't have out-of-the box support for RUN_TYPE=hybrid with MOM6-MARBL. If we want to start a new hybrid run using BGC from 374, one approach is to output all MARBL tracer values on a z grid in 374 and have the model interpolate that to the model grid in the new run. We're not writing out all MARBL tracers in 374. I have a modified diag_table.unresolved to add them.

We don't have all diagnostic output needed to understand the global carbon balance. My modified diag_table.unresolved has the missing terms added.

Please copy /glade/u/home/klindsay/hannay/diag_table.unresolved.374 to SourceMods/src.mom/diag_table.unresolved