readheader()

void datrw::isacstream::readheader ( )

private

Definition at line 102 of file isacstream.cc.

References datrw::sac::SACheader::b, datrw::sac::SACheader::cmpaz, datrw::sac::SACheader::cmpinc, DATRW_assert, DATRW_debug, datrw::sac::SACheader::delta, datrw::sac::SACheader::e, datrw::idatstream::free(), datrw::sac::SACheader::iftype, datrw::idatstream::info(), datrw::sac::SACheader::kcmpnm, datrw::sac::SACheader::khole, datrw::sac::SACheader::knetwk, datrw::sac::SACheader::kstnm, datrw::sac::SACheader::leven, datrw::idatstream::Mdebug, datrw::idatstream::Mis, datrw::idatstream::newtrace(), datrw::sac::SACheader::npts, datrw::sac::SACheader::nzhour, datrw::sac::SACheader::nzjday, datrw::sac::SACheader::nzmin, datrw::sac::SACheader::nzmsec, datrw::sac::SACheader::nzsec, datrw::sac::SACheader::nzyear, datrw::sac::read_sac_header(), datrw::idatstream::setinfo(), datrw::idatstream::setlast(), datrw::idatstream::settracefree(), datrw::idatstream::setwid2(), datrw::sac::SACheader::stdp, datrw::sac::SACheader::stel, datrw::sac::SACheader::stla, datrw::sac::SACheader::stlo, and datrw::idatstream::wid2().

Referenced by fseries(), and skipseries().

  {
    DATRW_debug(this->Mdebug, "isacstream::readheader()",
                  "start reading header");
    this->newtrace();
    sff::FREE free;
    sff::WID2 wid2;
    sff::INFO info;
    std::ostringstream oss;
    DATRW_debug(this->Mdebug, "isacstream::readheader()",
                  "calling datrw::sac::read_sac_header");
    datrw::sac::SACheader hd=datrw::sac::read_sac_header(this->Mis);
    DATRW_debug(this->Mdebug, "isacstream::readheader()",
                  "returned from datrw::sac::read_sac_header");

    oss << "location: ";
    if (hd.stla>0) { oss << hd.stla << "°N; "; }
    else { oss << -hd.stla << "°S; "; }
    if (hd.stlo>0) { oss << hd.stlo << "°E; "; }
    else { oss << -hd.stlo << "°W"; }
    free.append(oss.str());
    DATRW_debug(this->Mdebug, "isacstream::readheader()", oss.str());
    free.append("values from SAC binary header " 
                "that do not fit into SFF trace header");
    wid2.dt=hd.delta;
    oss.str("");
    oss << "initial index: " << hd.b;
    oss << "; final index: " << hd.e;
    free.append(oss.str());
    info.cs=sff::CS_spherical;
    info.cx=hd.stla;
    info.cy=hd.stlo;
    info.cz=hd.stel;
    oss.str("");
    oss << hd.stdp;
    free.append("station depth: " + oss.str() + " m");
    wid2.hang=hd.cmpaz;
    wid2.vang=hd.cmpinc;
    wid2.date=libtime::TAbsoluteTime(hd.nzyear, 1, 1, hd.nzhour,
                                     hd.nzmin, hd.nzsec, hd.nzmsec);
    wid2.date.setdoy(hd.nzjday);
    wid2.nsamples=hd.npts;
    wid2.station=std::string(hd.kstnm,8);
    wid2.channel=std::string(hd.kcmpnm,8);
    wid2.auxid=std::string(hd.khole,8);
    free.append("network name: " + std::string(hd.knetwk, 8));
    DATRW_assert(hd.iftype==1, "unexpected file type");
    DATRW_assert(hd.leven==1, "uneven sampling");
    this->setwid2(wid2);
    this->setinfo(info);
    this->settracefree(free);
    this->setlast();

    DATRW_debug(this->Mdebug, "isacstream::readheader()",
                  "finished; return from function");
  }

Here is the call graph for this function:

Here is the caller graph for this function: