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 void SetPathway(PathwayParameters p)
{
double PsiFactor = 1.0;
var j = new LeafTemperatureParameters()
{
TMin = 0.0,
TOpt = 30.0,
TMax = 45.0,
C = 0.911017958600129,
Beta = 1.0
};
var g = new LeafTemperatureParameters()
{
TMin = 0.0,
TOpt = 29.2338417788683,
TMax = 45.0,
C = 0.875790608584141,
Beta = 1.0
};
var rubiscoCarboxylation = new TemperatureResponseValues()
{
At25 = 17.52 * 273.422964228666,
Factor = 93720.0
};
var rubiscoOxygenation = new TemperatureResponseValues()
{
At25 = 1.34 * 165824.064155384,
Factor = 33600.0
};
var rubiscoCarboxylationToOxygenation = new TemperatureResponseValues()
{
At25 = 13.07 * 4.59217066521612,
Factor = 35713.1987127717
};
var pepc = new TemperatureResponseValues()
{
At25 = 75,
Factor = 36300
};
var rubiscoActivity = new TemperatureResponseValues()
{
Factor = 65330.0
};
var respiration = new TemperatureResponseValues()
{
Factor = 46390
};
var pepcActivity = new TemperatureResponseValues()
{
Factor = 57043.2677590512
};
p.PEPRegeneration = 400;
p.SpectralCorrectionFactor = 0.15;
p.PS2ActivityFraction = 0.1;
p.BundleSheathConductance = 0.003;
p.MaxRubiscoActivitySLNRatio = 1.1 * PsiFactor;
p.MaxElectronTransportSLNRatio = 1.9484 * PsiFactor;
p.RespirationSLNRatio = 0.0 * PsiFactor;
p.MaxPEPcActivitySLNRatio = 1.0 * PsiFactor;
p.MesophyllCO2ConductanceSLNRatio = 0.00412 * PsiFactor;
p.ExtraATPCost = 0.75;
p.IntercellularToAirCO2Ratio = 0.7;
p.RubiscoCarboxylation = rubiscoCarboxylation;
p.RubiscoOxygenation = rubiscoOxygenation;
p.RubiscoCarboxylationToOxygenation = rubiscoCarboxylationToOxygenation;
p.RubiscoActivity = rubiscoActivity;
p.PEPc = pepc;
p.PEPcActivity = pepcActivity;
p.Respiration = respiration;
p.ElectronTransportRateParams = j;
p.MesophyllCO2ConductanceParams = g;
p.MesophyllElectronTransportFraction = p.ExtraATPCost / (3.0 + p.ExtraATPCost);
p.FractionOfCyclicElectronFlow = 0.25 * p.ExtraATPCost;
p.ATPProductionElectronTransportFactor = (3.0 - p.FractionOfCyclicElectronFlow) / (4.0 * (1.0 - p.FractionOfCyclicElectronFlow));
}
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);
}
public static PathwayParameters SetUpPathway(
double jTMin,
double jTOpt,
double jTMax,
double jC,
double jBeta,
double gTMin,
double gTOpt,
double gTMax,
double gC,
double gBeta,
double KcAt25,
double KcFactor,
double KoAt25,
double KoFactor,
double VcVoAt25,
double VcVoFactor,
double KpAt25,
double KpFactor,
double VcFactor,
double RdFactor,
double VpFactor,
double pepRegeneration,
double spectralCorrectionFactor,
double ps2ActivityFraction,
double bundleSheathConductance,
double maxRubiscoActivitySLNRatio,
double maxElectronTransportSLNRatio,
double respirationSLNRatio,
double maxPEPcActivitySLNRatio,
double mesophyllCO2ConductanceSLNRatio,
double extraATPCost,
double intercellularToAirCO2Ratio
)
{
var j = new LeafTemperatureParameters
{
TMin = jTMin,
TOpt = jTOpt,
TMax = jTMax,
C = jC,
Beta = jBeta
};
var g = new LeafTemperatureParameters
{
TMin = gTMin,
TOpt = gTOpt,
TMax = gTMax,
C = gC,
Beta = gBeta,
};
var Kc = new TemperatureResponseValues
{
At25 = KcAt25,
Factor = KcFactor
};
var Ko = new TemperatureResponseValues
{
At25 = KoAt25,
Factor = KoFactor
};
var VcVo = new TemperatureResponseValues
{
At25 = VcVoAt25,
Factor = VcVoFactor
};
var Kp = new TemperatureResponseValues
{
At25 = KpAt25,
Factor = KpFactor
};
var Vc = new TemperatureResponseValues
{
Factor = VcFactor
};
var Rd = new TemperatureResponseValues
{
Factor = RdFactor
};
var Vp = new TemperatureResponseValues
{
Factor = VpFactor
};
var PP = new PathwayParameters
{
PEPRegeneration = pepRegeneration,
SpectralCorrectionFactor = spectralCorrectionFactor,
PS2ActivityFraction = ps2ActivityFraction,
BundleSheathConductance = bundleSheathConductance,
MaxRubiscoActivitySLNRatio = maxRubiscoActivitySLNRatio,
MaxElectronTransportSLNRatio = maxElectronTransportSLNRatio,
RespirationSLNRatio = respirationSLNRatio,
MaxPEPcActivitySLNRatio = maxPEPcActivitySLNRatio,
MesophyllCO2ConductanceSLNRatio = mesophyllCO2ConductanceSLNRatio,
ExtraATPCost = extraATPCost,
IntercellularToAirCO2Ratio = intercellularToAirCO2Ratio,
RubiscoCarboxylation = Kc,
RubiscoOxygenation = Ko,
RubiscoCarboxylationToOxygenation = VcVo,
PEPc = Kp,
RubiscoActivity = Vc,
Respiration = Rd,
PEPcActivity = Vp,
ElectronTransportRateParams = j,
MesophyllCO2ConductanceParams = g
};
PP.MesophyllElectronTransportFraction = PP.ExtraATPCost / (3.0 + PP.ExtraATPCost);
PP.FractionOfCyclicElectronFlow = 0.25 * PP.ExtraATPCost;
PP.ATPProductionElectronTransportFactor = (3.0 - PP.FractionOfCyclicElectronFlow) / (4.0 * (1.0 - PP.FractionOfCyclicElectronFlow));
return(PP);
}
// 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);
}
public static void SetPathway(PathwayParameters p)
{
double PsiFactor = 0.4;
var j = new LeafTemperatureParameters()
{
TMin = 0,
TOpt = 37.8649150880407,
TMax = 55,
C = 0.711229539802063,
Beta = 1
};
var g = new LeafTemperatureParameters()
{
TMin = 0,
TOpt = 42,
TMax = 55,
C = 0.462820450976839,
Beta = 1,
};
var rubiscoCarboxylation = new TemperatureResponseValues()
{
At25 = 1210,
Factor = 64200
};
var rubiscoOxygenation = new TemperatureResponseValues()
{
At25 = 292000,
Factor = 10500
};
var rubiscoCarboxylationToOxygenation = new TemperatureResponseValues()
{
At25 = 5.51328906454566,
Factor = 21265.4029552906
};
var pepc = new TemperatureResponseValues()
{
At25 = 75,
Factor = 36300
};
var rubiscoActivity = new TemperatureResponseValues()
{
Factor = 78000
};
var respiration = new TemperatureResponseValues()
{
Factor = 46390
};
var pepcActivity = new TemperatureResponseValues()
{
Factor = 57043.2677590512
};
p.PEPRegeneration = 120;
p.SpectralCorrectionFactor = 0.15;
p.PS2ActivityFraction = 0.1;
p.BundleSheathConductance = 0.003;
p.MaxRubiscoActivitySLNRatio = 0.465 * PsiFactor;
p.MaxElectronTransportSLNRatio = 2.7 * PsiFactor;
p.RespirationSLNRatio = 0.0 * PsiFactor;
p.MaxPEPcActivitySLNRatio = 1.55 * PsiFactor;
p.MesophyllCO2ConductanceSLNRatio = 0.0135 * PsiFactor;
p.ExtraATPCost = 2;
p.MesophyllElectronTransportFraction = 0.4;
p.IntercellularToAirCO2Ratio = 0.45;
p.RubiscoCarboxylation = rubiscoCarboxylation;
p.RubiscoOxygenation = rubiscoOxygenation;
p.RubiscoCarboxylationToOxygenation = rubiscoCarboxylationToOxygenation;
p.PEPc = pepc;
p.RubiscoActivity = rubiscoActivity;
p.Respiration = respiration;
p.PEPcActivity = pepcActivity;
p.ElectronTransportRateParams = j;
p.MesophyllCO2ConductanceParams = g;
}
