offset.cc

Go to the documentation of this file.

#define TF_OFFSET_CC_VERSION \
  "TF_OFFSET_CC   V1.2"

#include <cmath>
#include <sffxx.h>

namespace sff {

  double offset(const SRCE& srce, const INFO& info,
                const double& radius)
  {
    double retval;
    SFF_assert((srce.cs == info.cs),
                "ERROR (offset): inconsistent coordinate systems!");
    if (srce.cs == CS_spherical)
    {
      const double convfac=0.017453293;
      double rtet=(90.-info.cx)*convfac;
      double stet=(90.-srce.cx)*convfac;
      double rphi=info.cy*convfac;
      double sphi=srce.cy*convfac;
      double cosepi=cos(stet)*cos(rtet)+sin(stet)*sin(rtet)*cos(sphi-rphi);
      retval=acos(cosepi)*radius*1.e3;
    }
    else if (srce.cs == CS_cartesian)
    {
      double dx=srce.cx-info.cx;
      double dy=srce.cy-info.cy;
      retval=sqrt(dx*dx+dy*dy);
    }
    else
    {
      SFF_abort("ERROR (offset): unknown coordinate system!");
    }
    return(retval);
  } // double offset(const SRCE& srce, const INFO& info;
    //               const double& radius=6371.)

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

  double sourcedistance(const SRCE& srce, const INFO& info)
  {
    double retval;
    SFF_assert((srce.cs == info.cs),
                "ERROR (sourcedistance): inconsistent coordinate systems!");
    SFF_assert((srce.cs == CS_cartesian),
                "ERROR (sourcedistance): "
                "only provided in Cartesian coordinates!");
    if (srce.cs == CS_cartesian)
    {
      double dx=srce.cx-info.cx;
      double dy=srce.cy-info.cy;
      double dz=srce.cz-info.cz;
      retval=sqrt(dx*dx+dy*dy+dz*dz);
    }
    else
    {
      SFF_abort("ERROR (sourcedistance): unknown coordinate system!");
    }
    return(retval);
  } // double sourcedistance(const SRCE& srce, const INFO& info)

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

  double offsetdeg(const SRCE& srce, const INFO& info, 
                   const double& radius)
  {
    double retval=offset(srce, info, radius)*1.e-3*360./(2.*M_PI*radius);
    return(retval);
  } // double offsetdeg(const SRCE& srce, const INFO& info,
    //                  const double& radius=6371.)

} // namespace sff

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

/* 
 * Fortran function
 *
      character*1 cs,rs
      real c1,c2,c3,r1,r2,r3
cE
      real distance
      real rtet, stet, rphi, sphi
      real convfac, earthradius, cosepi
      parameter(convfac=0.017453293, earthradius=6371.)
c
      if (cs.ne.rs) stop 'ERROR (sffu_offset): different reference frames'
      if (cs.eq.'C') then
        distance=sqrt((c1-r1)**2+(c2-r2)**2+(c3-r3)**2)
      elseif (cs.eq.'S') then
c        print *,'source ',c1,c2,c3
c        print *,'receiver ',r1,r2,r3
        rtet=(90.-r1)*convfac 
        stet=(90.-c1)*convfac 
        rphi=r2*convfac
        sphi=c2*convfac
c        print *,'source ',stet,sphi
c        print *,'receiver ',rtet,rphi
        cosepi=cos(stet)*cos(rtet)+sin(stet)*sin(rtet)*cos(sphi-rphi)
        distance=acos(cosepi)*earthradius*1.e3
c        print *,'distance ',distance
      else
        stop 
     &'ERROR (sffu_offset): reference frame is neither cartesian nor spherical'
      endif
      sffu_offset=distance
      return
      end
cS
*/

/* ----- END OF offset.cc ----- */