fapid_sff_rtracefi.cc

Go to the documentation of this file.

#define TF_FAPID_SFF_RTRACEFI_CC_VERSION \
  "TF_FAPID_SFF_RTRACEFI_CC   V1.0   "


#include <fapidxx/fapidsff.h>
#include <fapidxx/fileunit.h>
#include <fapidxx/helper.h>
#include <fapidxx/error.h>

using namespace fapidxx;

int sff_rtracefi__(integer *lu, real *tanf, real *dt, char *wid2line, 
                   integer *nsamp, real *fdata, integer *idata, char *code,
                   logical *last, integer *nline, char *lines,
                   integer *lindim, integer *lenmax, char *cs, real *c1,
                   real *c2, real *c3, integer *nstack, integer *ierr,
                   ftnlen wid2line_len, ftnlen code_len,
                   ftnlen lines_len, ftnlen cs_len)
{
  int retval=0;
  *ierr=0;
  try {
    datrw::ianystream &is=istreammanager(static_cast<int>(*lu));
    sff::WID2 wid2;
    sff::INFO info;
    datrw::Tfseries iseries;
    is >> iseries;
    FAPIDXX_fuassert((static_cast<int>(iseries.size())<=(*nsamp)), *lu,
                     "sff_rtracefi__: too many samples");
    int nsamples=iseries.size();
    aff::LinearShape shape(0, nsamples-1, 0);
    datrw::Tfseries series(shape, aff::SharedHeap<real>(fdata, *nsamp));
    series.copyin(iseries);
    *last = is.last() ? 1 : 0;
    is >> wid2;
    std::string ocode("");
    if (is.hasfree()) { ocode.append("F"); }
    if (is.hasinfo()) { ocode.append("I"); }
    if (!is.last()) { ocode.append("D"); }
    if (is.hasinfo()) { is >> info; }
    fillfstring(ocode, code, code_len);
    *tanf=static_cast<real>(maketanf(wid2.date));
    *dt=static_cast<real>(wid2.dt);
    fillfstring(wid2.line(), wid2line, wid2line_len);
    *nsamp=series.size();
    char thecs=sff::coosysID(info.cs);
    fillfstring(std::string(&thecs, 1), cs, cs_len);
    *c1=static_cast<real>(info.cx);
    *c2=static_cast<real>(info.cy);
    *c3=static_cast<real>(info.cz);
    *nstack=info.nstacks;
    sff::FREE free;
    if (is.hasfree()) { is >> free; }
    freeblock(free, nline, lines, lindim, lenmax, lines_len);
  }
  catch(...) {
    *ierr=1;
  }
  return retval;
} // int sff_rtracefi__

/* ----- END OF fapid_sff_rtracefi.cc ----- */