Skip to content

Include chemical heating in mid-top chemistry compsets#1478

Open
fvitt wants to merge 9 commits intoESCOMP:cam_developmentfrom
fvitt:midtop_waphys
Open

Include chemical heating in mid-top chemistry compsets#1478
fvitt wants to merge 9 commits intoESCOMP:cam_developmentfrom
fvitt:midtop_waphys

Conversation

@fvitt
Copy link
Collaborator

@fvitt fvitt commented Feb 2, 2026

This includes chemical heating in mid-top compsets. This is done by including reactions needed for chemical heating in the middle atmosphere which is achieved via waccm physics.

Closes #1476

Here we adopt S1 for stratospheric chemistry used in low-top configurations, S2 for stratospheric chemistry used in mid-top configurations.

fvitt added 8 commits December 9, 2025 08:12
	modified:   src/chemistry/mozart/mo_jeuv.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/chem_mech.doc
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/chem_mech.in
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/chem_mods.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/m_rxt_id.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/m_spc_id.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_adjrxt.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_indprd.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_lin_matrix.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_lu_factor.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_lu_solve.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_nln_matrix.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_phtadj.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_prod_loss.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_rxt_rates_conv.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_setrxt.F90
	modified:   src/chemistry/pp_trop_strat_mam5_ts4/mo_sim_dat.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/chem_mech.doc
	modified:   src/chemistry/pp_waccm_t4ma_mam5/chem_mech.in
	modified:   src/chemistry/pp_waccm_t4ma_mam5/chem_mods.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/m_spc_id.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_indprd.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_lin_matrix.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_lu_factor.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_lu_solve.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_nln_matrix.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_prod_loss.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_rxt_rates_conv.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_setrxt.F90
	modified:   src/chemistry/pp_waccm_t4ma_mam5/mo_sim_dat.F90
	modified:   bld/namelist_files/namelist_defaults_cam.xml
	modified:   bld/namelist_files/use_cases/hist_trop_strat_t4s_cam7.xml
	modified:   src/chemistry/mozart/mo_waccm_hrates.F90
	modified:   bld/namelist_files/namelist_defaults_cam.xml
	renamed:    src/chemistry/pp_trop_strat_mam5_vbs -> src/chemistry/pp_trop_strat_mam5_t1s1
	renamed:    src/chemistry/pp_trop_strat_mam5_ts2 -> src/chemistry/pp_trop_strat_mam5_t2s1
	renamed:    src/chemistry/pp_trop_strat_mam5_ts4 -> src/chemistry/pp_trop_strat_mam5_t4s2
	new file:   src/chemistry/pp_trop_strat_mam5_t1s2

	modified:   bld/build-namelist
	modified:   bld/config_files/definition.xml
	modified:   bld/configure
	modified:   bld/namelist_files/namelist_defaults_cam.xml
	modified:   bld/namelist_files/use_cases/hist_trop_strat_nudged_cam6.xml
	modified:   cime_config/config_component.xml
	modified:   cime_config/config_compsets.xml
	modified:   bld/namelist_files/namelist_defaults_cam.xml
	modified:   bld/namelist_files/use_cases/hist_trop_strat_nudged_cam6.xml
	modified:   bld/namelist_files/use_cases/hist_trop_strat_t4s_cam7.xml
	modified:   cime_config/config_component.xml
	modified:   cime_config/config_compsets.xml
@fvitt fvitt added this to the CESM3 milestone Feb 2, 2026
@fvitt fvitt requested a review from tilmes February 2, 2026 22:40
@fvitt fvitt self-assigned this Feb 2, 2026
add_default($nl, 'megan_mapped_emisfctrs', 'val'=>'.false.');
}
if ($chem =~ /trop_strat_mam4_vbs/ or $chem =~ /trop_strat_mam5_vbs/) {
if ($chem =~ /trop_strat_mam4_vbs/ or $chem =~ /trop_strat_mam5_vbs/ or $chem =~ /trop_strat_mam5_t1s/) {
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It appears that trop_strat_mam5_vbs is being phased out/renamed and this should be removed.

Also, namelist_definition.xml still has reference to trop_strat_mam5_vbs as well

@cacraigucar
Copy link
Collaborator

@tilmes - I will need to have a review from you as well before Francis can proceed with tagging this PR. That or you can remove yourself as a reviewer

Copy link
Collaborator

@tilmes tilmes left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I have briefly check the code, and agree that we only change the L93 MT -cam7 physics case for now to T4S2 and T1S2, while all the other vert. resoluitons will not be changed.

	modified:   bld/build-namelist
	modified:   bld/namelist_files/namelist_definition.xml
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants