two changes in bgc submodules from pflotran v5.0.0

thorntonpe · Aug 15, 2024 · 16dbcf1 · 16dbcf1
1 parent fefaaa9
commit 16dbcf1
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/components/elm/src/external_models/emi/src/em/alquimia/ExternalModelAlquimiaMod.F90 b/components/elm/src/external_models/emi/src/em/alquimia/ExternalModelAlquimiaMod.F90
@@ -44,6 +44,8 @@ module ExternalModelAlquimiaMod
 
  type, public, extends(em_base_type) :: em_alquimia_type
 
+ integer :: natural_id
+
  integer :: index_l2e_col_dz
  integer :: index_l2e_col_zi
 
@@ -2217,7 +2219,7 @@ recursive subroutine run_onestep(this,dt,num_cuts,max_cuts)
 
  porosity=this%chem_state%porosity
  call this%chem%ReactionStepOperatorSplit(this%chem_engine, actual_dt, this%chem_properties, this%chem_state, &
- this%chem_aux_data, this%chem_status)
+ this%chem_aux_data, this%natural_id, this%chem_status)
  ! Reset porosity because Pflotran tends to mess it up
  this%chem_state%porosity=porosity
  ! write(iulog,*),'Converged =',this%chem_status%converged,"ncuts =",num_cuts
@@ -2362,7 +2364,7 @@ recursive subroutine run_column_onestep(this,dt,num_cuts,max_cuts, &
 
  porosity_tmp=this%chem_state%porosity
  call this%chem%ReactionStepOperatorSplit(this%chem_engine, actual_dt, this%chem_properties, this%chem_state, &
- this%chem_aux_data, this%chem_status)
+ this%chem_aux_data, this%natural_id, this%chem_status)
  ! Reset porosity because Pflotran tends to mess it up
  this%chem_state%porosity=porosity_tmp