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

224 Appendix
sp := firstSp;
WHILE sp NIL DO (* declare numbers *)
sp^.statArrN := notExistingStatArray;
DeclStatArray( sp^.statArrN, arrLen );
DeclDispMV( sp^.statArrN, fp, sp^.meanNr, fp, timeIsIndep );
IF withNr THEN
DeclStatArrayForOutput( sp^.statArrN, sp^.longDescNr, "#/ha", minValNr );
END;
sp := sp^.next;
END;
statArrLAI := notExistingStatArray;
DeclStatArray( statArrLAI, arrLen );
DeclDispMV( statArrLAI, fp, meanLAI, fp, timeIsIndep );
DeclStatArrayForOutput( statArrLAI, "Leaf area index", "m^2/m^2", 0.0 );
statArrTotBio := notExistingStatArray;
DeclStatArray( statArrTotBio, arrLen );
DeclDispMV( statArrTotBio, fp, meanTotBio, fp, timeIsIndep );
DeclStatArrayForOutput( statArrTotBio, "Total biomass", "t/ha", 0.0 );
statArrTotNr := notExistingStatArray;
DeclStatArray( statArrTotNr, arrLen );
DeclDispMV( statArrTotNr, fp, meanTotNr, fp, timeIsIndep );
DeclStatArrayForOutput( statArrTotNr, "Total number", "#/ha", 0.0 );
FOR i:= MIN(Litter) TO MAX(Litter) DO
statArrLitt[i] := notExistingStatArray;
DeclStatArray( statArrLitt[i], arrLen );
DeclDispMV( statArrLitt[i], fp, meanLitt[i], fp, timeIsIndep );
END;
DeclStatArrayForOutput( statArrLitt[leafFast], "Litterfall (foliage fast)", "t/ha", 0.0 );
DeclStatArrayForOutput( statArrLitt[leafMedium], "Litterfall (foliage medium)", "t/ha", 0.0 );
DeclStatArrayForOutput( statArrLitt[leafSlow], "Litterfall (foliage slow)", "t/ha", 0.0 );
DeclStatArrayForOutput( statArrLitt[twigs], "Litterfall (twigs)", "t/ha", 0.0 );
DeclStatArrayForOutput( statArrLitt[roots], "Litterfall (fine roots)", "t/ha", 0.0 );
DeclStatArrayForOutput( statArrLitt[wood], "Litterfall (wood)", "t/ha", 0.0 );
END DeclareStatArrays;
PROCEDURE DisplayStatArrays;
VAR i: Litter;
BEGIN
IF NOT ExperimentAborted() THEN
sp := firstSp;
WHILE sp NIL DO
DisplayArray( sp^.statArrB, TRUE, prob950 );
DisplayArray( sp^.statArrN, TRUE, prob950 );
sp := sp^.next;
END;
DisplayArray( statArrLAI, TRUE, prob950 );
DisplayArray( statArrTotBio, TRUE, prob950 );
DisplayArray( statArrTotNr, TRUE, prob950 );
FOR i:= MIN(Litter) TO MAX(Litter) DO
DisplayArray( statArrLitt[i],TRUE, prob950 );
END;
END;
END DisplayStatArrays;
(***********************************)
(* Procedures for model dynamics *)
(***********************************)
PROCEDURE Initialize;
VAR j, x0New, y0New, z0New, arrLen: INTEGER;
tend, c, h, er: REAL;
ok: BOOLEAN;
(*
PROCEDURE SetImmigPar( desc: ARRAY OF CHAR; kImmYrNew: REAL );
VAR found: BOOLEAN;
BEGIN
found := FALSE;
sp := firstSp;
WHILE (sp NIL) AND NOT found DO
IF CompareStrings( sp^.shortDescBio, desc ) = equal THEN
found := TRUE;
sp^.p.kImmYr := kImmYrNew;
END;
sp := sp^.next;
END; (* WHILE *)
IF NOT found THEN HALT END;
END SetImmigPar;
*)
PROCEDURE SetSpeciesParameter( VAR sensSp, modelSp: SpeciesPtr; nrRun: INTEGER );
VAR par: INTEGER;
BEGIN
par := nrRun MOD 14;
ResetAllSpeciesParameters;
IF (par = 1) AND (nrRun 1) THEN (* move to a new species *)
sensSp := sensSp^.next;
modelSp := modelSp^.next;
END;
(* now assign a new species parameter *)
IF par = 1 THEN
modelSp^.p.kA1 := sensSp^.p.kA1;
modelSp^.p.kA2 := sensSp^.p.kA2;
modelSp^.p.kC1 := sensSp^.p.kC1;
modelSp^.p.kC2 := sensSp^.p.kC2;
ELSIF par = 2 THEN
modelSp^.p.kDm := sensSp^.p.kDm;
modelSp^.p.kB2 := 2.0*(modelSp^.p.kHm - 137.0) / modelSp^.p.kDm;
modelSp^.p.kB3 := modelSp^.p.kB2 / 2.0 / modelSp^.p.kDm;
ELSIF par = 3 THEN