main.cpp

/*
*  main.cpp
*  ReGen 1.0.1
*
*  Created by Martin Köchy on Thu Feb 06 2003.
*
*/

#include <stdlib.h>
#include <fstream>
#include <iostream>
#include <sstream>
#include <stack>
#include "ReGen.h"


int main (int argc, char * const argv[]) {
  
  /*********** declarations ************************/
  
  int nP = 0; // number of supplied parameters
  int y; // length of time series, years
  
  ReGen* pTS; // pointer to data base of time series
  
  char answer = '?';
  
  std::string RFileName("ReGen time series.txt");
  std::string PList; 
  
  std::deque <float> manualP;
  RAINPARAMETERS* pP;
  
  do
  {
        /*********** set parameters ************************/
        if (argc==1)
        { 
          // ask for values
          std::cout << "This program generates stochastic time series of precipitation." << std::endl;
          std::cout << "It first determines the seasonal probability that a day is a rainy day." << std::endl;
          std::cout << "Then it determines the volume for that day using a negative exponential" << std::endl;
          std::cout << "distribution whose mean varies seasonally." << std::endl;
          std::cout << "Both seasonal functions have a Gaussian peak shape with three parameters." << std::endl;
          std::cout << "The location and width should be in month units starting with August = 0." << std::endl;
          std::cout << "These parameters can be supplied directly or the program can calculate" << std::endl;
          std::cout << "them if the mean annual precipitation is supplied." << std::endl;
          std::cout << "The parameters can be entered interactively after the start of the" << std::endl;
          std::cout << "program or added as a space-delimited list after the program" << std::endl;
          std::cout << "name when it is launched from the command line.\n" << std::endl;
          std::cout << "Parameters must be supplied in the following order:" << std::endl;
          std::cout << "1. Number of years (365 d) in the time series" << std::endl;
          std::cout << "2. rainy day probability" << std::endl;
          std::cout << "2.1. amplitude" << std::endl;
          std::cout << "2.2. location of maximum" << std::endl;
          std::cout << "2.3. width" << std::endl;
          
          std::cout << "3. rainy day volume" << std::endl;
          std::cout << "3.1. amplitude" << std::endl;
          std::cout << "3.2. location of maximum" << std::endl;
          std::cout << "3.3. width" << std::endl;
          std::cout << "4. optionally: a factor [-0.3,+0.3] to change the amplitude\n while maintaining the annual volume\n" << std::endl;
          std::cout << "Alternatively, the number of years, mean annual precipitation [90, 900], \nand the optional factor can be supplied.\n\n" << std::endl;
          
          float parameter = 0;
          char input [60];
          std::cout << "Enter parameters:" << std::endl;
          std::cin.getline(input,60);
          istringstream sPList(input);
          while (sPList >> parameter) // input funktioniert nicht!
                manualP.push_back(parameter);
          nP=manualP.size()-1;
          
          /*    std::cout << "For further details about how the time series are generated and how" << std::endl;
          std::cout << "internal parameters have been derived refer to the documentation" << std::endl;
          std::cout << "or http://martinkoechy.de/research/software/ ." << std::endl;
          */
        }
        else // if no interactive input read from command line arguments
        {
          for(int s=1; s<argc; s++)
                manualP.push_back(atof(argv[s]));
          nP = argc-1;
        }
        
        y = int(manualP[0]);
        
        if (nP==1)
        {
          // calculate parameters, standard scenario
          pP = new RAINPARAMETERS (int(manualP[1]));
        }
        else if (nP==2)
        {
          // calculate parameters, apply change to daily pattern
          pP = new RAINPARAMETERS (int(manualP[1]), manualP[2]);
        }
        else if (nP == 6)
        {
          // parameters provided, standard scenario
          pP = new RAINPARAMETERS (manualP[1], manualP[2], manualP[3], manualP[4], manualP[5], manualP[6]);
        }
        else if (nP == 7)
        {
          // parameters provided, apply change to daily pattern
          pP = new RAINPARAMETERS (manualP[1], manualP[2], manualP[3], manualP[4], manualP[5], manualP[6], manualP[7]);
        }
        else
        {
          if (argc < 2)
          { 
                // error message
                std::cerr << "Error: You supplied the wrong number of arguments." << std::endl;
                std::cerr << "Either supply the mean annual precipitation and optionally the" << std::endl;
                std::cerr << "fractional change to the amplitdue of daily rain volume," << std::endl;
                std::cerr << "or all six parameters and optionally the fractional change.\n" << std::endl;
                do
                {
                  std::cerr << "Quit [Q] or enter manually [M]?" << std::endl;
                  std::cin >> answer;
                }
                while (!(answer == 'Q' || answer == 'M' || answer == 'q' || answer == 'm'));
                if (answer=='q' || answer=='Q') exit(argc);
          }
        }
  } while (answer=='m');
  
  
  /*********** calculations ************************/
  
  pTS = new ReGen(pP);
  
  /*********** write output ************************/
  
  std::ofstream RainDocu;
  RainDocu.open(RFileName.c_str(), std::ios::app);
  RainDocu << "yr\tday\train" << std::endl;
  
  for (int j = 0; j<y; j++)
  {
        pTS->drawYear();
        for (int t=0; t<days_in_year; t++)
        {
          RainDocu << j << "\t" << t << "\t" << pTS->getRain(t) << std::endl;
        }
  }
  
  RainDocu.close();
  delete(pTS);
  delete (pP);
  return 0;
}

---