main()

int main	(	int	iargc,
		char *	argv[]
	)

Definition at line 50 of file gsexx_write_example.cc.

References examplefile, GSE2::waveform::TWID2::Fauxid, GSE2::waveform::TWID2::Fcalib, GSE2::waveform::TWID2::Fcalper, GSE2::waveform::TWID2::Fchannel, GSE2::waveform::TWID2::Fday, GSE2::waveform::TWID2::Fhang, GSE2::waveform::TWID2::Fhour, GSE2::waveform::TWID2::Finstype, GSE2::waveform::TWID2::Fminute, GSE2::waveform::TWID2::Fmonth, GSE2::waveform::TWID2::Fsamprate, GSE2::waveform::TWID2::Fsamps, GSE2::waveform::TWID2::Fstation, GSE2::waveform::TWID2::Fsubformat, GSE2::waveform::TWID2::Fvang, GSE2::waveform::TWID2::Fyear, GSEXX_WRITE_EXAMPLE_CVSID, GSEXX_WRITE_EXAMPLE_VERSION, GSE2::waveform::TDAT2sum::hot(), GSE2::waveform::TWID2::line(), and GSE2::waveform::SF_CM6.

{
  cout << GSEXX_WRITE_EXAMPLE_VERSION << endl;
  cout << GSEXX_WRITE_EXAMPLE_CVSID << endl;
  cout << endl;
  cout << "This program provides an example of basic writing to" << endl
       << "GSE data files. The program will write two traces" << endl
       << "of CM6 encoded data to a file called '" 
       << examplefile << "'." << endl;
  cout << endl;

  // open output file
  std::ofstream os(examplefile); 

  // write first time series
  {
    // number of samples to be written
    const int nsamples=2000; 

    // prepare WID2 line
    GSE2::waveform::TWID2 wid2line; 
    wid2line.Fyear=2006; 
    wid2line.Fmonth=1; 
    wid2line.Fday=25; 
    wid2line.Fhour=10; 
    wid2line.Fminute=30; 
    wid2line.Fseconds=12.456789; 
    wid2line.Fsamps=nsamples; 
    wid2line.Fsamprate=20.; 
    wid2line.Fstation="BFO"; 
    wid2line.Finstype="STS2"; 
    wid2line.Fchannel="BHZ"; 
    wid2line.Fauxid="any"; 
    wid2line.Fvang=0.; 
    wid2line.Fhang=90.; 
    wid2line.Fcalib=1.e5; 
    wid2line.Fcalper=10.; 
    wid2line.Fsubformat=GSE2::waveform::SF_CM6; 
      
    // report WID2 line
    std::cout << wid2line.line();

    // write WID2 line to data file
    os << wid2line.line();

    // create CM6 encoding data writer for appropriate number of samples
    GSE2::waveform::TDAT2writeCM6 writer(nsamples);

    // write data
    int i=0;
    while (writer.hot())
    {
      // calculate synthetic sample
      int sample=int(1.e5*std::sin(i*2.*3.141592653*5./(nsamples-35)));
      // output sample to data file
      os << writer(sample);
      i++;
    }
  }
    
  // write second time series
  {
    // number of samples to be written
    const int nsamples=2300; 

    // prepare WID2 line with the bare minimum of information required
    GSE2::waveform::TWID2 wid2line; 
    wid2line.Fsamps=nsamples; 
    wid2line.Fsamprate=1.; 
    wid2line.Fsubformat=GSE2::waveform::SF_CM6; 

    // report WID2 line
    std::cout << wid2line.line();

    // write WID2 line to data file
    os << wid2line.line();

    // create CM6 encoding data writer for appropriate number of samples
    GSE2::waveform::TDAT2writeCM6 writer(nsamples);

    // write data
    int i=0;
    while (writer.hot())
    {
      // calculate synthetic sample
      int sample=int(1.2e5*(std::sin(i*2.*3.141592653*5./(nsamples-35))
                          +std::cos(i*2.*3.141592653*5.5/(nsamples-35))));
      // output sample to data file
      os << writer(sample);
      i++;
    }
  }

  // that's it!
  cout << "data file is written..." << endl;
  
} // main

Here is the call graph for this function: