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

218 Appendix
IF rTVect[Dec] < uWiT THEN uWiT := rTVect[Dec] END;
END;
END StochWeatherOutput;
PROCEDURE Terminate;
BEGIN
IF ( (CurrentSimNr() MOD TRUNC(exp.nrRuns+0.5)) = 0 )
AND ExperimentRunning() AND (exp.type = manyRuns) THEN
DisplayStatArrays;
END;
EnableCommand( fMenu, chooseEnvCmd );
END Terminate;
PROCEDURE DeclareDDCorrectionTabF;
VAR T, corrDD: ARRAY [0..11] OF REAL;
BEGIN
T[0] := -20.0; corrDD[0] := 0.0;
T[1] := -10.0; corrDD[1] := 0.0;
T[2] := -5.0; corrDD[2] := 1.28;
T[3] := -2.5; corrDD[3] := 3.3;
T[4] := 0.0; corrDD[4] := 8.52;
T[5] := 2.5; corrDD[5] := 22.02;
T[6] := 5.24; corrDD[6] := 62.56;
T[7] := 7.5; corrDD[7] := 39.0;
T[8] := 10.0; corrDD[8] := 23.12;
T[9] := 12.5; corrDD[9] := 13.71;
T[10]:= 16.11; corrDD[10]:= 6.45;
T[11]:= 17.5; corrDD[11]:= 9.76;
DeclTabF( corrDDTab, T, corrDD, 12, FALSE, "DD correction",
"T", "corr", "T", "corr", -20.0, 30.0, 0.0, 70.0);
END DeclareDDCorrectionTabF;
PROCEDURE DeclStochWeatherObjects;
VAR i: INTEGER;
BEGIN
DeclMV( uAET, 0.0, 1000.0, "Actual evapotranspiration", "uAET" , "mm/y",
notOnFile, writeInTable, notInGraph);
DeclMV( uPET, 0.0, 1000.0, "Potential evapotranspiration", "uPET" , "mm/y",
notOnFile, notInTable, notInGraph);
DeclMV( uDD, 0.0, 3000.0, "Degree-days", "uDD" , "d*°C",
notOnFile, writeInTable, notInGraph);
DeclMV( uWiT, -30.0, 30.0, "Winter temperature", "uWiT" , "°C",
notOnFile, writeInTable, notInGraph);
DeclMV( uDrStr, 0.0, 1.0, "Drought stress", "uDrStr" , "--",
notOnFile, writeInTable, notInGraph);
DeclMV( meanAET, 0.0, 1000.0, "Average AET", "meanAET" , "mm/y",
notOnFile, notInTable, notInGraph);
DeclMV( meanDD, 0.0, 3000.0, "Average degree-days", "meanDD" , "d*°C",
notOnFile, notInTable, notInGraph);
DeclMV( meanWiT, 0.0, 1.0, "Average winter temperature", "meanWiT" , "°C",
notOnFile, notInTable, notInGraph);
DeclMV( meanDrStr, 0.0, 1.0, "Average drought stress", "meanDrStr" , "--",
notOnFile, notInTable, notInGraph);
DeclSoilMoistureParameters;
RemoveAllWeatherTabFuncs;
DeclareDDCorrectionTabF;
FOR i:= 1 TO 20 DO time[i] := FLOAT(i)*50.0; modifier[i] := 0.0 END;
DeclTabF( tSuTab, time, modifier, 20, TRUE, "Summer T",
"Year", "∆ (°C)", "year", "°C", 0.0, 5000.0, -10.0, 10.0);
DeclTabF( tWiTab, time, modifier, 20, TRUE, "Winter T",
"Year", "∆ (°C)", "year", "°C", 0.0, 5000.0, -10.0, 10.0);
DeclTabF( pSuTab, time, modifier, 20, TRUE, "Summer P",
"Year", "∆ (cm/month)", "year", "cm/month", 0.0, 5000.0, -5.0, 5.0);
DeclTabF( pWiTab, time, modifier, 20, TRUE, "Winter P",
"Year", "∆ (cm/month)", "year", "cm/month", 0.0, 5000.0, -5.0, 5.0);
FOR i:= 1 TO 20 DO modifier[i] := 1.0 END;
DeclTabF( tSDSuTab, time, modifier, 20, TRUE, "Summer SD(T)",
"Year", "Mult", "year", "--", 0.0, 5000.0, 0.5, 2.0);
DeclTabF( tSDWiTab, time, modifier, 20, TRUE, "Winter SD(T)",
"Year", "Mult", "year", "--", 0.0, 5000.0, 0.5, 2.0);
DeclTabF( pSDSuTab, time, modifier, 20, TRUE, "Summer SD(P)",
"Year", "Mult", "year", "--", 0.0, 5000.0, 0.5, 2.0);
DeclTabF( pSDWiTab, time, modifier, 20, TRUE, "Winter SD(P)",
"Year", "Mult", "year", "--", 0.0, 5000.0, 0.5, 2.0);
declEMod := stochE;
END DeclStochWeatherObjects;
PROCEDURE DeclStochEnvModel;
BEGIN
DeclM(fe, discreteTime, Initialize, NoInput, StochWeatherOutput, NoDynamic, Terminate,
DeclStochWeatherObjects, modDescrStoch, modIdent, NoAbout);
DeclMonitoringProc( Monitoring );
END DeclStochEnvModel;
PROCEDURE SkipFileHeader;
VAR ch: CHAR;
BEGIN
ReadChar( inF, ch );
WHILE NOT EOF(inF) AND (chEOL) DO ReadChar( inF, ch ) END;
ReadChar( inF, ch );
WHILE NOT EOF(inF) AND (chEOL) DO ReadChar( inF, ch ) END;
END SkipFileHeader;
(* Weather data from file *)
(****************************)