Skip to content

Commit

Permalink
GEOS-Chem species unit conversion updates for compatibility with 14.3
Browse files Browse the repository at this point in the history
GEOS-Chem unit variable is now integer and not string

Signed-off-by: Lizzie Lundgren <[email protected]>
  • Loading branch information
lizziel committed Feb 28, 2024
1 parent 666da3c commit bba809e
Showing 1 changed file with 5 additions and 3 deletions.
8 changes: 5 additions & 3 deletions src/chemistry/geoschem/chemistry.F90
Original file line number Diff line number Diff line change
Expand Up @@ -990,6 +990,7 @@ subroutine chem_init(phys_state, pbuf2d)
use TaggedDiagList_Mod, only : Init_TaggedDiagList, Print_TaggedDiagList
use Time_Mod, only : Accept_External_Date_Time
use Ucx_Mod, only : Init_Ucx
use Unitconv_Mod, only : MOLES_SPECIES_PER_MOLES_DRY_AIR
use Vdiff_Mod, only : Max_PblHt_For_Vdiff

TYPE(physics_state), INTENT(IN ) :: phys_state(BEGCHUNK:ENDCHUNK)
Expand Down Expand Up @@ -1420,7 +1421,7 @@ subroutine chem_init(phys_state, pbuf2d)
ENDIF

! Start with v/v dry (CAM standard)
State_Chm(I)%Spc_Units = 'v/v dry'
State_Chm(I)%Spc_Units = MOLES_SPECIES_PER_MOLES_DRY_AIR

ENDDO
Input_Opt%amIRoot = MasterProc
Expand Down Expand Up @@ -1864,7 +1865,8 @@ subroutine chem_timestep_tend( state, ptend, cam_in, cam_out, dT, pbuf, fh2o )
use Time_Mod, only : Accept_External_Date_Time
use Toms_Mod, only : Compute_Overhead_O3
use UCX_Mod, only : Set_H2O_Trac
use Unitconv_Mod, only : Convert_Spc_Units, KG_SPECIES_PER_KG_DRY_AIR, UNIT_STR
use Unitconv_Mod, only : Convert_Spc_Units, UNIT_STR
use Unitconv_Mod, only : KG_SPECIES_PER_KG_DRY_AIR
use Wetscav_Mod, only : Setup_Wetscav

REAL(r8), INTENT(IN) :: dT ! Time step
Expand Down Expand Up @@ -2100,7 +2102,7 @@ subroutine chem_timestep_tend( state, ptend, cam_in, cam_out, dT, pbuf, fh2o )

! 2. Copy tracers into State_Chm
! Data was received in kg/kg dry
State_Chm(LCHNK)%Spc_Units = 'kg/kg dry'
State_Chm(LCHNK)%Spc_Units = KG_SPECIES_PER_KG_DRY_AIR
! Initialize ALL State_Chm species data to zero, not just tracers
DO N = 1, State_Chm(LCHNK)%nSpecies
State_Chm(LCHNK)%Species(N)%Conc = 0.0e+0_fp
Expand Down

0 comments on commit bba809e

Please sign in to comment.