// This routine will typically be triggered after soil tests for a particular field have been updated
// This routine will recalculate for all crops in a field the nutrients all that are dependant on the soil tests
public void RecalcCropsSoilTestMessagesByField(string fieldName)
{
CalculateCropRequirementRemoval ccrr = new CalculateCropRequirementRemoval(_ud, _sd);
//iterate through the crops and update the crop requirements
List <FieldCrop> fieldCrops = _ud.GetFieldCrops(fieldName);
if (fieldCrops.Count > 0)
{
foreach (var _crop in fieldCrops)
{
CropType crpTyp = new CropType();
FieldCrop cf = _ud.GetFieldCrop(fieldName, _crop.id);
if (cf.cropId != null)
{
Crop cp = _sd.GetCrop(Convert.ToInt32(cf.cropId));
crpTyp = _sd.GetCropType(cp.croptypeid);
}
else
{
crpTyp.modifynitrogen = false;
}
CropRequirementRemoval crr = new CropRequirementRemoval();
ccrr.cropid = Convert.ToInt16(_crop.cropId);
ccrr.previousCropid = _crop.prevCropId;
ccrr.yield = _crop.yield;
ccrr.crudeProtien = _crop.crudeProtien;
ccrr.coverCropHarvested = _crop.coverCropHarvested;
ccrr.fieldName = fieldName;
crr = ccrr.GetCropRequirementRemoval();
if (!crpTyp.modifynitrogen)
{
cf.reqN = crr.N_Requirement;
}
cf.reqP2o5 = crr.P2O5_Requirement;
cf.reqK2o = crr.K2O_Requirement;
cf.remN = crr.N_Removal;
cf.remP2o5 = crr.P2O5_Removal;
cf.remK2o = crr.K2O_Removal;
_ud.UpdateFieldCrop(fieldName, cf);
}
}
}
public CropRequirementRemoval GetCropRequirementRemoval()
{
decimal n_removal = 0;
ConversionFactor _cf = _sd.GetConversionFactor();
CropRequirementRemoval crr = new CropRequirementRemoval();
Crop crop = _sd.GetCrop(cropid);
// For testing we're using Soil Test Kelowna of 65ppm, should actually by 100+ (say 101)
//
// Nutrient removal
// P2O5 = crop.cropremovalfactor_P2O5 * yield
// K2O = crop.cropremovalfactor_K2O * yield
// N = crop.cropremovalfactor_N * yield
// if crude_protien exists and is not default
// crop.cropremovalfactor_N = (crude protien / (0.625 * 0.5))
// N = crop.cropremovalfactor_N * yield
// else
// N = crop.cropremovalfactor_N * yield
//
// Note for Cover crops (only)
// if Cover Crop harvested
// don't change numbers
// if Cover crop not harvested
// set all removal amts to zero
if (!crudeProtien.HasValue || (crudeProtien.HasValue && crudeProtien.Value == 0))
{
decimal tmpDec;
if (decimal.TryParse(crop.cropremovalfactor_N.ToString(), out tmpDec))
{
n_removal = tmpDec * yield;
}
else
{
n_removal = 0;
}
}
else
{
n_removal = decimal.Divide(Convert.ToDecimal(crudeProtien), (_cf.n_protein_conversion * _cf.unit_conversion)) * yield;
}
crr.P2O5_Removal = Convert.ToInt32(crop.cropremovalfactor_P2O5 * yield);
crr.K2O_Removal = Convert.ToInt32(crop.cropremovalfactor_K2O * yield);
crr.N_Removal = Convert.ToInt32(n_removal);
if (coverCropHarvested.HasValue && coverCropHarvested.Value == false)
{
crr.P2O5_Removal = 0;
crr.K2O_Removal = 0;
crr.N_Removal = 0;
}
// Crop Requirement
// P205
// get region.soil_test_phospherous_region_cd
// get phosphorous_crop_group_region_cd = crop_stp_regioncd(cropid, region.soil_test_phospherous_region_cd)
// get stp_kelowna_range.id usign default ST Kelowna (65) between range_low and range_high
// get P2O5 recommedation = stp_recommend(stp_kelowna_range.id, region.soil_test_phospherous_region_cd, phosphorous_crop_group_region_cd)
//
// same as above for K2O
//
// For N
// Look up crop.n_recommcd
// if 0 or 1
// N = crop.n_recomm_lbperac
// if 2
// same as N removal
// if 3
// get default yield = cropyield(cropid, locationid)
// N = (yield / default yield) * crop.n_recomm_lbperac
int regionid = _ud.FarmDetails().farmRegion.Value;
Region region = _sd.GetRegion(regionid);
Field fld = _ud.GetFieldDetails(fieldName);
if (fld.soilTest == null)
{
fld.soilTest = new SoilTest();
DefaultSoilTest dt = _sd.GetDefaultSoilTest();
fld.soilTest.valNO3H = dt.nitrogen;
fld.soilTest.ValP = dt.phosphorous;
fld.soilTest.valK = dt.potassium;
fld.soilTest.valPH = dt.pH;
fld.soilTest.ConvertedKelownaK = dt.convertedKelownaK;
fld.soilTest.ConvertedKelownaP = dt.convertedKelownaP;
}
int _STP = fld.soilTest.ConvertedKelownaP;
if (_STP == 0)
{
_STP = _cf.defaultSoilTestKelownaP;
}
int _STK = fld.soilTest.ConvertedKelownaK;
if (fld.soilTest.ConvertedKelownaK == 0)
{
_STK = _cf.defaultSoilTestKelownaK;
}
// p2o5 recommend calculations
CropSTPRegionCd cropSTPRegionCd = _sd.GetCropSTPRegionCd(cropid, region.soil_test_phospherous_region_cd);
int? phosphorous_crop_group_region_cd = cropSTPRegionCd.phosphorous_crop_group_region_cd;
STPKelownaRange sTPKelownaRange = _sd.GetSTPKelownaRangeByPpm(_STP);
int stp_kelowna_range_id = sTPKelownaRange.id;
if (phosphorous_crop_group_region_cd == null)
{
crr.P2O5_Requirement = 0;
}
else
{
STPRecommend sTPRecommend = _sd.GetSTPRecommend(stp_kelowna_range_id, region.soil_test_phospherous_region_cd, Convert.ToInt16(phosphorous_crop_group_region_cd));
crr.P2O5_Requirement = Convert.ToInt32(Convert.ToDecimal(sTPRecommend.p2o5_recommend_kgperha) * _cf.kgperha_lbperac_conversion);
}
// k2o recommend calculations
CropSTKRegionCd cropSTKRegionCd = _sd.GetCropSTKRegionCd(cropid, region.soil_test_potassium_region_cd);
int? potassium_crop_group_region_cd = cropSTKRegionCd.potassium_crop_group_region_cd;
STKKelownaRange sTKKelownaRange = _sd.GetSTKKelownaRangeByPpm(_STK);
int stk_kelowna_range_id = sTKKelownaRange.id;
if (potassium_crop_group_region_cd == null)
{
crr.K2O_Requirement = 0;
}
else
{
STKRecommend sTKRecommend = _sd.GetSTKRecommend(stk_kelowna_range_id, region.soil_test_potassium_region_cd, Convert.ToInt16(potassium_crop_group_region_cd));
crr.K2O_Requirement = Convert.ToInt32(Convert.ToDecimal(sTKRecommend.k2o_recommend_kgperha) * _cf.kgperha_lbperac_conversion);
}
// n recommend calculations -note the excel n_recommd are zero based, the static data is 1 based
switch (crop.n_recommcd)
{
case 1:
crr.N_Requirement = Convert.ToInt16(crop.n_recomm_lbperac);
break;
case 2:
crr.N_Requirement = Convert.ToInt16(crop.n_recomm_lbperac);
break;
case 3:
crr.N_Requirement = crr.N_Removal;
break;
case 4:
CropYield cropYield = _sd.GetCropYield(cropid, region.locationid);
if (cropYield.amt != null)
{
crr.N_Requirement = Convert.ToInt16(decimal.Divide(yield, Convert.ToDecimal(cropYield.amt)) * crop.n_recomm_lbperac);
}
else
{
crr.N_Requirement = 0;
}
break;
}
// if a previous crop has been ploughed dowm account for the N in the field (passed in as a credit)
crr.N_Requirement = crr.N_Requirement - nCredit;
// only reduce to 0
crr.N_Requirement = (crr.N_Requirement < 0) ? 0 : crr.N_Requirement;
return(crr);
}
