src/Navaids/ils.hxx

   1 // ils.hxx -- navaid class
   2 //
   3 // Written by Curtis Olson, started April 2000.
   4 //
   5 // Copyright (C) 2000  Curtis L. Olson - curt@flightgear.org
   6 //
   7 // This program is free software; you can redistribute it and/or
   8 // modify it under the terms of the GNU General Public License as
   9 // published by the Free Software Foundation; either version 2 of the
  10 // License, or (at your option) any later version.
  11 //
  12 // This program is distributed in the hope that it will be useful, but
  13 // WITHOUT ANY WARRANTY; without even the implied warranty of
  14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 // General Public License for more details.
  16 //
  17 // You should have received a copy of the GNU General Public License
  18 // along with this program; if not, write to the Free Software
  19 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  20 //
  21 // $Id$
  22
  23
  24 #ifndef _FG_ILS_HXX
  25 #define _FG_ILS_HXX
  26
  27
  28 #include <simgear/compiler.h>
  29 #include <simgear/math/sg_geodesy.hxx>
  30 #include <simgear/misc/sgstream.hxx>
  31
  32 #ifdef SG_HAVE_STD_INCLUDES
  33 #  include <istream>
  34 #elif defined( SG_HAVE_NATIVE_SGI_COMPILERS )
  35 #  include <iostream.h>
  36 #elif defined( __BORLANDC__ )
  37 #  include <iostream>
  38 #else
  39 #  include <istream.h>
  40 #endif
  41
  42 #if ! defined( SG_HAVE_NATIVE_SGI_COMPILERS )
  43 SG_USING_STD(istream);
  44 #endif
  45
  46
  47 #define FG_ILS_DEFAULT_RANGE 18
  48
  49 class FGILS {
  50
  51     char ilstype;
  52     char ilstypename[6];
  53     char aptcode[5];
  54     char rwyno[4];
  55     int  locfreq;
  56     char locident[5];           // official ident
  57     double locheading;
  58     double loclat;
  59     double loclon;
  60     double x, y, z;
  61     bool has_gs;
  62     double gselev;
  63     double gsangle;
  64     double gslat;
  65     double gslon;
  66     double gs_x, gs_y, gs_z;
  67     bool has_dme;
  68     double dmelat;
  69     double dmelon;
  70     double dme_x, dme_y, dme_z;
  71     double omlat;
  72     double omlon;
  73     double mmlat;
  74     double mmlon;
  75     double imlat;
  76     double imlon;
  77
  78     // for failure modeling
  79     string trans_ident;         // transmitted ident
  80     bool loc_failed;            // localizer failed?
  81     bool gs_failed;             // glide slope failed?
  82     bool dme_failed;            // dme failed?
  83
  84 public:
  85
  86     inline FGILS(void) {}
  87     inline ~FGILS(void) {}
  88
  89     inline char get_ilstype() const { return ilstype; }
  90     inline char *get_ilstypename() { return ilstypename; }
  91     inline char *get_aptcode() { return aptcode; }
  92     inline char *get_rwyno() { return rwyno; }
  93     inline int get_locfreq() const { return locfreq; }
  94     inline char *get_locident() { return locident; }
  95     inline string get_trans_ident() { return trans_ident; }
  96     inline double get_locheading() const { return locheading; }
  97     inline double get_loclat() const { return loclat; }
  98     inline double get_loclon() const { return loclon; }
  99     inline double get_x() const { return x; }
 100     inline double get_y() const { return y; }
 101     inline double get_z() const { return z; }
 102     inline bool get_has_gs() const { return has_gs; }
 103     inline double get_gselev() const { return gselev; }
 104     inline double get_gsangle() const { return gsangle; }
 105     inline double get_gslat() const { return gslat; }
 106     inline double get_gslon() const { return gslon; }
 107     inline double get_gs_x() const { return gs_x; }
 108     inline double get_gs_y() const { return gs_y; }
 109     inline double get_gs_z() const { return gs_z; }
 110     inline bool get_has_dme() const { return has_dme; }
 111     inline double get_dmelat() const { return dmelat; }
 112     inline double get_dmelon() const { return dmelon; }
 113     inline double get_dme_x() const { return dme_x; }
 114     inline double get_dme_y() const { return dme_y; }
 115     inline double get_dme_z() const { return dme_z; }
 116     inline double get_omlat() const { return omlat; }
 117     inline double get_omlon() const { return omlon; }
 118     inline double get_mmlat() const { return mmlat; }
 119     inline double get_mmlon() const { return mmlon; }
 120     inline double get_imlat() const { return imlat; }
 121     inline double get_imlon() const { return imlon; }
 122
 123     friend istream& operator>> ( istream&, FGILS& );
 124 };
 125
 126
 127 inline istream&
 128 operator >> ( istream& in, FGILS& i )
 129 {
 130     double f;
 131     in >> i.ilstype;
 132
 133     if ( i.ilstype == '[' )
 134           return in;
 135
 136     in >> i.ilstypename >> i.aptcode >> i.rwyno
 137        >> f >> i.locident >> i.locheading >> i.loclat >> i.loclon
 138        >> i.gselev >> i.gsangle >> i.gslat >> i.gslon
 139        >> i.dmelat >> i.dmelon
 140        >> i.omlat >> i.omlon
 141        >> i.mmlat >> i.mmlon
 142        >> i.imlat >> i.imlon;
 143
 144
 145     i.locfreq = (int)(f*100.0 + 0.5);
 146
 147     // generate cartesian coordinates
 148     Point3D geod, cart;
 149
 150     geod = Point3D( i.loclon * SGD_DEGREES_TO_RADIANS, i.loclat * SGD_DEGREES_TO_RADIANS, i.gselev );
 151     cart = sgGeodToCart( geod );
 152     i.x = cart.x();
 153     i.y = cart.y();
 154     i.z = cart.z();
 155
 156     if ( i.gslon < SG_EPSILON && i.gslat < SG_EPSILON ) {
 157         i.has_gs = false;
 158     } else {
 159         i.has_gs = true;
 160
 161         geod = Point3D( i.gslon * SGD_DEGREES_TO_RADIANS, i.gslat * SGD_DEGREES_TO_RADIANS, i.gselev );
 162         cart = sgGeodToCart( geod );
 163         i.gs_x = cart.x();
 164         i.gs_y = cart.y();
 165         i.gs_z = cart.z();
 166         // cout << "gs = " << cart << endl;
 167     }
 168
 169     if ( i.dmelon < SG_EPSILON && i.dmelat < SG_EPSILON ) {
 170         i.has_dme = false;
 171     } else {
 172         i.has_dme = true;
 173
 174         geod = Point3D( i.dmelon * SGD_DEGREES_TO_RADIANS, i.dmelat * SGD_DEGREES_TO_RADIANS, i.gselev);
 175         cart = sgGeodToCart( geod );
 176         i.dme_x = cart.x();
 177         i.dme_y = cart.y();
 178         i.dme_z = cart.z();
 179         // cout << "dme = " << cart << endl;
 180     }
 181
 182     i.trans_ident = "I";
 183     i.trans_ident += i.locident;
 184     i.loc_failed = i.gs_failed = i.dme_failed = false;
 185
 186     // return in >> skipeol;
 187     return in;
 188 }
 189
 190
 191 #endif // _FG_ILS_HXX