arts-docs-2.4/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& x0,

                 const Numeric& y0,

                 const Numeric& z0,

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


void cycle_lat_lon(Numeric& lat, Numeric& lon);


#endif  // geodetic_h