b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321

[cecilehannay]

Purpose:
Repeat 320

• update ocean configuration (tx2_3v3),
• frozen runoff mapping file is a merge of two mapping files: 100 km spreading for the SH and 250 km spreading for the NH.
• Start from a G case at year 8.
• Add volcano bugfix as in 319.

---

Description:
Same as 320 (#304) except:

• change run type

./xmlchange RUN_TYPE=startup

• change runoff

./xmlchange ROF2OCN_ICE_RMAPNAME=/glade/work/gmarques/cesm/tx2_3/runoff_mapping/map_r05_to_tx2_3_nnsm_e100r100sh_e250r250nh_merged_260307.nc
./xmlchange ROF2OCN_LIQ_RMAPNAME=/glade/work/gmarques/cesm/tx2_3/runoff_mapping/map_jra_to_tx2_3_nnsm_e100r100_260306.nc

• change grid

./xmlchange ICE_DOMAIN_MESH=/glade/work/gmarques/cesm/tx2_3/mesh/ESMF_mesh_tx2_3v3_260305_cdf5.nc
./xmlchange OCN_DOMAIN_MESH=/glade/work/gmarques/cesm/tx2_3/mesh/ESMF_mesh_tx2_3v3_260305_cdf5.nc
./xmlchange MASK_MESH=/glade/work/gmarques/cesm/tx2_3/mesh/ESMF_mesh_tx2_3v3_260305_cdf5.nc

• add sourcemods for MAM5 volcanos:

/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/SourceMods/src.cam/
aerosol_optics_cam.F90  
phys_prop.F90

• add to user_nl_cam corrected rad_climate and extra variables to look at MAM5 volcanos

add to fincl1: 'CME', 'CMEIOUT', 'AODSSdn','CME', 'CMEIOUT', 'AODSSdn', 'EXTINCTdn','VOLC_MMR1','VOLC_MMR2','VOLC_MMR3','VOLC_MMR5','VOLC_SAD',
'VOLC_MASS5','BURDENSO4dn'

rad_climate = 'A:Q:H2O', 'N:O2:O2', 'A:CO2:CO2', 'N:ozone:O3', 'A:N2O:N2O', 'A:CH4:CH4', 'N:CFC11STAR:CFC11', 'A:CFC12:CFC12',          'M:mam4_mode1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode1_rrtmg_aeronetdust_sig1.6_dgnh.48_c140304.nc',           'M:mam4_mode2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode2_rrtmg_aitkendust_c141106.nc',          'M:mam4_mode3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode3_rrtmg_aeronetdust_c141106.nc',         'M:mam4_mode4:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode4_rrtmg_c130628.nc',           'N:VOLC_MMR1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode1_c210211.nc',           'N:VOLC_MMR2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode2_c210211.nc',           'N:VOLC_MMR3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.2_mode3_c210211.nc',          'N:VOLC_MMR5:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.2_mode5_c260228.nc'

• user_nl_cice

grid_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.nc"
grid_format = 'nc'
ice_ic      = "/glade/derecho/scratch/gmarques/archive/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.cice.r.0008-01-01-00000.nc"
kmt_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.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

• user_nl_mom

INPUTDIR = "INPUT"
GRID_FILE = "ocean_hgrid_250930.nc"
TOPO_FILE = "ocean_topo_tx2_3v3_260305.nc"
TOPO_EDITS_FILE = ""
TIDAL_ENERGY_FILE = "tidal_energy_tx2_3v3_conserve_260306.nc"
CHL_FILE = "seawifs-clim-1997-2010-tx2_3v3_260305.nc"
CHANNEL_LIST_FILE = "channels_tx2_3v3_260304.txt"
GEOTHERMAL_FILE = "geothermal_davies2013_tx2_3v3_260305.nc"
MLE%
CR_FILE = "mle_cr_labsea_mask_tx2_3v3_260306.nc"
%MLE

MEKE_GEOMETRIC_ALPHA = 0.07
SMAGORINSKY_KH = False
MAXTRUNC = 1.0E6

! Leith+E Settings
TAPER_LEITHY = True
LEITH_BI_CONST = 85.0
LEITHY_CK = 1.0
LEITHY_DEPTH = 1000.0
LEITHY_WIDTH = 500.0

! Background KH
USE_KH_BG_2D = True
KH_BG_2D_FILENAME = "KH_BG_2D_tx2_3v3_260306.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

KHTR_MIN = 50.

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_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.mom6.r.0008-01-01-00000.nc"

TS_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.mom6.r.0008-01-01-00000.nc"
VELOCITY_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.mom6.r.0008-01-01-00000.nc"

TEMP_IC_VAR = "Temp"
SALT_IC_VAR = "Salt"

• create INPUT with all required files

mkdir /glade/derecho/scratch/hannay/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/run/INPUT 
cd /glade/derecho/scratch/hannay/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/run/INPUT 
cp -L /glade/derecho/scratch/gmarques/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/run/INPUT/* . 

---

Case directory:

• Locally (if still available):
  /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321

• On github:
  https://github.com/NCAR/cesm_dev/tree/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321

---

Sandbox:

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

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

---

Diagnostics:

• Diags (if available):
  https://webext.cgd.ucar.edu/BLT1850/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/

---

Output:

• Raw output on derecho (if still available): /glade/derecho/scratch/hannay/archive/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/
• CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321

---

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

---

Extra details:

[cecilehannay]

Add extra variables for Gunter:

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

[cecilehannay]

The run crashed in MOM in year 30 with the error:

dec0970.hsn.de.hpc.ucar.edu 3433: WARNING from PE   617: Extreme surface sfc_state detected: i=  84 j= 477 lon=-283.993 lat=  72.281 x=-231.333 y\
=  87.552 D= 1.1179E+01 SSH= 4.7348E+00 SST=-3.0066E+00 SSS= 5.5745E+01 U-= 3.9000E-06 U+= 2.0746E-01 V-= 1.5171E-02 V+= 1.3288E-07
dec0970.hsn.de.hpc.ucar.edu 3433:
dec0970.hsn.de.hpc.ucar.edu 3433:
dec0970.hsn.de.hpc.ucar.edu 3433: FATAL from PE   617: There were a total of         1 locations detected with extreme surface values!
dec0970.hsn.de.hpc.ucar.edu 3433:
dec0970.hsn.de.hpc.ucar.edu 3433: cesm.exe           0000000004175949  mpp_mod_mp_mpp_er          72  mpp_util_mpi.inc
dec0970.hsn.de.hpc.ucar.edu 3433: cesm.exe           0000000002D9ABC8  mom_error_handler         191  MOM_error_handler.F90
dec0970.hsn.de.hpc.ucar.edu 3433: cesm.exe           00000000029EC871  mom_mp_extract_su        4414  MOM.F90
dec0970.hsn.de.hpc.ucar.edu 3433: cesm.exe           00000000029CF54E  mom_mp_step_mom_         1110  MOM.F90
dec0970.hsn.de.hpc.ucar.edu 3433: cesm.exe           0000000002986FD9  mom_ocean_model_n         658  mom_ocean_model_nuopc.F90
dec0970.hsn.de.hpc.ucar.edu 3433: cesm.exe           0000000002941E66  mom_cap_mod_mp_mo        1865  mom_cap.F90

[cecilehannay]

RESTOM In 321 and 322 seems weird:

[iangrooms]

Surface salinity reaches >55 PSU in that spot. In 316 the surface salinity stays below 15 PSU there.

There appears to be no river runoff into the Arctic in 321, which would explain the excessive salinity. I think that the liquid runoff map is using the JRA map instead of the r05 map, which is the problem. The JRA source map is on a 1/4 degree grid, while the r05 map is on a 1/2 degree grid, so the liquid runoff in 321 is all going to the wrong places. I think it should have been

./xmlchange ROF2OCN_LIQ_RMAPNAME=/glade/work/gmarques/cesm/tx2_3/runoff_mapping/map_r05_to_tx2_3_nnsm_e100r100_260306.nc

[gustavo-marques]

Yes, we are using the wrong file (JRA instead of r05) for ROF2OCN_LIQ_RMAPNAME. This is also the case for #306. Thanks for catching this, @iangrooms.
@cecilehannay, sorry for missing this when I gave you the instructions. We will need to redo these runs using ./xmlchange ROF2OCN_LIQ_RMAPNAME=/glade/work/gmarques/cesm/tx2_3/runoff_mapping/map_r05_to_tx2_3_nnsm_e100r100_260306.nc

[cecilehannay]

No worries. I will restart with that file.