ReadPeriodList.C

Go to the documentation of this file.

/*
# Copyright (C) 2019 Gabriele Vedovato
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see <https://www.gnu.org/licenses/>.
*/
 
 
#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;
}