AIRES library reference: xslant.

c     Routine directly callable from C/C++ (it could be necessary to
c     append an underscore, "_", to its name depending on the C/C++
c     compiler used).
      function xslant(xvertical, x0, cozenith, zcground)
c     Evaluating the slanted path (in g/cm2) along an axis starting at
c     altitude zcground, and inclination given by cozenith (cosine of
c     the zenith angle). The path corresponds to the "segment" whose
c     ends possess vertical depths xvertical and x0, respectively. The
c     entire segment must lie above the ground surface. If x0 = 0 then
c     the path up to the top of the atmosphere is evaluated. The
c     earth's curvature is taken into account to perform the
c     calculations.
c     The slanted path is defined via
c                    point2
c          Xslant = INTEGRAL'  rho(zv) dl
c                    point1
c     where rho(zv) is the density of the medium, which is supposed to
c     depend only on zv, the vertical altitude. The prime in the
c     integral indicates that it is done along the specified axis.
c     The atmospheric model must be properly initialized before
c     invoking this routine.
c     Written by: S. J. Sciutto, La Plata 1998, 2001, 2003.
c     Parameters:
c     ==========
c     xvertical....... (input, double precision) Vertical depth of
c                      the path's end point 1 (g/cm2).
c     x0.............. (input, double precision) Vertical depth of
c                      the path's end point 2 (g/cm2). If x0 is zero
c                      then the path is evaluated up to the top of the
c                      atmosphere.
c     cozenith........ (input, double precision) Cosine of the axis'
c                      zenith angle. Must be in the range (0, 1].
c     zcground........ (input, double precision) Central altitude of
c                      the ground level, which defines the intersection
c                      between the slanted axis and the z-axis.
c     Return value: (double precision) The corresponding slanted path,
c     ============  in g/cm2. The returned value is zero on error.

These pages are maintained by
Sergio J. Sciutto
Back to AIRES library index
Back to AIRES home [HOME]