On the Ecology of Mountainous Forests in a Changing Climate: A ...

Appendix 229
Swiss Federal Institute of Technology Zurich ETHZ
CH-8092 Zurich
Switzerland
Last revision of definition: 2.3.1992 hb
****************************************)
FROM FCPBase IMPORT SpeciesPtr;
PROCEDURE DeclSubMParameters;
(* declares ForClim-P parameters *)
PROCEDURE LitterInitialization( t0, kPatchSize: REAL );
(* initializes litter production and assigns accumulated litter to the
corresponding ForestBase variables;
t0 is simulation start time (used for initialization);
LitterInitialization is to be called in the Output proc
*)
PROCEDURE TreeDeathAndGrowth( VAR firstSp: SpeciesPtr; VAR cumLA: ARRAY OF REAL;
gLAI: REAL );
PROCEDURE TreeEstablishment( VAR firstSp: SpeciesPtr;
uWiT: REAL; uDD, gLAI, totalTrees, kPatchSize: REAL );
PROCEDURE UpdateGap( VAR firstSp: SpeciesPtr; uDD, uDrStr: REAL );
END FCPDynamic.
IMPLEMENTATION MODULE FCPDynamic;
(*
Implementation and Revisions:
============================
Author Date Description
------ ---- -----------
hb 18. 1.1991 First implementation (DM 2.02, MacMETH 2.6.2)
hb 30. 1.1991 PROCs TreeDeath and TreeGrowth merged (V0.1)
hb 26. 2.1991 Version 0.3 implemented
hb 21. 3.1991 Version 0.4 implemented
hb 28. 1.1992 More efficient establishment at the beginning of the
simulation introduced
hb 30. 1.1992 List management changed in TreeDeathAndGrowth
hb 11. 3.1992 Formulation of browsing changed
hb 13. 8.1992 New formulation for light growth factor (kShaw = [1…9])
hb 10.11.1992 Adaptation of limiting factor statistics to new equations
hb 1. 3.1993 Adaptation for new gH-diameter relationship in maximum
growth equation
hb 8. 6.1993 calculation of litter output fixed (twigs were too high)
*)
FROM ForestBase IMPORT Power, Rmax, uBrPr, Litter, uLitt;
FROM FCPBase IMPORT SpeciesPtr, CohortPtr, CreateCohort, DeleteCohort,
MergeCohorts, SpeciesType;
FROM FCPMon IMPORT WriteTreeRings;
FROM RandGen IMPORT U;
FROM SYSTEM IMPORT Exp;
FROM SimBase IMPORT DeclP, RTCType, CurrentTime, GetGlobSimPars;
VAR kInitDBH, kMinAbsInc, kMinRelInc, kLAtt, kDeathP, kSlowGrP, kSlowGrYrs,
kEstP, kEstNr, kTwig, kAFW, kRSR: REAL;
littF: ARRAY [1..3] OF REAL;
littW, littT, littR: REAL;
kFRT: ARRAY [deciduous..coniferous] OF REAL;
PROCEDURE DeclSubMParameters;
BEGIN
DeclP(kMinAbsInc, 0.03, 0.0, 1.0, noRtc, 'Min. abs. growth (vigorous tree)', 'kMinAbsInc', 'cm');
DeclP(kMinRelInc, 0.1, 0.0, 1.0, noRtc, 'Min. rel. growth (vigorous tree)', 'kMinRelInc', '%');
DeclP(kLAtt, 0.25, 0.0, 1.0, noRtc, 'Light attenuation coefficient', 'kLAtt', '--');
DeclP(kDeathP, 4.605, 0.0, 10.0, noRtc, 'Death probability coefficient', 'kDeathP', '--');
DeclP(kSlowGrP, 0.368, 0.0, 1.0, noRtc, 'Slow growth enhanced mortality', 'kSlowGrP', '--');
DeclP(kSlowGrYrs, 2.0, 0.0, 10.0, noRtc, 'Nr of slow growth years required', 'kSlowGrYrs', '#');
DeclP(kEstP, 0.1, 0.0, 1.0, noRtc, 'Probability of tree establishment', 'kEstP', '--');
DeclP(kEstNr, 0.006, 0.0, 1.0, noRtc, 'Max. establishment per species', 'kEstNr', '#/m2*yr');
DeclP(kInitDBH, 1.27, 0.0, 10.0, noRtc, 'DBH of new trees', 'kInitDBH', 'cm');
DeclP(kFRT[coniferous], kFRT[coniferous], 0.0, 10.0, noRtc, 'Foliage retention time (conifers)', 'kFRT', 'years');
DeclP(kTwig, 0.0025, 0.0, 0.1, noRtc, 'Twig litter production parameter', 'kTwig', 'kg/cm2');
DeclP(kAFW, 0.92, 0.5, 1.0, noRtc, 'Ash-free weight of litter', 'kAFW', '--');
DeclP(kRSR, 4.0, 0.0, 10.0, noRtc, 'Root:shoot ratio of litter', 'kRSR', '--');
END DeclSubMParameters;
PROCEDURE LitterInitialization( t0, kPatchSize: REAL );
VAR conv: REAL;
BEGIN
IF CurrentTime() = t0 THEN (* initialize variables *)
littW := 0.0; littT := 0.0; littR := 0.0;
littF[1] := 0.0; littF[2] := 0.0; littF[3] := 0.0;
END;
(* assign ForestBase variables, and reset internal variables to zero *)
conv := 10.0 / kPatchSize; (* conversion factor kg/gap -> t/ha *)
uLitt[leafFast] := littF[1]*conv; littF[1] := 0.0;
uLitt[leafMedium] := littF[2]*conv; littF[2] := 0.0;
uLitt[leafSlow] := littF[3]*conv; littF[3] := 0.0;
uLitt[wood] := littW*conv; littW := 0.0;
uLitt[twigs] := littT*conv; littT := 0.0;