/// <summary>
/// Init the universe for a new run.
/// </summary>
public void Init(SiteItem site, ManagementItem management)
{
Weather_ = new Weather(this);
Crop_ = new Crop(this);
Soil_ = new Soil(this);
thermalTimeWrapper_ = new ThermalTimeBiomaWrapper(this);
ShootTemperature_ = new ShootTemperature(this);
meteorologyWrapper_ = new MeteoBiomaWrapper(this);
Weather_.Init(site, management);
CurrentDate = management.FinalSowingDate.AddDays(-1);
Soil_.Init(Weather_.MinTemp(CurrentDate), Weather_.MaxTemp(CurrentDate), Weather_.GetMeanTemperatureBeforeSowing());
Crop_.Init(thermalTimeWrapper_.CumulTT, CurrentDate.AddDays(1));
CalcChangeinNFertilisation();
GrowthDay = 0;
CumAirTempFromSowing = 0;
CumMaxCanopyTempFromSowing = 0;
CumMaxAirTempFromSowing = 0;
if (management.Species == "Maize")
{
switchMaize = true;
}
else
{
switchMaize = false;
}
}
private void ShowDetails <T>(ItemDetails mgmtDetails, ManagementItem item) where T : IncomeGenerator
{
if (item.Item is Boat boat)
{
mgmtDetails.gameObject.SetActive(true);
mgmtDetails.UpdateDetails(item.Item.Descriptor, boat.CurrentHealth);
}
else if (item.Item is CrewMember)
{
mgmtDetails.gameObject.SetActive(true);
mgmtDetails.UpdateDetails(item.Item.Descriptor, 0f);
}
}
public void EstimateSowingWindows(SiteItem site, ManagementItem management)
{
if (site.SowingWindowType == 1) // JRC winter sowing window
{
EstimateSowingWindow_JRC_Winter(cache, site);
site.ESD_W = ESD_Winter(cache.Start.Value.Year);
site.LSD_W = LSD_Winter(cache.Start.Value.Year);
}
if (site.SowingWindowType == 2) // JRC spring sowing window
{
EstimateSowingWindow_JRC_Spring(cache, site);
site.ESD_S = ESD_Spring(cache.Start.Value.Year, site);
site.LSD_S = LSD_Spring(cache.Start.Value.Year, site);
}
if (site.SowingWindowType == 3) // SiriusQuality method (based on nominal sowing date)
{
EstimateSowingWindow_Sirius(cache, site, management);
site.ESD_Sir = ESD_Sirius(cache.Start.Value.Year, site);
site.LSD_Sir = LSD_Sirius(cache.Start.Value.Year, site);
}
}
public void DoStop()
{
varietyDef = null;
soilDef = null;
siteDef = null;
managementDef = null;
parametersDef = null;
runOptionsDef = null;
var nbRun = savedUniverses != null ? savedUniverses.Count : 0;
for (var i = 0; i < nbRun; ++i)
{
var universe = savedUniverses[i];
if (universe != null)
{
universe.Dispose();
savedUniverses[i] = null;
}
}
savedUniverses = null;
currentUniverse = null;
}
private void linkLabelName_LinkClicked(object sender, LinkLabelLinkClickedEventArgs e)
{
ManagementItem item = managementItemsBindingSource1.Current as ManagementItem;
if (item == null)
{
return;
}
using (var renameDialog = new RenameItem())
{
renameDialog.Item = item;
if (renameDialog.ShowDialog(this) == DialogResult.OK)
{
renameDialog.Apply(n => n.ManagementItem, (n, s) => n.ManagementItem = s,
m => m.ManagementItem, (m, s) => m.ManagementItem = s);
}
}
//update the list after the renaming
if (experimentNameBindingSource != null)
{
updateListBox1((string)experimentNameBindingSource.Current, item.Name);
}
}
///<summary>
///Update the run. This does the loop calculation over the year and calculate universe state.
///</summary>
///<param name="variety">The genotypeDef used in the run.</param>
///<param name="soil">The soilDef used in the run.</param>
///<param name="site">The siteDef used in the run.</param>
///<param name="management">The managementDef used in the run.</param>
///<param name="Inparameters">The parameter set used in the run.</param>
///<param name="runOptions">The run option set used in the run.</param>
public void Start(
CropParameterItem variety,
SoilItem soil,
SiteItem site,
ManagementItem management,
CropParameterItem Inparameters,
RunOptionItem runOptions)
{
this.DoStop();
// Initialisation
//Ph=new Phenology.Phenology(new Universe(this));
createParameters(variety, Inparameters);
this.varietyDef = variety;
this.soilDef = soil;
this.siteDef = site;
this.managementDef = management;
this.parametersDef = Inparameters;
this.runOptionsDef = runOptions;
if (variety != null &&
management != null &&
Inparameters != null &&
runOptions != null &&
site != null &&
soil != null)
{
}
else
{
throw new InvalidOperationException("Some input are null.");
}
savedUniverses = new RawData <Universe>(400);
//weather = new Weather(this);
currentUniverse = new Universe(this);
///<Behnam>
///<Comment>Estimating sowing date or using a fixed sowing date</Comment>
IsSowDateEstimate = management.IsSowDateEstimate;
SowingWindowType = site.SowingWindowType;
management.FinalSowingDate = management.SowingDate;
currentUniverse.Init(site, management);
currentUniverse.Weather_.EstimateSowingWindows(site, management);
currentUniverse.IsSowDateEstimate = IsSowDateEstimate;
double SkipDays = Math.Max(1, management.SkipDays);
SowingDate = management.SowingDate;
// UseActualBase = runOptionsDef.UseActualBase;
DoInteractions = runOptionsDef.DoInteractions;
InteractionsW = runOptionsDef.InteractionsW;
InteractionsN = runOptionsDef.InteractionsN;
InteractionsT = runOptionsDef.InteractionsT;
InteractionsWN = runOptionsDef.InteractionsWN;
InteractionsWT = runOptionsDef.InteractionsWT;
InteractionsNT = runOptionsDef.InteractionsNT;
InteractionsWNT = runOptionsDef.InteractionsWNT;
if (IsSowDateEstimate)
{
var year1 = SowingDate.Year;
if (SowingWindowType == 0) // Fixed sowing window
{
site.FinalMinSowingDate = currentUniverse.Weather_.ESD_Fixed(year1, site);
site.FinalMaxSowingDate = currentUniverse.Weather_.LSD_Fixed(year1, site);
}
else if (SowingWindowType == 1) // JRC winter sowing window
{
site.FinalMinSowingDate = currentUniverse.Weather_.ESD_Winter(year1);
site.FinalMaxSowingDate = currentUniverse.Weather_.LSD_Winter(year1);
}
else if (SowingWindowType == 2) // JRC spring sowing window
{
site.FinalMinSowingDate = currentUniverse.Weather_.ESD_Spring(year1, site);
site.FinalMaxSowingDate = currentUniverse.Weather_.LSD_Spring(year1, site);
}
else if (SowingWindowType == 3) // SiriusQuality method (based on nominal sowing date)
{
site.FinalMinSowingDate = currentUniverse.Weather_.ESD_Sirius(year1, site);
site.FinalMaxSowingDate = currentUniverse.Weather_.LSD_Sirius(year1, site);
}
flagReSow = true;
MinSowDate = site.FinalMinSowingDate.AddDays(-SkipDays);
MaxSowDate = site.FinalMaxSowingDate;
SowingDate = MinSowDate;
management.FinalSowingDate = MinSowDate;
}
else
{
flagReSow = false;
management.FinalSowingDate = SowingDate;
site.FinalMinSowingDate = site.MinSowingDate;
site.FinalMaxSowingDate = site.MaxSowingDate;
}
currentUniverse.Init(site, management);
SaveCurrentUniverse();
currentUniverse.CurrentDate = currentUniverse.CurrentDate.AddDays(1);
//currentUniverse.Crop_.Sow(currentUniverse.calculateDailyThermalTime_.CumulTT, currentUniverse.CurrentDate);
int numdSkipDays = 0;
int numdFreezing = 0;
int numdSoilTemp = 0;
int numdSoilMois = 0;
double numdPcp = 0;
///<Behnam (2015.11.27)>
///<Comment>To enable using soil depth without any limitations</Comment>
checkLay = Convert.ToInt16(Math.Ceiling(management.CheckDepth / 5)); // Soil layer to be ckecked
checkLay = Math.Min(checkLay, currentUniverse.Soil_.Layers.Count) - 1;
///</Behnam>
double CheckDaysTemp = management.CheckDaysTemp;
double CheckDaysPcp = management.CheckDaysPcp;
double TAveThr = management.TAveThr;
double TMinThr = management.TMinThr;
double SoilMoistThr = management.SoilMoistThr;
double CumPcpThr = management.CumPcpThr;
double SoilWorkabThr = management.SoilWorkabThr;
///<Relaxation variables>
double origTAveThr = management.TAveThr;
double origTMinThr = management.TMinThr;
double origSoilMoistThr = management.SoilMoistThr;
double origCumPcpThr = management.CumPcpThr;
double origSoilWorkabThr = management.SoilWorkabThr;
double poTAveThr = 1.0;
double poTMinThr = 1.0;
double poSoilMoistThr = 1.0;
double poCumPcpThr = 1.0;
double poSoilWorkabThr = 1.0;
double length = (int)(site.FinalMaxSowingDate - site.FinalMinSowingDate).TotalDays;
double TAveThrRelax = 0;
double TMinThrRelax = -4;
double SoilMoistThrRelax = 0;
double CumPcpThrRelax = 3;
double SoilWorkabThrRelax = 1.5;
///</Relaxation variables>
///<Relaxation curves>
poTAveThr = DecreasePower(TAveThr, TAveThrRelax, length);
poTMinThr = DecreasePower(TMinThr, TMinThrRelax, length);
poSoilMoistThr = DecreasePower(SoilMoistThr, SoilMoistThrRelax, length);
poCumPcpThr = DecreasePower(CumPcpThr, CumPcpThrRelax, length);
poSoilWorkabThr = IncreasePower(SoilWorkabThr, SoilWorkabThrRelax, length);
///</Relaxation curves>
bool oldLimitedWater = runOptionsDef.UnlimitedWater;
bool oldLimitedNitrogen = runOptionsDef.UnlimitedNitrogen;
bool oldLimitedTemperature = runOptionsDef.UnlimitedTemperature;
bool initSoilAfterSowing = false;
double minTempInitSoilAfterSowing = 0.0;
double maxTempInitSoilAfterSowing = 0.0;
double meanTempInitSoilAfterSowing = 0.0;
while (!currentUniverse.Crop_.IsEnd)
{
if (!IsSowDateEstimate && !flagReSow)
///<Comment>Non-stressed conditions are only applied after sowing</Comment>
{
if (DoInteractions)
{
///<Parallel runs>
double OutputTotalDM_W = 0;
double OutputTotalDM_N = 0;
double OutputTotalDM_T = 0;
double OutputTotalDM_WN = 0;
double OutputTotalDM_WT = 0;
double OutputTotalDM_NT = 0;
double OutputTotalDM_WNT = 0;
/// Behnam (2016.07.13): The runOld.RunOptionDef.UseActualBase was previously used. Pierre wanted it to be deleted.
/// At one stage, we decided to use Potential conditions as the base run. Now it is not working anymore.
if (true)
{
///<Use actual conditions as the base run>
if (InteractionsW)
{
OutputTotalDM_W = RunUnlimited(currentUniverse, true, false, false);
}
if (InteractionsN)
{
OutputTotalDM_N = RunUnlimited(currentUniverse, false, true, false);
}
if (InteractionsT)
{
OutputTotalDM_T = RunUnlimited(currentUniverse, false, false, true);
}
if (InteractionsWN)
{
OutputTotalDM_WN = RunUnlimited(currentUniverse, true, true, false);
}
if (InteractionsWT)
{
OutputTotalDM_WT = RunUnlimited(currentUniverse, true, false, true);
}
if (InteractionsNT)
{
OutputTotalDM_NT = RunUnlimited(currentUniverse, false, true, true);
}
if (InteractionsWNT)
{
OutputTotalDM_WNT = RunUnlimited(currentUniverse, true, true, true);
}
currentUniverse.isUnlimitedWater = false;
currentUniverse.isUnlimitedNitrogen = false;
currentUniverse.isUnlimitedTemperature = false;
currentUniverse.Crop_.areRootsToBeGrown = true;
if (initSoilAfterSowing)
{
currentUniverse.Soil_.Init(minTempInitSoilAfterSowing, maxTempInitSoilAfterSowing, meanTempInitSoilAfterSowing, false);
currentUniverse.ShootTemperature_.MinShootTemperature = 0;
currentUniverse.ShootTemperature_.MaxShootTemperature = 0;
currentUniverse.ShootTemperature_.MeanShootTemperature = 0;
currentUniverse.CumAirTempFromSowing = 0.0;
currentUniverse.CumMaxCanopyTempFromSowing = 0.0;
currentUniverse.CumMaxAirTempFromSowing = 0.0;
currentUniverse.PrevAppliedStagesN = 1;
currentUniverse.PrevAppliedStagesIrr = 1;
//currentUniverse.NFertChange = 1;
currentUniverse.CalcChangeinNFertilisation();
initSoilAfterSowing = false;
}
currentUniverse.RunDayStep();
}
else
{
///<Use non-stressed conditions as the base run>
if (InteractionsW)
{
OutputTotalDM_W = RunUnlimited(currentUniverse, false, true, true);
}
if (InteractionsN)
{
OutputTotalDM_N = RunUnlimited(currentUniverse, true, false, true);
}
if (InteractionsT)
{
OutputTotalDM_T = RunUnlimited(currentUniverse, true, true, false);
}
if (InteractionsWN)
{
OutputTotalDM_WN = RunUnlimited(currentUniverse, false, false, true);
}
if (InteractionsWT)
{
OutputTotalDM_WT = RunUnlimited(currentUniverse, false, true, false);
}
if (InteractionsNT)
{
OutputTotalDM_NT = RunUnlimited(currentUniverse, true, false, false);
}
if (InteractionsWNT)
{
OutputTotalDM_WNT = RunUnlimited(currentUniverse, false, false, false);
}
currentUniverse.isUnlimitedWater = true;
currentUniverse.isUnlimitedNitrogen = true;
currentUniverse.isUnlimitedTemperature = true;
currentUniverse.Crop_.areRootsToBeGrown = true;
if (initSoilAfterSowing)
{
currentUniverse.Soil_.Init(minTempInitSoilAfterSowing, maxTempInitSoilAfterSowing, meanTempInitSoilAfterSowing, false);
currentUniverse.ShootTemperature_.MinShootTemperature = 0;
currentUniverse.ShootTemperature_.MaxShootTemperature = 0;
currentUniverse.ShootTemperature_.MeanShootTemperature = 0;
currentUniverse.CumAirTempFromSowing = 0.0;
currentUniverse.CumMaxCanopyTempFromSowing = 0.0;
currentUniverse.CumMaxAirTempFromSowing = 0.0;
currentUniverse.PrevAppliedStagesN = 1;
currentUniverse.PrevAppliedStagesIrr = 1;
//currentUniverse.NFertChange = 1;
currentUniverse.CalcChangeinNFertilisation();
}
currentUniverse.RunDayStep();
}
currentUniverse.Crop_.OutputTotalDM_W = OutputTotalDM_W;
currentUniverse.Crop_.OutputTotalDM_N = OutputTotalDM_N;
currentUniverse.Crop_.OutputTotalDM_T = OutputTotalDM_T;
currentUniverse.Crop_.OutputTotalDM_WN = OutputTotalDM_WN;
currentUniverse.Crop_.OutputTotalDM_WT = OutputTotalDM_WT;
currentUniverse.Crop_.OutputTotalDM_NT = OutputTotalDM_NT;
currentUniverse.Crop_.OutputTotalDM_WNT = OutputTotalDM_WNT;
///</Parallel runs>
}
else
{
///<Normal run>
currentUniverse.isUnlimitedWater = oldLimitedWater;
currentUniverse.isUnlimitedNitrogen = oldLimitedNitrogen;
currentUniverse.isUnlimitedTemperature = oldLimitedTemperature;
currentUniverse.Crop_.areRootsToBeGrown = true;
if (initSoilAfterSowing)
{
currentUniverse.Soil_.Init(minTempInitSoilAfterSowing, maxTempInitSoilAfterSowing, meanTempInitSoilAfterSowing, false);
currentUniverse.ShootTemperature_.MinShootTemperature = 0;
currentUniverse.ShootTemperature_.MaxShootTemperature = 0;
currentUniverse.ShootTemperature_.MeanShootTemperature = 0;
currentUniverse.CumAirTempFromSowing = 0.0;
currentUniverse.CumMaxCanopyTempFromSowing = 0.0;
currentUniverse.CumMaxAirTempFromSowing = 0.0;
currentUniverse.PrevAppliedStagesN = 1;
currentUniverse.PrevAppliedStagesIrr = 1;
//currentUniverse.NFertChange = 1;
currentUniverse.CalcChangeinNFertilisation();
initSoilAfterSowing = false;
}
currentUniverse.RunDayStep();
///</Normal run>
}
}
else
{
///<Comment>Non-stressed conditions are only applied after sowing</Comment>
currentUniverse.isUnlimitedWater = false;
currentUniverse.isUnlimitedNitrogen = false;
currentUniverse.isUnlimitedTemperature = false;
currentUniverse.Crop_.areRootsToBeGrown = false;
currentUniverse.RunDayStep();
}
if (IsSowDateEstimate)
{
// Checking sowing conditions:
///<Skip days>
numdSkipDays += 1;
flagSkipDays = (numdSkipDays <= SkipDays);
///</Skip days>
flagMaxDate = true;
flagSoilMois = true;
flagCumPcp = true;
flagSoilTemp = true;
flagFreezing = true;
flagWorkabil = true;
if (!flagSkipDays)
{
if (management.DoRelax)
{
///<Relaxation>
///<Comment>The criteria are relaxed here using four exponential curves, if DoRelax=TRUE</Comment>
var n = numdSkipDays - SkipDays - 1;
TAveThr = ExpoDecrease(TAveThr, TAveThrRelax, origTAveThr, n, poTAveThr);
TMinThr = ExpoDecrease(TMinThr, TMinThrRelax, origTMinThr, n, poTMinThr);
SoilMoistThr = ExpoDecrease(SoilMoistThr, SoilMoistThrRelax, origSoilMoistThr, n, poSoilMoistThr);
CumPcpThr = ExpoDecrease(CumPcpThr, CumPcpThrRelax, origCumPcpThr, n, poCumPcpThr);
SoilWorkabThr = ExpoIncrease(SoilWorkabThr, SoilWorkabThrRelax, origSoilWorkabThr, n, poSoilWorkabThr);
///</Relaxation>
}
///<Average daily air temperature temperature>
numdSoilTemp = 0;
for (var i = 1; i <= CheckDaysTemp; ++i)
{
if (currentUniverse.Weather_.MeanTemp(currentUniverse.CurrentDate.AddDays(i)) >= TAveThr)
{
numdSoilTemp += 1;
}
}
flagSoilTemp = (numdSoilTemp < CheckDaysTemp);
///</Average daily air temperature temperature>
///<Minimum daily air temperature (frost risk)>
numdFreezing = 0;
for (var i = 1; i <= CheckDaysTemp; ++i)
{
if (currentUniverse.Weather_.MinTemp(currentUniverse.CurrentDate.AddDays(i)) >= TMinThr)
{
numdFreezing += 1;
}
}
flagFreezing = (numdFreezing < CheckDaysTemp);
///</Minimum daily air temperature (frost risk)>
///<Cumulative precipitation>
numdPcp = 0;
for (var i = 1; i <= CheckDaysPcp; ++i)
{
numdPcp += currentUniverse.Weather_.Rain(currentUniverse.CurrentDate.AddDays(i)) / MMwaterToGwater;
}
flagCumPcp = (numdPcp < CumPcpThr);
///</Cumulative precipitation>
///<Soil moisture>
if (currentUniverse.Soil_.Layers[checkLay].AvWater < SoilMoistThr * currentUniverse.Soil_.Layers[checkLay].MaxAvWater)
{
numdSoilMois = 0;
}
if (currentUniverse.Soil_.Layers[checkLay].AvWater >= SoilMoistThr * currentUniverse.Soil_.Layers[checkLay].MaxAvWater)
{
numdSoilMois += 1;
}
flagSoilMois = (numdSoilMois == 0);
///</Soil moisture>
///<Soil compaction (workability)>
flagWorkabil = ((currentUniverse.Soil_.Layers[checkLay].AvWater +
currentUniverse.Soil_.Layers[checkLay].ExWater) / currentUniverse.Soil_.Layers[checkLay].FcWater
> SoilWorkabThr);
///</Soil compaction (workability)>
}
///<Sowing window>
flagMaxDate = (currentUniverse.CurrentDate < MaxSowDate.AddDays(-1));
///</Sowing window>
flagReSow = flagMaxDate & (flagSkipDays || flagSoilMois || flagCumPcp ||
flagSoilTemp || flagFreezing || flagWorkabil);
}
if (IsSowDateEstimate) // If sowing date search is ON, at least one day before sowing must be simulated.
{
SowingDate = currentUniverse.CurrentDate.AddDays(1);
currentUniverse.Crop_.Dispose();
currentUniverse.Crop_ = new CropModel.Crop(currentUniverse);
//Debug
//currentUniverse.CumMaxCanopyTempFromSowing = 0.0;
//??currentUniverse.Crop_.RootLength = 0.0;
//??currentUniverse.Crop_.Tau = 1.0;
//??currentUniverse.Crop_.SumInternodesLength = 0.0;
//??currentUniverse.Crop_.PotentialWaterOnLeaves = 0.0;
//??currentUniverse.Crop_.RootLength = 0.0;
//currentUniverse.Crop_.LeafNumber = 0.0;
//currentUniverse.Crop_.Ntip = 0.0;
/*currentUniverse.ShootTemperature_.MinShootTemperature=0.0;
* currentUniverse.ShootTemperature_.MaxShootTemperature=0.0;
* //currentUniverse.Soil_.Init(currentUniverse.Weather_.MinTemp(currentUniverse.CurrentDate/*.AddDays(1)*//*)/*, currentUniverse.Weather_.MaxTemp(currentUniverse.CurrentDate/*.AddDays(1)*//*), (currentUniverse.Weather_.MinTemp(currentUniverse.CurrentDate/*.AddDays(1)*//*) + currentUniverse.Weather_.MaxTemp(currentUniverse.CurrentDate/*.AddDays(1)*//*))/2, false);*/
//currentUniverse.Soil_.energyBalanceBiomaWrapper_.Init(false);
//Debug /**/
//currentUniverse.Soil_.MaximumCanopyTemperature = currentUniverse.Weather_.MaxTemp(currentUniverse.CurrentDate/*.AddDays(1)*/);//??
//currentUniverse.Soil_.MinimumCanopyTemperature = currentUniverse.Weather_.MinTemp(currentUniverse.CurrentDate/*.AddDays(1)*/);//??
//currentUniverse.Soil_.SoilMinTemperature = currentUniverse.Weather_.MinTemp(currentUniverse.CurrentDate/*.AddDays(1)*/);//??
//currentUniverse.Soil_.SoilMaxTemperature = currentUniverse.Weather_.MaxTemp(currentUniverse.CurrentDate/*.AddDays(1)*/);//??;
//currentUniverse.ShootTemperature_.MinShootTemperature = 0;
//currentUniverse.ShootTemperature_.MaxShootTemperature = 0;
//currentUniverse.ShootTemperature_.MeanShootTemperature = 0;
//currentUniverse.Crop_.Universe_.thermalTimeWrapper_.Init();
initSoilAfterSowing = true;
minTempInitSoilAfterSowing = currentUniverse.Weather_.MinTemp(currentUniverse.CurrentDate /*.AddDays(1)*/);
maxTempInitSoilAfterSowing = currentUniverse.Weather_.MaxTemp(currentUniverse.CurrentDate /*.AddDays(1)*/);
// meanTempInitSoilAfterSowing = (currentUniverse.Weather_.MinTemp(currentUniverse.CurrentDate/*.AddDays(1)*/) + currentUniverse.Weather_.MaxTemp(currentUniverse.CurrentDate/*.AddDays(1)*/)) / 2;
//Debug /**/
currentUniverse.Crop_.Universe_.thermalTimeWrapper_.Init();
currentUniverse.Crop_.Init(currentUniverse.thermalTimeWrapper_.CumulTT, currentUniverse.CurrentDate.AddDays(1));
currentUniverse.CurrentDate = currentUniverse.CurrentDate.AddDays(1);
meanTempInitSoilAfterSowing = currentUniverse.Weather_.GetMeanTemperatureBeforeSowing(SowingDate);
}
if (IsSowDateEstimate & flagReSow) // At least one of the sowing criteria is not met.
{
IsSowDateEstimate = true; // Sowing date search continues;
SaveCurrentUniverse();
}
else if (IsSowDateEstimate & !flagReSow) // All of the sowing criteria are met.
{
IsSowDateEstimate = false; // Sowing date search ends;
currentUniverse.IsSowDateEstimate = false;
management.FinalSowingDate = SowingDate;
currentUniverse.CalcChangeinNFertilisation();
SaveCurrentUniverse();
}
else if (!IsSowDateEstimate) // Sowing date search is OFF or all of the sowing criteria are met.
{
SaveCurrentUniverse();
currentUniverse.CurrentDate = currentUniverse.CurrentDate.AddDays(1);
}
}
}
private void SellItem(ManagementItem item)
{
_playerController.SellItem(item);
_items.Remove(item);
Destroy(item.gameObject);
}
public void EstimateSowingWindow_Sirius(Cache cache, SiteItem site, ManagementItem management)
{
DateTime date1;
DateTime date2;
DateTime datex;
int daynx;
int ndays;
double tav10d;
double pcp3d;
int HalfSowWindow = (int)(InitSowWindow / 2);
List <int> esd = new List <int>();
List <int> lsd = new List <int>();
var TempThr = site.TempThr;
var TempSum = site.TempSum;
var PcpSum = site.PcpSum;
var CheckDaysPcp = site.CheckDaysPcp;
var CheckDaysTemp = site.CheckDaysTemp;
var year1 = cache.Start.Value.Year;
var year2 = cache.End.Value.Year;
var yy = management.SowingDate.AddDays(HalfSowWindow).Year - management.SowingDate.AddDays(-HalfSowWindow).Year;
date1 = new DateTime(year1, management.SowingDate.Month, management.SowingDate.Day);
esd_Sir = date1.DayOfYear;
lsd_Sir = date1.DayOfYear;
esd.Add(esd_Sir);
lsd.Add(lsd_Sir);
for (var year = year1; year <= year2; ++year)
{
date1 = new DateTime(year, management.SowingDate.Month, management.SowingDate.Day);
date2 = date1.AddDays(HalfSowWindow);
ndays = (int)(date2 - date1).TotalDays;
if (date2 <= cache.End.Value && date1 >= cache.Start.Value)
{
for (var dayn = 1; dayn <= ndays; ++dayn)
{
datex = date1.AddDays(dayn - 1);
daynx = datex.DayOfYear;
tav10d = cache.TemperatureMean(datex.AddDays(-1), CheckDaysTemp);
if (tav10d >= TempThr)
{
pcp3d = cache.Rain(datex.AddDays(-1), CheckDaysPcp);
if (pcp3d >= PcpSum)
{
var xx = datex.Year - year;
lsd.Add(xx * 365 + daynx);
break;
}
}
}
}
date2 = new DateTime(year, management.SowingDate.Month, management.SowingDate.Day);
date1 = date2.AddDays(-HalfSowWindow);
if (date1 < cache.Start.Value)
{
date1 = cache.Start.Value;
}
ndays = (int)(date2 - date1).TotalDays;
if (date2 <= cache.End.Value)
{
for (var dayn = 1; dayn <= ndays; ++dayn)
{
datex = date2.AddDays(-dayn + 1);
daynx = datex.DayOfYear;
tav10d = cache.TemperatureMean(datex.AddDays(-1), CheckDaysTemp);
if (tav10d >= TempThr)
{
pcp3d = cache.Rain(datex.AddDays(-1), CheckDaysPcp);
if (pcp3d >= PcpSum)
{
var xx = datex.Year - year;
esd.Add(xx * 365 + daynx);
break;
}
}
}
}
}
esd_Sir = (int)Percentile(esd, 0.25);
lsd_Sir = (int)Percentile(lsd, 0.75);
}
///<summary>
///Read the weather file and compute weekTemp, threeDayRain and weekRain.
///</summary>
public void Init(SiteItem site, ManagementItem management)
{
isUnlimitedTemperature = false;
var weatherFiles = site.WeatherFiles;
if (site.Format == WeatherFileFormat.YearJdayMinMaxRainRad)
{
format = NoWindVp;
}
else
{
if (site.Format == WeatherFileFormat.YearJdayMinMaxRainRadWindVp)
{
format = WindVp;
}
else
{
format = WindVpHourly;
}
}
var input = new WeatherInput {
Files = new List <string>(), Format = format, MaxConsecutiveMissingValuesAllowed = 5
};
foreach (var file in weatherFiles)
{
input.Files.Add(file.AbsoluteFile);
}
cache = Cache.Get(site.GeoPosition, input);
if (!cache.Start.HasValue)
{
throw new InvalidOperationException("Weather has no starting date.");
}
if (cache.Start.Value > management.SowingDate.AddDays(-(NumberOfFirstDays + 1)))
{
throw new InvalidOperationException("Sowing date - number of first days needed by DeepLayer " + management.SowingDate.ToString("u").Split()[0] + " - (" + NumberOfFirstDays + ")" + " is before the first daily weather data " + cache.Start.Value.ToString("u").Split()[0] + ".");
}
if (cache.Start.Value.AddDays(cache.Count) <= management.SowingDate)
{
throw new InvalidOperationException("Sowing date " + management.SowingDate.ToString("u").Split()[0] + " is after the last daily weather data " + cache.Start.Value.AddDays(cache.Count).ToString("u").Split()[0] + ".");
}
///<Behnam>
if (management.IsSowDateEstimate & site.SowingWindowType != 1 & site.SowingWindowType != 3 &
site.MinSowingDate > site.MaxSowingDate)
{
throw new InvalidOperationException("Minimum sowing date " +
site.MinSowingDate.ToString("u").Split()[0] +
" is after maximum sowing date " + site.MaxSowingDate.ToString("u").Split()[0] + ".");
}
if (management.IsSowDateEstimate & site.SowingWindowType != 1 & site.SowingWindowType != 3 &
site.MinSowingDate.AddDays(-management.SkipDays - NumberOfFirstDays - 1) < cache.Start.Value)
{
throw new InvalidOperationException("First date of simulation " +
site.MinSowingDate.AddDays(-management.SkipDays - NumberOfFirstDays - 1).ToString("u").Split()[0] +
" is before the first daily weather data " + cache.Start.Value.ToString("u").Split()[0] + ".");
}
if (management.IsSowDateEstimate & site.SowingWindowType != 1 & site.SowingWindowType != 3 &
site.MaxSowingDate >= cache.Start.Value.AddDays(cache.Count))
{
throw new InvalidOperationException("Maximum sowing date " + site.MaxSowingDate.ToString("u").Split()[0] +
" is after the last daily weather data " +
cache.Start.Value.AddDays(cache.Count).ToString("u").Split()[0] + ".");
}
}
