Skip to content

Conversation

@fvitt
Copy link
Collaborator

@fvitt fvitt commented Jan 7, 2026

Closes #1469

fvitt added 2 commits January 6, 2026 16:06
	new file:   cime_config/testdefs/testmods_dirs/cam/ts_slh_mam4_outfrq9s/shell_commands
	new file:   cime_config/testdefs/testmods_dirs/cam/ts_slh_mam4_outfrq9s/user_nl_cam
	new file:   cime_config/testdefs/testmods_dirs/cam/ts_slh_mam4_outfrq9s/user_nl_clm
	new file:   cime_config/testdefs/testmods_dirs/cam/wa_slh_mam4_outfrq9s/shell_commands
	new file:   cime_config/testdefs/testmods_dirs/cam/wa_slh_mam4_outfrq9s/user_nl_cam
	new file:   cime_config/testdefs/testmods_dirs/cam/wa_slh_mam4_outfrq9s/user_nl_clm
@fvitt fvitt self-assigned this Jan 7, 2026
@fvitt fvitt marked this pull request as draft January 8, 2026 00:09
@fvitt
Copy link
Collaborator Author

fvitt commented Jan 9, 2026

From @RafaPedroFernandez (CC'd @CAlbertoCuevas):

Please note that I have just made a push for 2 F90 routines in my following GIT:
https://github.com/RafaPedroFernandez/CAM/tree/slh_chemistry_n04_cam6_4_097_rpf

Which is also updated local in derecho:
/glade/u/home/rpfernan/cesm_versions_rpf/derecho_cesm/git_fetch/cesm3_cam6dev_francis/slh_chemistry_n04_cam6_4_097_rpf

Routines are
chemistry.F90:
For SLH, we must call iodine_emissions_srf before the loop to outfld( sflxnam(m), cam_in%cflx(:ncol,m), ncol,lchnk ) in order for SFHOI and SFI2 to be properly written in the output file. If we do not so, the SFHOI and SFI2 are zero in the output file but are somehow included in cam_in%cflx array that is passed to other routines

mo_gas_phase_chmdr.F90
For Iodine washout ... we have always used SAD_ICETROP, which is computed based on CLDICE. This CLDICE was before computed as only ice
cldice(:ncol,k) = q(:ncol,k,cldice_ndx)
and now considers also SNOW
cldice(:ncol,k) = q(:ncol,k,cldice_ndx) + q(:ncol,k,snow_ndx)

Including snow in the variable makes the iodine deposition too efficient, and then we created a new array (ONLYICE) just to compute iodine wet-depositions without affecting other model developments.

	modified:   src/chemistry/mozart/chemistry.F90
	modified:   src/chemistry/mozart/mo_gas_phase_chemdr.F90
@fvitt fvitt marked this pull request as ready for review January 13, 2026 23:11
	modified:   src/chemistry/mozart/mo_usrrxt.F90
Copy link
Collaborator

@cacraigucar cacraigucar left a comment

Choose a reason for hiding this comment

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

Partial review - @jimmielin will be reviewing the chemistry section of this PR. My review was on the setup side.

I noticed that there are testing use_cases setup, but I do not see any regression tests being added. Is this an oversight?

Comment on lines +2064 to +2065
<xs_short_file>atm/waccm/phot/xs_short_c221005mm_slh_c250930.nc</xs_short_file>
<xs_long_file >atm/waccm/phot/temp_prs_GT200nm_c221005mm_slh_c250930.nc</xs_long_file>
Copy link
Collaborator

Choose a reason for hiding this comment

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

These two file still need to be "rimport"ed

Comment on lines +6312 to +6340
<entry id="SSAhno3_ScalingFactor" type="real" category="cam_chem"
group="slh_nl" valid_values="" >
Scaling-factor for SLH Halocarbon surface emissions
Default: set by build-namelist
</entry>

<entry id="SSAn2o5_ScalingFactor" type="real" category="cam_chem"
group="slh_nl" valid_values="" >
Scaling-factor for SLH Halocarbon surface emissions
Default: set by build-namelist
</entry>

<entry id="LIQfraprx_ScalingFactor_I" type="real" category="cam_chem"
group="slh_nl" valid_values="" >
Scaling-factor for SLH Halocarbon surface emissions
Default: set by build-namelist
</entry>

<entry id="ICEfraprx_ScalingFactor_I" type="real" category="cam_chem"
group="slh_nl" valid_values="" >
Scaling-factor for SLH Halocarbon surface emissions
Default: set by build-namelist
</entry>

<entry id="ICEfraprx_ScalingFactor_Br" type="real" category="cam_chem"
group="slh_nl" valid_values="" >
Scaling-factor for SLH Halocarbon surface emissions
Default: set by build-namelist
</entry>
Copy link
Collaborator

Choose a reason for hiding this comment

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

Comments probably need updating (since they are all the same)

<value compset="HIST_CAM60%CFIRE">1995-01-01</value>
<value compset="RCP[2468]_CAM\d+">2005-01-01</value>
<value compset="_CAM.*%NUDG" >2010-01-01</value>
<value compset="_CAM.*%NUDG" >2000-01-01</value>
Copy link
Collaborator

Choose a reason for hiding this comment

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

Why is this PR changing the date for nudging runs?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

Status: No status

Development

Successfully merging this pull request may close these issues.

2 participants