ALMaSS/GooseManagement/elements_8h_source.html

//

// elements.h

//

/*

*******************************************************************************************************

Copyright (c) 2011, Christopher John Topping, University of Aarhus

All rights reserved.


Redistribution and use in source and binary forms, with or without modification, are permitted provided

that the following conditions are met:


Redistributions of source code must retain the above copyright notice, this list of conditions and the

following disclaimer.

Redistributions in binary form must reproduce the above copyright notice, this list of conditions and

the following disclaimer in the documentation and/or other materials provided with the distribution.


THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR

IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND

FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS

BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,

BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR

BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT

LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS

SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

********************************************************************************************************

*/


#ifndef ELEMENTS_H

#define ELEMENTS_H


#define FMDEBUG

#include <fstream>

extern class Landscape * g_landscape_p;


class APoint;


typedef unsigned int LE_Signal;


static CfgBool l_el_read_bug_percentage_file("ELEM_READ_BUG_PERCENTAGE_FILE", CFG_CUSTOM, false);

static CfgStr l_el_bug_percentage_file("ELEM_BUG_PERCENTAGE_FILE", CFG_CUSTOM,"bugpercents.txt");


#define LE_SIG_NO_FUNGICIDE   0x0001

#define LE_SIG_NO_INSECTICIDE 0x0002

#define LE_SIG_NO_SYNG_INSECT 0x0004

#define LE_SIG_NO_GROWTH_REG  0x0008

#define LE_SIG_NO_MOLLUSC     0x0010

#define LE_SIG_NO_HERBICIDE   0x0020


// Number of user defined storage locations for farm management

// states. For ridiculously active farm management plans this

// might have to be increased. Ie. bump this number if there is

// a chance of more than 10 farm events for an given polygon on

// any given day.

#define EL_MAX_USERSPACE 10 // Has to be a #define.


// Decay time for treatment generated tramlines (in days).

#define EL_TRAMLINE_DECAYTIME (g_el_tramline_decaytime_days.value())

extern CfgInt g_el_tramline_decaytime_days;

// Delay time before regrowth of weeds after herbicide treatment.

#define EL_HERBICIDE_DELAYTIME (g_el_herbicide_delaytime_days.value())

#define EL_STRIGLING_DELAYTIME (g_el_strigling_delaytime_days.value())

extern CfgInt g_el_herbicide_delaytime_days;

extern CfgInt g_el_strigling_delaytime_days;


// Will be declared in other .h files, but these are needed here to break circular class definitions.

class Farm;

extern class LE_TypeClass *g_letype;


class LE_TypeClass

{

 public:

  TTypesOfLandscapeElement TranslateEleTypes(int EleReference);

  TTypesOfVegetation       TranslateVegTypes(int VegReference);


  int BackTranslateEleTypes(TTypesOfLandscapeElement EleReference);

  int BackTranslateVegTypes(TTypesOfVegetation VegReference);

  int VegTypeToCurveNum(    TTypesOfVegetation VegReference);

};


class LE {

public:

        LE(void);

  virtual ~LE( void );

  int m_tried_to_do;

  int m_squares_in_map;

  bool m_skylarkscrapes;

  //double m_volegrazing;

  // /** \brief Variable used to reduce the grazing pressure */

  //double m_voleResregrowth;

  //double m_volegrazingDensity;

  //void AddVoleGrazing(int a_number) { m_volegrazing+=a_number; }

  //double GetVoleGrazing( void ) {return m_volegrazing; }

  //void CalcVoleGrazingDensity( void ) { m_volegrazing *= m_voleResregrowth; m_volegrazingDensity = (m_volegrazing/m_area); }

  //double GetVoleGrazingDensity( void ) {return m_volegrazingDensity; }

  long  m_user[ EL_MAX_USERSPACE ];

  bool  m_repeat_start;

  LE_Signal GetSignal( void ) { return m_signal_mask; }

  void      SetSignal( LE_Signal a_signal ) { m_signal_mask = a_signal; }

  int  GetSubType( void ) { return m_subtype; }

  void SetSubType( int a_subtype ) { m_subtype = a_subtype; }

  int  GetALMaSSEleType(void) { return m_almass_le_type; }

  void SetALMaSSEleType(int a_type) { m_almass_le_type = a_type; }

  virtual TTypesOfLandscapeElement GetElementType( void ) { return m_type; }

  void   SetElementType( int a_type ) { m_type = (TTypesOfLandscapeElement) a_type; }

  void   Trace( int a_value );

  void   ResetTrace( void );

  long   GetRunNum( void ) { return m_running; }

  void   BumpRunNum( void ) { m_running++; }

  long   GetMgtLoopDetectDate ( void ) { return m_management_loop_detect_date; }

  void   SetMgtLoopDetectDate ( long a_num ) { m_management_loop_detect_date = a_num; }

  long   GetMgtLoopDetectCount ( void ) { return m_management_loop_detect_count; }

  void   SetMgtLoopDetectCount ( long a_num ) { m_management_loop_detect_count = a_num; }

  void   Tick( void );

  virtual TTypesOfVegetation

         GetVegType( void ) { return tov_None; }

  int    GetVegStore( void ) { return m_vege_danger_store; }

  void   SetVegStore( int a_veg )

    { m_vege_danger_store = a_veg; }

  void ZeroVegAge() { m_vegage=0; }

  int GetVegAge() { return m_vegage; }

  virtual void  DoDevelopment( void );

  virtual void  RecalculateBugsNStuff( void ) { }

  virtual int GetVegPhase(void) { return vegphase_foobar; }

  virtual double GetVegCover(void) { return 0.0; }

  virtual double GetVegHeight( void ) { return 0.0; }

  virtual APoint GetCentroid();

  virtual int GetCentroidX() { return m_centroidx; }

  virtual int GetCentroidY() { return m_centroidy; }

  virtual void SetCentroid(int x, int y) { m_centroidx=x; m_centroidy=y; }

  virtual double GetDigestability( void ) { return 0.0; }

  virtual int GetVegDensity( void ) { return 0; }

  virtual bool GetSkScrapes( void ) { return false; }

  virtual bool IsMatureCereal() { return false; }

  virtual bool IsCereal() { return false; }

  virtual bool IsGrass() { return false; }

  virtual bool GetVegPatchy(void) { return false; }

  //virtual void SetVegPatchy( bool p ) { ;}

  virtual void SetVegPatchy( bool /* p */ ) { ;}

  virtual double GetDayDegrees( void ) { return m_ddegs; }

  virtual double GetLAGreen( void ) { return 0.0; }

  virtual double GetLATotal( void ) { return 0.0; }

  virtual double GetGreenBiomass( void ) { return 0.0; }

  virtual double GetDeadBiomass( void ) { return 0.0; }

  virtual double GetVegBiomass( void ) { return 0.0; }

  virtual double GetWeedBiomass( void ) { return 0.0; }

  virtual double GetInsectPop( void ) { return 0.0; }

  virtual void  SetInsectPop( double  /* i */ ) { ; }

  virtual void  Insecticide( double /* a_fraction */ ) { ; }

  virtual void  InsectMortality( double /* a_fraction */ ) { ; }

  virtual void  ReduceWeedBiomass( double /* a_fraction */ ) { ; }

  virtual void  ToggleCattleGrazing( void ) {;}

  virtual void  TogglePigGrazing( void ) {;}

  virtual void  SetVegBiomass( int /* a_biomass */ ) {;}

  virtual void  SetVegType( TTypesOfVegetation /* a_type */, TTypesOfVegetation /* a_weed_type */  ) {;}

  virtual void  SetVegHeight(double /* a_height */) { ; }

  virtual void  SetVegParameters(double /* a_height */, double /* a_LAtotal */, double /* a_LAgreen */, double /* a_WeedBiomass */) { ; }

  virtual void SetCropData(double, double, double, TTypesOfVegetation, double, int) { ; }

  virtual void SetCropDataAll(double,double,double, double,TTypesOfVegetation, double, double , int, double , bool, double) {;}

  virtual void  StoreLAItotal() {;}

  virtual void  SetGrowthPhase( int /* a_phase */ ) {;}

  virtual void  ForceGrowthInitialize( void ) {;}

  virtual void  ForceGrowthTest( void ) {;}

  virtual void  ForceGrowthDevelopment( void ) {;}

  virtual void  ZeroVeg( void ) {;}

  virtual void  ReduceVeg( double /* a_reduc */ ) {;}

  virtual void  ReduceVeg_Extended(double /* a_reduc */) { ; }

  virtual void  GrazeVeg_Extended(double /* a_reduc */) { ; }

  virtual double GetTrafficLoad(void) { return 0.0; }

  virtual TTypesOfVegetation GetPreviousCrop(int /* a_index */) { return tov_None;  }

  bool  GetPoison( void ) { return m_poison; }

  void  SetPoison( bool a_poison ) { m_poison = a_poison; }

  bool  GetHigh( void ) { return m_high; }

  int  GetCattleGrazing( void ) { return m_cattle_grazing; }

  bool  GetPigGrazing( void ) { return m_pig_grazing; }

  bool  HasTramlines( void ) { return ( m_tramlinesdecay>0 ); }

  bool  IsRecentlyMown( void ) { return ( m_mowndecay>0 ); }

  int   IsRecentlySprayed( void ) { return ( m_herbicidedelay ); }

  int   GetPoly( void ) { return m_poly; }

  int   GetMapIndex( void ) { return m_map_index; }

  void  SetMapIndex( int a_map_index ) { m_map_index = a_map_index; }

  int   GetOwnerFile( void ) { return m_owner_file; }

  int   GetOwnerIndex( void ) { return m_owner_index; }

  int   GetLastTreatment( int *a_index );

  Farm* GetOwner( void ) { return m_owner; }

  double GetArea( void ) { return m_area; }

  double GetGooseGrazingForage(GooseSpecies a_goose) { return m_goosegrazingforage[a_goose]; }

  double GetBirdSeed( void ) {

          return m_birdseedforage;

  }

  void SetBirdSeed( double a_forage ) {

          m_birdseedforage = a_forage;

  }

  double GetBirdMaize(void) {

          return m_birdmaizeforage;

  }

 /*

 double GetIsCereal(void) {

          if ( m_goosegrazingforage[gs_foobar] == 1) return true;

          return false;

  }

*/

  void SetBirdMaize( double a_forage ) {

          m_birdmaizeforage = a_forage;

  }

  int GetCountryDesignation( void ) { return m_countrydesignation; }

  void SetCountryDesignation( int a_designation ) { m_countrydesignation = a_designation; }

  long  GetOldDays( void ) { return m_olddays; }

  void  SetOwner( Farm *a_owner, int a_owner_num, int a_owner_index )

  {

    m_owner = a_owner;

    m_owner_file = a_owner_num;

    m_owner_index = a_owner_index;

  }

  void  SetHigh( bool a_high ) { m_high = a_high; }

  void  SetPoly( int a_poly ) { m_poly = a_poly; }

  void  SetMapValid( bool a_valid ) { m_is_in_map = a_valid; }

  bool  GetMapValid( void ) { return m_is_in_map; }

  void  SetArea( double a_area ) { m_area = a_area; }

  virtual void GrazeVegetation( double /* m_grams */, bool )

  {

          ;

  }

  virtual void GrazeVegetationTotal( double /* m_grams */ )

  {

          ;

  }

  int   GetSoilType( ) {

          return m_soiltype;

  }

  int   GetSoilTypeR( ) {

          switch (m_soiltype)

          {

                  case 0: return 3;

                  case 1: return 1;

                  case 2: return 3;

                  case 3: return 3;

                  case 4: return 0;

                  case 5: return 1;

                  case 6: return 0;

                  case 7: return 0;

                  case 8: return 1;

                  case 9: return 0;

                  case 10: return 0;

                  case 11: return 0;

                  case 12: return 0;

                  case 13: return 3;

                  default:

                          return 3;

          }

  }

  void  SetSoilType( int a_st ) {

          m_soiltype = a_st;

  }

  void SetLastSownVeg(TTypesOfVegetation a_tov) {

          m_lastsownvegtype = a_tov;

  }

  TTypesOfVegetation GetLastSownVeg() {

          return m_lastsownvegtype;

  }

  void  SetLastTreatment(int a_treatment);

  void SetSprayedToday(bool a_didit) { m_sprayedtoday = a_didit; }

  bool GetSprayedToday() { return m_sprayedtoday; }

  void  SetCopyTreatment(int a_treatment);

  void  SetOldDays( long a_days ) { m_olddays = a_days; }

  void  SetRotIndex( int a_index ) { m_rot_index = a_index; }


  int   GetRotIndex( void ) { return m_rot_index; }

  void  SetTramlinesDecay( int a_decaytime_days )

    { m_tramlinesdecay = a_decaytime_days; }

  void  SetHerbicideDelay( int a_decaytime_days )

    { m_herbicidedelay = a_decaytime_days; }

  void  SetMownDecay( int a_decaytime_days )

    { m_mowndecay = a_decaytime_days; }

  void  AddArea( double a_area_diff ) { m_area += a_area_diff; }

  LE*   GetBorder( void ) { return m_border; }

//  LE* GetUnsprayedMargin( void ) { return m_unsprayedmargin; }

  int   GetUnsprayedMarginPolyRef( void ) { return m_unsprayedmarginpolyref; }

  void  SetBorder( LE* a_border ) { m_border = a_border; }

//  void  SetUnsprayedMargin( LE* a_unsprayedmargin ) { m_unsprayedmargin = a_unsprayedmargin; }

  void  SetUnsprayedMarginPolyRef( int a_unsprayedmargin ) { m_unsprayedmarginpolyref = a_unsprayedmargin; }

  void  SetPesticideCell( int a_cell ) { m_PesticideGridCell = a_cell; }

  int   GetPesticideCell( ) { return m_PesticideGridCell; }

  void  SetValidXY( int a_valid_x, int a_valid_y )

    {

      m_valid_x = a_valid_x;

      m_valid_y = a_valid_y;

    }

  int   GetValidX( void ) { return m_valid_x; }

  int   GetValidY( void ) { return m_valid_y; }


  int   GetMaxX( void ) { return m_maxx; }

  int   GetMaxY( void ) { return m_maxy; }

  int   GetMinX( void ) { return m_minx; }

  int   GetMinY( void ) { return m_miny; }

  void  SetMaxX( int x ) { m_maxx=x; }

  void  SetMaxY( int y ) { m_maxy=y; }

  void  SetMinX( int x ) { m_minx=x; }

  void  SetMinY( int y ) { m_miny=y; }

  int   GetMDates(int a, int b) {return MDates[a][b];}

  void  SetMDates(int a, int b, int c) {MDates[a][b]=c;}

  int   GetMConstants(int a) {return MConsts[a];}

  void  SetMConstants(int a, int c) {MConsts[a]=c;}

  double GetVegDDegs() {

          return m_vegddegs; //

  }

  int GetOpenness( void ) { return m_openness; }

  void SetGooseNos(int a_number, int a_day) { m_gooseNos[a_day] = a_number; }

  void SetGooseSpNos(int a_number, int a_day, GooseSpecies a_goose) {

          m_gooseSpNos[a_day][a_goose] = a_number;

  }

  void SetGooseNosTimed(int a_number, int a_day) { m_gooseNosTimed[a_day] = a_number; }

  void SetGooseSpNosTimed(int a_number, int a_day, GooseSpecies a_goose) {

          m_gooseSpNosTimed[a_day][a_goose] = a_number;

  }

  void SetGooseRoostDist(int a_dist, GooseSpecies a_goose) { m_dist_to_closest_roost[a_goose] = a_dist; }

  int GetGooseNosToday();

  int GetGooseNosTodayTimed();

  int GetGooseNos( );

  int GetQuarryNos();

  int GetGooseSpNosToday(GooseSpecies a_goose);

  int GetGooseSpNosTodayTimed(GooseSpecies a_goose);

  int GetGooseRoostDist(GooseSpecies a_goose);

  void SetOpenness( int a_openness ) { m_openness=a_openness; }

  virtual void ResetDigestability() { ; }

  void DoCopy(const LE* a_Le);

  virtual void SetMaleNewtPresent(bool) { ; };

  virtual bool IsMaleNewtPresent() { return false; };

  void SetStubble(bool a_flag) { m_instubble = a_flag; }

  bool GetStubble() { return m_instubble; }

  /* *********************** Bee Helper Methods *********************** */

  bool GetOsmiaNest() {

          if ((m_currentOsmiaNests < m_maxOsmiaNests) && (g_rand_uni() < m_OsmiaNestProb)) return true;

          return false;

  }


  void ReleaseOsmiaNest() {

          m_currentOsmiaNests--;

  }


  void UpdateOsmiaNesting()

  {

          if (m_maxOsmiaNests <= 0) m_OsmiaNestProb = 0.0;

          else {

                  m_OsmiaNestProb = 1.0 - (m_currentOsmiaNests / m_maxOsmiaNests);

          }

  }


  void IncOsmiaNesting()

  {

          m_currentOsmiaNests++;

  }


  void SetMaxOsmiaNests(double a_noNests)

  {

          m_maxOsmiaNests = int(a_noNests * m_area);

          if (m_maxOsmiaNests < 1) m_maxOsmiaNests = 0;

  }

  /* ********************* End Bee Helper Methods ********************* */


  void SetPollenNectarData(int a_almasstype);

  virtual PollenNectarQuality GetPollen() { return m_pollenquality; }

  virtual PollenNectarQuality GetNectar() { return m_nectarquality; }

  virtual double GetTotalPollen() { return m_totalPollen; }

  virtual double GetTotalNectar() { return m_totalNectar; }

  void SetPollenNectarCurves(PollenNectarDevelopmentCurve* a_pollencurnve, PollenNectarDevelopmentCurve* a_nectarcurve) {

          m_pollencurve = a_pollencurnve;

          m_nectarcurve = a_nectarcurve;

  }

protected:

#ifdef FMDEBUG

  int   m_ptrace[256]; // Change bitmap in 'elements.cpp' if this

  int   m_pdates[256]; // needs to be increased!

  int   m_pindex;      // Points to the next *empty* location, so last

                       // program step is in the location below this...

#endif

  long  m_running;

  LE    *m_border;

  int    m_unsprayedmarginpolyref;

  // May be -1 if no unsprayed margin or this element is a margin!

  // Valid map coordinate inside this landscape element

  // before(!) adding the border.

  int   m_valid_x;

  int   m_valid_y;

  int m_centroidx;

  int m_centroidy;

  int   m_vege_danger_store;

  int   m_PesticideGridCell;

  // Whether this polygon is actually mentioned in the main map.

  bool  m_is_in_map;

  // Quick hack to discern different hedgebank types.

  int   m_subtype;

  int m_maxx;

  int m_minx;

  int m_maxy;

  int m_miny;

  long  m_olddays;

  long  m_management_loop_detect_date;

  long  m_management_loop_detect_count;

  Farm *m_owner;

  int m_countrydesignation;

  int   m_owner_file;

  int   m_owner_index;

  int   m_rot_index;

  bool  m_sprayedtoday;

  int   m_poly;

  // Copy of the index stored in the main map at run time.

  // Needed by the pesticide machinery.

  int   m_map_index;

  int   m_almass_le_type;

  int   m_farmfunc_tried_to_do;

  int   m_tramlinesdecay;

  int   m_mowndecay;

  int   m_herbicidedelay;

  TTypesOfLandscapeElement m_type;

  double m_area;

  int m_gooseNos[366];

  int m_gooseSpNos[366][gs_foobar];

  int m_gooseNosTimed[366];

  int m_gooseSpNosTimed[366][gs_foobar];

  bool m_instubble;

  double m_birdseedforage;

  double m_birdmaizeforage;

  double m_goosegrazingforage[ gs_foobar ];

  TTypesOfVegetation m_lastsownvegtype;

  int m_openness;

  double m_dist_to_closest_roost[gs_foobar];

  int m_vegage;

  int m_soiltype;

  int m_days_since_insecticide_spray;

  bool  m_high;

  bool  m_poison;

  int  m_cattle_grazing;

  int m_default_grazing_level;

  bool  m_pig_grazing;

  double m_vegddegs; // Total of day degrees for the crop/veg on this element.

  double m_yddegs; // Yesterdays total of day degrees for this element.

  double m_ddegs;  // Total of day degrees for this element when today

                  // is taken into consideration.

  vector <int> m_lasttreat;

  unsigned int m_lastindex;

  // Only one copy of these is stored for use by all landscape elements.

  static double m_monthly_traffic[];

  static double m_largeroad_load[];

  static double m_smallroad_load[];

  int MDates[2][25];

  int MConsts[10];

  LE_Signal m_signal_mask;

  double m_OsmiaNestProb;

  int m_maxOsmiaNests;

  int m_currentOsmiaNests;

  PollenNectarQuality m_pollenquality;

  PollenNectarQuality m_nectarquality;

  double m_totalPollen;

  double m_totalNectar;

  double m_PollenM2;

  double m_NectarM2;

  PollenNectarDevelopmentCurve* m_pollencurve;

  PollenNectarDevelopmentCurve* m_nectarcurve;

};


//---------------------------------------------------------------------------


struct CropActualValues {

        bool taken;

        TTypesOfVegetation tov_type;

        double biomass_at_harvest;

        bool harvested;

        double area;

        int no_herb_app; //number of herbicide applications

        int missed_herb_app;

        int no_fi_app;

        int missed_fi_app;

};


class VegElement : public LE {

public:

  VegElement( void );

  virtual ~VegElement( void ) { ; };

  virtual void  DoDevelopment( void );

  virtual double GetDigestability( void ) { return m_digestability; }

  virtual double GetLAGreen( void ) { return m_LAgreen; }

  virtual double GetLATotal( void ) { return m_LAtotal; }

  virtual double GetVegCover( void ) { return m_veg_cover; }

  virtual double GetVegHeight( void ) { return m_veg_height; }

  virtual double GetVegBiomass( void ) { return m_veg_biomass; }

  virtual double GetDeadBiomass( void ) { return m_dead_biomass; }

  virtual double GetGreenBiomass( void ) { return m_green_biomass; }

  virtual double GetWeedBiomass( void ) { return m_weed_biomass; }

  virtual int GetVegDensity( void ) { return m_veg_density; }

  virtual bool GetSkScrapes( void ) { return m_skylarkscrapes; }

  virtual bool GetVegPatchy( void ) { return m_veg_patchy; }

  virtual void SetVegPatchy( bool p ) { m_veg_patchy=p; }

  virtual double GetInsectPop( void ) { return m_insect_pop; }

  virtual void SetInsectPop( double insects ) { m_insect_pop=insects; }

  virtual TTypesOfVegetation GetVegType( void ) { return m_vege_type; }

  virtual void  Insecticide( double a_fraction )

  {

          m_insect_pop *= a_fraction;

      m_days_since_insecticide_spray=28;

  }

  virtual bool IsCereal();

  virtual bool IsMatureCereal();

  virtual bool IsMaize();

  virtual bool IsGrass();

  virtual bool IsGooseGrass();

  virtual void  InsectMortality(double a_fraction);

  virtual void  ReduceWeedBiomass( double a_fraction )

  {

                m_weed_biomass *= a_fraction;

  }

  virtual void  ToggleCattleGrazing(void) { if (m_cattle_grazing == 0) m_cattle_grazing = m_default_grazing_level; else m_cattle_grazing = 0;; }

  virtual void  TogglePigGrazing( void )        { m_pig_grazing = !m_pig_grazing; }

  virtual void GrazeVegetation( double a_grams, bool a_force );

  virtual void GrazeVegetationTotal( double a_grams);

  // Outdated! Use ReduceVeg() below instead. Total biomass is

  // (re)calculated from the total leaf area, it cannot just be

  // set or changed to an arbitrary value.

  //virtual void  SetVegBiomass( int a_veg_biomass )

  //  { m_veg_biomass = a_veg_biomass; }

  virtual void SetVegType(TTypesOfVegetation a_vege_type, TTypesOfVegetation a_weed_type);

  virtual void SetCropData(double,double,double,TTypesOfVegetation, double, int);

  virtual void SetCropDataAll(double,double,double, double,TTypesOfVegetation, double, double , int, double , bool, double);

  virtual void SetVegHeight(double a_veg_height) { m_veg_height = a_veg_height; }

  virtual void  SetVegParameters(double a_veg_height, double a_LAtotal, double a_LAgreen, double a_WeedBiomass)  {

      m_veg_height = a_veg_height;

      m_LAtotal    = a_LAtotal;

      m_LAgreen    = a_LAgreen;

      m_weed_biomass = a_WeedBiomass;

      //RecalculateBugsNStuff();

  }

  virtual void  SetGrowthPhase( int a_phase );

  virtual void  StoreLAItotal() { m_oldLAtotal = m_LAtotal; }

  virtual void  ForceGrowthInitialize( void );

  virtual void  ForceGrowthTest( void );

  virtual void  ForceGrowthDevelopment( void );

  virtual void  ZeroVeg( void );

  virtual void  ReduceVeg(double a_reduc);

  virtual void  ReduceVeg_Extended(double a_reduc);


  //added 09.04.13 - moved from Field

        //seems all get function not necessary

        void Set_taken(bool a_taken, int index){m_CropDataStorage[index].taken = a_taken;};

        bool Get_taken(int index){return m_CropDataStorage[index].taken;};

        void Set_tov_type(TTypesOfVegetation a_tov_type, int index){m_CropDataStorage[index].tov_type = a_tov_type;};

        //TTypesOfVegetation Get_tov_type (int index){return m_CropDataStorage[index].tov_type;};

        void Set_area_in_crop_data(double a_area){m_CropDataStorage[0].area = a_area;};

        void Set_biomass_at_harvest(double a_biomass, int index){m_CropDataStorage[index].biomass_at_harvest = a_biomass;};

        //double Get_biomass_at_harvest(int index){return m_CropDataStorage[index].biomass_at_harvest;};

        void Set_harvested(){m_CropDataStorage[0].harvested = true;};

        bool Get_harvested(int index){return m_CropDataStorage[index].harvested;};

        //void Add_no_herb_app(int index){m_CropDataStorage[index].no_herb_app++;};

        void Add_no_herb_app(){m_CropDataStorage[0].no_herb_app++;};

        //int Get_no_herb_app (int index){return m_CropDataStorage[index].no_herb_app;};

        //void Add_no_fi_app(int index){m_CropDataStorage[index].no_fi_app++;};

        void Add_no_fi_app(){m_CropDataStorage[0].no_fi_app++;};

        //int Get_no_fi_app (int index){return m_CropDataStorage[index].no_fi_app;};


        //void Add_missed_herb_app(int index){m_CropDataStorage[index].missed_herb_app++;}; //index not necessary - always save at position 0.

        //int Get_missed_herb_app (int index){return m_CropDataStorage[index].missed_herb_app;};

        //void Add_missed_fi_app(int index){m_CropDataStorage[index].missed_fi_app++;};

        //int Get_missed_fi_app (int index){return m_CropDataStorage[index].missed_fi_app;};

        void Add_missed_herb_app(){m_CropDataStorage[0].missed_herb_app++;};

        void Add_missed_fi_app(){m_CropDataStorage[0].missed_fi_app++;};


        //void Set_crop_index(int a_index){m_crop_index = a_index;}; //probably not needed

        //int Get_crop_index(void){return m_crop_index;};

        void Set_CropDataStorage(int index, CropActualValues a_struct){m_CropDataStorage[index] = a_struct;};

        CropActualValues Get_CropDataStorage (int index){return m_CropDataStorage[index];};


        void Clean_CropDataStorage(int index);

        //09.04

        virtual void ResetDigestability()

        {

                for (int i = 0; i<32; i++)  m_oldnewgrowth[i] = 0.0;

        }


        virtual int GetVegPhase() { return m_veg_phase; }


protected:

  TTypesOfVegetation m_vege_type;

  int   m_curve_num;

  int   m_weed_curve_num;

  int   m_veg_phase;

  double   m_growth_scaler;

  int   m_nutrient_status;

  double m_LAtotal;

  double m_oldLAtotal;

  double m_LAgreen;

  double m_insect_pop;

  double m_veg_biomass;

  double m_total_biomass;

  double m_total_biomass_old;

  double m_green_biomass;

  double m_dead_biomass;

  double m_weed_biomass;

  double m_veg_height;

  int m_veg_density;

  bool m_veg_patchy;

  double m_veg_cover;

  bool  m_forced_phase_shift;

  double m_digestability;

  double m_oldnewgrowth[32];

  int m_newoldgrowthindex;

  double m_newgrowthsum;

  double m_force_Weed;

  double m_force_LAgreen;

  double m_force_LAtotal;

  double m_force_veg_height;

  bool  m_force_growth;

  ifstream m_ifile;

  void  ReadBugPercentageFile( void );

  virtual void  RecalculateBugsNStuff( void );

  void CalcGooseForageResources();

  void  RandomVegStartValues( double *a_LAtotal,double *a_LAgreen,double *a_veg_height,double *a_weed_biomass );

  void  ForceGrowthSpringTest( void );

  void PollenNectarPhenologyCalculation();

  void ResetGeese(void);


  //added 09.04 .13 - moved from Field!

  CropActualValues m_CropDataStorage [2];

  int m_crop_index;

};


class Field : public VegElement

{

public:

  Field( void );

  virtual void DoDevelopment( void );

  TTypesOfVegetation GetPreviousCrop(int a_index);


protected:


};


class PermPasture : public VegElement

{

public:

  PermPasture( void );

};


class PermPastureLowYield : public VegElement

{

 public:

  PermPastureLowYield( void );

};


class PermPastureTussocky : public VegElement

{

 public:

  PermPastureTussocky( void );

};


class PermanentSetaside : public VegElement

{

public:

  PermanentSetaside( void );

};


class Hedges : public VegElement

{

public:

  Hedges( void );

  virtual void DoDevelopment( void ) {

    VegElement::DoDevelopment();

    m_insect_pop=0;

  }

};


class HedgeBank : public VegElement

{

public:

  HedgeBank( void );

  virtual void DoDevelopment( void );

};


class BeetleBank : public VegElement

{

public:

  BeetleBank( void );

  virtual void DoDevelopment( void );

};


class RoadsideVerge : public VegElement

{

  // Attributes

  void Cutting(int a_today);

  long m_DateCut;


 public:

  RoadsideVerge( void );

  // Is needed due to special development for this element

  // type, cutting of the roadside verges.

  virtual void DoDevelopment( void );

};


class WaterBufferZone : public VegElement

{

        void ResetingVeg(int a_today);

public:

        WaterBufferZone(void);

        // Is needed due to special development for this element

        // We want to reset the vegetation every year

        virtual void DoDevelopment(void);

};


class RoadsideSlope : public RoadsideVerge

{

 public:

  RoadsideSlope( void );

  virtual void DoDevelopment( void );

};


class GreenElement : public VegElement

{

public:

        GreenElement( void );

        virtual ~GreenElement( void ) { ; }

  virtual void DoDevelopment( void );

};


class Scrub : public GreenElement

{

public:

  Scrub( void );

  virtual void DoDevelopment( void ) {

    VegElement::DoDevelopment();

    m_insect_pop=0;

  }

};


class Marsh : public GreenElement

{

public:

  Marsh( void );

};


class Saltmarsh : public Marsh

{

public:

  Saltmarsh( void );

};


class Heath : public GreenElement

{

public:

  Heath( void );

};


class RiversidePlants : public GreenElement

{

public:

  RiversidePlants( void );

};


class NaturalGrassDry : public GreenElement

{

public:

  NaturalGrassDry( void );

  virtual void DoDevelopment();

};


class NaturalGrassWet : public GreenElement

{

public:

  NaturalGrassWet( void );

  virtual void DoDevelopment();

};


class Wasteland : public GreenElement

{

public:

        Wasteland( void );

        virtual void DoDevelopment();

};


class FieldBoundary : public GreenElement

{

public:

  FieldBoundary( void );

};


class UnsprayedFieldMargin : public GreenElement

{

public:

  UnsprayedFieldMargin( void );

  virtual void DoDevelopment( void );

};


class RiversideTrees : public GreenElement

{

public:

  RiversideTrees( void );

};


class Vildtager : public GreenElement

{

public:

        Vildtager( void );

};


class Railway : public GreenElement

{

public:

  Railway( void );

};


//---------------------------------------------------------------------------


class ForestElement : public VegElement

{

public:

  ForestElement( void );

  virtual void DoDevelopment( void ) {

    VegElement::DoDevelopment();

    m_insect_pop=0;

  }

};


class DeciduousForest : public ForestElement

{

public:

  DeciduousForest( void );

};


class Copse : public ForestElement

{

public:

  Copse( void );

};


class IndividualTree : public ForestElement

{

public:

        IndividualTree(void);

};


class ConiferousForest : public ForestElement

{

public:

  ConiferousForest( void );

};


class WoodlandMargin : public ForestElement

{

public:

        WoodlandMargin(void);

};


class MixedForest : public ForestElement

{

public:

  MixedForest( void );

};


class YoungForest : public ForestElement

{

public:

  YoungForest( void );

};


class WoodyEnergyCrop : public ForestElement

{

public:

        WoodyEnergyCrop( void );

};


class Orchard : public GreenElement

{

  void Cutting(int a_today);

  long m_DateCut;


public:

  Orchard( void );

  // Is needed due to special development for this element

  // type, cutting of the roadside verges.

  virtual void  DoDevelopment( void );

};


class OrchardBand : public GreenElement

{

  long m_LastSprayed;


public:

  OrchardBand( void );

  // Is needed due to special development for this element

  // type, cutting of the roadside verges.

  virtual void  DoDevelopment( void );

};


class MownGrass : public GreenElement

{

  void Cutting(int a_today);

  long m_DateCut;


public:

  MownGrass( void );

  virtual void  DoDevelopment( void );

};


class NonVegElement : public LE

{

public:

  NonVegElement( void );

};


class ChameleonLE : public NonVegElement

{

public:

        ChameleonLE(void);

};


class Garden : public NonVegElement

{

public:

        Garden(void);

};


class Building : public NonVegElement

{

public:

  Building( void );

};


class Pylon : public NonVegElement

{

public:

        Pylon(void);

};


class WindTurbine : public NonVegElement

{

public:

        WindTurbine(void);

};


class PlantNursery : public GreenElement

{

public:

        PlantNursery(void);

};


class Urban : public NonVegElement

{

public:

  Urban( void );

};


class AmenityGrass : public NonVegElement

{

public:

  AmenityGrass( void );

};


class Parkland  : public NonVegElement

{

public:

  Parkland( void );

};


class UrbanNoVeg : public NonVegElement

{

public:

  UrbanNoVeg( void );

};


class UrbanVeg : public NonVegElement

{

public:

        UrbanVeg(void);

};


class UrbanPark : public NonVegElement

{

public:

  UrbanPark( void );

};


class Suburban : public NonVegElement

{

public:

  Suburban( void );

};


class RuralResidential : public NonVegElement

{

public:

  RuralResidential( void );

};


class BuiltUpWithParkland : public NonVegElement

{

public:

  BuiltUpWithParkland( void );

};


class SandDune : public NonVegElement

{

public:

  SandDune( void );

};


class BareRock : public NonVegElement

{

public:

  BareRock( void );

};


class StoneWall : public NonVegElement

{

public:

  StoneWall( void );

};


class PitDisused : public NonVegElement

{

public:

  PitDisused( void );

};


class Saltwater : public NonVegElement

{

public:

  Saltwater( void );

};


class Freshwater : public NonVegElement {

public:

        Freshwater(void);

};


class FishFarm : public NonVegElement {

public:

        FishFarm(void);

};


class Fence : public NonVegElement

{

public:

        Fence(void);

};


class Pond : public Freshwater {

protected:

        double m_LarvalFood;

        double m_LarvalFoodScaler;

        double m_pondpesticide;

        bool m_MaleNewtPresent;

        double m_pondquality;

public:

        Pond( void );

        virtual void  DoDevelopment(void);

        void CalcPondPesticide();

        void CalcLarvalFood();

        bool SubtractLarvalFood(double a_food);

        double SupplyPondPesticide() { return m_pondpesticide; }

        virtual void SetMaleNewtPresent(bool a_ispresent) { m_MaleNewtPresent = a_ispresent; };

        bool IsMaleNewtPresent() { return m_MaleNewtPresent; };

};


class RiverBed : public NonVegElement

{

public:

        RiverBed(void);

};


class River : public NonVegElement

{

public:

        River(void);

};


class Canal : public NonVegElement

{

public:

        Canal(void);

};


class RefuseSite : public NonVegElement

{

public:

        RefuseSite(void);

};


class DrainageDitch : public NonVegElement

{

public:

        DrainageDitch(void);

};


class Coast : public NonVegElement

{

public:

  Coast( void );

};


class ActivePit : public NonVegElement

{

public:

  ActivePit( void );

};


class LargeRoad : public NonVegElement

{

public:

  virtual double GetTrafficLoad( void );

  LargeRoad( void );

};


class SmallRoad : public NonVegElement

{

public:

  virtual double GetTrafficLoad( void );

  SmallRoad( void );

};


class Track : public NonVegElement

{

public:

  Track( void );

};


class MetalledPath : public NonVegElement

{

public:

  MetalledPath( void );

};


class Carpark : public NonVegElement

{

public:

  Carpark( void );

};


class Churchyard : public NonVegElement

{

public:

  Churchyard( void );

};


class Stream : public NonVegElement

{

public:

  Stream( void );

};


class HeritageSite : public NonVegElement

{

public:

  HeritageSite( void );

};


#endif // ELEMENTS_H