Aphid_Population_Manager.h

Go to the documentation of this file.

/*
*******************************************************************************************************
Copyright (c) 2021, Xiaodong Duan, Aarhus University
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 Aphid_Population_ManagerH
#define Aphid_Population_ManagerH
//---------------------------------------------------------------------------
 
//---------------------------------------------------------------------------
 
#include <Landscape/ls.h>
#include "../SubPopulation/SubPopulation_Population_Manager.h"
#include "../SubPopulation/SubPopulation.h"
#include <vector>
class Aphid;
class Aphid_Population_Manager;
class Landscape;
class Population_Manager;
class SubPopulation_Population_Manager;
class SubPopulation;
 
class struct_Aphid
{
 public:
  int x;
  int y;
  int w;
  int h;
  int species;
  Landscape* L;
  Aphid_Population_Manager * NPM;
  bool empty_flag;
  double* starting_suitability;
  double* starting_popu_density;
  int index_x;
  int index_y;
  TTypesOfLandscapeElement winter_landscape_host;
  TTypesOfLandscapeElement summer_landscape_host;
  bool farm_flag;
  double random_biomass_weight;
};
 
enum TTypesOfAphidDevelopmentSeason: unsigned
{
  toAphidSpringDev=0,
  toAphidLateSpringDev,
  toAphidSummerDev,
  toAphidFallDev,
  toAphidHibernate
};
 
class Aphid_Population_Manager : public SubPopulation_Population_Manager
{
   string m_aphid_specie_name;
   double m_aphid_egg_mortality;
   double m_aphid_egg_mortality_sd;
   double m_aphid_nymph_mortality_para_a;
   double m_aphid_nymph_mortality_para_b;
   double m_aphid_nymph_mortality_para_c;
   double m_aphid_adult_mortality_para_a;
   double m_aphid_adult_mortality_para_b;
   double m_aphid_adult_mortality_para_c;
   double m_aphid_adult_longevity_para_a;
   double m_aphid_adult_longevity_para_b;
   double m_aphid_egg_starting_degree_days;
   double m_aphid_egg_hatch_chance_para_a;
   double m_aphid_egg_hatch_chance_para_b;
   double m_aphid_egg_hatch_chance_threshold;
   double m_aphid_egg_hatch_delay_degree_days;
   double m_aphid_egg_hatch_chance_para_a_after_threshold;
   double m_aphid_egg_hatch_chance_para_b_after_threshold;
   double m_aphid_egg_hatch_chance_turing_degree_days;
   double m_aphid_egg_hatch_chance_max_degree_days;
   double m_aphid_nymph_with_wings_development_day_degrees;
   double m_aphid_nymph_with_wings_development_day_degrees_sd;
   double m_aphid_nymph_without_wings_development_day_degrees;
   double m_aphid_nymph_without_wings_development_day_degrees_sd;
   double m_aphid_egg_production_para_a;
   double m_aphid_egg_production_para_b;
   double m_aphid_egg_production_para_c;
   double m_aphid_egg_production_para_sd;
   double m_aphid_nymph_production_para_a1;
   double m_aphid_nymph_production_para_b1;
   double m_aphid_nymph_production_para_a2;
   double m_aphid_nymph_production_para_b2;
   double m_aphid_nymph_production_changing_add;
   double m_aphid_nymph_production_temp_weight_para_a;
   double m_aphid_nymph_production_temp_weight_para_b;
   double m_aphid_nymph_production_temp_weight_para_c;
   double m_aphid_alate_proportion_para_a;
   double m_aphid_alate_proportion_para_b;
   double m_aphid_alate_proportion_para_c;
   double m_aphid_male_reproduction_day_length;
   double m_aphid_female_reproduction_day_length;
   double m_aphid_winter_host_day_length;
   double m_aphid_summer_host_day_length;
   vector <double> m_aphid_adult_mortality_inflection_point_dd;
   vector <double> m_aphid_adult_mortality_vec;
   double m_aphid_sexual_reproduction_average_add;
   double m_aphid_sexual_reproduction_add_sd;
 
   TTypesOfAphidDevelopmentSeason m_current_developtype;
   int (Aphid_Population_Manager::*m_next_life_func_pointer)(int current_stage);
   void (Aphid_Population_Manager::*m_update_mortality_array_pointer)();
   int (Aphid_Population_Manager::*m_offspring_func_pointer)(int current_stage, double *offspring_num, double a_age, double a_density, double a_growth_stage, double* a_propotion, bool winter_host_flag);
   bool (Aphid_Population_Manager::*m_is_ready_next_liff_func_pointer)(int a_life_stage);
   int m_max_adult_dd;
   vector <double> m_nymph_per_day;
   vector <double> m_egg_per_day; 
   double m_nymph_produced_num_weight;
   bool m_egg_reproduction_flag;
   double m_accu_ddeg_egg_reproduction;
   bool m_first_egg_hatch_flag;
   double m_max_day_degree_adult;
   #ifdef __APHID_CALIBRATION
      FILE* m_aphid_calib_crop_output;
      FILE* m_aphid_calib_tole_output;
      FILE* m_aphid_calib_land_output;
   #endif
 
public:
   Aphid_Population_Manager(Landscape* L, string a_aphid_specie_name = "English grain aphid");
   virtual ~Aphid_Population_Manager();
   Aphid* CreateObjects(TAnimal *pvo, struct_Aphid* data, int number);
   virtual int calNextStage(int current_stage, double density = 1);
   virtual void updateDevelopmentSeason();
   virtual unsigned supplyDevelopmentSeason () {return m_current_developtype;}
   virtual void initialisePopulation();
   int calNextStageShared(int current_stage);
   void initialiseSimWithEggs(void);
   virtual void updateMortalityArray(void);
   void updateMortalityArrayShared(void);
   virtual int calOffspringStage(int current_stage, double *offspring_num=NULL, double a_age=1, double a_density=-1, double a_growth_stage=0, double* a_propotion=NULL, bool winter_host_flag=false);
   int calOffspringStageEGandPEA(int current_stage, double *offspring_num=NULL, double a_age=1, double a_density=-1, double a_growth_stage=0, double* a_propotion=NULL, bool winter_host_flag=false);
   int calOffspringStageBBandPP(int current_stage, double *offspring_num=NULL, double a_age=1, double a_density=-1, double a_growth_stage=0, double* a_propotion=NULL, bool winter_host_flag=false);
   virtual bool isEnoughNextLifeStage(int a_life_stage);
   bool isEnoughNextLifeStageShared(int a_life_stage);
   double supplyNymphNumWeight(void) {return m_nymph_produced_num_weight;}
   virtual void doSpeciesLastThing();
   virtual void writeCalibrationFiles(void);
   virtual double calHatchChance();
   #ifdef __APHID_CALIBRATION
      bool openAphidCalibFiles(void);
      void closeAphidCalibFiles(void);
      void writeAphidCalibData(void);
      void writeRowAphidCalibData(FILE* target_file, CfgArray_Int poly_id_list);
   #endif
 
   double SupplyAphidDensity(int a_x, int a_y);
   void KillAphid(int a_x, int a_y, double a_num);
   
};
 
#endif