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

Appendix 217
END; (* IF *)
prevDay := currDay;
currDay := prevDay + kDays;
END; (* FOR *)
uAET := uAET * 10.0; (* convert AET & PET from cm to mm *)
uPET := uPET * 10.0;
uDrStr := (uPET - uAET) / uPET;
END SoilMoisture;
(*****************************)
(* Monitoring & statistics *)
(*****************************)
PROCEDURE DeclareStatArrays( arrLen : INTEGER );
BEGIN
uAETStatArray := notExistingStatArray;
DeclStatArray( uAETStatArray, arrLen );
DeclDispMV( uAETStatArray, fe, meanAET, fe, timeIsIndep );
DeclStatArrayForOutput( uAETStatArray, "Actual evapotranspiration", "mm/yr", 0.0 );
uDDStatArray := notExistingStatArray;
DeclStatArray( uDDStatArray, arrLen );
DeclDispMV( uDDStatArray, fe, meanDD, fe, timeIsIndep );
DeclStatArrayForOutput( uDDStatArray, "Degree-days", "°C*days", 0.0 );
uDrStrStatArray := notExistingStatArray;
DeclStatArray( uDrStrStatArray, arrLen );
DeclDispMV( uDrStrStatArray, fe, meanDrStr, fe, timeIsIndep );
DeclStatArrayForOutput( uDrStrStatArray, "Drought stress", "%/100", 0.0 );
uWiTStatArray := notExistingStatArray;
DeclStatArray( uWiTStatArray, arrLen );
DeclDispMV( uWiTStatArray, fe, meanWiT, fe, timeIsIndep );
DeclStatArrayForOutput( uWiTStatArray, "Min. winter temperature", "°C", -100.0 );
END DeclareStatArrays;
PROCEDURE DisplayStatArrays;
BEGIN
IF NOT ExperimentAborted() THEN
DisplayArray( uAETStatArray, TRUE, prob950 );
DisplayArray( uDDStatArray, TRUE, prob950 );
DisplayArray( uDrStrStatArray, TRUE, prob950 );
DisplayArray( uWiTStatArray, 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() 0.0 THEN kPMod := 1.0 + site.kSlAsp*0.125;
ELSE kPMod := 1.0 + site.kSlAsp*0.063;
END;
IF ExperimentRunning() AND (CurrentSimNr() = 1 ) THEN
IF exp.type = manyRuns THEN
GetGlobSimPars(tzero, tend, h, er, c , hm);
arrLen := TRUNC((tend-tzero)/hm+1.5);
DeclareStatArrays( arrLen );
ELSIF exp.type = estimEquil THEN
DeclareStatArrays( 1 );
ELSE
(* do nothing *)
END;
END;
END Initialize;
PROCEDURE StochWeatherOutput;
BEGIN
WITH site DO
CorrWeatherGenerator( rTVect, rPVect );
DegreeDays( uDD, rTVect );
SoilMoisture( uDrStr, uAET, rTVect, rPVect );
uWiT := rTVect[Jan];
IF rTVect[Feb] < uWiT THEN uWiT := rTVect[Feb] END;