#include <weather.h>

Public Member Functions
	Weather (const char *a_weatherdatafile)

	~Weather (void)

double	GetDDDegs (long a_date)

double	GetGlobalRadiation (long a_date)

double	GetGlobalRadiation (void)

double	GetTemp (long a_date)
	Get the temperature on a particular date. More...

double	GetTemp (void)
	Get the temperature today. More...

double	GetHumidity (void)
	Get the humidity score today. More...

double	GetMeanTemp (long a_date, unsigned int a_period)

double	GetRain (long a_date)

double	GetRain (void)

double	GetWind (long a_date)

double	GetWind (void)

int	GetWindDirection (void)

bool	GetSnow (long a_date)

bool	GetSnow (void)

double	GetSnowDepth (void)
	Get the current snow depth. More...

bool	Raining (void)

double	GetRainPeriod (long a_date, unsigned int a_period)

double	GetWindPeriod (long a_date, unsigned int a_period)

double	GetTempPeriod (long a_date, unsigned int a_period)

void	Tick (void)

Private Attributes
vector< double >	m_rain

vector< double >	m_wind

vector< double >	m_winddir

vector< double >	m_temp

double	m_temptoday

double	m_raintoday

double	m_windtoday

int	m_winddirtoday

bool	m_snowtoday

bool	m_rainingtoday

double	m_insolation

double	m_humiditytoday

bool	m_wind_valid

bool	m_winddir_valid

long	m_datemodulus

double	m_snowdepth
	The snow depth in cm. More...

Detailed Description

Definition at line 385 of file weather.h.

Constructor & Destructor Documentation

◆ Weather()

Weather::Weather ( const char * a_weatherdatafile )

Definition at line 167 of file weather.cpp.

 {
   int NoDays, Day, Month, Year, FirstYear=0, LastYear=0;
   double Temp, Rain, Wind;
  
   ifstream inFile(a_weatherdatafile);
   if (!inFile) {
     g_msg->Warn(WARN_FILE, "Weather::Weather(): Unable to open file",
         a_weatherdatafile );
     exit(1);
   }
  
   string buf;
   std::getline(inFile,buf);
   NoDays=std::stoi(buf);
  
  
   m_rain.resize( NoDays );
   m_wind.resize( NoDays );
   m_temp.resize( NoDays );
  
   bool storing = false; // Whether we are storing weather data.
   unsigned int index = 0;
  
   g_date->SetLastYear( 0 );
  
  
   for ( int i=0; i<NoDays; i++) {
       inFile >> Year >> Month >> Day >> Temp >> Wind >> Rain;
  
     if ( Month == 2 && Day == 29 ) {
       // Skip leap days.
       continue;
     }
  
     if ( Month == 1 && Day == 1 && !storing &&
      (
       l_weather_starting_year.value() == 0 ||
       l_weather_starting_year.value() == Year
       )
      ) {
       // Commence storing of data from Jan. 1st of the first
       // year requested.
       storing = true;
       g_date->SetFirstYear( Year );
       FirstYear = Year;
       g_date->Reset();
     }
  
  
     if ( storing ) {
             //if (Temp > 0.0 )
       //    DDegs += Temp;
       m_rain[ index ] = Rain;
       m_wind[ index ] = Wind;
       m_temp[ index ] = Temp;
       index++;
     }
  
     //cout << "i: " << i << " index: " << index << " Month: " << Month << " Year: " << Year <<"\n";
  
     if ( (Month == 12) && (Day == 31) && (storing) ) {
         cout << "LastYear: " << Year << "\n";
       // Found yet another full year worth of weather data.
       g_date->SetLastYear( Year );
       LastYear = Year;
     }
   }
  
   // Did we find at least one full year worth of data?
   if ( g_date->GetLastYear() == 0 ) {
     // Nope...
     g_msg->Warn(WARN_FILE, "Weather::Weather(): Weather data file did",
         "not contain at least one year of data!" );
     exit(1);
   }
   //fclose( ifile );
   m_datemodulus = (LastYear - FirstYear + 1)*365;
   m_snowdepth = 0;
   Tick();
 }

References g_date, g_msg, Calendar::GetLastYear(), l_weather_starting_year, m_datemodulus, m_rain, m_snowdepth, m_temp, m_wind, Calendar::Reset(), Calendar::SetFirstYear(), Calendar::SetLastYear(), Tick(), CfgInt::value(), MapErrorMsg::Warn(), and WARN_FILE.

◆ ~Weather()

Weather::~Weather ( void )

Definition at line 249 of file weather.cpp.

250 {

251 }

Member Function Documentation

◆ GetDDDegs()

double Weather::GetDDDegs ( long a_date )

Definition at line 155 of file weather.cpp.

 {
   double temp = m_temp[ a_date%m_datemodulus ];
   if ( temp < 0.0 ) {
     temp = 0.0;
   }
  
   return temp;
 }

References m_datemodulus, and m_temp.

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

◆ GetGlobalRadiation() [1/2]

double Weather::GetGlobalRadiation ( long a_date )

inline

Definition at line 414 of file weather.h.

414 { return (double)c_insolation[a_date % m_datemodulus]; }

References c_insolation.

Referenced by Landscape::SupplyGlobalRadiation().

◆ GetGlobalRadiation() [2/2]

double Weather::GetGlobalRadiation ( void )

inline

Definition at line 415 of file weather.h.

415 { return m_insolation; }

◆ GetHumidity()

double Weather::GetHumidity ( void )

inline

Get the humidity score today.

Definition at line 421 of file weather.h.

421 { return m_humiditytoday; }

Referenced by Landscape::SupplyHumidity().

◆ GetMeanTemp()

double Weather::GetMeanTemp	(	long	a_date,
		unsigned int	a_period
	)

Definition at line 120 of file weather.cpp.

 {
   double sum = 0.0;
  
   for ( unsigned int i=0; i<a_period; i++ )
     sum += GetTemp( a_date - i );
  
   return sum/(double)a_period;
 }

References GetTemp().

Referenced by Landscape::SupplyMeanTemp().

◆ GetRain() [1/2]

double Weather::GetRain ( long a_date )

inline

Definition at line 423 of file weather.h.

423 { return m_rain[a_date %m_datemodulus]; }

Referenced by Landscape::SupplyRain().

◆ GetRain() [2/2]

double Weather::GetRain ( void )

inline

Definition at line 424 of file weather.h.

424 { return m_raintoday; }

Referenced by GetRainPeriod(), and Tick().

◆ GetRainPeriod()

double Weather::GetRainPeriod	(	long	a_date,
		unsigned int	a_period
	)

Definition at line 253 of file weather.cpp.

 {
   double sum = 0.0;
  
   for ( unsigned int i=0; i<a_period; i++ )
     sum += GetRain( a_date - i );
   return sum;
 }

References GetRain().

Referenced by Farm::BulbHarvest(), Farm::BurnStrawStubble(), Farm::CutToHay(), Farm::Harvest(), OptimisingFarm::Harvest(), Farm::HarvestLong(), Farm::HayBailing(), Farm::HayTurning(), Farm::RowCultivation(), Farm::StrawChopping(), Farm::Strigling(), Farm::StriglingSow(), Farm::StubbleHarrowing(), and Landscape::SupplyRainPeriod().

◆ GetSnow() [1/2]

bool Weather::GetSnow ( long a_date )

Definition at line 132 of file weather.cpp.

 {
   if ( a_date == g_date->Date() ) {
     return m_snowtoday;
   }
  
   int weatherday = a_date % m_datemodulus;
   bool snow = false;
  
   if ( m_temp[ weatherday ] < 0 &&
        m_rain[ weatherday ] > 1 &&
        random(100) < 50 ) {
     snow = true;
   }
  
   //  if ( ((dayinyear<90) || (dayinyear>330)) && (rand()%4==1) )
   //  snow = true;
  
   return snow;
 }

References Calendar::Date(), g_date, m_datemodulus, m_rain, m_snowtoday, and m_temp.

Referenced by Landscape::SupplySnowcover().

◆ GetSnow() [2/2]

bool Weather::GetSnow ( void )

inline

Definition at line 429 of file weather.h.

429 { return m_snowtoday; }

◆ GetSnowDepth()

double Weather::GetSnowDepth ( void )

inline

Get the current snow depth.

Definition at line 431 of file weather.h.

431 { return m_snowdepth; }

Referenced by Landscape::SupplySnowDepth().

◆ GetTemp() [1/2]

double Weather::GetTemp ( long a_date )

inline

Get the temperature on a particular date.

Definition at line 417 of file weather.h.

417 { return m_temp[a_date % m_datemodulus]; }

Referenced by Farm::FA_Manure(), Farm::FA_Sludge(), Farm::FA_Slurry(), Farm::FP_Manure(), Farm::FP_Sludge(), Farm::FP_Slurry(), Landscape::SupplyTemp(), and Pesticide::Tick().

◆ GetTemp() [2/2]

double Weather::GetTemp ( void )

inline

Get the temperature today.

Definition at line 419 of file weather.h.

419 { return m_temptoday; }

Referenced by GetMeanTemp(), and GetTempPeriod().

◆ GetTempPeriod()

double Weather::GetTempPeriod	(	long	a_date,
		unsigned int	a_period
	)

Sums the temperature for the period from a_date back a_period days.

Parameters

[in]	a_date	the day to start summing degrees
[in]	a_period	the number of days períod to sum

Returns: The sum of day degrees

Definition at line 272 of file weather.cpp.

 {
   double sum = 0.0;
  
   for ( unsigned int i=0; i<a_period; i++ )
     sum += GetTemp( a_date - i );
  
   return sum;
 }

References GetTemp().

Referenced by Landscape::SupplyTempPeriod(), and Tick().

◆ GetWind() [1/2]

double Weather::GetWind ( long a_date )

inline

Definition at line 425 of file weather.h.

425 { return m_wind[a_date % m_datemodulus]; }

Referenced by Farm::FungicideTreat(), OptimisingFarm::FungicideTreat(), Farm::GrowthRegulator(), Farm::HerbicideTreat(), OptimisingFarm::HerbicideTreat(), Farm::InsecticideTreat(), OptimisingFarm::InsecticideTreat(), Farm::Molluscicide(), Farm::ProductApplication(), and Landscape::SupplyWind().

◆ GetWind() [2/2]

double Weather::GetWind ( void )

inline

Definition at line 426 of file weather.h.

426 { return m_windtoday; }

Referenced by GetWindPeriod().

◆ GetWindDirection()

int Weather::GetWindDirection ( void )

inline

Definition at line 427 of file weather.h.

427 { return m_winddirtoday; }

Referenced by Landscape::SupplyWindDirection().

◆ GetWindPeriod()

double Weather::GetWindPeriod	(	long	a_date,
		unsigned int	a_period
	)

Definition at line 262 of file weather.cpp.

 {
   double sum = 0.0;
  
   for ( unsigned int i=0; i<a_period; i++ )
     sum += GetWind( a_date - i );
  
   return sum;
 }

References GetWind().

Referenced by Landscape::SupplyWindPeriod().

◆ Raining()

bool Weather::Raining ( void )

inline

Definition at line 432 of file weather.h.

432 { return m_rainingtoday; }

Referenced by Farm::AutumnHarrow(), Farm::AutumnPlough(), Farm::AutumnRoll(), Farm::AutumnSow(), Farm::BedForming(), Farm::Biocide(), Farm::CutOrch(), Farm::CutToSilage(), Farm::CutWeeds(), Farm::DeepPlough(), Farm::FA_AmmoniumSulphate(), Farm::FA_Calcium(), Farm::FA_GreenManure(), Farm::FA_ManganeseSulphate(), Farm::FA_Manure(), Farm::FA_NPK(), Farm::FA_NPKS(), Farm::FA_PK(), Farm::FA_RSM(), Farm::FA_Sludge(), Farm::FA_Slurry(), Farm::FlowerCutting(), Farm::FP_AmmoniumSulphate(), Farm::FP_Calcium(), Farm::FP_GreenManure(), Farm::FP_LiquidNH3(), Farm::FP_ManganeseSulphate(), Farm::FP_Manure(), Farm::FP_NPK(), Farm::FP_NPKS(), Farm::FP_PK(), Farm::FP_RSM(), Farm::FP_Sludge(), Farm::FP_Slurry(), Farm::FungicideTreat(), OptimisingFarm::FungicideTreat(), Farm::GrowthRegulator(), Farm::HeavyCultivatorAggregate(), Farm::HerbicideTreat(), OptimisingFarm::HerbicideTreat(), Farm::HillingUp(), Farm::InsecticideTreat(), OptimisingFarm::InsecticideTreat(), Farm::Molluscicide(), Farm::PreseedingCultivator(), Farm::PreseedingCultivatorSow(), Farm::ProductApplication(), Farm::ShallowHarrow(), Farm::SpringHarrow(), Farm::SpringPlough(), Farm::SpringRoll(), Farm::SpringSow(), Farm::SpringSowWithFerti(), Farm::StrawCovering(), Farm::StrawRemoval(), Farm::StriglingHill(), Farm::StubbleCultivatorHeavy(), Farm::StubblePlough(), Farm::Swathing(), Farm::Water(), and Farm::WinterPlough().

◆ Tick()

void Weather::Tick ( void )

The complicated calculation of weatherday is due to the potential to run off either end of the rain data, especially since this is called on Day0

Humidity calculation. This is based on an assumption that humidity is relative to rainfall in the days before. The humidity score is the mean of rainfall over the last 5 days assuming 50% loss each day. Therefore rainfall of 10,0,0,2,0 mm in the previous 5 days would give a humidity of 0.625+0+0+1+0 = 1.625/5 = 0.325 But we also need to calculate temperature into this. Lets assume <10 degrees has no effect, 10-20 and this figure is halved, 20-40 and it is half again.

datemodulus is added to prevent negative overun on the rain or temperature data array, important since this is called on Day0

Wind directions is based on probabilities. Can also use the WindDirections array which allows montly probabilities

Definition at line 58 of file weather.cpp.

                          {
     int weatherday = (g_date->Date() + m_datemodulus) % m_datemodulus;
     double snowtempthreshold = -1.0;
     if (m_snowdepth > 0.0 && m_temp[ weatherday ] < snowtempthreshold) {
         m_snowdepth -= 0.1;  // We decay the snow depth by 0.1 cm per day when temperatures are sub zero
     }
     else if (m_snowdepth > 0.0) // If temperatures are above 0.0, we decay snow 1 cm per degree C
     {
         m_snowdepth -= m_temp[ weatherday ];
         if (m_snowdepth < 0.0) m_snowdepth = 0.0;
     }
  
     if (m_temp[ weatherday ] < snowtempthreshold && m_rain[ weatherday ] > 0.0) {
         m_snowdepth += m_rain[ weatherday ];  // rain is in mm and snow in cm, so we get the conversion for free.
     }
     if (m_temp[ weatherday ] > snowtempthreshold && m_rain[ weatherday ] > 0.0) {
         m_snowdepth -= m_rain[ weatherday ];  // rain is in mm and snow in cm, so we get the conversion for free.
         if (m_snowdepth < 0.0) m_snowdepth = 0.0;
     }
     if (m_snowdepth > 0.0) {
         m_snowtoday = true;
     }
     else {
         m_snowtoday = false;
     }
  
     m_insolation = c_insolation[ g_date->DayInYear() ];
     m_temptoday = m_temp[ weatherday ];
     m_raintoday = m_rain[ weatherday ];
     m_windtoday = m_wind[ weatherday ];
  
     if (m_raintoday > 0.001)
         m_rainingtoday = true;
     else
         m_rainingtoday = false;
     double rainscale = 0;
     for (int i = 0; i < 5; i++) {
         rainscale += GetRain( (weatherday + m_datemodulus) - (i + 1) ) * pow( 0.5, i );
     }
     double temp = GetTempPeriod( weatherday + m_datemodulus - 1, 5 ) / 5.0;
     if (temp > 10.0) rainscale *= 0.5;
     if (temp > 20.0) rainscale *= 0.5;
     if (temp > 30.0) rainscale *= 0.5;
     m_humiditytoday = rainscale;
     double chance = g_rand_uni();
     if (chance < 0.1) m_winddirtoday = 2; // South
     else if (chance < 0.3) m_winddirtoday = 1; // East
     else if (chance < 0.5) m_winddirtoday = 0; // North
     else m_winddirtoday = 3; // West
 }

References c_insolation, Calendar::Date(), Calendar::DayInYear(), g_date, g_rand_uni, GetRain(), GetTempPeriod(), m_datemodulus, m_humiditytoday, m_insolation, m_rain, m_rainingtoday, m_raintoday, m_snowdepth, m_snowtoday, m_temp, m_temptoday, m_wind, m_winddirtoday, and m_windtoday.

Referenced by Landscape::Tick(), and Weather().

Member Data Documentation

◆ m_datemodulus

long Weather::m_datemodulus

private

Definition at line 404 of file weather.h.

Referenced by GetDDDegs(), GetSnow(), Tick(), and Weather().

◆ m_humiditytoday

double Weather::m_humiditytoday

private

Definition at line 399 of file weather.h.

Referenced by Tick().

◆ m_insolation

double Weather::m_insolation

private

Definition at line 398 of file weather.h.

Referenced by Tick().

◆ m_rain

vector<double> Weather::m_rain

private

Definition at line 387 of file weather.h.

Referenced by GetSnow(), Tick(), and Weather().

◆ m_rainingtoday

bool Weather::m_rainingtoday

private

Definition at line 397 of file weather.h.

Referenced by Tick().

◆ m_raintoday

double Weather::m_raintoday

private

Definition at line 393 of file weather.h.

Referenced by Tick().

◆ m_snowdepth

double Weather::m_snowdepth

private

The snow depth in cm.

Definition at line 407 of file weather.h.

Referenced by Tick(), and Weather().

◆ m_snowtoday

bool Weather::m_snowtoday

private

Definition at line 396 of file weather.h.

Referenced by GetSnow(), and Tick().

◆ m_temp

vector<double> Weather::m_temp

private

Definition at line 390 of file weather.h.

Referenced by GetDDDegs(), GetSnow(), Tick(), and Weather().

◆ m_temptoday

double Weather::m_temptoday

private

Definition at line 392 of file weather.h.

Referenced by Tick().

◆ m_wind

vector<double> Weather::m_wind

private

Definition at line 388 of file weather.h.

Referenced by Tick(), and Weather().

◆ m_wind_valid

bool Weather::m_wind_valid

private

Definition at line 401 of file weather.h.

◆ m_winddir

vector<double> Weather::m_winddir

private

Definition at line 389 of file weather.h.

◆ m_winddir_valid

bool Weather::m_winddir_valid

private

Definition at line 402 of file weather.h.

◆ m_winddirtoday

int Weather::m_winddirtoday

private

Definition at line 395 of file weather.h.

Referenced by Tick().

◆ m_windtoday

double Weather::m_windtoday

private

Definition at line 394 of file weather.h.

Referenced by Tick().

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

C:/Repo/ALMaSS_all/Landscape/weather.h
C:/Repo/ALMaSS_all/Landscape/weather.cpp

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Weather()

◆ ~Weather()

Member Function Documentation

◆ GetDDDegs()

◆ GetGlobalRadiation() [1/2]

◆ GetGlobalRadiation() [2/2]

◆ GetHumidity()

◆ GetMeanTemp()

◆ GetRain() [1/2]

◆ GetRain() [2/2]

◆ GetRainPeriod()

◆ GetSnow() [1/2]

◆ GetSnow() [2/2]

◆ GetSnowDepth()

◆ GetTemp() [1/2]

◆ GetTemp() [2/2]

◆ GetTempPeriod()

◆ GetWind() [1/2]

◆ GetWind() [2/2]

◆ GetWindDirection()

◆ GetWindPeriod()

◆ Raining()

◆ Tick()

Member Data Documentation

◆ m_datemodulus

◆ m_humiditytoday

◆ m_insolation

◆ m_rain

◆ m_rainingtoday

◆ m_raintoday

◆ m_snowdepth

◆ m_snowtoday

◆ m_temp

◆ m_temptoday

◆ m_wind

◆ m_wind_valid

◆ m_winddir

◆ m_winddir_valid

◆ m_winddirtoday

◆ m_windtoday