FarmMaps/CarbonLib
13
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

aw3567: added 10cm depth measurement
All checks were successful
FarmMaps.Develop/CarbonLib/pipeline/head This commit looks good
Details

Browse Source
This commit is contained in:
Francisco Salas
2022-02-25 11:41:32 +01:00
parent 1d9d213514
commit 597f42001b
7 changed files with 85 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

75
CarbonCalculation.cs
View File

@@ -139,7 +139,7 @@ namespace FarmMapsAPI.Carbon
Area = cropField.Area,
HistoricalCropData = mappedCropData,
SoilType = MapSoil(cropField.SoilType),
SoilProperty = new SoilProperty { Clay_Content = cropField.SoilProperty.Clay_Content, OM_Const = cropField.SoilProperty.OM_Const },
SoilProperty = new SoilProperty { Clay_Content = cropField.SoilProperty.Clay_Content, OM_Const = cropField.SoilProperty.OM_Const, Depth10 = cropField.SoilProperty.Depth10 },
GreenManures = greenManures
});
foreach (var crops in cropField.HistoricalCropData)
@@ -242,6 +242,7 @@ namespace FarmMapsAPI.Carbon
// cropAreaSheet (external input)
var cropAreaField = GetCropAreaField(farmInput);
var year = GetYears(farmInputSet, farmInput);
//////////////////////////////////////////////
// precalculations
@@ -251,7 +252,7 @@ namespace FarmMapsAPI.Carbon
//soilPropSoilType
var soilPropSoilType = GetSoilPropSoilType(soilTypeData, soilPropertiesField);
//soilPropertiesField (again)
//soilPropertiesField (again) --> remove this
foreach (var item in soilPropertiesField)
{
var soilProps = soilPropSoilType.Keys.Where(x => x.Item1 == item.Key).ToList();
@@ -263,21 +264,64 @@ namespace FarmMapsAPI.Carbon
item.Value.OM_Corr = sumCorr;
}
//var oM_corr_SoilType = new Dictionary<string, List<ColumnValue>>();
var oM_corr_SoilType = new Dictionary<(string,string), double>();
foreach (var item in soilPropertiesField)
{
var soilOfItem = soilTypeData[item.Key].FirstOrDefault();
var initialSoilKey = historicalCropdata.Keys.FirstOrDefault(x => x.Item1 == item.Key && x.Item2 == year.First());
var initialSoil = historicalCropdata[initialSoilKey].FirstOrDefault();
var soilPropertyDepth10 = item.Value.Depth10;
if (initialSoil == "Temp_gras" && soilPropertyDepth10)
{
oM_corr_SoilType.Add((item.Key, soilOfItem), item.Value.OM_Cont * 0.97);
}
else if (cropSet.Grass_nat.Contains(initialSoil) && soilPropertyDepth10 && (soilOfItem == "Klei" || soilOfItem == "Zand"))
{
if (soilOfItem == "Klei")
{
oM_corr_SoilType.Add((item.Key, soilOfItem), item.Value.OM_Cont * 0.67);
}
else
{
oM_corr_SoilType.Add((item.Key, soilOfItem), item.Value.OM_Cont * 0.81);
}
}
else
{
oM_corr_SoilType.Add((item.Key, soilOfItem), item.Value.OM_Cont);
}
}
var oM_Corr = new Dictionary<string, double>();
foreach (var item in oM_corr_SoilType)
{
oM_Corr.Add(item.Key.Item1, item.Value);
}
// soilBulkDensities
Dictionary<string, SoilBulkDensityType> soilBulkDensities = GetSoilBulkDensities(soilPropertiesField, soilTypeData, soilPropSoilType);
Dictionary<string, SoilBulkDensityType> soilBulkDensities = GetSoilBulkDensities(soilPropertiesField, soilTypeData, soilPropSoilType, oM_corr_SoilType);
foreach (var soilBulkDensity in soilBulkDensities)
{
if (soilPropSoilType.ContainsKey((soilBulkDensity.Key, "Klei")))
{
soilPropSoilType[(soilBulkDensity.Key, "Klei")].Bulk_Density = (1.0 / soilBulkDensity.Value.Clay);
soilBulkDensity.Value.Clay = (1.0 / soilBulkDensity.Value.Clay);
}
if (soilPropSoilType.ContainsKey((soilBulkDensity.Key, "Zand")))
{
soilPropSoilType[(soilBulkDensity.Key, "Zand")].Bulk_Density = (1.0 / soilBulkDensity.Value.Sand);
soilBulkDensity.Value.Sand = (1.0 / soilBulkDensity.Value.Sand);
}
}
var bulkDens = new Dictionary<string, double>();
foreach (var soilBulkDensity in soilBulkDensities)
{
bulkDens.Add(soilBulkDensity.Key, (soilBulkDensity.Value.Clay + soilBulkDensity.Value.Sand));
}
foreach (var item in soilPropertiesField)
{
var soilProps = soilPropSoilType.Keys.Where(x => x.Item1 == item.Key).ToList();
@@ -295,8 +339,7 @@ namespace FarmMapsAPI.Carbon
// loop rothc
//////////////////////////////////////////
var soc_Field = new Dictionary<string, List<YearColumnValue>>();
var year = GetYears(farmInputSet, farmInput);
var soc_Field = new Dictionary<string, List<YearColumnValue>>();
var cropAreaField1 = new Dictionary<(string, int), List<ColumnValue>>();
var carbonInputs = new GamsThreeKeyParameter();
@@ -798,7 +841,7 @@ namespace FarmMapsAPI.Carbon
foreach (var field in cropAreaField1.Where(x => x.Key.Item2 == loopYear))
{
var valueList = new List<ColumnValue>();
var result = parameterSet.LayerDepth * soilPropertiesField[field.Key.Item1].Bulk_Density * soilPropertiesField[field.Key.Item1].OM_Corr * 0.54;
var result = parameterSet.LayerDepth * bulkDens[field.Key.Item1] * oM_Corr[field.Key.Item1] * 0.54;
foreach (var crop in field.Value)
{
valueList.Add(new ColumnValue(crop.Key, result));
@@ -1544,18 +1587,18 @@ namespace FarmMapsAPI.Carbon
summary.Initial = sum;
}
foreach (var item in soilPropertiesField)
foreach(var item in bulkDens)
{
var result = item.Value.Bulk_Density;
var result = item.Value;
foreach (var lYear in allYears)
{
soc_Summary[(item.Key, lYear)].Add(new ColumnValue("Bulk_Density", result));
}
}
foreach (var item in soilPropertiesField)
foreach (var item in bulkDens)
{
var result = item.Value.Bulk_Density;
var result = item.Value;
foreach (var lYear in allYears)
{
var summaries = cropFieldSummary.Single(x => x.Name == item.Key);
@@ -1729,7 +1772,7 @@ namespace FarmMapsAPI.Carbon
var soilPropertiesField = new Dictionary<string, SoilPropertyType>();
foreach (var field in farmInput.CropFields)
{
soilPropertiesField.Add(field.Name, new SoilPropertyType { OM_Cont = field.SoilProperty.OM_Const.Value, Clay_Content = field.SoilProperty.Clay_Content.Value, OM_Corr = field.SoilProperty.OM_Const.Value <= 20 ? field.SoilProperty.OM_Const.Value : 0 });
soilPropertiesField.Add(field.Name, new SoilPropertyType { OM_Cont = field.SoilProperty.OM_Const.Value, Clay_Content = field.SoilProperty.Clay_Content.Value, OM_Corr = field.SoilProperty.OM_Const.Value <= 20 ? field.SoilProperty.OM_Const.Value : 0, Depth10 = field.SoilProperty.Depth10.Value });
}
return soilPropertiesField;
@@ -1772,7 +1815,7 @@ namespace FarmMapsAPI.Carbon
return cropAreaField;
}
private Dictionary<string, SoilBulkDensityType> GetSoilBulkDensities(Dictionary<string, SoilPropertyType> soilPropertiesField, Dictionary<string, List<string>> soilTypeData, Dictionary<(string, string), SoilPropertyType> soilPropSoilType)
private Dictionary<string, SoilBulkDensityType> GetSoilBulkDensities(Dictionary<string, SoilPropertyType> soilPropertiesField, Dictionary<string, List<string>> soilTypeData, Dictionary<(string, string), SoilPropertyType> soilPropSoilType, Dictionary<(string, string), double> oM_corr_SoilType)
{
var soilBulkDensities = new Dictionary<string, SoilBulkDensityType>();
foreach (var item in soilPropertiesField)
@@ -1781,15 +1824,17 @@ namespace FarmMapsAPI.Carbon
if (soilTypeData[item.Key].Any(x => x == "Klei"))
{
var soilBulkDensityKey = soilPropSoilType.Keys.Single(x => x.Item1 == item.Key && x.Item2 == "Klei");
var clay = (0.6117 + (0.003601 * item.Value.Clay_Content)) +
(0.002172 * soilPropSoilType[soilBulkDensityKey].OM_Corr * soilPropSoilType[soilBulkDensityKey].OM_Corr) +
(0.01715 * Math.Log(soilPropSoilType[soilBulkDensityKey].OM_Corr));
(0.002172 * oM_corr_SoilType[(item.Key, "Klei")] * oM_corr_SoilType[(item.Key, "Klei")]) +
(0.01715 * Math.Log(oM_corr_SoilType[(item.Key, "Klei")]));
soilBulkDensities[item.Key].Clay = clay;
}
if (soilTypeData[item.Key].Any(x => x == "Zand"))
{
var soilBulkDensityKey = soilPropSoilType.Keys.Single(x => x.Item1 == item.Key && x.Item2 == "Zand");
soilBulkDensities[item.Key].Sand = 0.667 + (0.021 * soilPropSoilType[soilBulkDensityKey].OM_Corr);
soilBulkDensities[item.Key].Sand = 0.667 + (0.021 * oM_corr_SoilType[(item.Key, "Zand")]);
}
}
return soilBulkDensities;