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

Appendix 241
CONST lem = 3;
VAR ef: FormFrame;
ok: BOOLEAN;
cl: INTEGER;
button: ARRAY [MIN(Litter)..MAX(Litter)] OF RadioButtonID;
allTypesB, selectedbutton: RadioButtonID;
BEGIN
cl := 2;
WriteLabel(cl, lem-2, "Select litter type :"); INC(cl);
DefineRadioButtonSet(selectedbutton);
IF NOT allTypes THEN selectedbutton := button[type] END;
RadioButton( allTypesB, cl, lem, "All types"); INC(cl,2);
RadioButton( button[leafFast], cl, lem, "Fast decaying foliage"); INC(cl);
RadioButton( button[leafMedium], cl, lem, "Medium decaying foliage"); INC(cl);
RadioButton( button[leafSlow], cl, lem, "Slowly decaying foliage"); INC(cl);
RadioButton( button[roots], cl, lem, "Roots"); INC(cl);
RadioButton( button[twigs], cl, lem, "Twigs"); INC(cl);
RadioButton( button[wood], cl, lem, "Wood"); INC(cl);
INC(cl);
CheckBox( cl, lem, "Every year", everyYear );
INC(cl);
ef.x := 0; ef.y := -1; (* display entry form in middle of screen *)
ef.lines := cl+1; ef.columns := 30;
UseEntryForm(ef,ok);
IF ok THEN
IF selectedbutton = allTypesB THEN allTypes := TRUE;
ELSE allTypes := FALSE;
END;
IF selectedbutton = button[leafFast] THEN type := leafFast;
ELSIF selectedbutton = button[leafMedium] THEN type := leafMedium;
ELSIF selectedbutton = button[leafSlow] THEN type := leafSlow;
ELSIF selectedbutton = button[roots] THEN type := roots;
ELSIF selectedbutton = button[twigs] THEN type := twigs;
ELSIF selectedbutton = button[wood]
END; (* IF *)
END; (* IF *)
END LitterInput;
THEN type := wood;
(********************************)
(* Procedure ModelObjects *)
(********************************)
PROCEDURE ModelObjects;
VAR i: Litter;
BEGIN
DeclSV( HOM, HOMNew, 0.0, 0.0, 100.0,
'Humus organic matter', 'HOM', 't/ha');
DeclSV( HN, HNNew, 0.0, 0.0, 10.0,
'Humus nitrogen', 'HN', 't/ha');
DeclMV( uAvN, 0.0, 500.0, 'Available nitrogen','uAvN','kg/ha',
notOnFile,writeInTable,notInGraph);
DeclMV( litterM, 0.0, 500.0, 'Litter organic matter','LOM','t/ha',
notOnFile,writeInTable,notInGraph);
DeclMV( litterN, 0.0, 50.0, 'Litter nitrogen','LN','t/ha',
notOnFile,writeInTable,notInGraph);
DeclMV( HOM, 0.0, 500.0, 'Humus organic matter','HOM','t/ha',
notOnFile,writeInTable,notInGraph);
DeclMV( HN, 0.0, 50.0, 'Humus nitrogen','HN','t/ha',
notOnFile,writeInTable,notInGraph);
DeclMV( gLNC, 0.0, 0.05, 'Litter N:C ratio','gLNC','--',
notOnFile,notInTable,notInGraph);
DeclMV( humCN, 0.0, 100.0, 'Humus C:N ratio','gHCN','--',
notOnFile,notInTable,notInGraph);
DeclMV( totCO2, 0.0, 25.0, 'total C(O2) emission','CO2(tot)','t/ha*yr',
notOnFile,notInTable,notInGraph);
DeclMV( totSOM, 0.0, 500.0, 'total soil organic matter','totSOM','t/ha',
notOnFile,writeInTable,notInGraph);
DeclMV( meanAvN, 0.0, 500.0, 'average available nitrogen','meanAvN','t/ha',
notOnFile,notInTable,notInGraph);
DeclMV( meanLitM, 0.0, 500.0, 'average litter organic matter','meanLOM','t/ha',
notOnFile,notInTable,notInGraph);
DeclMV( meanHumM, 0.0, 500.0, 'average humus organic matter','meanHOM ','t/ha',
notOnFile,notInTable,notInGraph);
DeclP( kNAtm, 5.0, 0.0, 100.0, rtc,
'Atmospheric N input', 'kNAtm', 'kg/ha*yr');
DeclP( kAET, 1200.0, 0.0, 2000.0, rtc,
'AET multiplier parameter', 'kAET', 'mm/yr');
DeclP( kMin, 0.035, 0.0, 1.0, rtc,
'Humus decay in absence of litter', 'kMin', '%');
DeclP( kLeach, 0.16, 0.0, 1.0, rtc,
'Leaching from leaf litter', 'kLeach', '%');
DeclP( decMlt, 1.0, 0.0, 1.0, rtc,
'Decay multiplier', 'decMlt', '--');
DeclP( k1, 0.9804, 0.0, 10.0, rtc,
'Regression parameter', 'k1', '--');
DeclP( k2, 0.09352, 0.0, 1.0, rtc,
'Regression parameter', 'k2', '--');
DeclP( k3, -0.4956, -1.0, 1.0, rtc,
'Regression parameter', 'k3', '--');
DeclP( k4, 0.00193, 0.0, 1.0, rtc,
'Regression parameter', 'k4', '--');
DeclP( k5, 0.0079702, 0.0, 1.0, rtc,
'Regression parameter', 'k5', '--');
DeclP( k6, -1.3173E-4, -1.0, 0.0, rtc,
'Regression parameter', 'k6', '--');
DeclP( kNC, 0.005, 0.0, 0.1, rtc,
'N immobilized per unit weight loss', 'kNC', '--');
FOR i := MIN(Litter) TO MAX(Litter) DO
DeclP( kInitN[i], kInitN[i], 0.0, 0.05, rtc,
'Initial nitrogen percentage', 'kInitN', '%/100');
DeclP( kCritN[i], kCritN[i], 0.0, 0.1, rtc,
'Critical N percentage', 'kCritN', '%/100');
END;
END ModelObjects;