test_dat2cm6()

void GSE2::tests::test_dat2cm6

(

)

Definition at line 163 of file gsexx_tests.cc.

References GSE2::waveform::TWID2::Fauxid, GSE2::waveform::TWID2::Fchannel, GSE2::waveform::TWID2::Fday, GSE2::waveform::TWID2::Finstype, GSE2::waveform::TWID2::Fmonth, GSE2::waveform::TWID2::Fsamprate, GSE2::waveform::TWID2::Fsamps, GSE2::waveform::TWID2::Fstation, GSE2::waveform::TWID2::Fyear, GSE2::waveform::TDAT2sum::hot(), GSE2::waveform::TWID2::line(), and GSE2::waveform::TWID2::read().

Referenced by main().

{
  std::cout << std::endl;
  std::cout << "Test the DAT2 class with CM6 subformat."
    << std::endl;
  std::cout << "1. write waveform"
    << std::endl;

  const int msamples=2000;
  int data[msamples];
  int indata[msamples];
  for (int i=0; i<msamples; i++)
  {
    data[i]=int(1.e5*std::sin(i*2.*3.141592653*5./(msamples-35)));
    indata[i]=0;
  }

  GSE2::waveform::TWID2 wid2line;
  wid2line.Fsamps=msamples;
  wid2line.Fsamprate=1.;
  wid2line.Fstation="BFO";
  wid2line.Finstype="STS-2";
  wid2line.Fchannel="BHZ";
  std::cout << wid2line.line();

  GSE2::waveform::TDAT2writeCM6 writer(msamples);
  int i=0;
  while (writer.hot())
  {
    if (i >= msamples)
    {
      std::cerr << "ERROR: missed last sample!" << std::endl;
      abort();
    }
    std::cout << writer(data[i]);
    i++;
  }
 
  // write to file
  {
    std::ofstream os("junk.dat");
    os << wid2line.line();
    GSE2::waveform::TDAT2writeCM6 fwriter(msamples);
    int i=0;
    while (fwriter.hot())
    {
      if (i >= msamples)
      {
        std::cerr << "ERROR: missed last sample!" << std::endl;
        abort();
      }
      os << fwriter(data[i]);
      i++;
    }
  }

  {
    char c;
    std::cerr << "manipulate junk.dat - if you like to..." << std::endl;
    std::cin.get(c);
  }

  GSE2::waveform::TWID2 newwid2line;
  newwid2line.Fstation="STA";
  newwid2line.Fchannel="CHA";
  newwid2line.Fauxid="AUX";
  newwid2line.Fyear=1;
  newwid2line.Fmonth=2;
  newwid2line.Fday=3;

  // read from file
  {
    std::ifstream is("junk.dat");
    GSE2::waveform::TDAT2readCM6 freader(msamples);
    newwid2line.read(is);
    int i=0;
    while (freader.hot())
    {
      if (i >= msamples)
      {
        std::cerr << "ERROR: missed last sample!" << std::endl;
        abort();
      }
      indata[i] =freader(is);
      i++;
    }
  }
 
  // write to second file
  {
    std::ofstream os("junk2.dat");
    os << newwid2line.line();
    GSE2::waveform::TDAT2writeCM6 fwriter(msamples);
    int i=0;
    while (fwriter.hot())
    {
      if (i >= msamples)
      {
        std::cerr << "ERROR: missed last sample!" << std::endl;
        abort();
      }
      os << fwriter(indata[i]);
      i++;
    }
  }
      
} // test_dat2cm6

Here is the call graph for this function:

Here is the caller graph for this function: