b.e30_alpha09a.B1850C_MTso.ne30_t233_wgx3.346

[gustavo-marques]

Purpose:
Evaluate the new CESM tag (alpha09a) in conjunction with the deep convection (ZM) bug fix.

---

Description:
Same as 345 (#337), but using alpha09a.

---

Case directory:

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

• On github:
  https://github.com/NCAR/cesm_dev/tree/b.e30_alpha09a.B1850C_MTso.ne30_t233_wgx3.346

---

Sandbox:

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

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

---

Diagnostics:

• Diags (if available):
  https://webext.cgd.ucar.edu/BLT1850/b.e30_alpha09a.B1850C_MTso.ne30_t233_wgx3.346/

---

Output:

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

---

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_alpha09a.B1850C_MTso.ne30_t233_wgx3.346

• Create new case

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

• 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 OCN_DIAG_MODE

./xmlchange OCN_DIAG_MODE=spinup

• change to CAM_CONFIG_OPTS

./xmlchange CAM_CONFIG_OPTS='-pcols 9 -phys cam7 -model_top mt'

• put the case in the database

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

• start from 341 at year 107

./xmlchange RUN_REFCASE=b.e30_alpha08o.B1850C_MTso.ne30_t232_wgx3.341
./xmlchange RUN_REFDATE=0107-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'

flbc_file  = '/glade/u/home/jzhan166/scripts/LBC/clean/LBC_17500116-20221216_CMIP7_0p5degLat_OCSupdate_c250523.nc'

clubb_gamma_coef                =  0.298     
clubb_gamma_coefb               =  0.298   

• user_nl_clm

reseed_dead_plants = .true.
use_excess_ice_streams=.false. 

• user_nl_cpl

histaux_l2x1yrg = .true.

• user_nl_cice

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

• SourceMods
  There are SourceMods in cam.
  Make sure to copy these *F90 files

/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha08o.B1850C_MTso.ne30_t232_wgx3.345/SourceMods/src.cam/ 
air_composition.F90  
physpkg.F90 
zm_conv_momtran.F90 

[dabail10]

The volume does seem higher in 346 compared to 341.

https://webext.cgd.ucar.edu/BLT1850/b.e30_alpha09a.B1850C_MTso.ne30_t233_wgx3.346/ice/html/infrastructure/index.html

[iangrooms]

Two ocean changes in alpha09a that might be responsible for differences compared to 345:

BBL_EFFIC 0.2 [345] -> 0.01 [346]
MLE_DENSITY_DIFF 0.02 [345] -> 0.03 [346]

Both changes are moving to recommended values in the literature. They were tested in #322 (run 336 has these changes, compared to run 335 without). There may be other differences between 345 and 346 that are relevant; these are simply the changes in MOM6 parameters.

[iangrooms]

ACC transport in 336 is very similar to and slightly better than in 335 over about 30 years, so it seems these two ocean parameter changes are probably not responsible for the ACC transport reduction going from 345 to 346.

SST is very similar in 335 and 336, but maybe a little warmer in 336, which reduces a cold bias. It's plausible that the MLE_DENSITY_DIFF parameter change is warming SSTs in alpha09a (346). Having a greater MLE density difference means that the MLD computed by the MLE scheme will be deeper. Since the strength of the submesoscale restratification scales with the MLD, this would lead to stronger restratification which is consistent with warmer SSTs.

[gustavo-marques]

Thanks, @iangrooms!
I will run a clone of 346 with MLE_DENSITY_DIFF = 0.02 to see if that helps.