#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 "../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 "../Bembidion/bembidion_all.h"
#include "../Hare/hare_all.h"
#include "../Spider/spider_all.h"
#include "../Spider/SpiderPopulationManager.h"
#include "../MarshFritillary/MarshFritillaryHdrs.h"
#include "../GooseManagement/GooseMemoryMap.h"
#include "../GooseManagement/Goose_Base.h"
#include "../BatchALMaSS/CurveClasses.h"
#include "../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"

Macros
#define	_CRT_SECURE_NO_DEPRECATE

Functions
bool	BeginningOfMonth ()

static CfgInt	cfg_DayInMonth ("PRB_DAYINMONTH", CFG_CUSTOM, 1)

static CfgBool	cfg_UseEasyPop ("VOLE_USEEASYPOP", 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)

static CfgBool	cfg_VoleCatastrophe_on ("VOLE_CATASTROPHE_ON", CFG_CUSTOM, false)

void	CloseDownSim ()

void	CreateLandscape ()

bool	CreatePopulationManager ()

void	delay (int secs)

void	DumpVegAreaData (int a_day)

void	FloatToDouble (double &d, float f)

void	GetProbeInput_ini ()

void	ImpactProbeReport (int a_time)

int	main ()

void	PredProbeReportDay0 (int a_time)

void	ProbeReport (int a_time)

int	random (int a_range)

bool	ReadBatchINI ()

void	RunTheSim ()

void	SpeciesSpecificActions ()

void	SpeciesSpecificReporting ()

Variables
CfgBool	cfg_dumpvegjan

CfgStr	cfg_dumpvegjanfile

CfgBool	cfg_dumpvegjune

CfgStr	cfg_dumpvegjunefile

CfgBool	cfg_fixed_random_sequence

Landscape *	g_ALandscape

ALMaSS_MathFuncs	g_AlmassMathFuncs

Population_Manager *	g_AManager

char *	g_files [100]

Hunter_Population_Manager *	g_Hunter_Population_Manager

unsigned	g_NoOfPredProbes

int	g_NoProbes

PopulationManagerList	g_PopulationManagerList

TPredator_Population_Manager *	g_PredatorManager

char *	g_Predfiles [100]

char	g_PredResultsDir [255]

boost::variate_generator< base_generator_type &, boost::uniform_real<> >	g_rand_uni

const double	g_randmaxp = RAND_MAX + 1.0

char	g_ResultsDir [255]

int	g_Species

int	g_Steps

int	g_time

int	g_torun

int	g_Year

Macro Definition Documentation

#define _CRT_SECURE_NO_DEPRECATE

Function Documentation

bool BeginningOfMonth ( )

Referenced by Population_Manager::Catastrophe(), Population_Manager::ImpactProbeReport(), and Population_Manager::ProbeReport().

static CfgInt cfg_DayInMonth	(	"PRB_DAYINMONTH"	,
		CFG_CUSTOM	,
		1
	)

static

static CfgBool cfg_UseEasyPop	(	"VOLE_USEEASYPOP"	,
		CFG_CUSTOM	,
		false
	)

static

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

static

static CfgInt cfg_VoleCatastrophe_mortality	(	"VOLE_CATASTROPHE_M"	,
		CFG_CUSTOM	,
		90
	)

static

static CfgBool cfg_VoleCatastrophe_on	(	"VOLE_CATASTROPHE_ON"	,
		CFG_CUSTOM	,
		false
	)

static

void CloseDownSim ( )

References probe_data::CloseFile(), g_ALandscape, g_NoProbes, g_PredatorManager, Landscape::SimulationClosingActions(), and Population_Manager::TheProbe.

Referenced by main().

                     {
     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 ];
       delete g_AManager;
       if ( g_Species == 1 ) {
         delete g_PredatorManager;
       }
     }
         delete g_ALandscape;
     delete g_cfg;
 }

void CreateLandscape ( )

Referenced by main().

                        {
   // Create the landscape
   g_ALandscape = new Landscape( "TIALMaSSConfig.cfg", "ErrorFile.txt" );
 }

bool CreatePopulationManager ( )

References g_Hunter_Population_Manager, g_NoOfPredProbes, g_PredatorManager, Population_Manager::OpenTheBreedingPairsProbe(), Population_Manager::OpenTheBreedingSuccessProbe(), Population_Manager::OpenTheFledgelingProbe(), Population_Manager::SetNoProbes(), PopulationManagerList::SetPopulation(), TOP_Beetle, TOP_Goose, TOP_Hare, TOP_Hunters, TOP_MarshFritillary, TOP_Newt, TOP_Partridge, TOP_Predators, TOP_Rabbit, TOP_RoeDeer, TOP_Skylark, TOP_Spider, and TOP_Vole.

Referenced by main().

                                {
   // SET UP THE ANIMAL POPULATION
   // THE LANDSCAPE MUST BE SETUP BEFORE THE CALL HERE
     if ( g_Species == 0 ) {
                 Skylark_Population_Manager * skMan = new Skylark_Population_Manager( g_ALandscape );
                 g_AManager = skMan;
                 g_PopulationManagerList.SetPopulation(g_AManager, TOP_Skylark);
                 g_AManager->OpenTheBreedingPairsProbe();
                 g_AManager->OpenTheBreedingSuccessProbe();
                 g_AManager->OpenTheFledgelingProbe();
         }
 
     if ( g_Species == 1 ) {
                 Vole_Population_Manager * vMan = new Vole_Population_Manager( g_ALandscape );
                 g_AManager = vMan;
             g_PopulationManagerList.SetPopulation(g_AManager, TOP_Vole);
                 g_PredatorManager = new TPredator_Population_Manager( g_ALandscape, vMan );
                 g_PredatorManager->SetNoProbes( g_NoOfPredProbes);
                 g_PopulationManagerList.SetPopulation(g_PredatorManager, TOP_Predators);
         }
   if ( g_Species == 2 ) {
     Spider_Population_Manager * spMan = new Spider_Population_Manager( g_ALandscape );
         g_PopulationManagerList.SetPopulation(g_AManager, TOP_Spider);
     g_AManager = spMan;
   }
 
   if ( g_Species == 3 ) {
     Bembidion_Population_Manager * bMan = new Bembidion_Population_Manager( g_ALandscape );
         g_PopulationManagerList.SetPopulation(g_AManager, TOP_Beetle);
     g_AManager = bMan;
   }
 
   if ( g_Species == 4 ) { // Hare
     THare_Population_Manager * hMan = new THare_Population_Manager( g_ALandscape );
         g_PopulationManagerList.SetPopulation(g_AManager, TOP_Hare);
     g_AManager = hMan;
   }
 
   if ( g_Species == 5 ) {
     Partridge_Population_Manager * pMan = new Partridge_Population_Manager( g_ALandscape );
     g_AManager = pMan;
         g_PopulationManagerList.SetPopulation(g_AManager, TOP_Partridge);
   }
 
   if ( g_Species == 6 ) {
     Goose_Population_Manager * gMan = new Goose_Population_Manager( g_ALandscape );
     g_AManager = gMan;
         g_PopulationManagerList.SetPopulation(g_AManager, TOP_Goose);
     g_Hunter_Population_Manager = new Hunter_Population_Manager( g_ALandscape );
         g_PopulationManagerList.SetPopulation(g_Hunter_Population_Manager, TOP_Hunters);
   }
 
   if ( g_Species == 7 ) {
     MarshFritillary_Population_Manager * mMan = new  MarshFritillary_Population_Manager( g_ALandscape );
     g_AManager = mMan;
         g_PopulationManagerList.SetPopulation(g_AManager, TOP_MarshFritillary);
   }
 #ifdef __dormouse  //LADA
   if (g_Species == 8) {
           Dormouse_Population_Manager * dMan = new Dormouse_Population_Manager(g_ALandscape);
           g_AManager = dMan;
           g_PopulationManagerList.SetPopulation(g_AManager, TOP_Dormouse);
   }
 #endif // __dormouse
   if (g_Species == 9)
   {
           RoeDeer_Population_Manager * rMan = new RoeDeer_Population_Manager( g_ALandscape );
           g_AManager = rMan;
           g_PopulationManagerList.SetPopulation( g_AManager, TOP_RoeDeer );
   }
 
   if (g_Species == 10)
   {
           Rabbit_Population_Manager * rMan = new Rabbit_Population_Manager(g_ALandscape);
           g_AManager = rMan;
           g_PopulationManagerList.SetPopulation(g_AManager, TOP_Rabbit);
   }
 
   if (g_Species == 11)
   {
           Newt_Population_Manager * rMan = new Newt_Population_Manager(g_ALandscape);
           g_AManager = rMan;
           g_PopulationManagerList.SetPopulation(g_AManager, TOP_Newt);
   }
 
   g_AManager->SetNoProbes(g_NoProbes);
   return true;
 }

void delay ( int secs )

References FloatToDouble().

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

void DumpVegAreaData ( int a_day )

Referenced by Landscape::GetVegArea().

void FloatToDouble	(	double &	d,
		float	f
	)

Referenced by delay(), and VegElement::ReadBugPercentageFile().

                                        {
     char * num = 0;
     num = new char[_CVTBUFSIZE];
     errno_t err = _gcvt_s(num, _CVTBUFSIZE,f,8);
     if (err!=0) {
       assert(0);
     }
     d = atof(num);
         delete [] num;
 }

void GetProbeInput_ini ( )

References g_files, g_NoOfPredProbes, g_NoProbes, g_PredatorManager, g_Predfiles, g_ResultsDir, probe_data::OpenFile(), Population_Manager::ProbeFileInput(), probe_data::SetFile(), Population_Manager::TheProbe, and Landscape::Warn().

Referenced by main().

                          {
   char Nme[ 511 ];
  ofstream * AFile=NULL;
   for ( int NProbes = 0; NProbes < g_NoProbes; NProbes++ ) {
     // Must read the probe from a file
     g_AManager->TheProbe[ NProbes ] = new probe_data;
     g_AManager->ProbeFileInput( ( char * ) g_files[ NProbes ], NProbes );
     char NoProbesString[ 255 ];
     sprintf( NoProbesString, "%sProbe.res", g_ResultsDir );
     if ( NProbes == 0 ) {
       AFile = g_AManager->TheProbe[ NProbes ]->OpenFile( NoProbesString );
       if ( !AFile ) {
         g_ALandscape->Warn( "BatchALMSS - cannot open Probe File", NULL );
         exit( 1 );
       }
         } else
       g_AManager->TheProbe[ NProbes ]->SetFile( AFile );
   }
   if ( g_Species == 1 ) {
     for ( int NProbes = 0; NProbes < ( int )g_NoOfPredProbes; NProbes++ ) {
       g_PredatorManager->TheProbe[ NProbes ] = new probe_data;
       g_PredatorManager->ProbeFileInput( ( char * ) g_Predfiles[ NProbes ], NProbes );
       sprintf( Nme, "%sPredProbe%d.res", g_ResultsDir, NProbes + 1 );
       //      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 );
     }
   }
   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];
   }
 }

void ImpactProbeReport ( int a_time )

Referenced by Population_Manager::Catastrophe(), and Population_Manager::SpeciesSpecificReporting().

int main ( )

References cfg_fixed_random_sequence, CloseDownSim(), CreateLandscape(), CreatePopulationManager(), GetProbeInput_ini(), ReadBatchINI(), RunTheSim(), and Landscape::SetThePopManager().

            {
 #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 (cfg_fixed_random_sequence.value()) srand( 0 ); else srand( (int)time( NULL ) );
         CreateLandscape();
         printf( "Landscape Created\n" ); // So now set up a Population Manager
         if (!ReadBatchINI()) {
                 char ch;
                 cout << "Problem with ini file";
                 cin >> ch;
         }
         if (!CreatePopulationManager()) return false; else
                 printf( "Population Created\n" ); // Now got to get the probe files read in
         g_ALandscape->SetThePopManager( g_AManager );
         GetProbeInput_ini();
         // Ready to go
         RunTheSim();
         CloseDownSim();
         return 0;
 }

void PredProbeReportDay0 ( int a_time )

void ProbeReport ( int a_time )

Referenced by Population_Manager::Catastrophe(), and Population_Manager::SpeciesSpecificReporting().

int random ( int a_range )

References g_rand_uni.

Referenced by Landscape::AddBeetleBanks(), APoint::APoint(), Landscape::BorderScan(), RoadsideVerge::DoDevelopment(), Landscape::Landscape(), VegElement::RandomVegStartValues(), Landscape::ReadPolys2(), Population_Manager::Shuffle(), and UnsprayedFieldMargin::UnsprayedFieldMargin().

                         {
   /* Want to raise exception on this?
          if ( a_range <= 0 )
          return 0;
         */
         //int result = (int)(((double) rand() / g_randmaxp ) * a_range);
         return (int) (g_rand_uni()*a_range);
         //return result;
 }

bool ReadBatchINI ( )

References g_files, g_NoOfPredProbes, g_NoProbes, g_Predfiles, g_PredResultsDir, and g_ResultsDir.

Referenced by main().

                     {
   // Must read the TIBatch.INI
   // Read the INI file
   FILE * Fi = NULL;
   char answer = 'X';
   Fi=fopen("BatchALMaSS.ini", "r" );
   while ( Fi==NULL ) {
       // Issue and error warning
       cout << "INI File Missing: ";
       cout << "BatchALMaSS.ini" << "\n";
           cout << "Try Again?";
       cin >> answer;
     if ( answer != 'Y' ) exit(0);
         Fi = fopen("BatchALMaSS.ini", "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 );
   // 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 );
   }
   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 );
   switch (g_Species) {
           case 0:
                   printf("Running Skylarks\n");
                   break;
           case 1:
                   printf("Running Voles\n");
                   break;
           case 2:
                   printf("Running Spiders\n");
                   break;
           case 3:
                   printf("Running Beetles\n");
                   break;
           case 4:
                   printf("Running Hares\n");
                   break;
           case 5:
                   printf("Running Partridges\n");
                   break;
           case 6:
                   printf("Running Goose Management\n");
                   break;
           case 7:
                   printf("Running Marsh Fritillary\n");
                   break;
           case 8:
                   printf("Running Dormouse\n");
                   break;
           case 9:
                   printf("Running Roe Deer\n");
                   break;
           case 10:
                   printf("Running Rabbits\n");
                   break;
           case 11:
                   printf("Running Newts\n");
                   break;
   }
   return true;
 }

void RunTheSim ( )

References Landscape::DumpVegAreaData(), g_Hunter_Population_Manager, g_PredatorManager, g_torun, Population_Manager::Run(), Population_Manager::SpeciesSpecificReporting(), Landscape::SupplyDayInMonth(), Landscape::SupplyMonth(), and Landscape::TurnTheWorld().

Referenced by main().

                  {
   for ( int i = 0; i < g_torun; i++ ) {
     g_ALandscape->TurnTheWorld();
     // Update the Date
     g_time++;
 
         printf("%d\r",g_time);
     int day = g_ALandscape->SupplyDayInMonth();
     int month = g_ALandscape->SupplyMonth();
     if ( ( day == 1 ) && ( month == 1 ) ) g_Year++;
 
         if (g_Species == 1) g_PredatorManager->Run(1);
         if (g_Species == 6)
         {
                 for (int tenmin = 0; tenmin<144; tenmin++)
                 {
                         g_AManager->Run(1); // Goose Model
                         g_Hunter_Population_Manager->Run(1);
                         if (g_date->TickMinute10()) g_date->TickHour();
                 }
         }
         else g_AManager->Run(1);
 
     char str[255];
         strcpy(str,g_AManager->SpeciesSpecificReporting(g_Species,g_time));
 
     g_ALandscape->DumpVegAreaData( g_time );
   }
 }