DK_EnergyCrop_Perm Class Reference

DK_EnergyCrop_Perm class
. More...

#include <DK_EnergyCrop_Perm.h>

Inheritance diagram for DK_EnergyCrop_Perm:

Public Member Functions
virtual bool	Do (Farm *a_farm, LE *a_field, FarmEvent *a_ev)
	The one and only method for a crop management plan. All farm actions go through here. More...
	DK_EnergyCrop_Perm (TTypesOfVegetation a_tov, TTypesOfCrops a_toc, Landscape *a_L)
void	SetUpFarmCategoryInformation ()
Public Member Functions inherited from Crop
virtual	~Crop ()
	Crop (TTypesOfVegetation a_tov, TTypesOfCrops a_toc, Landscape *a_L)
int	GetFirstDate (void)
void	ChooseNextCrop (int a_no_dates)
	Chooses the next crop to grow in a field. More...
int	GetCropClassification ()
void	SetCropClassification (int a_classification)
void	GetCropType (TTypesOfCrops a_toc)
TTypesOfCrops	GetCropType ()
FarmManagementCategory	LookUpManagementCategory (int a_todo)
bool	GetForceSpringOK ()

Additional Inherited Members
Public Attributes inherited from Crop
TTypesOfVegetation	m_tov
Protected Member Functions inherited from Crop
void	SimpleEvent (long a_date, int a_todo, bool a_lock)
	Adds an event to this crop management. More...
void	SimpleEvent_ (long a_date, int a_todo, bool a_lock, Farm *a_farm, LE *a_field)
	Adds an event to this crop management without relying on member variables. More...
bool	StartUpCrop (int a_spring, std::vector< std::vector< int >> a_flexdates, int a_todo)
	Holds the translation between the farm operation enum for each cropand the farm management category associated with it More...
bool	AphidDamage (LE *a_field)
	Compares aphid numbers per m2 with a threshold to return true if threshold is exceeded. More...
Protected Attributes inherited from Crop
Farm *	m_farm
LE *	m_field
FarmEvent *	m_ev
int	m_first_date
int	m_count
int	m_last_date
int	m_ddegstoharvest
int	m_base_elements_no
Landscape *	m_OurLandscape
bool	m_forcespringpossible = false
	Used to signal that the crop can be forced to start in spring. More...
TTypesOfCrops	m_toc
	The Crop type in terms of the TTypesOfCrops list (smaller list than tov, no country designation) More...
int	m_CropClassification
	Contains information on whether this is a winter crop, spring crop, or catch crop that straddles the year boundary (0,1,2) More...
vector< FarmManagementCategory >	m_ManagementCategories
	Holds the translation between the farm operation enum for each crop and the farm management category associated with it. More...
Static Protected Attributes inherited from Crop
static int	m_date_modifier = 0
	Holds a value that shifts test pesticide use by this many days in crops modified to use it. More...

Detailed Description

DK_EnergyCrop_Perm class
.

See DK_EnergyCrop.h::DK_EnergyCropToDo for a complete list of all possible events triggered codes by the management plan. When triggered these events are handled by Farm and are available as information for other objects such as animal and bird models.

Constructor & Destructor Documentation

DK_EnergyCrop_Perm()

DK_EnergyCrop_Perm::DK_EnergyCrop_Perm ( TTypesOfVegetation  a_tov, TTypesOfCrops  a_toc, Landscape *  a_L  )

inline

                                                                                     : Crop(a_tov, a_toc, a_L)
   {
        // When we start it off, the first possible date for a farm operation is ...
        // This information is used by other crops when they decide how much post processing of 
        // the management is allowed after harvest before the next crop starts.
        m_first_date=g_date->DayInYear( 1,12 );
        SetUpFarmCategoryInformation();
   }

References Calendar::DayInYear(), g_date, Crop::m_first_date, and SetUpFarmCategoryInformation().

Member Function Documentation

Do()

bool DK_EnergyCrop_Perm::Do ( Farm *  a_farm, LE *  a_field, FarmEvent *  a_ev  )

virtual

The one and only method for a crop management plan. All farm actions go through here.

Called every time something is done to the crop by the farmer in the first instance it is always called with a_ev->todo set to start, but susequently will be called whenever the farmer wants to carry out a new operation.
This method details all the management and relationships between operations necessary to grow and ALMaSS crop.

Reimplemented from Crop.

{
    m_farm = a_farm;
    m_field = a_field;
    m_ev = a_ev;
    bool done = false; // The boolean value done indicates when we are totally finished with this plan (i.e. it is set to true).
    int d1 = 0;
    TTypesOfVegetation l_tov = tov_DKEnergyCrop_Perm;
    // Depending what event has occured jump to the correct bit of code
    switch (m_ev->m_todo)
    {
    case dk_ec_start:
    {
        DK_EC_YEARS_AFTER_EST = false;
        // dk_ec_start just sets up all the starting conditions and reference dates that are needed to start a dk_bfp1
 
        a_field->ClearManagementActionSum();
        m_last_date = g_date->DayInYear(15, 10); // Should match the last flexdate below
           //Create a 2d array of 1 plus the number of operations you use. Change only 4+1 to what you need in the line below
        std::vector<std::vector<int>> flexdates(3 + 1, std::vector<int>(2, 0));
        // Set up the date management stuff
                // Start and stop dates for all events after harvest
        flexdates[0][1] = g_date->DayInYear(30, 9); // last possible day of harvest
        // Now these are done in pairs, start & end for each operation. If its not used then -1
        flexdates[1][0] = -1; // This date will be moved back as far as necessary and potentially to flexdates 1 (start op 1)
        flexdates[1][1] = g_date->DayInYear(15, 10); // This date will be moved back as far as necessary and potentially to flexdates 1 (end op 1) herbi
        flexdates[2][0] = -1; // This date will be moved back as far as necessary and potentially to flexdates 2 (start op 2)
        flexdates[2][1] = g_date->DayInYear(15, 10); // This date will be moved back as far as necessary and potentially to flexdates 2 (end op 2) ferti
        flexdates[3][0] = -1; // This date will be moved back as far as necessary and potentially to flexdates 3 (start op 3)
        flexdates[3][1] = g_date->DayInYear(15, 10); // This date will be moved back as far as necessary and potentially to flexdates 3 (end op 3) ferti
 
        // Below if this is a spring crop use 365, otherwise first parameter is always 0, second parameter is fixed, and the third is the start up operation in the first year
        int isSpring = 365;
        if (StartUpCrop(isSpring, flexdates, int(dk_ec_deep_plough_sand))) break;
 
        // End single block date checking code. Please see next line comment as well.
        // Reinit d1 to first possible starting date.
    // Here we queue up the first event
        if ((DK_EC_YEARS_AFTER_EST + g_date->GetYearNumber()) % 1 == 1)
        {
            d1 = g_date->OldDays() + g_date->DayInYear(25, 2) + isSpring;
            SimpleEvent(d1, dk_ec_deep_plough_sand, false);
        }
        else if ((DK_EC_YEARS_AFTER_EST + g_date->GetYearNumber()) % 1 == 0)
        {
            d1 = g_date->OldDays() + g_date->DayInYear(1, 3) + isSpring;
            SimpleEvent(d1, dk_ec_cut_trees, false);
        }
        break;
    }
    break;
 
    // LKM: This is the first real farm operation - done if much seed weed or setaside field
    case dk_ec_deep_plough_sand:
        if (m_field->GetSoilType() == tos_Sand || m_field->GetSoilType() == tos_LoamySand || m_field->GetSoilType() == tos_SandyLoam || m_field->GetSoilType() == tos_SandyClayLoam) // on sandy soils (NL ZAND & LOSS)
        {
            if (a_field->GetVegBiomass() > 0)
            {
                if (m_ev->m_lock || m_farm->DoIt_prob(1.0)) {
                    if (!m_farm->DeepPlough(m_field, 0.0, g_date->DayInYear(10, 3) - g_date->DayInYear())) {
                        SimpleEvent(g_date->Date() + 1, dk_ec_deep_plough_sand, true);
                        break;
                    }
                }
            }
            SimpleEvent(g_date->Date(), dk_ec_roll1_sand, false);
            break;
        }
        else SimpleEvent(g_date->OldDays() + g_date->DayInYear(25, 2), dk_ec_deep_plough_clay, false);
        break;
 
    case dk_ec_roll1_sand:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->SpringRoll(m_field, 0.0, g_date->DayInYear(12, 3) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_roll1_sand, true);
                break;
            }
        }
        SimpleEvent(g_date->Date() + 1, dk_ec_slurry1_s_sand, false);
        break;
 
    case dk_ec_slurry1_s_sand:
        if (m_farm->IsStockFarmer()) { //Stock Farmer
            if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
            {
                if (!a_farm->FA_Slurry(m_field, 0.0, g_date->DayInYear(14, 3) - g_date->DayInYear())) {
                    SimpleEvent(g_date->Date() + 1, dk_ec_slurry1_s_sand, true);
                    break;
                }
            }
            SimpleEvent(g_date->Date() + 1, dk_ec_planting_sand, false);
            break;
        }
        else
        SimpleEvent(g_date->Date(), dk_ec_npk1_p_sand, false);
        break;
 
    case dk_ec_npk1_p_sand:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->FP_NPK(m_field, 0.0, g_date->DayInYear(14, 3) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_npk1_p_sand, true);
                break;
            }
        }
        SimpleEvent(g_date->Date() +1, dk_ec_planting_sand, false);
        break;
 
    case dk_ec_planting_sand:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->SpringSow(m_field, 0.0, g_date->DayInYear(15, 3) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_planting_sand, true);
                break;
            }
        }
        SimpleEvent(g_date->Date() + 1, dk_ec_roll2, false);
        break;
 
    case dk_ec_deep_plough_clay:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->DeepPlough(m_field, 0.0, g_date->DayInYear(25, 4) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_deep_plough_clay, true);
                break;
            }
        }
        SimpleEvent(g_date->Date(), dk_ec_roll1_clay, false);
        break;
 
    case dk_ec_roll1_clay:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->SpringRoll(m_field, 0.0, g_date->DayInYear(27, 4) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_roll1_clay, true);
                break;
            }
        }
        SimpleEvent(g_date->Date() + 1, dk_ec_slurry1_s_clay, false);
        break;
 
    case dk_ec_slurry1_s_clay:
        if (m_farm->IsStockFarmer()) { //Stock Farmer
            if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
            {
                if (!a_farm->FA_Slurry(m_field, 0.0, g_date->DayInYear(29, 4) - g_date->DayInYear())) {
                    SimpleEvent(g_date->Date() + 1, dk_ec_slurry1_s_clay, true);
                    break;
                }
            }
            SimpleEvent(g_date->Date() + 1, dk_ec_planting_clay, false);
            break;
        }
        else
            SimpleEvent(g_date->Date(), dk_ec_npk1_p_clay, false);
        break;
 
    case dk_ec_npk1_p_clay:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->FP_NPK(m_field, 0.0, g_date->DayInYear(29, 4) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_npk1_p_clay, true);
                break;
            }
        }
        SimpleEvent(g_date->Date() + 1, dk_ec_planting_clay, false);
        break;
 
    case dk_ec_planting_clay:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->SpringSow(m_field, 0.0, g_date->DayInYear(31, 4) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_planting_clay, true);
                break;
            }
        }
        SimpleEvent(g_date->Date() + 1, dk_ec_roll2, false);
        break;
 
    case dk_ec_roll2:
        if (!a_farm->SpringRoll(m_field, 0.0, g_date->DayInYear(1, 5) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_roll2, true);
            break;
        }
 
        // if time then herbicide before May, otherwise after May (not when budding in ~May)
        d1 = g_date->Date();
        if (d1 < g_date->OldDays() + g_date->DayInYear(30, 4)) {
            SimpleEvent(g_date->Date()+1, dk_ec_herbicide1, false);
            break;
        }
        else SimpleEvent(g_date->OldDays() + g_date->DayInYear(1, 6), dk_ec_herbicide1, false);
        break;
 
    case dk_ec_herbicide1:
        if (!a_farm->HerbicideTreat(m_field, 0.0, g_date->DayInYear(30, 6) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_herbicide1, true);
            break;
        }
 
        SimpleEvent(g_date->Date()+1, dk_ec_shallowharrow1, false);
        break;
    case dk_ec_shallowharrow1:      
        if (!a_farm->ShallowHarrow(m_field, 0.0, g_date->DayInYear(1, 8) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_shallowharrow1, true);
            break;
        }
        SimpleEvent(g_date->Date()+1, dk_ec_shallowharrow2, false);
        break;
    case dk_ec_shallowharrow2: 
        if (!a_farm->ShallowHarrow(m_field, 0.0, g_date->DayInYear(2, 8) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_shallowharrow2, true);
            break;
        }
 
        // After first growth season
        SimpleEvent(g_date->OldDays() + g_date->DayInYear(1, 3)+365, dk_ec_cut_trees, false);
        break;
 
    case dk_ec_cut_trees:
        if (!a_farm->Pruning(m_field, 0.0, g_date->DayInYear(15, 3) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_cut_trees, true);
            break;
        }
 
        // here comes a fork of parallel events:
        SimpleEvent(g_date->Date()+1, dk_ec_herbicide2, false); // herbi thread
        SimpleEvent(g_date->Date() + 1, dk_ec_slurry2_s, false); // fertilizer thread
        SimpleEvent(g_date->Date() + 1, dk_ec_shallowharrow3, false); // cultivation thread - main thread
        break;
 
    case dk_ec_herbicide2:
        if (!a_farm->HerbicideTreat(m_field, 0.0, g_date->DayInYear(31, 3) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_herbicide2, true);
            break;
        }
        // if time then herbicide before May, otherwise after May (not when budding in ~May)
        d1 = g_date->Date();
        if (d1 < g_date->OldDays() + g_date->DayInYear(30, 4)) {
            SimpleEvent(g_date->Date() + 1, dk_ec_herbicide3, false);
            break;
        }
        else SimpleEvent(g_date->OldDays() + g_date->DayInYear(1, 6), dk_ec_herbicide3, false);
        break;
 
    case dk_ec_herbicide3:
        if (!a_farm->HerbicideTreat(m_field, 0.0, g_date->DayInYear(30, 6) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_herbicide3, true);
            break;
        }
        break; // end of thread
 
    case dk_ec_slurry2_s:
        if (m_farm->IsStockFarmer()) { //Stock Farmer
            if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
            {
                if (!a_farm->FA_Slurry(m_field, 0.0, g_date->DayInYear(31, 3) - g_date->DayInYear())) {
                    SimpleEvent(g_date->Date() + 1, dk_ec_slurry2_s, true);
                    break;
                }
            }
            break; //end of thread
        }
        else
        SimpleEvent(g_date->Date(), dk_ec_npk2_p, false);
        break;
 
    case dk_ec_npk2_p:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->FP_NPK(m_field, 0.0, g_date->DayInYear(31, 3) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_npk2_p, true);
                break;
            }
        }
        break; //end of thread
 
    case dk_ec_shallowharrow3:
        if (!a_farm->ShallowHarrow(m_field, 0.0, g_date->DayInYear(1, 8) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_shallowharrow3, true);
            break;
        }
        SimpleEvent(g_date->Date() + 1, dk_ec_shallowharrow4, false);
        break;
 
    case dk_ec_shallowharrow4:
        if (!a_farm->ShallowHarrow(m_field, 0.0, g_date->DayInYear(2, 8) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_shallowharrow4, true);
            break;
        }
        SimpleEvent(g_date->OldDays() + g_date->DayInYear(1, 8), dk_ec_harvest, false);
        break;
 
    case dk_ec_harvest:
        if (!a_farm->Harvest(m_field, 0.0, m_field->GetMDates(1, 0) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_harvest, true);
            break;
        }
 
        //fork with parallel events:
        SimpleEvent(g_date->Date() + 1, dk_ec_herbicide4, false); // herbi thread
        SimpleEvent(g_date->Date() + 1, dk_ec_slurry3_s, false); // fertilizer thread - main thread
        break;
 
    case dk_ec_herbicide4:
        if (!a_farm->HerbicideTreat(m_field, 0.0, m_field->GetMDates(1, 1) - g_date->DayInYear())) {
            SimpleEvent(g_date->Date() + 1, dk_ec_herbicide4, true);
            break;
        }
        break; // end of thread
 
    case dk_ec_slurry3_s:
        if (m_farm->IsStockFarmer()) { //Stock Farmer
            if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
            {
                if (!a_farm->FA_Slurry(m_field, 0.0, m_field->GetMDates(1, 2) - g_date->DayInYear())) {
                    SimpleEvent(g_date->Date() + 1, dk_ec_slurry3_s, true);
                    break;
                }
            }
            done = true;
            break;
        }
        else
        SimpleEvent(g_date->Date(), dk_ec_npk3_p, false);
        break;
 
    case dk_ec_npk3_p:
        if (m_ev->m_lock || m_farm->DoIt_prob(1.00))
        {
            if (!a_farm->FP_NPK(m_field, 0.0, m_field->GetMDates(1, 3) - g_date->DayInYear())) {
                SimpleEvent(g_date->Date() + 1, dk_ec_npk3_p, true);
                break;
            }
        }
        done = true;
        break;
        // So we are done, and somewhere else the farmer will queue up the start event of the next crop
        // END of MAIN THREAD
        break;
        default:
        g_msg->Warn(WARN_BUG, "DK_EnergyCrop_Perm::Do(): "
            "Unknown event type! ", "");
        exit(1);
    }
    return done;
}

References LE::ClearManagementActionSum(), Calendar::Date(), Calendar::DayInYear(), Farm::DeepPlough(), dk_ec_cut_trees, dk_ec_deep_plough_clay, dk_ec_deep_plough_sand, dk_ec_harvest, dk_ec_herbicide1, dk_ec_herbicide2, dk_ec_herbicide3, dk_ec_herbicide4, dk_ec_npk1_p_clay, dk_ec_npk1_p_sand, dk_ec_npk2_p, dk_ec_npk3_p, dk_ec_planting_clay, dk_ec_planting_sand, dk_ec_roll1_clay, dk_ec_roll1_sand, dk_ec_roll2, dk_ec_shallowharrow1, dk_ec_shallowharrow2, dk_ec_shallowharrow3, dk_ec_shallowharrow4, dk_ec_slurry1_s_clay, dk_ec_slurry1_s_sand, dk_ec_slurry2_s, dk_ec_slurry3_s, dk_ec_start, DK_EC_YEARS_AFTER_EST, Farm::DoIt_prob(), Farm::FA_Slurry(), Farm::FP_NPK(), g_date, g_msg, LE::GetMDates(), LE::GetSoilType(), LE::GetVegBiomass(), Calendar::GetYearNumber(), Farm::Harvest(), Farm::HerbicideTreat(), Farm::IsStockFarmer(), Crop::m_ev, Crop::m_farm, Crop::m_field, Crop::m_last_date, FarmEvent::m_lock, FarmEvent::m_todo, Calendar::OldDays(), Farm::Pruning(), Farm::ShallowHarrow(), Crop::SimpleEvent(), Farm::SpringRoll(), Farm::SpringSow(), Crop::StartUpCrop(), tos_LoamySand, tos_Sand, tos_SandyClayLoam, tos_SandyLoam, tov_DKEnergyCrop_Perm, MapErrorMsg::Warn(), and WARN_BUG.

SetUpFarmCategoryInformation()

void DK_EnergyCrop_Perm::SetUpFarmCategoryInformation ( )

inline

                                       {
       const int elements = 2 + (dk_ec_foobar - DK_EC_BASE);
       m_base_elements_no = DK_EC_BASE - 2;
 
       FarmManagementCategory catlist[elements] =
       {
            fmc_Others, // zero element unused but must be here
            fmc_Others, // dk_ec_start = 1, // Compulsory, must always be 1 (one).
            fmc_Others, // dk_ec_sleep_all_day = // DK_EC_BASE,
            fmc_Cultivation,    // dk_ec_deep_plough_sand,
            fmc_Cultivation,    //dk_ec_roll1_sand,
            fmc_Fertilizer, // dk_ec_slurry1_s_sand,
            fmc_Fertilizer, // dk_ec_npk1_p_sand,
            fmc_Cultivation,    // dk_ec_deep_plough_clay,
            fmc_Cultivation,    //dk_ec_roll1_clay,
            fmc_Fertilizer, // dk_ec_slurry1_s_clay,
            fmc_Fertilizer, // dk_ec_npk1_p_clay,
            fmc_Others,     //dk_ec_planting_sand,
            fmc_Others,     //dk_ec_planting_clay,
            fmc_Cultivation,    //dk_ec_roll2,
            fmc_Herbicide,  // dk_ec_herbicide1,
            fmc_Cultivation,    // dk_ec_shallow_harrow1,
            fmc_Cultivation,    // dk_ec_shallow_harrow2,
            fmc_Cutting,    // dk_ec_cut_trees, 
            fmc_Herbicide,  // dk_ec_herbicide2,
            fmc_Fertilizer, // dk_ec_slurry2_s,
            fmc_Fertilizer, // dk_ec_npk2_p,
            fmc_Herbicide,  // dk_ec_herbicide3,
            fmc_Cultivation,    // dk_ec_shallow_harrow3,
            fmc_Cultivation,    // dk_ec_shallow_harrow4,
            fmc_Harvest,    // dk_ec_harvest,
            fmc_Herbicide,  // dk_ec_herbicide4,
            fmc_Fertilizer, // dk_ec_slurry3_s,
            fmc_Fertilizer, // dk_ec_npk3_p,
                // no foobar entry
 
       };
       // Iterate over the catlist elements and copy them to vector
       copy(begin(catlist), end(catlist), back_inserter(m_ManagementCategories));
   }

References DK_EC_BASE, dk_ec_foobar, fmc_Cultivation, fmc_Cutting, fmc_Fertilizer, fmc_Harvest, fmc_Herbicide, fmc_Others, Crop::m_base_elements_no, and Crop::m_ManagementCategories.

Referenced by DK_EnergyCrop_Perm().

---

The documentation for this class was generated from the following files:

• /builds/ALMaSS/ALMaSS_MIDox/Landscape/cropprogs/DK_EnergyCrop_Perm.h
• /builds/ALMaSS/ALMaSS_MIDox/Landscape/cropprogs/DK_EnergyCrop_Perm.cpp