widXio.cc

Go to the documentation of this file.

#define TF_WIDXIO_CC_VERSION \
  "TF_WIDXIO_CC   V1.1"

#include <sffxx.h>
#include <sstream>

namespace sff {

  namespace helper {

    int nsignificantdigits(double v, const bool& debug=false)
    {
      const double epsilon=1.e-30;
      if (debug) 
      { 
        std::cerr << "DEBUG (nsignificantdigits): "
          << "entered function for value:" << v << std::endl; 
      }
      if (v<0) { v = -v; }
      int n=0;
      if (v<epsilon) 
      {
        n=2;
      }
      else
      {
        double basefactor=pow(10.,-floor(log10(v)));
        double v1,v2;
        do
        { 
          double factor=basefactor*pow(10.,n);
          v1=v*factor;
          v2=floor(v*factor);
          n++;
          if (debug) 
          { 
            std::cerr << "DEBUG (nsignificantdigits): "
              << "basefactor: " << basefactor << " "
              << "factor: " << factor << " "
              << "n: " << n << " "
              << "v1: " << v1 << " "
              << "v2: " << v2 << " "
              << std::endl; 
          }
        }
        while (v1 != v2);
      }
      return (n);
    } // int nsignificantdigits(const double& v)

  } // namespace helper

  /*----------------------------------------------------------------------*/

  std::string WIDXline(const sff::WID2& wid2, const bool& debug)
  {
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): entered function" <<
                                      std::endl; 
    }
    std::ostringstream os;
    os << sff::WIDXID;
    os << " " << wid2.date.hierarchicalstring();
    os << " ";
    os.width(5);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.station.substr(0,5);
    os << " ";
    os.width(3);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.channel.substr(0,3);
    os << " ";
    os.width(4);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.auxid.substr(0,4);
    os << " ";
    os.width(6);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.instype.substr(0,6);
    os << " " << wid2.nsamples;
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): partial result:" <<
                                      os.str() << std::endl; 
    }
    double rate=1./wid2.dt;
    double value=wid2.dt;
    if (sff::helper::nsignificantdigits(rate) < 4)
    {
      os << " r ";
      value=rate;
    }
    else
    {
      os << " i ";
    }
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): partial result:" <<
                                      os.str() << std::endl; 
    }
    os.setf(std::ios_base::scientific, std::ios_base::floatfield);
    int ndigits=sff::helper::nsignificantdigits(value, debug);
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): ndigits:" <<
                                      ndigits << std::endl; 
    }
    os.precision(ndigits);
    os.width(ndigits+6);
    os << value;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.calib, debug));
    os << wid2.calib;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.calper, debug));
    os << wid2.calper;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.hang, debug));
    os << wid2.hang;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.vang, debug));
    os << wid2.vang;
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): final result:" <<
                                      os.str() << std::endl; 
    }
    return(os.str());
  } // std::string WIDXline(const sff::WID2& wid2)

  /*----------------------------------------------------------------------*/

  sff::WID2 WIDXline(const std::string& line)
  {
    sff::WID2 retval;
    SFF_assert(line.substr(0,4) == WIDXID, 
               "wrong line ID");
    retval.date=libtime::TAbsoluteTime(line.substr(5,26));
    retval.station=line.substr(32,5);
    retval.channel=line.substr(38,3);
    retval.auxid=line.substr(42,4);
    retval.instype=line.substr(47,6);
    std::istringstream is(line.substr(54));
    is >> retval.nsamples;
    std::string type;
    is >> type;
    is >> retval.dt;
    if (type == "r") { retval.dt = 1./retval.dt; }
    is >> retval.calib;
    is >> retval.calper;
    is >> retval.hang;
    is >> retval.vang;
    return(retval);
  } // sff::WID2 WIDXline(const std::string& line)

} // namespace sff

/* ----- END OF widXio.cc ----- */