ALMaSS_CmdLine.cpp File Reference

#include "assert.h"
#include <cstdlib>
#include <cstring>
#include <stdio.h>
#include <stdlib.h>
#include <errno.h>
#include <iostream>
#include <fstream>
#include <time.h>
#include <cmath>
#include <vector>
#include <list>
#include <string>
#include "../BatchALMaSS/ALMaSS_Setup.h"
#include "../ALMaSSDefines.h"
#include "../Landscape/ls.h"
#include "../BatchALMaSS/MovementMap.h"
#include "../BatchALMaSS/BinaryMapBase.h"
#include "../BatchALMaSS/PopulationManager.h"
#include "../Beetles/Beetle_BaseClasses.h"
#include "../Beetles/PoecilusCupreus_All.h"
#include "../Beetles/Ladybird_All.h"
#include "../Beetles/Bembidion_All.h"
#include "../Vole/GeneticMaterial.h"
#include "../Skylark/skylarks_all.h"
#include "../Partridge/Partridge_All.h"
#include "../Partridge/Partridge_Population_Manager.h"
#include "../Vole/vole_all.h"
#include "../Vole/VolePopulationManager.h"
#include "../Vole/Predators.h"
#include "../Hare/Hare_All.h"
#include "../Spiders/Spider_BaseClasses.h"
#include "../Spiders/Erigone.h"
#include "../Spiders/Erigone_Population_Manager.h"
#include "../Spiders/Oedothorax.h"
#include "../Spiders/Oedothorax_Population_Manager.h"
#include "../GooseManagement/GooseMemoryMap.h"
#include "../GooseManagement/Goose_Base.h"
#include "../BatchALMaSS/CurveClasses.h"
#include "../People/Hunters/Hunters_all.h"
#include "../GooseManagement/Goose_Population_Manager.h"
#include "../RoeDeer/Roe_all.h"
#include "../RoeDeer/Roe_pop_manager.h"
#include "../Rabbit/Rabbit.h"
#include "../Rabbit/Rabbit_Population_Manager.h"
#include "../Newt/Newt.h"
#include "../Newt/Newt_Population_Manager.h"
#include "../Osmia/Osmia.h"
#include "../Osmia/Osmia_Population_Manager.h"
#include "../HoneyBee/HoneyBee.h"
#include "../HoneyBee/HoneyBee_Colony.h"
#include "../SubPopulation/SubPopulation.h"
#include "../SubPopulation/SubPopulation_Population_Manager.h"
#include "../OliveMoth/OliveMoth.h"
#include "../OliveMoth/OliveMoth_Population_Manager.h"
#include "../Lacewing/Lacewing.h"
#include "../Lacewing/Lacewing_Population_Manager.h"
#include "../Aphid/Aphid.h"
#include "../Aphid/Aphid_Population_Manager.h"
#include "../Bombus/Bombus.h"
#include "../Bombus/Bombus_Population_Manager.h"
#include "../People/Hunters/Disturbers_All.h"
#include "../Common/shared.h"

Macros
#define	_CRT_SECURE_NO_DEPRECATE

Functions
void	RunTheSim ()
void	CloseDownSim ()
bool	ReadBatchINI (const char *filename)
void	SetOutputFiles ()
void	GetProbeInput_ini (std::shared_ptr< Population_Manager_Base >, std::string filename)
bool	CreatePopulationManager ()
Landscape *	CreateLandscape ()
void	delay (int secs)
int	main (int argc, char *argv[])

Variables
const double	g_randmaxp = RAND_MAX + 1.0
static CfgInt	cfg_DayInMonth ("PRB_DAYINMONTH", CFG_CUSTOM, 1)
static CfgInt	cfg_DisturbersNo ("DISTURBERS_NO", CFG_CUSTOM, 1, 0, 1000)
static CfgBool	cfg_UseEasyPop ("VOLE_USEEASYPOP", CFG_CUSTOM, false)
static CfgBool	cfg_VoleCatastrophe_on ("VOLE_CATASTROPHE_ON", CFG_CUSTOM, false)
static CfgInt	cfg_VoleCatastrophe_interval ("VOLE_CATASTROPHE_I", CFG_CUSTOM, 365 *5)
static CfgInt	cfg_VoleCatastrophe_mortality ("VOLE_CATASTROPHE_M", CFG_CUSTOM, 90)
CfgStr	cfg_IniFileName ("INI_FILE", CFG_CUSTOM, "BatchALMaSS.ini")
CfgStr	cfg_AORProbeName ("AOR_PROBE_NAME", CFG_CUSTOM, "AOR_Probe.txt")
CfgStr	cfg_ProbeName ("PROBE_NAME", CFG_CUSTOM, "Probe.res")
CfgBool	cfg_OldProbeNames ("OLDSTYLEPROBENAMES", CFG_CUSTOM, true)
CfgBool	cfg_dumpvegjan
CfgStr	cfg_dumpvegjanfile
CfgBool	cfg_dumpvegjune
CfgStr	cfg_dumpvegjunefile
static CfgInt	cfg_tracelevel ("G_TRACELEVEL", CFG_CUSTOM, 4)
CfgBool	cfg_AorOutput_used
PopulationManagerList	g_PopulationManagerList
CfgInt	cfg_LadybirdAphidPreySpecies
	The number of aphid species in the simulation. More...
CfgStr	cfg_AphidSpeciesName
	This variable decides which aphid species to run. The default one is the English grain aphid. More...
CfgInt	cfg_LadybirdAphidPreyIndex
	The aphid species name. More...
Landscape *	g_ALandscape
std::shared_ptr< Population_Manager_Base >	g_AManager
std::shared_ptr< Hunter_Population_Manager >	g_Hunter_Population_Manager
std::shared_ptr< Disturber_Population_Manager >	g_Disturber_Population_Manager
std::shared_ptr< TPredator_Population_Manager >	g_PredatorManager
int	g_torun
char *	g_files [100]
char *	g_Predfiles [100]
char	g_ResultsDir [255]
char	g_PredResultsDir [255]
int	g_time
int	g_Year
int	g_Steps
int	g_NoProbes
unsigned	g_NoOfPredProbes
ALMaSS_MathFuncs	g_AlmassMathFuncs
TTypesOfPopulation	g_Species

Macro Definition Documentation

_CRT_SECURE_NO_DEPRECATE

#define _CRT_SECURE_NO_DEPRECATE

Function Documentation

CloseDownSim()

void CloseDownSim

(

)

                    {
    if (g_AManager)
    {
        g_ALandscape->SimulationClosingActions();
        // Close the probe file
        if (g_AManager->TheProbe[0] != NULL) g_AManager->TheProbe[0]->CloseFile();
        // delete all probes
        for (int i = 0; i < g_NoProbes; i++) delete g_AManager->TheProbe[i];
    }
    /*
     * the following delete statements will all be removed once we move to smart pointers entirely
     * For now, we can remove the tests, since removing a null pointer has no effect
     * */
 
    delete g_ALandscape;
    delete g_msg;
    delete g_crops;
    delete g_letype;
    delete g_weather;
    delete g_date;
    delete g_pest;
    delete g_nectarpollen;
 
    for ( int NProbes = 0; NProbes < g_NoProbes; NProbes++ ) {
        delete [] g_files[NProbes];
    }
    for ( int NProbes = 0; NProbes < ( int )g_NoOfPredProbes; NProbes++ ) {
        delete [] g_Predfiles[NProbes];
    }
}

References g_ALandscape, g_AManager, g_crops, g_date, g_files, g_letype, g_msg, g_nectarpollen, g_NoOfPredProbes, g_NoProbes, g_pest, g_Predfiles, g_weather, and Landscape::SimulationClosingActions().

Referenced by main().

CreateLandscape()

Landscape * CreateLandscape

(

)

                             {
    g_ALandscape = new Landscape();
    return g_ALandscape;
}

References g_ALandscape.

CreatePopulationManager()

bool CreatePopulationManager

(

)

                               {
  // SET UP THE ANIMAL POPULATION
  // THE LANDSCAPE MUST BE SETUP BEFORE THE CALL HERE
 
    if ( g_Species == TOP_NoSpecies ) {
        g_AManager = std::make_shared<None_Population_Manager>(g_ALandscape);
        // Do not set in g_PopulationManagerList
    }
    else if ( g_Species == TOP_Skylark ) {
        omp_set_num_threads(1);
        g_AManager = std::make_shared<Skylark_Population_Manager>(g_ALandscape);
        g_PopulationManagerList.SetPopulation(g_AManager, TOP_Skylark);
        std::static_pointer_cast<Skylark_Population_Manager>(g_AManager)->OpenTheBreedingPairsProbe();
        std::static_pointer_cast<Skylark_Population_Manager>(g_AManager)->OpenTheBreedingSuccessProbe();
        std::static_pointer_cast<Skylark_Population_Manager>(g_AManager)->OpenTheFledgelingProbe();
    }
    else if ( g_Species == TOP_Vole ) {
        omp_set_num_threads(1);
        g_AManager = std::make_shared<Vole_Population_Manager>(g_ALandscape);
        g_PopulationManagerList.SetPopulation(g_AManager, TOP_Vole);
        g_PredatorManager = std::make_shared<TPredator_Population_Manager> ( g_ALandscape, static_cast<Vole_Population_Manager*>(g_AManager.get()) );
        g_PredatorManager->SetNoProbesAndSpeciesSpecificFunctions( g_NoOfPredProbes);
        g_PopulationManagerList.SetPopulation(g_PredatorManager, TOP_Predators);
    }
  else if ( g_Species == TOP_Erigone ) {
    g_AManager = std::make_shared<Erigone_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Erigone);
 
  }
  else if (g_Species == TOP_Oedothorax) {
        g_AManager = std::make_shared<Oedothorax_Population_Manager>(g_ALandscape);
        g_PopulationManagerList.SetPopulation(g_AManager, TOP_Oedothorax);
    }
  else if ( g_Species == TOP_Bembidion ) {
    g_AManager = std::make_shared<Bembidion_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Bembidion);
 
  }
  else if (g_Species == TOP_PoecilusCupreus) {
        omp_set_num_threads(1);
        //std::shared_ptr<AphidsD_Population_Manager> apMan ( new AphidsD_Population_Manager(g_ALandscape));
        //g_PopulationManagerList.SetPopulation(apMan, TOP_AphidsD);
        g_AManager = std::make_shared<PoecilusCupreus_Population_Manager>(g_ALandscape);
        //std::shared_ptr<PoecilusCupreus_Population_Manager> lbMan(new PoecilusCupreus_Population_Manager(g_ALandscape)) ;
        //lbMan->Init();
        //g_AManager = lbMan;
        g_PopulationManagerList.SetPopulation(g_AManager, TOP_PoecilusCupreus);
 
    }
else if ( g_Species == TOP_Hare ) { // Hare
    omp_set_num_threads(1);
    g_AManager = std::make_shared<THare_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Hare);
 
  }
  else if ( g_Species == TOP_Partridge ) {
    omp_set_num_threads(1);
    g_AManager = std::make_shared<Partridge_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Partridge);
  }
  else if ( g_Species == TOP_Goose ) {
    omp_set_num_threads(1);
    g_AManager = std::make_shared<Goose_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Goose);
    g_Hunter_Population_Manager = std::make_shared<Hunter_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_Hunter_Population_Manager, TOP_Hunters);
    CfgInt cfg_DisturbersSetNo( "DISTURBERS_SETS_NO", CFG_CUSTOM, 1  );
    if (cfg_DisturbersSetNo.value()>0) {
        cout << "Disturbers are ON. Initiating." << endl;
        g_Disturber_Population_Manager = std::make_shared<Disturber_Population_Manager>(cfg_DisturbersNo.value(), g_ALandscape, TOP_Goose);
        g_PopulationManagerList.SetPopulation(g_Disturber_Population_Manager, TOP_Disturbers);
        g_Disturber_Population_Manager->Init();
    }
  }
  else if (g_Species == TOP_RoeDeer)
  {
    omp_set_num_threads(1);
    g_AManager =std::make_shared<RoeDeer_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation( g_AManager, TOP_RoeDeer );
  }
  else if (g_Species == TOP_Rabbit)
  {
    omp_set_num_threads(1);
    g_AManager = std::make_shared<Rabbit_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Rabbit);
  }
  else if (g_Species == TOP_Newt)
  {
    omp_set_num_threads(1);
    g_AManager = std::make_shared<Newt_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Newt);
  }
  else if (g_Species == TOP_Osmia)
  {
      g_AManager = std::make_shared<Osmia_Population_Manager>(g_ALandscape);
      g_PopulationManagerList.SetPopulation(g_AManager, TOP_Osmia);
  }
  else if (g_Species == TOP_ApisRAM)
  {
    omp_set_num_threads(1);
      g_AManager = std::make_shared<HoneyBee_Colony>(g_ALandscape);
      g_PopulationManagerList.SetPopulation(g_AManager, TOP_ApisRAM);
  }
  else if (g_Species == TOP_OliveMoth)
  {
      g_AManager = std::make_shared<OliveMoth_Population_Manager>(g_ALandscape);
      g_PopulationManagerList.SetPopulation(g_AManager, TOP_OliveMoth);
  }
 
  else if (g_Species == TOP_Lacewing)
  {
      g_AManager = std::make_shared<Lacewing_Population_Manager>(g_ALandscape);
      g_PopulationManagerList.SetPopulation(g_AManager, TOP_Lacewing);
  }
 
  else if (g_Species == TOP_Aphid)
  {
      g_AManager = std::make_shared<Aphid_Population_Manager>(g_ALandscape, cfg_AphidSpeciesName.value());
      g_PopulationManagerList.SetPopulation(g_AManager, TOP_Aphid);
  }
  
  else if (g_Species == TOP_Ladybird)
  {
    Aphid_Population_Manager** aphid_manager_ptr;
    if(cfg_LadybirdAphidPreySpecies.value() > 0) {
        string aphid_species_name_list[4] = {"English grain aphid", "Pea aphid",  "Aphis fabae", "Peach potato aphid"};
        aphid_manager_ptr = new Aphid_Population_Manager*[cfg_LadybirdAphidPreySpecies.value()];
        for(int i = 0; i < cfg_LadybirdAphidPreySpecies.value(); i++){
            int aphid_index = 0;
            if(cfg_LadybirdAphidPreySpecies.value() <=1 ){
                aphid_index = cfg_LadybirdAphidPreyIndex.value();
                }
            else{
                aphid_index = i;
            }
            cfg_AphidSpeciesName.set(aphid_species_name_list[aphid_index]);
            std::shared_ptr<Aphid_Population_Manager> apMan ( new Aphid_Population_Manager(g_ALandscape, aphid_species_name_list[aphid_index]));
            g_PopulationManagerList.SetPopulation(apMan, TOP_Aphid);
            aphid_manager_ptr[i] = apMan.get();
        }         
    }
    std::shared_ptr<Ladybird_Population_Manager> lbMan(new Ladybird_Population_Manager(g_ALandscape, 7)) ;
    g_AManager = lbMan;
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Ladybird);
    if(cfg_LadybirdAphidPreySpecies.value() > 0){
        for(int i = 0; i < cfg_LadybirdAphidPreySpecies.value(); i++){
            lbMan->SetAphidPopulationManager(aphid_manager_ptr[i]);
        }
        delete[] aphid_manager_ptr;
    }
  }
 
  else if (g_Species == TOP_Bombus)
  {
    omp_set_num_threads(1);
    g_AManager = std::make_shared<Bombus_Population_Manager>(g_ALandscape);
    g_PopulationManagerList.SetPopulation(g_AManager, TOP_Bombus);
  }
 
  else {
      g_msg->Warn("Population_Manager::CreatePopulationManager()  Unknown species number> ", int(g_Species));
      exit(0);
  }
 
  g_AManager->SetNoProbesAndSpeciesSpecificFunctions(g_NoProbes);
  #ifdef __FLOWERTESTING
    g_ALandscape->SetSpeciesFunctions(TOP_ApisRAM);
  #endif
  return true;
}

References cfg_AphidSpeciesName, CFG_CUSTOM, cfg_DisturbersNo, cfg_LadybirdAphidPreyIndex, cfg_LadybirdAphidPreySpecies, g_ALandscape, g_AManager, g_Disturber_Population_Manager, g_Hunter_Population_Manager, g_msg, g_NoOfPredProbes, g_NoProbes, g_PopulationManagerList, g_PredatorManager, g_Species, CfgStr::set(), PopulationManagerList::SetPopulation(), Landscape::SetSpeciesFunctions(), TOP_Aphid, TOP_ApisRAM, TOP_Bembidion, TOP_Bombus, TOP_Disturbers, TOP_Erigone, TOP_Goose, TOP_Hare, TOP_Hunters, TOP_Lacewing, TOP_Ladybird, TOP_Newt, TOP_NoSpecies, TOP_Oedothorax, TOP_OliveMoth, TOP_Osmia, TOP_Partridge, TOP_PoecilusCupreus, TOP_Predators, TOP_Rabbit, TOP_RoeDeer, TOP_Skylark, TOP_Vole, CfgInt::value(), CfgStr::value(), and MapErrorMsg::Warn().

Referenced by main().

delay()

void delay

(

int

secs

)

                     {
    time_t start_time, cur_time;
    time(&start_time);
    do
    {
       time(&cur_time);
     } while((cur_time - start_time) < secs);
}

GetProbeInput_ini()

void GetProbeInput_ini	(	std::shared_ptr< Population_Manager_Base >	p_PM,
		std::string	filename
	)

                                                                                         {
    ofstream * AFile=NULL;
 
    string NoProbes_String;
 
 
    NoProbes_String.append(g_ResultsDir);
 
    for ( int NProbes = 0; NProbes < g_NoProbes; NProbes++ ) {
        // Must read the probe from a file
        p_PM->TheProbe[ NProbes ] = new Probe_Data;
        p_PM->ProbeFileInput( ( char * ) g_files[ NProbes ], NProbes );
 
 
        NoProbes_String.append(ProbeName);
        if ( NProbes == 0 ) {
            AFile = p_PM->TheProbe[ NProbes ]->OpenFile( NoProbes_String );
            if ( !AFile ) {
                g_msg->Warn( "BatchALMSS - cannot open Probe File", "" );
                exit( 1 );
            }
        } else
            p_PM->TheProbe[ NProbes ]->SetFile( AFile );
    }
    if ( p_PM->SupplySimulationName() == "Vole Predator" ) {
        for ( int NProbes = 0; NProbes < ( int )g_NoOfPredProbes; NProbes++ ) {
            g_PredatorManager->TheProbe[ NProbes ] = new Probe_Data;
            g_PredatorManager->ProbeFileInput( ( char * ) g_Predfiles[ NProbes ], NProbes );
            string nme = g_ResultsDir;
            nme.append("PredProbe");
            nme.append(to_string(NProbes + 1));
            nme.append(".res");
            //      strcpy( Nme, g_PredResultsDir );
            if ( NProbes == 0 ) {
                AFile = g_PredatorManager->TheProbe[ NProbes ]->OpenFile( nme );
                if ( !AFile ) {
                    g_ALandscape->Warn( "BatchALMSS - cannot open Probe File", nme );
                    exit( 1 );
                }
            } else
                g_PredatorManager->TheProbe[ NProbes ]->SetFile( AFile );
        }
    }
 
}

References g_ALandscape, g_files, g_msg, g_NoOfPredProbes, g_NoProbes, g_PredatorManager, g_Predfiles, g_ResultsDir, Probe_Data::OpenFile(), MapErrorMsg::Warn(), and Landscape::Warn().

Referenced by SetOutputFiles().

main()

int main	(	int	argc,
		char *	argv[]
	)

{
#ifdef _DEBUG
    // Get the current state of the flag
    // and store it in a temporary variable
    int tmpFlag = _CrtSetDbgFlag( _CRTDBG_REPORT_FLAG );
 
    // Turn On (OR) - Keep freed memory blocks in the
    // heap's linked list and mark them as freed
    tmpFlag |= _CRTDBG_DELAY_FREE_MEM_DF;
 
    // Turn Off (AND) - prevent _CrtCheckMemory from
    // being called at every allocation request
    tmpFlag &= ~_CRTDBG_CHECK_ALWAYS_DF;
    // Set the new state for the flag
    _CrtSetDbgFlag( tmpFlag );
#endif
 
#if defined(_PROFILE)
    Profiler_Init();
#endif
        
    Sig_Init();
    
    const char* ConfigurationFile;
    std::string ErrorFileName;
    switch (argc)
    {
        case 3:
            cout << "Running with config file: " << argv[1] << "\n";
            ConfigurationFile = argv[1];
            cout << "Error file: " << argv[2] << "\n";
            ErrorFileName = argv[2];
            break;
        case 2:
            cout << "Running with config file: " << argv[1] << "\n";
            ConfigurationFile = argv[1];
            ErrorFileName = "ErrorFile.txt";
            break;
        default:
            ConfigurationFile = "TIALMaSSConfig.cfg";
            ErrorFileName = "ErrorFile.txt";
            break;
 
    }
 
    // Must come first. Used by the configurator below.
    g_msg = CreateErrorMsg();
    g_msg->SetFilename(ErrorFileName);
    g_msg->SetWarnLevel(WARN_MSG);
    g_cfg->ReadSymbols(ConfigurationFile);
    // Configurator instantiation is automatic.
    g_cfg = CreateConfigurator();
 
    //g_cfg->DumpAllSymbolsAndExit( "allsymbols.cfg" );
    // Read the BatchALMaSS.ini to be sure g_Species is set before its needed
    if (!ReadBatchINI(cfg_IniFileName.value())) {
        char ch;
        cout << "Problem with ini file";
        cin >> ch;
    }
 
    g_init_random_seed_fnc();
 
    g_date = CreateCalendar();
    cout << "Created Calendar Object" << endl;
 
    g_weather = CreateWeather();
    cout << "Created Weather Object " << endl;
 
    g_letype =  CreateLETypeClass();
    cout << "Created LE_TypeClass Object" << endl;
 
    g_crops = CreatePlantGrowthData();
    cout << "Created PlantGrowthData Object " << endl;
 
    g_nectarpollen = CreatePollenNectarDevelopmentData();
    cout << "Created PollenNectarDevelopmentData Object" << endl;
    cout << "Creating Landscape Object ..." << endl;
 
    g_ALandscape = new Landscape();
    cout << "Landscape Object Done" << endl;
 
    /* Now latitude and longitude given, or using DK as default. Calculate daylength */
    g_date->CreateDaylength(g_ALandscape->SupplyLatitude(), g_ALandscape->SupplyLongitude(), g_ALandscape->SupplyTimezone() );
 
    g_pest = CreatePesticide(g_ALandscape);
    cout << "Created Pesticide Object" << endl;
 
 
 
    g_ALandscape->SetThePopManagerList(&g_PopulationManagerList);
    if (!CreatePopulationManager())
        return false;
    else
        if (g_Species != -1)
        {
            cout << "Population Manager Created" << endl; // Now got to get the probe files read in
        }
        /* Landscape module is still operating with raw pointers (therefore we give the reference to g_PopulationManagerList):
        * this will be changed in the next stage of operation
        * However the pointers inside the g_PopulationManagerList are smart, this is not a problem, since all reference operations on the
        * smart pointers are the same as on the raw ones
        * */
 
    g_ALandscape->RunHiddenYear();
    SetOutputFiles();
    // Ready to go
    RunTheSim();
    CloseDownSim();
    
#if defined(_PROFILE)
    Profiler_End();
#endif
    
    return 0;
}

References cfg_IniFileName, CloseDownSim(), CreateCalendar(), CreateConfigurator(), Calendar::CreateDaylength(), CreateErrorMsg(), CreateLETypeClass(), CreatePesticide(), CreatePlantGrowthData(), CreatePollenNectarDevelopmentData(), CreatePopulationManager(), CreateWeather(), g_ALandscape, g_cfg, g_crops, g_date, g_init_random_seed_fnc(), g_letype, g_msg, g_nectarpollen, g_pest, g_PopulationManagerList, g_Species, g_weather, ReadBatchINI(), Landscape::RunHiddenYear(), RunTheSim(), MapErrorMsg::SetFilename(), SetOutputFiles(), Landscape::SetThePopManagerList(), MapErrorMsg::SetWarnLevel(), Landscape::SupplyLatitude(), Landscape::SupplyLongitude(), Landscape::SupplyTimezone(), CfgStr::value(), and WARN_MSG.

ReadBatchINI()

bool ReadBatchINI

(

const char *

filename

)

                                        {
  // Must read the TIBatch.INI
  // Read the INI file
  FILE * Fi = NULL;
  char answer = 'X';
  std::cout<<"Batch INI file: "<<filename<<endl;
  Fi=fopen(filename, "r" );
  while ( Fi==NULL ) {
      // Issue and error warning
      cout << "INI File Missing: ";
      cout << filename << endl;
      cout << "Try Again?";
      cin >> answer;
    if ( answer != 'Y' ) exit(0);
    Fi = fopen(filename, "r" );
  }
  char Data[ 255 ];
  fscanf( Fi, "%d\r", & g_NoProbes );
  for ( int i = 0; i < g_NoProbes; i++ ) {
    fscanf( Fi, "%s\r", Data);
    g_files[ i ] = new char[ 255 ];
    strcpy( g_files[ i ], Data );
  }
  fscanf( Fi, "%s\r", Data );
  strcpy( g_ResultsDir, Data );
  fscanf( Fi, "%d\r", & g_torun );
  g_torun *= 365; // the number of years and multiplied by 365 to get days
  fscanf( Fi, "%d\r", & g_Species );
  fclose( Fi );
    // Read the PredBatch.INI file
  if ( g_Species == 1) {
      FILE * Fi2 = NULL;
      answer = 'X';
      while ( !Fi2 ) {
        Fi2=fopen("VoleToxPreds.ini", "r" );
        if ( !Fi2 ) {
          // Issue and error warning
          cout << "Predator Batch File Missing: VoleToxPreds.INI";
          cout << "Try Again  (Y/N) ? ";
          cin >> answer;
          if ( answer != 'Y' ) return false;
        }
      }
      char Data2[ 255 ];
      fscanf( Fi2, "%d\r", & g_NoOfPredProbes );
      for ( int i = 0; i < ( int )g_NoOfPredProbes; i++ ) {
        fscanf( Fi2, "%s\r", Data2 );
        g_Predfiles[ i ] = new char[ 255 ];
        strcpy( g_Predfiles[ i ], Data2 );
      }
      fscanf( Fi2, "%s\r", Data2 );
      strcpy( g_PredResultsDir, Data2 );
      fclose( Fi2 );
  }
  switch (g_Species) {
  case TOP_NoSpecies:
      cout << "Running Landscape only,  animal species not selected" << endl;
      break;
  case TOP_Skylark:
      cout << "Running Skylarks" << endl;
      break;
  case TOP_Vole:
      cout << "Running Voles" << endl;
      break;
  case TOP_Erigone:
      cout << "Running Erigone" << endl;
      break;
  case TOP_Bembidion:
      cout << "Running Bembidon" << endl;
      break;
  case TOP_PoecilusCupreus:
      cout << "Running Poecilus" << endl;
      break;
  case TOP_Hare:
      cout << "Running Hares" << endl;
      break;
  case TOP_Partridge:
      cout << "Running Partridges" << endl;
      break;
  case TOP_Goose:
      cout << "Running Goose Management" << endl;
      break;
  case TOP_RoeDeer:
      cout << "Running Roe Deer" << endl;
      break;
  case TOP_Rabbit:
      cout << "Running Rabbits" << endl;
      break;
  case TOP_Newt:
      cout << "Running Newts" << endl;
      break;
  case TOP_Osmia:
      cout << "Running Osmia" << endl;
      break;
  case TOP_ApisRAM:
      cout <<"Running Honey Bee" <<endl;
      break;
  case TOP_Oedothorax:
      cout << "Running Oedothorax" << endl;
      break;
  case TOP_OliveMoth:
      cout << "Running Olive Moth" << endl;
      break;
  case TOP_Lacewing:
      cout << "Running Lacewing" << endl;
      break;
  case TOP_Aphid:
      cout << "Running Aphid" << endl;
      break;
  case TOP_Ladybird:
      cout << "Running Ladybird" << endl;
      break;
  case TOP_Bombus:
      cout << "Running Bombus" << endl;
      break;
  }
  return true;
}

References g_files, g_NoOfPredProbes, g_NoProbes, g_Predfiles, g_PredResultsDir, g_ResultsDir, g_Species, g_torun, TOP_Aphid, TOP_ApisRAM, TOP_Bembidion, TOP_Bombus, TOP_Erigone, TOP_Goose, TOP_Hare, TOP_Lacewing, TOP_Ladybird, TOP_Newt, TOP_NoSpecies, TOP_Oedothorax, TOP_OliveMoth, TOP_Osmia, TOP_Partridge, TOP_PoecilusCupreus, TOP_Rabbit, TOP_RoeDeer, TOP_Skylark, and TOP_Vole.

Referenced by main().

RunTheSim()

void RunTheSim

(

)

{
    for (int i = 0; i < g_torun; i++)
    {
        Sig_Update();
 
        g_ALandscape->TurnTheWorld();
        // Update the Date
        g_time++;
        int day = g_ALandscape->SupplyDayInMonth();
        int month = g_ALandscape->SupplyMonth();
    
        if ((day == 1) && (month == 1)){
             g_Year++;
             cout << "Year; " << g_Year << '\n';
        }
        cout << g_time << '\n';
 
        if (g_Species == TOP_Vole) 
            g_PredatorManager->Run(1);
 //       if (g_Species== TOP_Ladybird) g_PopulationManagerList.GetPopulation(TOP_AphidsD)->Run(1);
        if ((g_Species == TOP_Goose) || (g_Species == TOP_RoeDeer) || (g_Species == TOP_ApisRAM) || (g_Species == TOP_Bombus))
        {
            for (int tenmin = 0; tenmin < 144; tenmin++)
            {
                for (int p = 0; p < TOP_foobar; p++)
                {
                    int popu_num = g_PopulationManagerList.GetPopulationManagerNum(p);
                    if (popu_num < 1) continue;
                    for (int i = 0; i < popu_num; i++){
                        if (g_PopulationManagerList.GetPopulation(p) != nullptr){
                            g_PopulationManagerList.GetPopulation(p)->Run(1);
                            g_PopulationManagerList.GetPopulation(p)->SetStepCounterInADay(tenmin);
                        }
                    }
                }
 
                if (g_date->TickMinute10()) 
                    g_date->TickHour();
            }
        }
        else {
            for (int p = 0; p < TOP_foobar; p++){
                int popu_num = g_PopulationManagerList.GetPopulationManagerNum(p);
                if (popu_num < 1) continue;
                for (int i = 0; i < popu_num; i++){
                    if (g_PopulationManagerList.GetPopulation(p, i) != nullptr){
                        g_PopulationManagerList.GetPopulation(p, i)->Run(1);
                    }
                }
            }
        }
        //for 1 day-time step, always 0
        g_AManager->SetStepCounterInADay(0);
        if (g_Species != TOP_Aphid && g_Species != TOP_OliveMoth ){
            char str[255];
            strcpy(str, g_AManager->SpeciesSpecificReporting(g_Species, g_time));
        }
        g_ALandscape->DumpVegAreaData(g_time);
    }
}

References Landscape::DumpVegAreaData(), g_ALandscape, g_AManager, g_date, g_PopulationManagerList, g_PredatorManager, g_Species, g_time, g_torun, g_Year, PopulationManagerList::GetPopulation(), PopulationManagerList::GetPopulationManagerNum(), Population_Manager_Base::Run(), Population_Manager_Base::SetStepCounterInADay(), Landscape::SupplyDayInMonth(), Landscape::SupplyMonth(), Calendar::TickHour(), Calendar::TickMinute10(), TOP_Aphid, TOP_ApisRAM, TOP_Bombus, TOP_foobar, TOP_Goose, TOP_OliveMoth, TOP_RoeDeer, TOP_Vole, and Landscape::TurnTheWorld().

Referenced by main().

SetOutputFiles()

void SetOutputFiles

(

)

                     {
    if (cfg_OldProbeNames.value()){
        // Use the regular naming of the Probe.res and AOR_Probe (from the config)
        // If there is no PM set there is nothing to report so we only need to cover the case of the only PM
        for (int i = 0; i<TOP_foobar; i++){
            auto PM = g_PopulationManagerList.GetPopulation_smart(i);
            if (PM!= nullptr){
                if (cfg_AorOutput_used.value()) {
                    PM->OpenTheAOROutputProbe(cfg_AORProbeName.value());
                }
                GetProbeInput_ini(PM, cfg_ProbeName.value());
            }
        }
    }
    else{
    // We create the probes with naming format of simulationname_Probe.res and simulationname_AOR_Probe.txt
        for (int i = 0; i<TOP_foobar; i++){
            auto PM = g_PopulationManagerList.GetPopulation_smart(i);
            if (PM!= nullptr){
                if (cfg_AorOutput_used.value()) {
                    PM->OpenTheAOROutputProbe(PM->SupplySimulationName() + "_" + cfg_AORProbeName.value());
                }
                GetProbeInput_ini(PM, PM->SupplySimulationName()+"_"+cfg_ProbeName.value());
            }
 
        }
    }
}

References cfg_AorOutput_used, cfg_AORProbeName, cfg_OldProbeNames, cfg_ProbeName, g_PopulationManagerList, PopulationManagerList::GetPopulation_smart(), GetProbeInput_ini(), TOP_foobar, CfgBool::value(), and CfgStr::value().

Referenced by main().

Variable Documentation

cfg_AorOutput_used

CfgBool cfg_AorOutput_used

Referenced by SetOutputFiles().

cfg_AORProbeName

CfgStr cfg_AORProbeName("AOR_PROBE_NAME", CFG_CUSTOM, "AOR_Probe.txt")

Referenced by SetOutputFiles().

cfg_AphidSpeciesName

CfgStr cfg_AphidSpeciesName

This variable decides which aphid species to run. The default one is the English grain aphid.

Referenced by CreatePopulationManager().

cfg_DayInMonth

CfgInt cfg_DayInMonth("PRB_DAYINMONTH", CFG_CUSTOM, 1)

static

cfg_DisturbersNo

CfgInt cfg_DisturbersNo("DISTURBERS_NO", CFG_CUSTOM, 1, 0, 1000)

static

Referenced by CreatePopulationManager().

cfg_dumpvegjan

CfgBool cfg_dumpvegjan

Referenced by Landscape::DumpVegAreaData(), and Landscape::Landscape().

cfg_dumpvegjanfile

CfgStr cfg_dumpvegjanfile

Referenced by Landscape::DumpVegAreaData(), and Landscape::Landscape().

cfg_dumpvegjune

CfgBool cfg_dumpvegjune

Referenced by Landscape::DumpVegAreaData(), and Landscape::Landscape().

cfg_dumpvegjunefile

CfgStr cfg_dumpvegjunefile

Referenced by Landscape::DumpVegAreaData(), and Landscape::Landscape().

cfg_IniFileName

CfgStr cfg_IniFileName("INI_FILE", CFG_CUSTOM, "BatchALMaSS.ini")

Referenced by main().

cfg_LadybirdAphidPreyIndex

CfgInt cfg_LadybirdAphidPreyIndex

The aphid species name.

Referenced by CreatePopulationManager().

cfg_LadybirdAphidPreySpecies

CfgInt cfg_LadybirdAphidPreySpecies

The number of aphid species in the simulation.

Referenced by CreatePopulationManager().

cfg_OldProbeNames

CfgBool cfg_OldProbeNames("OLDSTYLEPROBENAMES", CFG_CUSTOM, true)

Referenced by SetOutputFiles().

cfg_ProbeName

CfgStr cfg_ProbeName("PROBE_NAME", CFG_CUSTOM, "Probe.res")

Referenced by SetOutputFiles().

cfg_tracelevel

CfgInt cfg_tracelevel("G_TRACELEVEL", CFG_CUSTOM, 4)

static

cfg_UseEasyPop

CfgBool cfg_UseEasyPop("VOLE_USEEASYPOP", CFG_CUSTOM, false)

static

cfg_VoleCatastrophe_interval

CfgInt cfg_VoleCatastrophe_interval("VOLE_CATASTROPHE_I", CFG_CUSTOM, 365 *5)

static

cfg_VoleCatastrophe_mortality

CfgInt cfg_VoleCatastrophe_mortality("VOLE_CATASTROPHE_M", CFG_CUSTOM, 90)

static

cfg_VoleCatastrophe_on

CfgBool cfg_VoleCatastrophe_on("VOLE_CATASTROPHE_ON", CFG_CUSTOM, false)

static

g_ALandscape

Landscape* g_ALandscape

Referenced by CloseDownSim(), CreateLandscape(), CreatePopulationManager(), GetProbeInput_ini(), main(), and RunTheSim().

g_AlmassMathFuncs

ALMaSS_MathFuncs g_AlmassMathFuncs

g_AManager

std::shared_ptr<Population_Manager_Base> g_AManager

Referenced by CloseDownSim(), CreatePopulationManager(), Field::GetAphidDensity(), and RunTheSim().

g_Disturber_Population_Manager

std::shared_ptr<Disturber_Population_Manager> g_Disturber_Population_Manager

Referenced by CreatePopulationManager().

g_files

char* g_files[100]

Referenced by CloseDownSim(), GetProbeInput_ini(), and ReadBatchINI().

g_Hunter_Population_Manager

std::shared_ptr<Hunter_Population_Manager> g_Hunter_Population_Manager

Referenced by CreatePopulationManager().

g_NoOfPredProbes

unsigned g_NoOfPredProbes

Referenced by CloseDownSim(), CreatePopulationManager(), GetProbeInput_ini(), and ReadBatchINI().

g_NoProbes

int g_NoProbes

Referenced by CloseDownSim(), CreatePopulationManager(), GetProbeInput_ini(), and ReadBatchINI().

g_PopulationManagerList

PopulationManagerList g_PopulationManagerList

Referenced by CreatePopulationManager(), main(), RunTheSim(), and SetOutputFiles().

g_PredatorManager

std::shared_ptr<TPredator_Population_Manager> g_PredatorManager

Referenced by CreatePopulationManager(), GetProbeInput_ini(), and RunTheSim().

g_Predfiles

char* g_Predfiles[100]

Referenced by CloseDownSim(), GetProbeInput_ini(), and ReadBatchINI().

g_PredResultsDir

char g_PredResultsDir[255]

Referenced by ReadBatchINI().

g_randmaxp

const double g_randmaxp = RAND_MAX + 1.0

g_ResultsDir

char g_ResultsDir[255]

Referenced by GetProbeInput_ini(), and ReadBatchINI().

g_Species

TTypesOfPopulation g_Species

Referenced by CreatePopulationManager(), main(), ReadBatchINI(), and RunTheSim().

g_Steps

int g_Steps

g_time

int g_time

THESE MUST BE IN THE CONFIGURATION FILES TO OVERRIDE THE DEFAULT BEETLE VALUES

BEETLE_LARVALSTAGESNUM (int) = 4 BEETLE_MAXDAILYMOVEMENT (int) = 500 BEETLE_DEVELOPMENTTHRESHOLDS (array double) = 11.5 13.8 13.6 13.6 13.9 12.9 12.4 0 0 0 BEETLE_DEVELOPMENTCONSTANTS (array double) = 42 22.8 20.4 23.1 38 63.6 151.9 0 0 0

Referenced by RunTheSim().

g_torun

int g_torun

Referenced by ReadBatchINI(), and RunTheSim().

g_Year

int g_Year

Referenced by RunTheSim().