#include <plants.h>

Public Member Functions
double	GetLAtotalDiff (double a_ddegs, double a_yddegs, int a_plant, int a_phase)

double	GetLAgreenDiff (double a_ddegs, double a_yddegs, int a_plant, int a_phase)

double	GetHeightDiff (double a_ddegs, double a_yddegs, int a_plant, int a_phase)

double	GetStartValue (int a_veg_type, int a_phase, int a_type)

bool	StartValid (int a_veg_type, int a_phase)

int	GetNumCrops ()

int	VegTypeToCurveNum (TTypesOfVegetation VegReference)

double	GetWeedPercent (TTypesOfVegetation a_letype)

double	GetBugPercentA (TTypesOfVegetation a_letype)

double	GetBugPercentB (TTypesOfVegetation a_letype)

double	GetBugPercentC (TTypesOfVegetation a_letype)

double	GetBugPercentD (TTypesOfVegetation a_letype)

bool	GetNutStatus (int a_plant_num)

bool	GetNutStatusExt (int a_plant)

	CropData (const char *a_cropcurvefile)

	~CropData ()

Private Member Functions
double	FindDiff (double a_ddegs, double a_yddegs, int a_plant, int a_phase, int a_type)

unsigned int	FindCropNum (const char *a_cropcurvefile)

void	SetVegNum (unsigned int a_i, const char *a_cropcurvefile)

void	MakeBugPercentArray (void)

void	ReadBugPercentageFile (void)

Private Attributes
vector< CropGrowth * >	m_growth

vector< int >	m_numbers

int	m_num_crops

FILE *	m_ifile

double *	m_weed_percent

double *	m_bug_percent_a

double *	m_bug_percent_b

double *	m_bug_percent_c

int *	m_bug_percent_d

Detailed Description

Definition at line 80 of file plants.h.

Constructor & Destructor Documentation

CropData::CropData ( const char * a_cropcurvefile )

Definition at line 174 of file plants.cpp.

References FindCropNum(), FloatToDouble(), m_growth, m_ifile, MaxNoInflections, and SetVegNum().

 {
   // Just finds out how many veg curves there are.
   unsigned int NoPlants = FindCropNum( a_vegcurvefile );
 
   for ( unsigned int i=0; i<NoPlants; i++) {
     CropGrowth* temp;
     temp = new CropGrowth;
     m_growth[ i ] = temp;
     SetVegNum( i, a_vegcurvefile );
 
     for (unsigned int j=0; j<5; j++) { // for each growth phase
       // 'Local' index into crop growth curves.
       int lk = 0;
       for (unsigned int k=0; k<MaxNoInflections; k++) {
     // for each inflection point
         int entry;
         fscanf( m_ifile, "%d", &entry );
     float f1=0,f2=0,f0=0;
     if ( entry == -1 ) {
           // Crop start data.
           m_growth[ i ]->m_start_valid[j] = true;
           fscanf( m_ifile, "%g %g %g",&f1,&f0,&f2);
           FloatToDouble(m_growth[ i ]->m_start[j][1],f1);
           FloatToDouble(m_growth[ i ]->m_start[j][0],f0);
           FloatToDouble(m_growth[ i ]->m_start[j][2],f2);
     } else {
       // Add inflection point to normal growth curves.
       m_growth[ i ]->m_dds[j][lk] = (double)entry;
       fscanf( m_ifile, "%g %g %g",&f1,&f0,&f2);
       FloatToDouble(m_growth[ i ]->m_slopes[j][1][lk],f1);
       FloatToDouble(m_growth[ i ]->m_slopes[j][0][lk],f0);
       FloatToDouble(m_growth[ i ]->m_slopes[j][2][lk],f2);
       lk++;
     }
       } // MaxNoInflections
     } // Growth Phases
   } // NoPlants
   fclose( m_ifile );
 }

CropData::~CropData ( )

Definition at line 215 of file plants.cpp.

References m_growth.

 {
   for ( unsigned int i=0; i<m_growth.size(); i++ )
     delete m_growth[i];
 }

Member Function Documentation

unsigned int CropData::FindCropNum ( const char * a_cropcurvefile )

private

Definition at line 124 of file plants.cpp.

References g_msg, m_growth, m_ifile, m_num_crops, m_numbers, MapErrorMsg::Warn(), and WARN_FILE.

Referenced by CropData().

 {
   int NoPlants;
 
   m_numbers.resize(201);
   for ( unsigned int i=0; i<201; i++) {
     m_numbers[ i ] = -1;
   }
 
   m_ifile = fopen(a_cropcurvefile, "r" );
   if (!m_ifile){
     g_msg->Warn(WARN_FILE, "CropData::CropData: Unable to open file",
         a_cropcurvefile );
     exit(1);
   }
 
   fscanf( m_ifile , "%d", &NoPlants );  // How many tables to read in
   m_growth.resize( NoPlants );
   m_num_crops = NoPlants;
 
   return NoPlants;
 }

double CropData::FindDiff	(	double	a_ddegs,
		double	a_yddegs,
		int	a_plant,
		int	a_phase,
		int	a_type
	)

private

Definition at line 41 of file plants.cpp.

References m_growth, m_numbers, and MaxNoInflections.

Referenced by GetHeightDiff(), GetLAgreenDiff(), and GetLAtotalDiff().

 {
   // Check for valid plant number at runtime?
   // This is broken for growth curves where one can risk passing
   // more than a single inflection point in the growth curve in a
   // single day...
 
   int index = m_numbers[ a_plant ];
   unsigned int oldindex=0, newindex=0;
 
   if ( m_growth[ index ]->m_dds[ a_phase ][ 0 ] == 99999 ) {
     return 0.0;
   }
 
   for ( unsigned int i=0; i<MaxNoInflections; i++ ) {
     // In other words: If the current value for summed day degrees
     // is smaller than the X position of the *next* inflection
     // point, then we are in the correct interval.
     if ( m_growth[ index ]->m_dds[ a_phase ][ i+1 ] > a_ddegs ) {
       newindex = i;
       break;
       //      return m_growth[ index ]->m_slopes[ a_phase ][ a_type ][i];
     }
   }
 
   for ( unsigned int i=0; i<MaxNoInflections; i++ ) {
     if ( m_growth[ index ]->m_dds[ a_phase ][ i+1 ] > a_yddegs ) {
       oldindex = i;
       break;
       //      return m_growth[ index ]->m_slopes[ a_phase ][ a_type ][i];
     }
   }
 
   double diff;
 
   if ( newindex > oldindex ) {
     // We have passed an inflection point between today and yesterday.
     // First add the increment from yesterdays day degree sum up to
     // the inflection point.
     double dddif =
       m_growth[ index ]->m_dds[ a_phase ][ newindex ] - a_yddegs;
     diff =
       m_growth[ index ]->m_slopes[ a_phase ][ a_type ][oldindex]*
       dddif;
 
     // Then from the inflection point up to today.
     dddif = a_ddegs -
       m_growth[ index ]->m_dds[ a_phase ][ newindex ];
     diff +=
       m_growth[ index ]->m_slopes[ a_phase ][ a_type ][ newindex ]*
       dddif;
   } else {
     // No inflection point passed.
     diff = m_growth[ index ]->m_slopes[ a_phase ][ a_type ][ newindex ] *
       (a_ddegs - a_yddegs);
   }
   return diff;
 }

double CropData::GetBugPercentA ( TTypesOfVegetation a_letype )

inline

Definition at line 129 of file plants.h.

                                                        {
     return m_bug_percent_a[ a_letype ];
   }

double CropData::GetBugPercentB ( TTypesOfVegetation a_letype )

inline

Definition at line 134 of file plants.h.

                                                        {
     return m_bug_percent_b[ a_letype ];
   }

double CropData::GetBugPercentC ( TTypesOfVegetation a_letype )

inline

Definition at line 139 of file plants.h.

                                                        {
     return m_bug_percent_c[ a_letype ];
   }

double CropData::GetBugPercentD ( TTypesOfVegetation a_letype )

inline

Definition at line 144 of file plants.h.

                                                        {
     return (double) m_bug_percent_d[ a_letype ];
   }

double CropData::GetHeightDiff	(	double	a_ddegs,
		double	a_yddegs,
		int	a_plant,
		int	a_phase
	)

Definition at line 116 of file plants.cpp.

References FindDiff().

Referenced by VegElement::DoDevelopment(), and VegElement::ForceGrowthSpringTest().

 {
   return FindDiff( a_ddegs , a_yddegs, a_plant, a_phase, 2 );
 }

double CropData::GetLAgreenDiff	(	double	a_ddegs,
		double	a_yddegs,
		int	a_plant,
		int	a_phase
	)

Definition at line 102 of file plants.cpp.

References FindDiff().

Referenced by VegElement::DoDevelopment(), and VegElement::ForceGrowthSpringTest().

 {
   return FindDiff( a_ddegs, a_yddegs, a_plant, a_phase, 0 );
 }

double CropData::GetLAtotalDiff	(	double	a_ddegs,
		double	a_yddegs,
		int	a_plant,
		int	a_phase
	)

Definition at line 109 of file plants.cpp.

References FindDiff().

Referenced by VegElement::DoDevelopment(), UnsprayedFieldMargin::DoDevelopment(), and VegElement::ForceGrowthSpringTest().

 {
   return FindDiff( a_ddegs, a_yddegs, a_plant, a_phase, 1 );
 }

int CropData::GetNumCrops ( )

inline

Definition at line 113 of file plants.h.

Referenced by CropRotation::GetFirstCrop().

                     {
     return m_num_crops;
   }

bool CropData::GetNutStatus ( int a_plant_num )

inline

Definition at line 148 of file plants.h.

Referenced by CropRotation::GetFirstCrop().

                                        {
     return m_growth[ a_plant_num ]->m_lownut;
   }

bool CropData::GetNutStatusExt ( int a_plant )

inline

Definition at line 152 of file plants.h.

                                       {
     return m_growth[ m_numbers[ a_plant ]]->m_lownut;
   }

double CropData::GetStartValue	(	int	a_veg_type,
		int	a_phase,
		int	a_type
	)

inline

Definition at line 105 of file plants.h.

Referenced by VegElement::ForceGrowthInitialize(), and VegElement::SetGrowthPhase().

                                                                   {
     return m_growth[ m_numbers[ a_veg_type ]]->m_start[ a_phase ] [ a_type ];
   }

double CropData::GetWeedPercent ( TTypesOfVegetation a_letype )

inline

Definition at line 124 of file plants.h.

                                                        {
     return m_weed_percent[ a_letype ];
   }

void CropData::MakeBugPercentArray ( void )

private

void CropData::ReadBugPercentageFile ( void )

private

void CropData::SetVegNum	(	unsigned int	a_i,
		const char *	a_cropcurvefile
	)

private

Definition at line 149 of file plants.cpp.

References g_msg, m_growth, m_ifile, m_numbers, MapErrorMsg::Warn(), and WARN_FILE.

Referenced by CropData().

 {
   int ThisPlant;
 
   // Find out what crop and what nutrient status
   fscanf( m_ifile, "%d", &ThisPlant);
 
   // Check if valid plant number (from the file).
   if ( ThisPlant < 0 || ThisPlant > 200 ) {
     g_msg->Warn(WARN_FILE, "CropData::FindCropNum(): Illegal plant number"
         " specified in", a_cropcurvefile );
     exit(1);
   }
 
   m_numbers[ ThisPlant ] = a_i;
 
   // if greater than 100 then it is low nutrient
   if ( ThisPlant > 100 ) {
     m_growth[ a_i ]->m_lownut = true;
   } else {
     m_growth[ a_i ]->m_lownut = false;
   }
 }

bool CropData::StartValid	(	int	a_veg_type,
		int	a_phase
	)

Definition at line 364 of file plants.cpp.

References m_growth, and m_numbers.

Referenced by VegElement::ForceGrowthInitialize(), and VegElement::SetGrowthPhase().

 {
   int a=m_numbers[ a_veg_type ];
   CropGrowth* p=m_growth[a];
   return p-> m_start_valid[ a_phase ];
 }

int CropData::VegTypeToCurveNum ( TTypesOfVegetation VegReference )

Definition at line 231 of file plants.cpp.

References g_msg, tov_AgroChemIndustryCereal, tov_BroadBeans, tov_Carrots, tov_CloverGrassGrazed1, tov_CloverGrassGrazed2, tov_FieldPeas, tov_FieldPeasSilage, tov_FieldPeasStrigling, tov_FodderBeet, tov_FodderGrass, tov_Heath, tov_Lawn, tov_Maize, tov_MaizeSilage, tov_MaizeStrigling, tov_NaturalGrass, tov_NoGrowth, tov_Oats, tov_OBarleyPeaCloverGrass, tov_OCarrots, tov_OCloverGrassGrazed1, tov_OCloverGrassGrazed2, tov_OCloverGrassSilage1, tov_OFieldPeas, tov_OFieldPeasSilage, tov_OFodderBeet, tov_OGrazingPigs, tov_OMaizeSilage, tov_OOats, tov_OPermanentGrassGrazed, tov_OPotatoes, tov_OrchardCrop, tov_OSBarleySilage, tov_OSeedGrass1, tov_OSeedGrass2, tov_OSetaside, tov_OSpringBarley, tov_OSpringBarleyClover, tov_OSpringBarleyExt, tov_OSpringBarleyGrass, tov_OSpringBarleyPigs, tov_OTriticale, tov_OWinterBarley, tov_OWinterBarleyExt, tov_OWinterRape, tov_OWinterRye, tov_OWinterWheat, tov_OWinterWheatUndersown, tov_PermanentGrassGrazed, tov_PermanentGrassLowYield, tov_PermanentGrassTussocky, tov_PermanentSetaside, tov_Potatoes, tov_PotatoesIndustry, tov_SeedGrass1, tov_SeedGrass2, tov_Setaside, tov_SpringBarley, tov_SpringBarleyCloverGrass, tov_SpringBarleyCloverGrassStrigling, tov_SpringBarleyPeaCloverGrassStrigling, tov_SpringBarleyPTreatment, tov_SpringBarleySeed, tov_SpringBarleySilage, tov_SpringBarleySKManagement, tov_SpringBarleySpr, tov_SpringBarleyStrigling, tov_SpringBarleyStriglingCulm, tov_SpringBarleyStriglingSingle, tov_SpringRape, tov_SugarBeet, tov_Triticale, tov_Wasteland, tov_WinterBarley, tov_WinterBarleyStrigling, tov_WinterRape, tov_WinterRapeStrigling, tov_WinterRye, tov_WinterRyeStrigling, tov_WinterWheat, tov_WinterWheatShort, tov_WinterWheatStrigling, tov_WinterWheatStriglingCulm, tov_WinterWheatStriglingSingle, tov_WWheatPControl, tov_WWheatPToxicControl, tov_WWheatPTreatment, tov_YoungForest, MapErrorMsg::Warn(), and WARN_FILE.

Referenced by BeetleBank::BeetleBank(), FieldBoundary::FieldBoundary(), Heath::Heath(), HedgeBank::HedgeBank(), Hedges::Hedges(), Marsh::Marsh(), MownGrass::MownGrass(), NaturalGrassDry::NaturalGrassDry(), NaturalGrassWet::NaturalGrassWet(), Orchard::Orchard(), OrchardBand::OrchardBand(), PermanentSetaside::PermanentSetaside(), VegElement::SetVegType(), UnsprayedFieldMargin::UnsprayedFieldMargin(), VegElement::VegElement(), Wasteland::Wasteland(), and YoungForest::YoungForest().

 {
   char error_num[20];
 
   switch (VegReference)
     {
     case tov_OSpringBarleyPigs:
     case tov_OSpringBarley:
     case tov_OSpringBarleyExt:
         return 101;
     case tov_SpringBarleySilage:
     case tov_SpringBarley:
     case tov_SpringBarleySpr:
     case tov_SpringBarleyPTreatment:
     case tov_SpringBarleySKManagement:
     case tov_SpringBarleyStrigling:
     case tov_SpringBarleyStriglingSingle:
     case tov_SpringBarleyStriglingCulm:
       return 1;
     case tov_WinterBarley:
     case tov_WinterBarleyStrigling:
       return 2;
     case tov_OWinterBarley:
     case tov_OWinterBarleyExt:
       return 102;
     case tov_WinterWheat:
     case tov_WinterWheatShort:
     case tov_WinterWheatStrigling:
     case tov_WinterWheatStriglingCulm:
     case tov_WinterWheatStriglingSingle:
     case tov_AgroChemIndustryCereal:
     case tov_WWheatPControl:
     case tov_WWheatPToxicControl:
     case tov_WWheatPTreatment:
       return 4;
     case tov_OWinterWheat:
     case tov_OWinterWheatUndersown:
         return 104;
     case tov_WinterRye:
     case tov_WinterRyeStrigling:
       return 5;
     case tov_OWinterRye:
     return 105;
     case tov_Oats:
       return 6;
     case tov_OOats:
       return 106;
     case tov_Maize:
     case tov_MaizeSilage:
     case tov_MaizeStrigling:
       return 8;
     case tov_OMaizeSilage:
       return 108;
     case tov_SpringBarleyCloverGrass:
     case tov_SpringBarleyCloverGrassStrigling:
     case tov_SpringBarleySeed:
     case tov_SpringBarleyPeaCloverGrassStrigling:
       return 13;
     case tov_OBarleyPeaCloverGrass:
     case tov_OSBarleySilage:
     case tov_OSpringBarleyGrass:
     case tov_OSpringBarleyClover:
       return 113;
     case tov_WinterRape:
     case tov_WinterRapeStrigling:
       return 22;
     case tov_OWinterRape:
       return 22;
     case tov_PermanentGrassLowYield:
     case tov_PermanentGrassTussocky:
         return 25;
     case tov_PermanentGrassGrazed:
     case tov_OPermanentGrassGrazed:
       return 26;
     case tov_SeedGrass1:
     case tov_SeedGrass2:
     case tov_OSeedGrass1:
     case tov_OSeedGrass2:
         return 27;
     case tov_FodderGrass:
     case tov_CloverGrassGrazed1:
     case tov_CloverGrassGrazed2:
     case tov_OCloverGrassGrazed1:
     case tov_OCloverGrassGrazed2:
     case tov_OCloverGrassSilage1:
     case tov_OGrazingPigs:
         return 29;
     case tov_OFieldPeas:
     case tov_OFieldPeasSilage:
     case tov_FieldPeas:
     case tov_FieldPeasStrigling:
     case tov_FieldPeasSilage:
     case tov_BroadBeans:
       return 30;
     case tov_Carrots:
       return 41;
     case tov_OCarrots:
       return 141;
     case tov_Potatoes: return 50;
     case tov_OPotatoes: return 150;
     case tov_PotatoesIndustry: return 50;
     case tov_SugarBeet:
     case tov_FodderBeet:
     case tov_OFodderBeet:
         return 60;
       // Special growth mode for green but unused elements.
       // tov_PermanentSetaside Does not change growth phase no matter
       // how hard one tries to do just that.
     case tov_PermanentSetaside: return 92;
     case tov_Heath:
     case tov_Setaside:
     case tov_OSetaside: return 112;
     case tov_OrchardCrop:
     case tov_YoungForest:
     case tov_NaturalGrass:
     case tov_Wasteland: return 90;
     case tov_NoGrowth: return 91;
     case tov_Lawn: return 94;
       //case tov_SpringWheat: return 3;
     case tov_OTriticale:
     case tov_Triticale: return 7;
     case tov_SpringRape: return 21;
 
     default: // No matching code so we need an error message of some kind
       sprintf( error_num, "%d", VegReference );
       g_msg->Warn( WARN_FILE,
            "CropData::VegTypeToCurveNum(): Unknown vegetation type:",
            error_num );
       exit( 1 );
     }
 }

Member Data Documentation

double* CropData::m_bug_percent_a

private

Definition at line 86 of file plants.h.

double* CropData::m_bug_percent_b

private

Definition at line 87 of file plants.h.

double* CropData::m_bug_percent_c

private

Definition at line 88 of file plants.h.

int* CropData::m_bug_percent_d

private

Definition at line 89 of file plants.h.

vector< CropGrowth * > CropData::m_growth

private

Definition at line 81 of file plants.h.

Referenced by CropData(), FindCropNum(), FindDiff(), SetVegNum(), StartValid(), and ~CropData().

FILE* CropData::m_ifile

private

Definition at line 84 of file plants.h.

Referenced by CropData(), FindCropNum(), and SetVegNum().

int CropData::m_num_crops

private

Definition at line 83 of file plants.h.

Referenced by FindCropNum().

vector< int > CropData::m_numbers

private

Definition at line 82 of file plants.h.

Referenced by FindCropNum(), FindDiff(), SetVegNum(), and StartValid().

double* CropData::m_weed_percent

private

Definition at line 85 of file plants.h.

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

D:/CJT/Repos/ALMaSS_all/Landscape/plants.h
D:/CJT/Repos/ALMaSS_all/Landscape/plants.cpp

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation