offset()

double sff::offset	(	const SRCE &	srce,
		const INFO &	info,
		const double &	radius
	)

return offset in meters

Definition at line 45 of file offset.cc.

References sff::SRCE::cs, sff::INFO::cs, CS_cartesian, CS_spherical, sff::SRCE::cx, sff::INFO::cx, sff::SRCE::cy, sff::INFO::cy, SFF_abort, and SFF_assert.

Referenced by offsetdeg().

  {
    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;

Here is the caller graph for this function: