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

232 Appendix
Module FCPGrFact
DEFINITION MODULE FCPGrFact;
(****************************************
Module FCPGrFact (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 Provides growth factors for the ForClim-P model
Programming
• Design
H. Bugmann 18.1.1991
• Implementation
H. Bugmann 18.1.1991
Swiss Federal Institute of Technology Zurich ETHZ
CH-8092 Zurich
Switzerland
Last revision of definition: 14.1.1993 hb
****************************************)
FROM FCPBase IMPORT SpeciesPtr;
CONST maxHeight = 600; (* in units of 10 cm, i.e. 60 m *)
PROCEDURE SoilNitrogenGrowthFactor( VAR sp: SpeciesPtr; uAvN: REAL );
PROCEDURE DegreeDayGrowthFactor( VAR sp: SpeciesPtr; uDD: REAL );
PROCEDURE DroughtGrowthFactor( VAR sp: SpeciesPtr; uDrStr: REAL );
PROCEDURE StandCharacteristics( VAR totalTrees, totalBiomass, gLAI: REAL;
VAR cumLA: ARRAY OF REAL; firstSp: SpeciesPtr );
END FCPGrFact.
IMPLEMENTATION MODULE FCPGrFact;
(*
*)
Implementation and Revisions:
============================
Author Date Description
------ ---- -----------
hb 21. 3.1991 First implementation (V0.1, DM 2.02, MacMETH 2.6.2)
hb 13. 8.1992 Adaptation for V2.0 (growth factors removed)
hb 6.11.1992 Calculation of leaf weight changed (Burger data)
hb 14. 1.1993 adapted for SNGF
FROM ForestBase
FROM FCPBase
FROM SYSTEM
IMPORT kPatchSize, Rmax;
IMPORT SpeciesPtr, CohortPtr;
IMPORT Exp, Ln, Sqrt;
VAR kN1, kN2: ARRAY [1..3] OF REAL;
PROCEDURE SoilNitrogenGrowthFactor( VAR sp: SpeciesPtr; uAvN: REAL );
(* Mitchell & Chandler (1939), Aber et al. (1979), Pastor & Post (1985) *)
VAR tol: INTEGER;
BEGIN
tol := sp^.p.kNTol;
sp^.gSNGF := Rmax( 1.0 - Exp( kN1[tol]*(uAvN-kN2[tol]) ) , 0.0 );
END SoilNitrogenGrowthFactor;
PROCEDURE DegreeDayGrowthFactor( VAR sp: SpeciesPtr; uDD: REAL );
(* Hellmers (1962), Botkin et al. (1972), Shugart (1984), Woodward (1987) *)
BEGIN
WITH sp^ DO
gDDGF := 4.0*(uDD-p.kDDMin)*(p.kDDMax-uDD)
/ ((p.kDDMax-p.kDDMin)*(p.kDDMax-p.kDDMin));
gDDGF := Rmax( 0.0, gDDGF );
END;
END DegreeDayGrowthFactor;
PROCEDURE DroughtGrowthFactor( VAR sp: SpeciesPtr; uDrStress: REAL );
(* Bassett (1964), Prentice & Helmisaari (1991) *)
VAR d: REAL;