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

Appendix 225
modelSp^.p.kHm := sensSp^.p.kHm;
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 = 4 THEN
modelSp^.p.kAm := sensSp^.p.kAm;
ELSIF par = 5 THEN
modelSp^.p.kG := sensSp^.p.kG;
ELSIF par = 6 THEN
modelSp^.p.kDDMin := sensSp^.p.kDDMin;
ELSIF par = 7 THEN
modelSp^.p.kDDMax := sensSp^.p.kDDMax;
ELSIF par = 8 THEN
modelSp^.p.kWiT := sensSp^.p.kWiT;
ELSIF par = 9 THEN
modelSp^.p.kDrT := sensSp^.p.kDrT;
ELSIF par = 10 THEN
modelSp^.p.kNTol := sensSp^.p.kNTol;
ELSIF par = 11 THEN
modelSp^.p.kBrow := sensSp^.p.kBrow;
ELSIF par = 12 THEN
modelSp^.p.kLy := sensSp^.p.kLy;
ELSIF par = 13 THEN
modelSp^.p.kLa := sensSp^.p.kLa;
ELSE (* par 14 *)
modelSp^.p.kLQ := sensSp^.p.kLQ;
END;
END SetSpeciesParameter;
BEGIN
DisableCommand( fMenu, statsCmd );
DisableCommand( fMenu, newSpeciesCmd );
DisableCommand( fMenu, chooseSiteCmd );
DisableCommand( fMenu, chooseWeatherCmd );
DisableCommand( fMenu, treeRingCmd );
AssignAllometricParameters;
IF CurrentSimNr() = 1 THEN
sp := firstSp;
WHILE sp NIL DO
FOR j := 1 TO 7 DO sp^.limCounter[j] := 0 END;
sp^.sumALGF := 0.0; sp^.sumDDGF := 0.0;
sp^.sumSMGF := 0.0; sp^.sumSNGF := 0.0; sp^.sumGF := 0.0;
sp := sp^.next;
END; (* WHILE *)
END; (* IF *)
GetGlobSimPars(tzero, tend, h, er, c , hm);
IF ExperimentRunning() AND (CurrentSimNr() = 1 ) THEN
IF exp.type = manyRuns THEN
arrLen := TRUNC((tend-tzero)/hm+1.5);
DeclareStatArrays( arrLen );
ELSIF exp.type = estimEquil THEN
DeclareStatArrays( 1 );
(* for sensitivity analysis of species parameters only: *)
IF exp.doSensAnalysis THEN
(*. ReadSpeciesFile( sensFirstSp, "", ok ); .*)
sensSp := sensFirstSp;
modelSp := firstSp;
END;
ELSE
(* do nothing *)
END;
END;
IF ExperimentRunning() AND (exp.type = estimEquil)
AND exp.doSensAnalysis THEN (* experiment for parameter sensitivity *)
(* execute this procedure only if the experiment is for the analysis
of parameter sensitivity of tree species! *)
(*. SetSpeciesParameter( sensSp, modelSp, CurrentSimNr() ); .*)
END;
IF exp.readState THEN
AssignString( exp.inFileName, pathAndFilename );
CreateStateFileName( fn );
Concat( pathAndFilename, fn );
ReadStateFromFile( ok, x0New, y0New, z0New, pathAndFilename );
IF ok AND (x0New undefSeed) THEN
SetRandomNumberSeeds( x0New, y0New, z0New );
END;
END;
END Initialize;
PROCEDURE Output;
BEGIN
LitterInitialization( tzero, kPatchSize );
UpdateGap( firstSp, uDD, uDrStr );
END Output;
PROCEDURE Input;
BEGIN
StandCharacteristics( totalTrees, totalBiomass, gLAI, cumLA, firstSp );
sp := firstSp;
WHILE sp NIL DO
DegreeDayGrowthFactor( sp, uDD );
DroughtGrowthFactor( sp, uDrStr);
SoilNitrogenGrowthFactor( sp, uAvN );
sp := sp^.next;
END; (* WHILE *)
END Input;
PROCEDURE Dynamic;
BEGIN
IF U() < kPDist THEN (* catastrophic disturbance occurs, all trees die *)
sp := firstSp;
WHILE sp NIL DO