public static void SetCanopy(CanopyParameters c)
{
c.Type = CanopyType.CCM;
c.AirCO2 = 370;
c.CurvatureFactor = 0.7;
c.DiffusivitySolubilityRatio = 0.047;
c.AirO2 = 210000;
c.DiffuseExtCoeff = 0.78;
c.DiffuseExtCoeffNIR = 0.8;
c.DiffuseReflectionCoeff = 0.036;
c.DiffuseReflectionCoeffNIR = 0.389;
c.LeafAngle = 60;
c.LeafScatteringCoeff = 0.15;
c.LeafScatteringCoeffNIR = 0.8;
c.LeafWidth = 0.05;
c.SLNRatioTop = 1.3;
c.MinimumN = 14;
c.Windspeed = 1.5;
c.WindSpeedExtinction = 1.5;
}
public void OnInitialised()
{
CP = Classic.SetUpCanopy(
CanopyType.C4, // Canopy type
363, // CO2 partial pressure
0.7, // Curvature factor
0.047, // Diffusivity-solubility ratio
210000, // O2 partial pressure
0.78, // Diffuse extinction coefficient
0.8, // Diffuse extinction coefficient NIR
0.036, // Diffuse reflection coefficient
0.389, // Diffuse reflection coefficient NIR
60, // Leaf angle
0.15, // Leaf scattering coefficient
0.8, // Leaf scattering coefficient NIR
0.15, // Leaf width
1.3, // SLN ratio at canopy top
14, // Minimum Nitrogen
1.5, // Wind speed
1.5); // Wind speed extinction
PP = Classic.SetUpPathway(
0, // jTMin
37.8649150880407, // jTOpt
55, // jTMax
0.711229539802063, // jC
1, // jBeta
0, // gTMin
42, // gTOpt
55, // gTMax
0.462820450976839, // gC
1, // gBeta
1210, // KcAt25
64200, // KcFactor
292000, // KoAt25
10500, // KoFactor
5.51328906454566, // VcVoAt25
21265.4029552906, // VcVoFactor
75, // KpAt25
36300, // KpFactor
78000, // VcFactor
46390, // RdFactor
57043.2677590512, // VpFactor
120, // pepRegeneration
0.15, // spectralCorrectionFactor
0.1, // ps2ActivityFraction
0.003, // bundleSheathConductance
0.465 * PsiFactor, // maxRubiscoActivitySLNRatio
2.7 * PsiFactor, // maxElectronTransportSLNRatio
0.0 * PsiFactor, // respirationSLNRatio
1.55 * PsiFactor, // maxPEPcActivitySLNRatio
0.0135 * PsiFactor, // mesophyllCO2ConductanceSLNRatio
2, // extraATPCost
0.45); // intercellularToAirCO2Ratio
//Set the LAI trigger
MyPaddock.Set("laiTrigger", LAITrigger);
}
public static CanopyParameters SetUpCanopy(
CanopyType type,
double airCO2,
double curvatureFactor,
double diffusivitySolubilityRatio,
double airO2,
double diffuseExtCoeff,
double diffuseExtCoeffNIR,
double diffuseReflectionCoeff,
double diffuseReflectionCoeffNIR,
double leafAngle,
double leafScatteringCoeff,
double leafScatteringCoeffNIR,
double leafWidth,
double slnRatioTop,
double minimumN,
double windspeed,
double windSpeedExtinction
)
{
var CP = new CanopyParameters
{
Type = type,
AirCO2 = airCO2,
CurvatureFactor = curvatureFactor,
DiffusivitySolubilityRatio = diffusivitySolubilityRatio,
AirO2 = airO2,
DiffuseExtCoeff = diffuseExtCoeff,
DiffuseExtCoeffNIR = diffuseExtCoeffNIR,
DiffuseReflectionCoeff = diffuseReflectionCoeff,
DiffuseReflectionCoeffNIR = diffuseReflectionCoeffNIR,
LeafAngle = leafAngle,
LeafScatteringCoeff = leafScatteringCoeff,
LeafScatteringCoeffNIR = leafScatteringCoeffNIR,
LeafWidth = leafWidth,
SLNRatioTop = slnRatioTop,
MinimumN = minimumN,
Windspeed = windspeed,
WindSpeedExtinction = windSpeedExtinction
};
return(CP);
}
public void OnInitialised()
{
/* IMPORTANT - Do NOT change the order of these values */
CP = Classic.SetUpCanopy(
CanopyType.C3, // Canopy type
370, // CO2 partial pressure
0.7, // Empirical curvature factor
0.047, // Diffusivity-solubility ratio
210000, // O2 partial pressure
0.78, // PAR diffuse extinction coefficient
0.8, // NIR diffuse extinction coefficient
0.036, // PAR diffuse reflection coefficient
0.389, // NIR diffuse reflection coefficient
60, // Leaf angle
0.15, // PAR leaf scattering coefficient
0.8, // NIR leaf scattering coefficient
0.05, // Leaf width
1.3, // SLN ratio at canopy top
14, // Minimum structural nitrogen
1.5, // Wind speed
1.5); // Wind speed profile distribution coefficient
PP = Classic.SetUpPathway(
0, // Electron transport minimum temperature
30.0, // Electron transport optimum temperature
45.0, // Electron transport maximum temperature
0.911017958600129, // Electron transport scaling constant
1, // Electron transport Beta value
6048.95289, //mesophyll conductance factor
273.422964228666, // Kc25 - Michaelis Menten constant of Rubisco carboxylation at 25 degrees C
93720, // KcFactor
165824.064155384, // Ko25 - Michaelis Menten constant of Rubisco oxygenation at 25 degrees C
33600, // KoFactor
4.59217066521612, // VcVo25 - Rubisco carboxylation to oxygenation at 25 degrees C
35713.19871277176, // VcVoFactor
0.0, // Kp25 - Michaelis Menten constant of PEPc activity at 25 degrees C (Unused in C3)
0.0, // KpFactor (Unused in C3)
65330, // VcFactor
46390, // RdFactor
57043.2677590512, // VpFactor
0.0, // PEPc regeneration (Unused in C3)
0.15, // Spectral correction factor
0.1, // Photosystem II activity fraction
0.003, // Bundle sheath CO2 conductance
1.1 * PsiFactor, // Max Rubisco activity to SLN ratio
1.85 * PsiFactor, // Max electron transport to SLN ratio
0.0 * PsiFactor, // Respiration to SLN ratio
1.0 * PsiFactor, // Max PEPc activity to SLN ratio
0.00412 * PsiFactor, // Mesophyll CO2 conductance to SLN ratio
0.75, // Extra ATP cost
0.7); // Intercellular CO2 to air CO2 ratio
//Set the LAI trigger
MyPaddock.Set("DCaPSTTriggerLAI", LAITrigger);
MyPaddock.Get("PsModelName1", out PsModelName1);
}
// The following event handler will be called once at the beginning of the simulation
[EventHandler] public void OnInitialised()
{
string path = "IntervalValues.csv";
stream = new FileStream(path, FileMode.Create);
writer = new StreamWriter(stream);
/* Do NOT change the order of these values */
CP = Classic.SetUpCanopy(
CanopyType.CCM, // Canopy type
370, // CO2 partial pressure
0.7, // Empirical curvature factor
0.047, // Diffusivity-solubility ratio (Used in CCM)
210000, // O2 partial pressure
0.78, // PAR diffuse extinction coefficient
0.8, // NIR diffuse extinction coefficient
0.036, // PAR diffuse reflection coefficient
0.389, // NIR diffuse reflection coefficient
60, // Leaf angle
0.15, // PAR leaf scattering coefficient
0.8, // NIR leaf scattering coefficient
0.05, // Leaf width
1.3, // SLN ratio at canopy top
14, // Minimum structural nitrogen
1.5, // Wind speed
1.5); // Wind speed profile distribution coefficient
PP = Classic.SetUpPathway(
0, // Electron transport minimum temperature
30.0, // Electron transport optimum temperature
45.0, // Electron transport maximum temperature
0.911017958600129, // Electron transport scaling constant
1, // Electron transport Beta value
0, // Mesophyll conductance minimum temperature
29.2338417788683, // Mesophyll conductance optimum temperature
45, // Mesophyll conductance maximum temperature
0.875790608584141, // Mesophyll conductance scaling constant
1, // Mesophyll conductance Beta value
17.52 * 273.422964228666, // Kc25 - Michaelis Menten constant of Rubisco carboxylation at 25 degrees C (Changed in CCM)*********************
93720, // KcFactor
1.34 * 165824.064155384, // Ko25 - Michaelis Menten constant of Rubisco oxygenation at 25 degrees C (Changed in CCM)*********************
33600, // KoFactor
13.07 * 4.59217066521612, // VcVo25 - Rubisco carboxylation to oxygenation at 25 degrees C (Changed in CCM)*********************
35713.19871277176, // VcVoFactor
75, // Kp25 - Michaelis Menten constant of PEPc activity at 25 degrees C (Used in CCM)
36300, // KpFactor (Used in CCM)
65330, // VcFactor
46390, // RdFactor
57043.2677590512, // VpFactor (Used in CCM)
400.0, // PEPc regeneration (Used in CCM)
0.15, // Spectral correction factor
0.1, // Photosystem II activity fraction (used in CCM)
0.003, // Bundle sheath CO2 conductance (Used in CCM)********************
1.1 * PsiFactor, // Max Rubisco activity to SLN ratio
1.9484 * PsiFactor, // Max electron transport to SLN ratio (Changed in CCM)********************
0.0 * PsiFactor, // Respiration to SLN ratio
1 * PsiFactor, // Max PEPc activity to SLN ratio (Changed in CCM)*********************0.373684157583268
0.00412 * PsiFactor, // Mesophyll CO2 conductance to SLN ratio
0.75, // Extra ATP cost (Changed in CCM)*********************
0.7); // Intercellular CO2 to air CO2 ratio
//Set the LAI trigger
MyPaddock.Set("DCaPSTTriggerLAI", LAITrigger);
MyPaddock.Get("PsModelName1", out PsModelName1);
}