ConvertChunkList.C

Go to the documentation of this file.

#define CHUNK_MAX_SIZE     100
 
// this macro converts Chunk_List.lst format 
// O3    1  1238166018  1238787954
// to Chunk_List.txt format
// O3    1  1238166018  1238787954 7.16   days  Apr 01 15:00 - Apr 08 19:45
 
// input  file = Chunk_List.lst
// output file = Chunk_List.txt
 
int ConvertChunkList(TString ifile) {
 
  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++;
  }
  in.clear(ios::goodbit);
  in.seekg(0, ios::beg);
  if(isize==0)  {cout << "Error : File " << ifile << " is empty" << endl;exit(1);}
  if(isize>CHUNK_MAX_SIZE) {cout << "Error : File " << ifile << " > " << CHUNK_MAX_SIZE << endl;exit(1);}
 
  char   run[CHUNK_MAX_SIZE][256];
  int    chunk[CHUNK_MAX_SIZE];
  double start[CHUNK_MAX_SIZE];
  double stop[CHUNK_MAX_SIZE];
 
  cout.precision(10); 
 
  TString ofile = ifile;
  ofile.ReplaceAll(".lst",".txt");
 
  ofstream out;
  out.open(ofile.Data());
  if (!out.good()) {cout << "Error Opening File : " << ofile << endl;exit(1);}
 
  int k=0;
  while(true) {
    // IN  -> O3     1  1238166018  1238787954
    in >> run[k] >> chunk[k] >> start[k] >> stop[k];
    if(!in.good()) break;
    if(run[k][0]=='#') continue;
//    cout << run[k] << "\t" << chunk[k] << "\t" << start[k] << "\t" << stop[k] << endl; 
    double days = (stop[k]-start[k])/(24.*3600.);
    wat::Time sdate(start[k]);
    wat::Time edate(stop[k]);
 
    TString month[12] = {"Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec"};
 
    char sout[1024];
    // OUT -> O3     1  1238166018  1238787954 7.16   days  Apr 01 15:00 - Apr 08 19:45
    sprintf(sout,"%s    %d  %d  %d %.2f  days    %s %02d %02d:%02d - %s %02d %02d:%02d",run[k],chunk[k],int(start[k]),int(stop[k]),days,
            month[sdate.GetMonth()-1].Data(),sdate.GetDay(),sdate.GetHour(),sdate.GetMinute(),
            month[edate.GetMonth()-1].Data(),edate.GetDay(),edate.GetHour(),edate.GetMinute());
 
    cout << sout << endl;
    out << sout << endl;
 
    k++;
  }
  in.close();
  out.close();
 
  exit(0);
}