-
Notifications
You must be signed in to change notification settings - Fork 29
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Refactoring: move cohort type to its own file. Functionalize in prepa…
…ration for OOP
- Loading branch information
marcadella
committed
Jan 9, 2025
1 parent
08591f8
commit a03df9f
Showing
7 changed files
with
222 additions
and
328 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,160 @@ | ||
| module md_cohort | ||
| !//////////////////////////////////////////////////////////////// | ||
| ! Module containing BiomeE state variable and parameter | ||
| ! definitions. | ||
| ! Code adopted from BiomeE https://doi.org/10.5281/zenodo.7125963. | ||
| !---------------------------------------------------------------- | ||
| use md_interface_biomee, only: myinterface, spec_data_type, MAX_LEVELS | ||
| use md_params_core, only: pi | ||
| use md_classdefs | ||
|
|
||
| ! define data types and constants | ||
| implicit none | ||
| !=============== Public types =========================================================== | ||
| public :: cohort_type | ||
|
|
||
| !=============== Public subroutines ===================================================== | ||
| public :: rootarea, W_supply, NSNmax | ||
|
|
||
| integer, public, parameter :: LEAF_OFF = 0 | ||
| integer, public, parameter :: LEAF_ON = 1 | ||
| integer, public, parameter :: PT_C3 = 0 ! physiology types | ||
| integer, public, parameter :: PT_C4 = 1 ! physiology types | ||
|
|
||
| !=============== Cohort level data type ============================================================= | ||
| type :: cohort_type | ||
|
|
||
| !===== Metadata | ||
| integer :: ccID = -1 ! cohort ID | ||
| integer :: layer = 1 ! the layer of this cohort (numbered from top, top layer=1) | ||
| real :: layerfrac = 0.0 ! fraction of layer area occupied by this cohort | ||
| integer :: firstlayer = 0 ! 0 = never been in the first layer; 1 = at least one year in first layer | ||
|
|
||
| !===== Population structure | ||
| real :: nindivs = 1.0 ! density of vegetation, tree/m2 | ||
| real :: age = 0.0 ! age of cohort, years | ||
| real :: topyear = 0.0 ! the years that a cohort is in top layer | ||
| real :: dbh = 0.0 ! diameter at breast height, m | ||
| real :: height = 0.0 ! vegetation height, m | ||
| real :: crownarea = 1.0 ! crown area, m2 tree-1 | ||
| real :: leafarea = 0.0 ! total area of leaves, m2 tree-1 | ||
|
|
||
| !===== Biological prognostic variables | ||
| integer :: species = 1 ! vegetation species | ||
| real :: gdd = 0.0 ! growing degree-day (phenology) | ||
| integer :: status = LEAF_OFF ! growth status of plant | ||
| real :: leaf_age = 0.0 ! leaf age (years) | ||
|
|
||
| !===== Organic pools | ||
| type(orgpool) :: pleaf ! leaf biomass, kg tree-1 | ||
| type(orgpool) :: proot ! root biomass, kg tree-1 | ||
| type(orgpool) :: psapw ! sapwood biomass, kg tree-1 | ||
| type(orgpool) :: pwood ! heartwood (non-living) biomass, kg tree-1 | ||
| type(orgpool) :: pseed ! biomass put aside for future progeny, kg tree-1 | ||
| type(orgpool) :: plabl ! labile pool, temporary storage of N and C, kg tree-1 | ||
|
|
||
| !===== Fast step fluxes, kg timestep-1 tree-1 | ||
| type(common_fluxes) :: fast_fluxes | ||
|
|
||
| real :: resl = 0.0 ! leaf respiration, kg C timestep-1 tree-1 | ||
| real :: resr = 0.0 ! root respiration, kg C timestep-1 tree-1 | ||
|
|
||
| !===== Daily fluxes, kg day-1 tree-1 | ||
| type(common_fluxes) :: daily_fluxes | ||
|
|
||
| !===== Annual fluxes, kg yr-1 tree-1 | ||
| type(common_fluxes) :: annual_fluxes | ||
| real :: NPPleaf = 0.0 ! C allocated to leaf, kg C yr-1 tree-1 | ||
| real :: NPProot = 0.0 ! C allocated to root, kg C yr-1 tree-1 | ||
| real :: NPPwood = 0.0 ! C allocated to wood, kg C yr-1 tree-1 | ||
|
|
||
| !===== Annual fluxes due to tree death, kg m-2 yr-1 | ||
| real :: n_deadtrees = 0.0 ! plant to soil N flux due to mortality (kg N m-2 yr-1) | ||
| real :: c_deadtrees = 0.0 ! plant to soil C flux due to mortality (kg C m-2 yr-1) | ||
| real :: m_turnover = 0.0 ! C turnover due to mortality and tissue turnover (kg C m-2 yr-1) | ||
| real :: deathrate = 0.0 ! Deathrate (0 to 1) | ||
|
|
||
| !===== Nitrogen model related variables | ||
| real :: bl_max = 0.0 ! Max. leaf biomass, kg C tree-1 | ||
| real :: br_max = 0.0 ! Max. fine root biomass, kg C tree-1 | ||
|
|
||
| !===== Water uptake-related variables | ||
| real :: WupL(MAX_LEVELS) = 0.0 ! normalized vertical distribution of uptake | ||
|
|
||
| !===== Photosynthesis variables | ||
| real :: An_op = 0.0 ! mol C/(m2 of leaf per year) | ||
| real :: An_cl = 0.0 ! mol C/(m2 of leaf per year) | ||
| real :: C_growth = 0.0 ! Carbon gain since last growth, kg C day-1 tree-1 | ||
| real :: N_growth = 0.0 ! Nitrogen used for plant tissue growth, kg N day-1 tree-1 | ||
| real :: resg = 0.0 ! growth respiration, kg C day-1 tree-1 | ||
|
|
||
| !===== Memory variables used for computing deltas | ||
| real :: DBH_ys = 0.0 ! DBH at the begining of a year (growing season) | ||
| real :: BA_ys = 0.0 ! Basal area at the beginning og a year | ||
|
|
||
| end type cohort_type | ||
|
|
||
| contains | ||
| function NSNmax(cohort) result(res) | ||
| real :: res | ||
| type(cohort_type) :: cohort | ||
|
|
||
| ! Local variable | ||
| type(spec_data_type) :: sp | ||
|
|
||
| sp = myinterface%params_species(cohort%species) | ||
|
|
||
| res = sp%fNSNmax * & | ||
| (cohort%bl_max / (sp%CNleaf0 * sp%leafLS) + cohort%br_max / sp%CNroot0) | ||
| end function NSNmax | ||
|
|
||
| function W_supply(cohort) result(res) | ||
| ! potential water uptake rate per unit time per tree | ||
| real :: res | ||
| type(cohort_type) :: cohort | ||
|
|
||
| res = sum(cohort%WupL(:)) | ||
| end function W_supply | ||
|
|
||
| function rootareaL(cohort, level) result(res) | ||
| ! Root length per layer, m of root/m | ||
| real :: res | ||
| integer :: level | ||
| type(cohort_type) :: cohort | ||
|
|
||
| res = rootarea(cohort) * myinterface%params_species(cohort%species)%root_frac(level) | ||
| end function rootareaL | ||
|
|
||
| function rootarea(cohort) result(res) | ||
| ! total fine root area per tree | ||
| real :: res | ||
| type(cohort_type) :: cohort | ||
|
|
||
| res = cohort%proot%c%c12 * myinterface%params_species(cohort%species)%SRA | ||
| end function rootarea | ||
|
|
||
| function lai(cohort) result(res) | ||
| ! Leaf area index: surface of leaves per m2 of crown | ||
| real :: res | ||
| type(cohort_type) :: cohort | ||
|
|
||
| res = cohort%leafarea / cohort%crownarea !(cohort%crownarea *(1.0-sp%internal_gap_frac)) | ||
| end function lai | ||
|
|
||
| function basal_area(cohort) result(res) | ||
| ! Tree basal area, m2 tree-1 | ||
| real :: res | ||
| type(cohort_type) :: cohort | ||
|
|
||
| res = pi/4 * cohort%dbh * cohort%dbh | ||
| end function basal_area | ||
|
|
||
| function volume(cohort) result(res) | ||
| ! Tree basal volume, m3 tree-1 | ||
| real :: res | ||
| type(cohort_type) :: cohort | ||
|
|
||
| res = (cohort%psapw%c%c12 + cohort%pwood%c%c12) / myinterface%params_species(cohort%species)%rho_wood | ||
| end function volume | ||
|
|
||
| end module md_cohort |
Oops, something went wrong.