arts-docs-2.2/doxygen/geodetic_8h_source.html

/* Copyright (C) 2012

   Patrick Eriksson <Patrick.Eriksson@chalmers.se>


   This program is free software; you can redistribute it and/or modify it

   under the terms of the GNU General Public License as published by the

   Free Software Foundation; either version 2, or (at your option) any

   later version.


   This program is distributed in the hope that it will be useful,

   but WITHOUT ANY WARRANTY; without even the implied warranty of

   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

   GNU General Public License for more details.


   You should have received a copy of the GNU General Public License

   along with this program; if not, write to the Free Software

   Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,

   USA. */


/*===========================================================================

  === File description

  ===========================================================================*/


#ifndef geodetic_h

#define geodetic_h


#include "interpolation.h"

#include "matpackI.h"


// 2D:


void cart2pol(

            Numeric&   r,

            Numeric&   lat,

      const Numeric&   x,

      const Numeric&   z,

      const Numeric&   lat0,

      const Numeric&   za0 );


void cart2poslos(

             Numeric&   r,

             Numeric&   lat,

             Numeric&   za,

       const Numeric&   x,

       const Numeric&   z,

       const Numeric&   dx,

       const Numeric&   dz,

       const Numeric&   ppc,

       const Numeric&   lat0,

       const Numeric&   za0 );


void distance2D(

            Numeric&   l,

      const Numeric&   r1,

      const Numeric&   lat1,

      const Numeric&   r2,

      const Numeric&   lat2 );


/*

void geomtanpoint2d(

             Numeric&    r_tan,

             Numeric&    lat_tan,

     ConstVectorView    refellipsoid,

       const Numeric&    r,

       const Numeric&    lat,

       const Numeric&    za );

*/


void line_circle_intersect(

         Numeric&   x,

         Numeric&   z,

   const Numeric&   xl,

   const Numeric&   zl,

   const Numeric&   dx,

   const Numeric&   dz,

   const Numeric&   xc,

   const Numeric&   zc,

   const Numeric&   r );


void pol2cart(

            Numeric&   x,

            Numeric&   z,

      const Numeric&   r,

      const Numeric&   lat );


void poslos2cart(

             Numeric&   x,

             Numeric&   z,

             Numeric&   dx,

             Numeric&   dz,

       const Numeric&   r,

       const Numeric&   lat,

       const Numeric&   za );


// 3D:


void cart2poslos(

             Numeric&   r,

             Numeric&   lat,

             Numeric&   lon,

             Numeric&   za,

             Numeric&   aa,

       const Numeric&   x,

       const Numeric&   y,

       const Numeric&   z,

       const Numeric&   dx,

       const Numeric&   dy,

       const Numeric&   dz,

       const Numeric&   ppc,

       const Numeric&   lat0,

       const Numeric&   lon0,

       const Numeric&   za0,

       const Numeric&   aa0 );


void cart2sph(

             Numeric&   r,

             Numeric&   lat,

             Numeric&   lon,

       const Numeric&   x,

       const Numeric&   y,

       const Numeric&   z,

       const Numeric&   lat0,

       const Numeric&   lon0,

       const Numeric&   za0,

       const Numeric&   aa0 );


void distance3D(

            Numeric&   l,

      const Numeric&   r1,

      const Numeric&   lat1,

      const Numeric&   lon1,

      const Numeric&   r2,

      const Numeric&   lat2,

      const Numeric&   lon2 );


void geompath_tanpos_3d(

             Numeric&    r_tan,

             Numeric&    lat_tan,

             Numeric&    lon_tan,

             Numeric&    l_tan,

       const Numeric&    r,

       const Numeric&    lat,

       const Numeric&    lon,

       const Numeric&    za,

       const Numeric&    aa,

       const Numeric&    ppc );


/*

void geomtanpoint(

             Numeric&    r_tan,

             Numeric&    lat_tan,

             Numeric&    lon_tan,

     ConstVectorView    refellipsoid,

       const Numeric&    r,

       const Numeric&    lat,

       const Numeric&    lon,

       const Numeric&    za,

       const Numeric&    aa );

*/


void latlon_at_aa(

         Numeric&   lat2,

         Numeric&   lon2,

   const Numeric&   lat1,

   const Numeric&   lon1,

   const Numeric&   aa,

   const Numeric&   ddeg );


void line_sphere_intersect(

         Numeric&   x,

         Numeric&   y,

         Numeric&   z,

   const Numeric&   xl,

   const Numeric&   yl,

   const Numeric&   zl,

   const Numeric&   dx,

   const Numeric&   dy,

   const Numeric&   dz,

   const Numeric&   xc,

   const Numeric&   yc,

   const Numeric&   zc,

   const Numeric&   r );


void los2xyz(

         Numeric&   za,

         Numeric&   aa,

   const Numeric&   r1,

   const Numeric&   lat1,

   const Numeric&   lon1,

   const Numeric&   x1,

   const Numeric&   y1,

   const Numeric&   z1,

   const Numeric&   x2,

   const Numeric&   y2,

   const Numeric&   z2 );


void poslos2cart(

              Numeric&   x,

              Numeric&   y,

              Numeric&   z,

              Numeric&   dx,

              Numeric&   dy,

              Numeric&   dz,

        const Numeric&   r,

        const Numeric&   lat,

        const Numeric&   lon,

        const Numeric&   za,

        const Numeric&   aa );


Numeric pos2refell_r(

       const Index&     atmosphere_dim,

       ConstVectorView  refellipsoid,

       ConstVectorView  lat_grid,

       ConstVectorView  lon_grid,

       ConstVectorView  rte_pos );


Numeric refell2r(

       ConstVectorView  refellipsoid,

       const Numeric&   lat );


Numeric refell2d(

       ConstVectorView  refellipsoid,

       ConstVectorView  lat_grid,

       const GridPos    gp );


Numeric sphdist(

   const Numeric&   lat1,

   const Numeric&   lon1,

   const Numeric&   lat2,

   const Numeric&   lon2 );


void sph2cart(

            Numeric&   x,

            Numeric&   y,

            Numeric&   z,

      const Numeric&   r,

      const Numeric&   lat,

      const Numeric&   lon );


// coord transform


void lon_shiftgrid(

            Vector&    longrid_out,

      ConstVectorView  longrid_in,

      const Numeric    lon );


#endif  // geodetic_h