ReadPeriodList.C

Go to the documentation of this file.

#define PERIOD_MAX_SIZE 10

/*
  TString run[PERIOD_MAX_SIZE];
  double  gps_start[PERIOD_MAX_SIZE];
  TString date_start[PERIOD_MAX_SIZE];
  double  gps_stop[PERIOD_MAX_SIZE];
  TString date_stop[PERIOD_MAX_SIZE];
  double  interval_day[PERIOD_MAX_SIZE];
  int     obs_time_sec[PERIOD_MAX_SIZE];
  double  obs_time_day[PERIOD_MAX_SIZE];

  int nPeriod = ReadPeriodList(ifile, run, gps_start, date_start, gps_stop, date_stop,
                                  interval_day, obs_time_sec, obs_time_day); 
  
*/

int ReadPeriodList(TString ifile, TString* run, double* gps_start, TString* date_start, double* gps_stop, TString* date_stop,
                                  double* interval_day, int* obs_time_sec, double* obs_time_day) {

  CWB::Toolbox::checkFile(ifile);

  // Open chunk list
  ifstream in;
  in.open(ifile.Data(),ios::in);
  if (!in.good()) {cout << "Error Opening File : " << ifile << endl;exit(1);}

  int isize=0;
  char str[1024];
  int fpos=0;
  while(true) {
    in.getline(str,1024);
    if (!in.good()) break;
    if(str[0] != '#') isize++;
  }
  cout << "size " << isize << endl;
  in.clear(ios::goodbit);
  in.seekg(0, ios::beg);
  if(isize==0)  {cout << "Error : File " << ifile << " is empty" << endl;exit(1);}
  if(isize>PERIOD_MAX_SIZE) {cout << "Error : File " << ifile << " > " << PERIOD_MAX_SIZE << endl;exit(1);}
 
  char srun[1024];
  char sdate_start[1024];
  char sdate_stop[1024];

  int k=0;
  while(true) {
    fpos=in.tellg();
    in.getline(str,1024);
    if (!in.good()) break;
    if(str[0] == '#' || str[0]=='\0') continue;
    in.seekg(fpos, ios::beg);
    in >> srun >> gps_start[k] >> sdate_start >> gps_stop[k] >> sdate_stop >> interval_day[k] >> obs_time_sec[k] >> obs_time_day[k];
    if(!in.good()) break;
    run[k]=srun;
    date_start[k]=sdate_start;
    date_stop[k]=sdate_stop;
    cout << "\t" << srun << "\t" << gps_start[k] << "\t" << sdate_start << "\t" << gps_stop[k] 
         << "\t" << sdate_stop << "\t" << obs_time_sec[k] << "\t" << obs_time_day[k] << endl;
    k++;
    if(k>=PERIOD_MAX_SIZE) {
      cout << "WARNING: period list exceed PERIOD_MAX_SIZE = " << PERIOD_MAX_SIZE << endl; 
      break;
    }
  }
  in.close();

  return k;
}