#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 ( const char * a_weatherdatafile )

Definition at line 167 of file weather.cpp.

References FloatToDouble(), 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.

 {
   FILE *ifile;
   int Magic, FileFormat, NoDays, Day, Month, Year, FirstYear=0, LastYear=0;
   float fTemp, fRain, fWind;
   double Temp, Rain, Wind;
   //char degree[20];
   //double DDegs = 0.0;
   ifile = fopen(a_weatherdatafile, "r" );
   if (!ifile) {
     g_msg->Warn(WARN_FILE, "Weather::Weather(): Unable to open file",
         a_weatherdatafile );
     exit(1);
   }
   fscanf( ifile, "%d", &Magic );
 
   if ( Magic != -1 ) {
     // Assume old file format.
     NoDays = Magic;
     FileFormat = 1;
   } else {
     fscanf( ifile, "%d %d", &FileFormat, &NoDays );
   }
 
   /*
   switch( FileFormat ) {
   case 1:
     break;
   default:
     g_msg->Warn(WARN_FILE, "Weather::Weather(): "
         "Unknown file format specified on line 1: ",
         a_weatherdatafile);
     g_msg->Warn(WARN_FILE, "First line must be: ",
         "-1 <file format specifier> <days worth of weather data>" );
     exit(1);
   }
   */
 
   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++) {
     // Read in the next data line.
     fscanf( ifile, "%d %d %d %g %g %g",
         &Year, &Month, &Day, &fTemp, &fWind, &fRain );
     FloatToDouble(Temp, fTemp);
     FloatToDouble(Rain, fRain);
     FloatToDouble(Wind, fWind);
     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++;
     }
 
     if ( Month == 12 && Day == 31 && storing ) {
       // Found yet another full year worth of weather data.
       g_date->SetLastYear( Year );
       LastYear = Year;
       //sprintf( degree, "%f", DDegs );
       //g_msg->Warn( WARN_FILE,
       //   "Sum of day degrees:",
       //     degree );
       //DDegs = 0.0;
     }
   }
 
   // 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();
 }

Weather::~Weather ( void )

Definition at line 274 of file weather.cpp.

275 {

276 }

Member Function Documentation

double Weather::GetDDDegs ( long a_date )

Definition at line 155 of file weather.cpp.

References m_datemodulus, and m_temp.

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

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

double Weather::GetGlobalRadiation ( long a_date )

inline

Definition at line 414 of file weather.h.

References c_insolation.

Referenced by Landscape::SupplyGlobalRadiation().

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

Weather::m_datemodulus

long m_datemodulus

Definition: weather.h:404

c_insolation

const double c_insolation[365]

Definition: weather.h:16

double Weather::GetGlobalRadiation ( void )

inline

Definition at line 415 of file weather.h.

415 { return m_insolation; }

Weather::m_insolation

double m_insolation

Definition: weather.h:398

double Weather::GetHumidity ( void )

inline

Get the humidity score today.

Definition at line 421 of file weather.h.

Referenced by Landscape::SupplyHumidity().

421 { return m_humiditytoday; }

Weather::m_humiditytoday

double m_humiditytoday

Definition: weather.h:399

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

Definition at line 120 of file weather.cpp.

References GetTemp().

Referenced by Landscape::SupplyMeanTemp().

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

double Weather::GetRain ( long a_date )

inline

Definition at line 423 of file weather.h.

Referenced by Landscape::SupplyRain().

423 { return m_rain[a_date %m_datemodulus]; }

Weather::m_rain

vector< double > m_rain

Definition: weather.h:387

Weather::m_datemodulus

long m_datemodulus

Definition: weather.h:404

double Weather::GetRain ( void )

inline

Definition at line 424 of file weather.h.

Referenced by GetRainPeriod(), and Tick().

424 { return m_raintoday; }

Weather::m_raintoday

double m_raintoday

Definition: weather.h:393

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

Definition at line 278 of file weather.cpp.

References GetRain().

Referenced by 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().

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

bool Weather::GetSnow ( long a_date )

Definition at line 132 of file weather.cpp.

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

Referenced by Landscape::SupplySnowcover().

 {
   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;
 }

bool Weather::GetSnow ( void )

inline

Definition at line 429 of file weather.h.

429 { return m_snowtoday; }

Weather::m_snowtoday

bool m_snowtoday

Definition: weather.h:396

double Weather::GetSnowDepth ( void )

inline

Get the current snow depth.

Definition at line 431 of file weather.h.

Referenced by Landscape::SupplySnowDepth().

431 { return m_snowdepth; }

Weather::m_snowdepth

double m_snowdepth

The snow depth in cm.

Definition: weather.h:407

double Weather::GetTemp ( long a_date )

inline

Get the temperature on a particular date.

Definition at line 417 of file weather.h.

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

417 { return m_temp[a_date % m_datemodulus]; }

Weather::m_temp

vector< double > m_temp

Definition: weather.h:390

Weather::m_datemodulus

long m_datemodulus

Definition: weather.h:404

double Weather::GetTemp ( void )

inline

Get the temperature today.

Definition at line 419 of file weather.h.

Referenced by GetMeanTemp(), and GetTempPeriod().

419 { return m_temptoday; }

Weather::m_temptoday

double m_temptoday

Definition: weather.h:392

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 297 of file weather.cpp.

References GetTemp().

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

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

double Weather::GetWind ( long a_date )

inline

Definition at line 425 of file weather.h.

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

425 { return m_wind[a_date % m_datemodulus]; }

Weather::m_datemodulus

long m_datemodulus

Definition: weather.h:404

Weather::m_wind

vector< double > m_wind

Definition: weather.h:388

double Weather::GetWind ( void )

inline

Definition at line 426 of file weather.h.

Referenced by GetWindPeriod().

426 { return m_windtoday; }

Weather::m_windtoday

double m_windtoday

Definition: weather.h:394

int Weather::GetWindDirection ( void )

inline

Definition at line 427 of file weather.h.

Referenced by Landscape::SupplyWindDirection().

427 { return m_winddirtoday; }

Weather::m_winddirtoday

int m_winddirtoday

Definition: weather.h:395

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

Definition at line 287 of file weather.cpp.

References GetWind().

Referenced by Landscape::SupplyWindPeriod().

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

bool Weather::Raining ( void )

inline

Definition at line 432 of file weather.h.

Referenced by Farm::AutumnHarrow(), Farm::AutumnPlough(), Farm::AutumnRoll(), Farm::AutumnSow(), Farm::CutOrch(), Farm::CutToSilage(), Farm::CutWeeds(), Farm::DeepPlough(), Farm::FA_AmmoniumSulphate(), Farm::FA_GreenManure(), Farm::FA_Manure(), Farm::FA_NPK(), Farm::FA_PK(), Farm::FA_Sludge(), Farm::FA_Slurry(), Farm::FP_GreenManure(), Farm::FP_LiquidNH3(), Farm::FP_ManganeseSulphate(), Farm::FP_Manure(), Farm::FP_NPK(), Farm::FP_NPKS(), Farm::FP_PK(), Farm::FP_Sludge(), Farm::FP_Slurry(), Farm::FungicideTreat(), OptimisingFarm::FungicideTreat(), Farm::GrowthRegulator(), Farm::HerbicideTreat(), OptimisingFarm::HerbicideTreat(), Farm::HillingUp(), Farm::InsecticideTreat(), OptimisingFarm::InsecticideTreat(), Farm::Molluscicide(), Farm::ProductApplication0(), Farm::ProductApplication1(), Farm::SpringHarrow(), Farm::SpringPlough(), Farm::SpringRoll(), Farm::SpringSow(), Farm::Swathing(), Farm::SynInsecticideTreat(), Farm::Water(), and Farm::WinterPlough().

432 { return m_rainingtoday; }

Weather::m_rainingtoday

bool m_rainingtoday

Definition: weather.h:397

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.

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().

                          {
     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
 }

Member Data Documentation

long Weather::m_datemodulus

private

Definition at line 404 of file weather.h.

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

double Weather::m_humiditytoday

private

Definition at line 399 of file weather.h.

Referenced by Tick().

double Weather::m_insolation

private

Definition at line 398 of file weather.h.

Referenced by Tick().

vector<double> Weather::m_rain

private

Definition at line 387 of file weather.h.

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

bool Weather::m_rainingtoday

private

Definition at line 397 of file weather.h.

Referenced by Tick().

double Weather::m_raintoday

private

Definition at line 393 of file weather.h.

Referenced by Tick().

double Weather::m_snowdepth

private

The snow depth in cm.

Definition at line 407 of file weather.h.

Referenced by Tick(), and Weather().

bool Weather::m_snowtoday

private

Definition at line 396 of file weather.h.

Referenced by GetSnow(), and Tick().

vector<double> Weather::m_temp

private

Definition at line 390 of file weather.h.

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

double Weather::m_temptoday

private

Definition at line 392 of file weather.h.

Referenced by Tick().

vector<double> Weather::m_wind

private

Definition at line 388 of file weather.h.

Referenced by Tick(), and Weather().

bool Weather::m_wind_valid

private

Definition at line 401 of file weather.h.

vector<double> Weather::m_winddir

private

Definition at line 389 of file weather.h.

bool Weather::m_winddir_valid

private

Definition at line 402 of file weather.h.

int Weather::m_winddirtoday

private

Definition at line 395 of file weather.h.

Referenced by Tick().

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:

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation