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

Appendix 237
Purpose Simulation model for soil organic mass & nitrogen dynamics
Remarks Based on the model by Pastor & Post (1985, 1986)
Implementation and Revisions:
============================
Author Date Description
------ ---- -----------
hb 5. 2.1992 First implementation (DM 2.2, MacMETH 3.2)
hb 7. 2.1992 Major bugs fixed; rewritten with mass as state var
hb 9. 2.1992 Rewritten (dynamic lists)
hb 19. 2.1992 Implementation finished
hb 23. 3.1992 Model renamed to FORCLIM-S
hb 23.12.1992 adapted for 6 litter types (usage within ForClim)
reverted to mass dynamics instead of mass
hb 5. 3.1993 adapted for output of statistical data
****************************************)
FROM SimMaster
FROM SimBase
FROM DMSystem
FROM DMMenus
IMPORT CurrentSimNr, ExperimentRunning, ExperimentAborted;
IMPORT DeclM, IntegrationMethod, DeclSV, RTCType, DeclP,
StashFiling, DeclMV, Tabulation, NoAbout, Graphing,
GetGlobSimPars, RemoveM, MDeclared, CurrentTime;
IMPORT CurrentDMLevel, InstallTermProc;
IMPORT Command, AccessStatus, Marking, Separator, InstallCommand,
DisableCommand, EnableCommand, InstallSeparator,
RemoveSeparatorAtCommand, RemoveCommand, SeparatorPosition;
FROM DMEntryForms IMPORT FormFrame, WriteLabel, RadioButtonID,
DefineRadioButtonSet, RadioButton, CheckBox, UseEntryForm;
FROM DMMessages IMPORT Warn, Abort;
FROM DMStorage IMPORT Allocate, Deallocate;
FROM ForestBase IMPORT fs, Litter, uLitt, uAvN, uAET, DeclConstSoil,
RemoveConstSoil, fMenu, exp, ExperimentType,
DeclStatArrayForOutput, DeclMonitoringProc,
RemoveMonitoringProc, Rmax, Rmin;
FROM StochStat IMPORT StatArray, Prob2Tail, DeclStatArray, notExistingStatArray,
DeclDispMV, DisplayArray, PutValue;
FROM SimGraphUtils IMPORT timeIsIndep;
CONST
modIdent = "ForClim-S";
modDescr = "ForClim-S: Soil C/N turnover model";
TYPE
LitterPtr = POINTER TO LitterNode;
LitterNode = RECORD
LOM, (* litter organic matter *)
LOMNew,
LOMinit, (* initial LOM *)
LN, (* litter nitrogen *)
LNNew,
gLign : REAL; (* lignin content *)
type: Litter; (* litter type *)
next: LitterPtr; (* pointer for list management *)
END;
VAR
startupLevel: CARDINAL;
first, prev: LitterPtr;
type: Litter;
statArrAvN, statArrLitM, statArrHumM : StatArray;
kInitN, kCritN: ARRAY [ MIN(Litter) .. MAX(Litter) ] OF REAL;
kNC, kLignA, kLignB,
gNMR, decMlt, kAET, kMin,
k1, k2, k3, k4, k5, k6, mLoss, kLeach, litterM, litterN, totSOM,
litterCO2, gLNC, gAETM, humusNMin, HOM, HOMNew, humCN,
HN, HNNew, gImmob, totCO2, totNimmob, totNMin, kNAtm,
meanAvN, meanLitM, meanHumM,
tzero, hm: REAL;
ok, allTypes, everyYear: BOOLEAN;
fcsCmd, litterCmd: Command;
PROCEDURE EmptyProc; BEGIN END EmptyProc;
(************************************)
(* Procedures for list management *)
(************************************)
PROCEDURE CreateLitterCohort( VAR first: LitterPtr; m: REAL;
type: Litter );
VAR litter: LitterPtr;
BEGIN
Allocate( litter, SIZE(LitterNode) );
IF litter = NIL THEN Abort( "", "Insufficient memory!", "" ) END;
litter^.next := first;
litter^.LOM := m;
litter^.LOMNew := m;
litter^.LOMinit := m;