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

Appendix 239
END Immobilization;
PROCEDURE Mineralization;
VAR humMDelta: REAL;
BEGIN
IF HOM 0.0 THEN
gAETM := Rmin( uAET / (kAET-uAET), 1.0 );
IF gLNC 0.0 THEN
(* new formulation based on data in Pastor et al. (1984), avoids occurrence
of mineralization values of ±∞ (pole in Pastor & Post formulation *)
humusNMin := Rmax( k5 + k6/gLNC, 0.0005 ) * decMlt * gAETM * HOM;
ELSE
humusNMin := HN * kMin * decMlt * gAETM;
END; (* IF *)
HNNew := HN - humusNMin;
humMDelta := HOM * humusNMin / HN;
IF HOM < humMDelta THEN HALT END;
HOMNew := HOM - humMDelta;
totCO2 := totCO2 + humMDelta;
totNMin := totNMin + humusNMin;
ELSE
HOMNew := 0.0;
HNNew := 0.0;
END; (* IF *)
END Mineralization;
(*****************************)
(* Monitoring & statistics *)
(*****************************)
PROCEDURE DeclareStatArrays( arrLen : INTEGER );
BEGIN
statArrAvN := notExistingStatArray;
DeclStatArray( statArrAvN, arrLen );
DeclDispMV( statArrAvN, fs, meanAvN, fs, timeIsIndep );
DeclStatArrayForOutput( statArrAvN, "Available nitrogen", "kg/ha", 0.0 );
statArrLitM := notExistingStatArray;
DeclStatArray( statArrLitM, arrLen );
DeclDispMV( statArrLitM, fs, meanLitM, fs, timeIsIndep );
DeclStatArrayForOutput( statArrLitM, "Litter mass", "t/ha", 0.0 );
statArrHumM := notExistingStatArray;
DeclStatArray( statArrHumM, arrLen );
DeclDispMV( statArrHumM, fs, meanHumM, fs, timeIsIndep );
DeclStatArrayForOutput( statArrHumM, "Humus mass", "t/ha", 0.0 );
END DeclareStatArrays;
PROCEDURE DisplayStatArrays;
BEGIN
IF NOT ExperimentAborted() THEN
DisplayArray( statArrAvN, TRUE, prob950 );
DisplayArray( statArrLitM, TRUE, prob950 );
DisplayArray( statArrHumM, TRUE, prob950 );
END;
END DisplayStatArrays;
PROCEDURE Monitoring;
VAR t: REAL;
index: INTEGER;
BEGIN
IF ExperimentRunning() THEN
IF exp.type = estimEquil THEN
t := 0.0;
index := 1;
ELSE
t := CurrentTime();
index := TRUNC( (t-tzero)/hm + 0.5 ) + 1;
END;
IF NOT ((exp.type = estimEquil) AND (CurrentTime()