Skip to content

allocate space to disturbance_rates in the restart file #1588

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
rgknox wants to merge 1 commit into
base: main
Choose a base branch
from
+11 −4
Open

allocate space to disturbance_rates in the restart file #1588

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
10 changes: 8 additions & 2 deletions main/FatesInterfaceMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,6 +38,7 @@ module FatesInterfaceMod
use FatesConstantsMod , only : TRS_regeneration
use FatesConstantsMod , only : g_per_kg
use FatesConstantsMod , only : n_landuse_cats
use FatesConstantsMod , only : n_dist_types
use FatesConstantsMod , only : primaryland
use FatesConstantsMod , only : secondaryland
use FatesConstantsMod , only : n_crop_lu_types
Expand Down Expand Up @@ -946,12 +947,17 @@ subroutine SetFatesGlobalElements2(use_fates)
end if

fates_maxElementsPerSite = max(fates_maxPatchesPerSite * fates_maxElementsPerPatch, &
numWatermem*numpft, num_vegtemp_mem, num_elements, nlevsclass*numpft*n_term_mort_types)
numWatermem*numpft, num_vegtemp_mem, num_elements*ncwd, num_elements*numpft, &
nlevsclass*numpft*n_term_mort_types)

if(hlm_use_planthydro==itrue)then
fates_maxElementsPerSite = max(fates_maxElementsPerSite, nshell*nlevsoi_hyd_max )
end if


! Need enough restart elements to accomodate sites(s)%disturbance_rates
fates_maxElementsPerSite = max(fates_maxElementsPerSite,n_landuse_cats*n_landuse_cats*n_dist_types)
fates_maxElementsPerSite = max(fates_maxElementsPerSite,n_landuse_cats*numpft)


! Set the maximum number of nutrient aquisition competitors per site
! This is used to set array sizes for the boundary conditions.
Expand Down
5 changes: 3 additions & 2 deletions main/FatesRestartInterfaceMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1764,9 +1764,10 @@ subroutine define_restart_vars(this, initialize_variables)
long_name='24-hour patch veg temp', &
units='K', initialize=initialize_variables,ivar=ivar, index = ir_tveg24_pa)

call this%DefineRMeanRestartVar(vname='fates_disturbance_rates',vtype=cohort_r8, &
call this%set_restart_var(vname='fates_disturbance_rates',vtype=cohort_r8, &
long_name='disturbance rates by donor land-use type, receiver land-use type, and disturbance type', &
units='1/day', initialize=initialize_variables,ivar=ivar, index = ir_disturbance_rates_siluludi)
units='1/day', flushval = flushzero, &
hlms='CLM:ALM',initialize=initialize_variables, ivar=ivar, index = ir_disturbance_rates_siluludi)

if ( hlm_regeneration_model == TRS_regeneration ) then

Expand Down