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

Appendix 219
PROCEDURE InitializeFWM;
BEGIN
DisableCommand( fMenu, chooseEnvCmd );
currWater := site.kFC; (* initial conditions in January of the first year *)
IF site.kSlAsp > 0.0 THEN kPMod := 1.0 + site.kSlAsp*0.125;
ELSE kPMod := 1.0 + site.kSlAsp*0.063;
END;
Lookup( inF, fileName, FALSE );
SkipFileHeader;
END InitializeFWM;
PROCEDURE FileWeatherOutput;
VAR t: REAL;
i: Month;
PROCEDURE GetWeatherDataFromFile;
BEGIN
GetReal( inF, t );
FOR i:= Jan TO Dec DO
GetReal( inF, rTVect[i] );
END;
FOR i:= Jan TO Dec DO
GetReal( inF, rPVect[i] );
END;
SkipGap( inF );
IF EOF(inF) THEN
Reset( inF );
SkipFileHeader;
END; (* IF *)
END GetWeatherDataFromFile;
BEGIN
GetWeatherDataFromFile;
DegreeDays( uDD, rTVect );
SoilMoisture( uDrStr, uAET, rTVect, rPVect );
uWiT := rTVect[Jan];
IF rTVect[Feb] < uWiT THEN uWiT := rTVect[Feb] END;
IF rTVect[Dec] < uWiT THEN uWiT := rTVect[Dec] END;
END FileWeatherOutput;
PROCEDURE TerminateFWM;
BEGIN
Close( inF );
EnableCommand( fMenu, chooseEnvCmd );
END TerminateFWM;
PROCEDURE DeclFileWeatherObjects;
BEGIN
RemoveAllWeatherTabFuncs;
DeclareDDCorrectionTabF;
DeclMV( uDD, 0.0, 2500.0, "Degree-days", "uDD" , "d*°C",
notOnFile, writeInTable, notInGraph);
DeclMV( uDrStr, 0.0, 1.0, "Drought stress", "uDrStr" , "-",
notOnFile, writeInTable, notInGraph);
DeclMV( uWiT, -30.0, 30.0, "Winter temperature", "uWiT" , "°C",
notOnFile, writeInTable, notInGraph);
DeclSoilMoistureParameters;
declEMod := fileW;
END DeclFileWeatherObjects;
PROCEDURE DeclFileWeatherModel;
BEGIN
IF declEMod = stochE THEN
RemoveMonitoringProc( Monitoring );
END;
DeclM(fe, discreteTime, InitializeFWM, NoInput, FileWeatherOutput, NoDynamic,
TerminateFWM, DeclFileWeatherObjects, modDescrFileW, modIdent, NoAbout);
END DeclFileWeatherModel;
(* Bioclimatic data from file *)
(********************************)
PROCEDURE TestTheBioclimaticFile( VAR allOK: BOOLEAN );
VAR x: REAL;
i: INTEGER;
PROCEDURE TestNum;
BEGIN
IF NOT legalNum AND allOK THEN
allOK := FALSE;
Warn( "The file you specified contains illegal numbers!",
"The stochastic environment model will be declared instead", "" );
END;
END TestNum;
BEGIN
allOK := (inF.res = done);
IF allOK THEN
SkipFileHeader;
WHILE NOT EOF(inF) DO
FOR i:= 1 TO 4 DO GetReal( inF, x ); TestNum END;
SkipGap( inF );
END; (* WHILE *)
ELSE
Warn( "The file you specified could not be found or opened!",
"The stochastic environment model will be declared instead", "" );
END; (* IF *)
END TestTheBioclimaticFile;