b.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.355

[gustavo-marques]

Purpose:
Check that applying a taper (depth = 900, power = 4) to both MLDs and BLDs used in the Bodner scheme avoids negative overturing circulation in the Nordic Seas and large oscillations in the ACC. Starting from year 101 of run 334.

Here is the taper function used:

---

Description:
Same as 354, but using Bonder with tapered MLDs and BLDs via srcMods.

---

Case directory:

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

• On github:
  https://github.com/NCAR/cesm_dev/tree/b.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.355

---

Sandbox:

• Locally (if still available): /glade/work/gmarques/cesm.sandboxes/cesm3_0_alpha09b

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

---

Diagnostics:

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

---

Output:

• Raw output on derecho (if still available): /glade/derecho/scratch/gmarques/archive/b.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.355/
• CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.355

---

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

Extra details:

Here is the instructions to reproduce b.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.355

• Create new case

./create_newcase --compset B1850C_MTso --res ne30pg3_t233_wg37 
--case /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.355 

• copy pe layout from previous case

/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha08o.B1850C_MTso.ne30_t232_wgx3.334/env_mach_pes.xml 

• change to CAM_CONFIG_OPTS

./xmlchange CAM_CONFIG_OPTS='-pcols 9' --amend

• put the case in the database

./xmlchange CASE_GIT_REPOSITORY=git@github.com:NCAR/cesm_dev.git

• start from 334 at year 101

./xmlchange RUN_REFCASE=b.e30_alpha08o.B1850C_MTso.ne30_t232_wgx3.334
./xmlchange RUN_REFDATE=0101-01-01
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=true
./xmlchange RUN_REFDIR=cesm2_init

• user_nl_cam

mfilt               =       0,       5,     20,      40,      12,      120,      1,   1
nhtfrq              =       0,     -24,    -24,      -3,       0,       -2,      0,  -8760
ndens               =       2,       2,      2,       2,       2,        1,      2,   1
interpolate_output  =  .true.,  .true., .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

reseed_dead_plants = .true.
use_excess_ice_streams=.false. 

• user_nl_mom

DIAG_COORD_DEF_RHO2 = "FILE:ocean_rho2_190917.nc,interfaces=rho2"
MLE_TAPER_FN_DEPTH = 900.
MLE_TAPER_FN_POWER = 4

• 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

• user_nl_cice

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

• SourceMods
  SourceMods/src.mom/MOM_mixed_layer_restrat.F90