test/src/test_pressure.cpp

   1 /* Copyright_License {
   2
   3   XCSoar Glide Computer - http://www.xcsoar.org/
   4   Copyright (C) 2000-2013 The XCSoar Project
   5   A detailed list of copyright holders can be found in the file "AUTHORS".
   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
   9   as published by the Free Software Foundation; either version 2
  10   of the License, or (at your option) any later version.
  11
  12   This program is distributed in the hope that it will be useful,
  13   but WITHOUT ANY WARRANTY; without even the implied warranty of
  14   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15   GNU 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., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
  20 }
  21 */
  22
  23
  24 #include "test_debug.hpp"
  25 #include "Atmosphere/AirDensity.hpp"
  26 #include "Atmosphere/Pressure.hpp"
  27
  28 static bool
  29 test_find_qnh()
  30 {
  31   AtmosphericPressure sp = AtmosphericPressure::Standard().QNHAltitudeToStaticPressure(fixed(100));
  32   AtmosphericPressure pres =
  33     AtmosphericPressure::FindQNHFromPressure(sp, fixed(100));
  34   return fabs(pres.GetHectoPascal() - fixed(1013.25)) < fixed(0.01);
  35 }
  36
  37 static bool
  38 test_find_qnh2()
  39 {
  40   AtmosphericPressure sp = AtmosphericPressure::Standard().QNHAltitudeToStaticPressure(fixed(100));
  41   AtmosphericPressure pres =
  42     AtmosphericPressure::FindQNHFromPressure(sp, fixed(120));
  43   if (verbose) {
  44     printf("%g\n",FIXED_DOUBLE(pres.GetHectoPascal()));
  45   }
  46   return fabs(pres.GetHectoPascal() - fixed(1015.6)) < fixed(0.1);
  47   // example, QNH=1014, ps=100203
  48   // alt= 100
  49   // alt_known = 120
  50   // qnh= 1016
  51 }
  52
  53 static bool
  54 test_qnh_to_static()
  55 {
  56   AtmosphericPressure pres = AtmosphericPressure::Standard();
  57   fixed p0 = pres.QNHAltitudeToStaticPressure(fixed(0)).GetPascal();
  58   if (verbose) {
  59     printf("%g\n",FIXED_DOUBLE(p0));
  60   }
  61   return fabs(p0-fixed(101325))<fixed(0.1);
  62 }
  63
  64 static bool
  65 test_qnh_round()
  66 {
  67   AtmosphericPressure sp = AtmosphericPressure::Standard().QNHAltitudeToStaticPressure(fixed(100));
  68   AtmosphericPressure pres =
  69     AtmosphericPressure::FindQNHFromPressure(sp, fixed(120));
  70   fixed h0 = pres.PressureAltitudeToQNHAltitude(fixed(100));
  71   if (verbose) {
  72     printf("%g\n",FIXED_DOUBLE(h0));
  73   }
  74   return fabs(h0-fixed(120))<fixed(1);
  75 }
  76
  77 static bool
  78 test_qnh_round2()
  79 {
  80   AtmosphericPressure sp = AtmosphericPressure::Standard().QNHAltitudeToStaticPressure(fixed(100));
  81   AtmosphericPressure pres =
  82     AtmosphericPressure::FindQNHFromPressure(sp, fixed(120));
  83   fixed h0 = pres.StaticPressureToQNHAltitude(pres);
  84   if (verbose) {
  85     printf("%g %g\n", FIXED_DOUBLE(pres.GetPascal()), FIXED_DOUBLE(h0));
  86   }
  87   return fabs(h0)<fixed(1);
  88 }
  89
  90 static bool
  91 test_isa_pressure(const fixed alt, const fixed prat)
  92 {
  93   const AtmosphericPressure pres = AtmosphericPressure::Standard();
  94   fixed p0 = pres.QNHAltitudeToStaticPressure(alt).GetPascal();
  95   if (verbose) {
  96     printf("%g\n",FIXED_DOUBLE(p0));
  97   }
  98   return fabs(p0/fixed(101325)-prat)<fixed(0.001);
  99 }
 100
 101 static bool
 102 test_isa_density(const fixed alt, const fixed prat)
 103 {
 104   fixed p0 = AirDensity(alt);
 105   if (verbose) {
 106     printf("%g\n",FIXED_DOUBLE(p0));
 107   }
 108   return fabs(p0/fixed(1.225)-prat)<fixed(0.001);
 109 }
 110
 111 int main(int argc, char** argv)
 112 {
 113
 114   if (!ParseArgs(argc,argv)) {
 115     return 0;
 116   }
 117
 118   plan_tests(9);
 119
 120   ok(test_find_qnh(),"find qnh 0-0",0);
 121   ok(test_find_qnh2(),"find qnh 100-120",0);
 122
 123   ok(test_qnh_to_static(),"qnh to static",0);
 124   ok(test_qnh_round(),"qnh round trip",0);
 125
 126   ok(test_qnh_round2(),"qnh round 2",0);
 127
 128   ok(test_isa_pressure(fixed(1524), fixed(0.8320)), "isa pressure at 1524m",0);
 129   ok(test_isa_pressure(fixed(6096), fixed(0.4594)), "isa pressure at 6096m",0);
 130
 131   ok(test_isa_density(fixed(1524), fixed(0.8617)), "isa density at 1524m",0);
 132   ok(test_isa_density(fixed(6096), fixed(0.5328)), "isa density at 6096m",0);
 133
 134   return exit_status();
 135 }