Skip to content

Commit

Permalink
Bug fix: avoid double-counting of HCl uptake on alkaline sea salt aer…
Browse files Browse the repository at this point in the history
…osols

This fix corrects the issue first reported by Viral Shah in
geoschem/geos-chem #1034.  We only need to do HCL uptake when
SO2 chemistry is called with when FullRun=.TRUE..  The quick fix was
to add FullRun to the IF block where the computations occur.

NOTE: This issue does not occur in the new KPP sulfur chemistry code:
KPP/fullchem/fullchem_SulfurChemFuncs.F90.

Signed-off-by: Bob Yantosca <yantosca@seas.harvard.edu>
  • Loading branch information
yantosca committed Dec 2, 2021
1 parent 4fe42a4 commit 841e216
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion GeosCore/sulfate_mod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -4705,7 +4705,7 @@ SUBROUTINE SEASALT_CHEM( I, J, L, &
ENDIF

! HCl lost [eq/timestep] converted back to [v/v/timestep]
IF ( id_HCl > 0 ) THEN
IF ( FullRun .and. id_HCl > 0 ) THEN
HCl_ss = ( TITR_HCl * AIRMW / AD(I,J,L) ) / 1000.0_fp
Spc(I,J,L,id_HCl) = MAX( HCl_vv - HCl_ss, MINDAT )
ENDIF
Expand Down

0 comments on commit 841e216

Please sign in to comment.