From 0276c54d84da1b7d77c864cb1c27601b727e54d5 Mon Sep 17 00:00:00 2001 From: Lizzie Lundgren Date: Tue, 23 Jul 2024 10:56:39 -0400 Subject: [PATCH] Uncomment setting Ca2, K, and Mg cation values passed to HETP Previously these values were zero due to an issue in ISORROPIA. With the update to HETP we can now enable these cations. Signed-off-by: Lizzie Lundgren --- CHANGELOG.md | 1 + GeosCore/aerosol_thermodynamics_mod.F90 | 34 ++++++++----------------- 2 files changed, 12 insertions(+), 23 deletions(-) diff --git a/CHANGELOG.md b/CHANGELOG.md index 53f17f99c..557b60c16 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -11,6 +11,7 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/), ### Changed - Copy values from `State_Chm%KPP_AbsTol` to `ATOL` and `State_Chm%KPP_RelTol` to `RTOL` for fullchem and Hg simulations +- Change previously zero Ca2, K, and Mg cation values passed to HETP to scaled SALA species concentrations ### Fixed - In `Headers/roundoff_mod.F90`, first cast and then only round off if `places > 0` diff --git a/GeosCore/aerosol_thermodynamics_mod.F90 b/GeosCore/aerosol_thermodynamics_mod.F90 index a5a1559e8..7dcb2b58c 100644 --- a/GeosCore/aerosol_thermodynamics_mod.F90 +++ b/GeosCore/aerosol_thermodynamics_mod.F90 @@ -673,29 +673,17 @@ SUBROUTINE DO_ATE( Input_Opt, State_Chm, State_Diag, & ! Total Cl- [mole/m3] TCL = ACL + GCL - !==================================================================== - ! NOTE: As of 11/2007, ISORROPIAII does not conserve mass when Ca,K,Mg - ! are non-zero. If you would like to consider Ca, K, Mg from seasalt - ! and dust, then ISORROPIA II ISRP4F routines must be debugged. - ! (hotp, bmy, 2/1/10) - ! This still applies in HETP (ewl, 4/22/2024) - ! - ! ! Total Ca2+ (1.16% by weight of seasalt) [mole/m3] - ! TCA = Spc(id_SALA)%Conc(I,J,L) * 0.0116e+0_fp * 1.d3 / - ! ( 40.08e+0_fp * VOL ) - ! - ! ! Total K+ (1.1% by weight of seasalt) [mole/m3] - ! TK = Spc(id_SALA)%Conc(I,J,L) * 0.0110e+0_fp * 1.d3 / - ! ( 39.102e+0_fp * VOL ) - ! - ! ! Total Mg+ (3.69% by weight of seasalt) [mole/m3] - ! TMG = Spc(id_SALA)%Conc(I,J,L) * 0.0369e+0_fp * 1.d3 / - ! ( 24.312e+0_fp * VOL ) - !==================================================================== - !! Set Ca, K, Mg to zero for time being (hotp, bmy, 2/1/10) - !TCA = 0e+0_fp - !TK = 0e+0_fp - !TMG = 0e+0_fp + ! Total Ca2+ (1.16% by weight of seasalt) [mole/m3] + TCA = Spc(id_SALA)%Conc(I,J,L) * 0.0116e+0_fp * 1.d3 / & + ( 40.08e+0_fp * VOL ) + + ! Total K+ (1.1% by weight of seasalt) [mole/m3] + TK = Spc(id_SALA)%Conc(I,J,L) * 0.0110e+0_fp * 1.d3 / & + ( 39.102e+0_fp * VOL ) + + ! Total Mg+ (3.69% by weight of seasalt) [mole/m3] + TMG = Spc(id_SALA)%Conc(I,J,L) * 0.0369e+0_fp * 1.d3 / & + ( 24.312e+0_fp * VOL ) ! Compute gas-phase NO3 IF ( id_HNO3 > 0 ) THEN