Add fates cross-grid seed dispersal mechanism

E3SM-Project · Jul 24, 2023 · 845ddfe · 845ddfe
1 parent c648dc2
commit 845ddfe
Show file tree

Hide file tree

Showing 4 changed files with 204 additions and 23 deletions.
diff --git a/components/elm/src/data_types/GridcellType.F90 b/components/elm/src/data_types/GridcellType.F90
@@ -126,7 +126,7 @@ subroutine grc_pp_init(this, begg, endg)
     allocate(this%MaxElevation (begg:endg)) ; this%MaxElevation (:) = nan
 
     allocate(this%landunit_indices(1:max_lunit, begg:endg)); this%landunit_indices(:,:) = ispval
-	
+
    ! allocate(this%topounit_indices (begg:endg,1:max_topounits)) ; this%topounit_indices (:,:) = ispval
 
   end subroutine grc_pp_init

diff --git a/components/elm/src/main/elm_driver.F90 b/components/elm/src/main/elm_driver.F90
@@ -1368,6 +1368,9 @@ subroutine elm_drv(doalb, nextsw_cday, declinp1, declin, rstwr, nlend, rdate)
     end do
     !$OMP END PARALLEL DO
 
+    ! Pass fates seed dispersal information to all nodes
+    if (use_fates .and. is_beg_curr_day()) call alm_fates%WrapSeedGlobal()
+
     ! ============================================================================
     ! Determine gridcell averaged properties to send to atm
     ! ============================================================================

diff --git a/components/elm/src/main/elm_initializeMod.F90 b/components/elm/src/main/elm_initializeMod.F90
@@ -84,7 +84,7 @@ subroutine initialize1( )
     use dynSubgridControlMod      , only: dynSubgridControl_init
     use filterMod                 , only: allocFilters
     use reweightMod               , only: reweight_wrapup
-    use topounit_varcon           , only: max_topounits, has_topounit, topounit_varcon_init    
+    use topounit_varcon           , only: max_topounits, has_topounit, topounit_varcon_init
     use elm_varctl                , only: use_top_solar_rad
     !
     ! !LOCAL VARIABLES: