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

Appendix 233
BEGIN
WITH sp^ DO
d := Rmax( uDrStress, p.kDrT );
gSMGF := Sqrt( 1.0 - uDrStress/d );
END;
END DroughtGrowthFactor;
PROCEDURE StandCharacteristics( VAR totalTrees, totalBiomass, gLAI: REAL;
VAR cumLA: ARRAY OF REAL; firstSp: SpeciesPtr );
CONST kB1 = 137.0; (* diameter at breast height *)
VAR sp: SpeciesPtr;
cohort: CohortPtr;
gFolA: REAL;
i, iHT: INTEGER;
BEGIN
totalBiomass := 0.0;
totalTrees := 0.0;
(* initialize canopy leaf biomass profile *)
FOR i:= 0 TO maxHeight-1 DO cumLA[i] := 0.0 END;
(* calculate stemwood & foliage biomass and leaf area *)
sp := firstSp;
WHILE sp NIL DO
WITH sp^ DO
biomass := 0.0;
nrTrees := 0.0;
cohort := firstCohort;
WHILE cohort NIL DO
nrTrees := nrTrees + FLOAT(cohort^.nrTrees);
WITH cohort^ DO
(* dry stemwood biomass after Sollins et al. (1973), Burger (1945-53),
Woods et al. (1991) *)
gSBio := 0.12 * Exp( 2.4*Ln(D) );
(* dry foliage biomass, leaf area and tree height after Burger (1945-53) *)
gFolW := p.kA1 * Exp( p.kA2*Ln(D) ) * p.kC1;
gFolA := p.kC2 * gFolW / p.kC1;
gH := kB1 + p.kB2*D - p.kB3*D*D;
(* sum leaf area for all trees of similar height *)
iHT := TRUNC( gH - kB1 ) DIV 10;
cumLA[iHT] := cumLA[iHT] + gFolA*FLOAT(nrTrees);
(* calculate species-specific biomass *)
biomass := biomass + (gSBio+gFolW)*FLOAT(nrTrees);
END; (* WITH *)
cohort := cohort^.next;
END; (* WHILE *)
totalBiomass := totalBiomass + biomass;
totalTrees := totalTrees + nrTrees;
sp := sp^.next;
END; (* WITH *)
END; (* WHILE *)
(* calculate cumulative leaf area index and leaf area index at forest floor*)
FOR i:= maxHeight - 2 TO 0 BY -1 DO
cumLA[i] := cumLA[i] + cumLA [i+1];
cumLA[i+1] := cumLA[i+1] / kPatchSize;
END; (* FOR *)
cumLA[0] := cumLA[0] / kPatchSize;
gLAI := cumLA[0];
END StandCharacteristics;
BEGIN
(* initialize soil nitrogen multiplier parameters (Aber et al. 1979 as modified
by Pastor & Post 1985) *)
kN1[1] := -0.016; kN2[1] := 2.245;
kN1[2] := -0.022; kN2[2] := 30.605;
kN1[3] := -0.016; kN2[3] := 43.973;
END FCPGrFact.
Definition module FCPMon
DEFINITION MODULE FCPMon;
(****************************************
Module FCPMon (Version 2.4)
Copyright ©1994 by Harald Bugmann and Swiss
Federal Institute of Technology Zürich ETHZ
Version written for:
'Dialog Machine' DM_V2.2 (User interface)
MacMETH_V3.2.1 (1-Pass Modula-2 implementation)
ModelWorks_V2.2 (Modelling & Simulation)
Purpose Monitoring of tree-ring chronologies and of tree frequency
distributions in the ForClim-P model; animation of tree growth
Programming
• Design
H. Bugmann 19.4.1991
• Implementation
H. Bugmann 19.4.1991
Swiss Federal Institute of Technology Zurich ETHZ