2011-12-28 05:32:21 -04:00
|
|
|
/// -*- tab-width: 4; Mode: C++; c-basic-offset: 4; indent-tabs-mode: nil -*-
|
2011-10-13 11:22:03 -03:00
|
|
|
|
2015-08-11 03:28:42 -03:00
|
|
|
#include "AP_Baro.h"
|
|
|
|
#include <AP_HAL/AP_HAL.h>
|
2014-01-05 01:36:12 -04:00
|
|
|
|
2012-10-11 14:53:21 -03:00
|
|
|
extern const AP_HAL::HAL& hal;
|
2011-10-13 11:22:03 -03:00
|
|
|
|
2014-10-19 16:22:51 -03:00
|
|
|
AP_Baro_HIL::AP_Baro_HIL(AP_Baro &baro) :
|
|
|
|
AP_Baro_Backend(baro)
|
2011-10-13 11:22:03 -03:00
|
|
|
{
|
2014-10-19 16:22:51 -03:00
|
|
|
_instance = _frontend.register_sensor();
|
2011-10-13 11:22:03 -03:00
|
|
|
}
|
|
|
|
|
2013-07-20 10:00:41 -03:00
|
|
|
// ==========================================================================
|
|
|
|
// based on tables.cpp from http://www.pdas.com/atmosdownload.html
|
|
|
|
|
|
|
|
/*
|
2014-10-19 16:22:51 -03:00
|
|
|
Compute the temperature, density, and pressure in the standard atmosphere
|
|
|
|
Correct to 20 km. Only approximate thereafter.
|
2013-07-20 10:00:41 -03:00
|
|
|
*/
|
2014-10-19 16:22:51 -03:00
|
|
|
void AP_Baro::SimpleAtmosphere(
|
2013-07-20 10:00:41 -03:00
|
|
|
const float alt, // geometric altitude, km.
|
|
|
|
float& sigma, // density/sea-level standard density
|
|
|
|
float& delta, // pressure/sea-level standard pressure
|
|
|
|
float& theta) // temperature/sea-level standard temperature
|
|
|
|
{
|
2014-10-19 16:22:51 -03:00
|
|
|
const float REARTH = 6369.0f; // radius of the Earth (km)
|
|
|
|
const float GMR = 34.163195f; // gas constant
|
|
|
|
float h=alt*REARTH/(alt+REARTH); // geometric to geopotential altitude
|
|
|
|
|
|
|
|
if (h < 11.0f) {
|
|
|
|
// Troposphere
|
|
|
|
theta=(288.15f-6.5f*h)/288.15f;
|
|
|
|
delta=powf(theta, GMR/6.5f);
|
|
|
|
} else {
|
|
|
|
// Stratosphere
|
|
|
|
theta=216.65f/288.15f;
|
|
|
|
delta=0.2233611f*expf(-GMR*(h-11.0f)/216.65f);
|
2013-07-20 10:00:41 -03:00
|
|
|
}
|
|
|
|
|
2014-10-19 16:22:51 -03:00
|
|
|
sigma = delta/theta;
|
2013-07-20 10:00:41 -03:00
|
|
|
}
|
|
|
|
|
|
|
|
|
2014-10-19 16:22:51 -03:00
|
|
|
/*
|
|
|
|
convert an altitude in meters above sea level to a presssure and temperature
|
|
|
|
*/
|
|
|
|
void AP_Baro::setHIL(float altitude_msl)
|
2011-10-13 11:22:03 -03:00
|
|
|
{
|
2013-07-20 10:00:41 -03:00
|
|
|
float sigma, delta, theta;
|
2013-07-20 01:37:47 -03:00
|
|
|
const float p0 = 101325;
|
2013-05-02 02:26:38 -03:00
|
|
|
|
2014-07-08 07:26:54 -03:00
|
|
|
SimpleAtmosphere(altitude_msl*0.001f, sigma, delta, theta);
|
2013-07-20 10:00:41 -03:00
|
|
|
float p = p0 * delta;
|
2014-07-08 07:26:54 -03:00
|
|
|
float T = 303.16f * theta - 273.16f; // Assume 30 degrees at sea level - converted to degrees Kelvin
|
2013-05-02 02:26:38 -03:00
|
|
|
|
2014-10-19 16:22:51 -03:00
|
|
|
setHIL(0, p, T);
|
2011-10-13 11:22:03 -03:00
|
|
|
}
|
2011-12-11 00:30:24 -04:00
|
|
|
|
2014-10-19 16:22:51 -03:00
|
|
|
/*
|
|
|
|
set HIL pressure and temperature for an instance
|
|
|
|
*/
|
|
|
|
void AP_Baro::setHIL(uint8_t instance, float pressure, float temperature)
|
|
|
|
{
|
|
|
|
if (instance >= _num_sensors) {
|
|
|
|
// invalid
|
|
|
|
return;
|
|
|
|
}
|
|
|
|
sensors[instance].pressure = pressure;
|
|
|
|
sensors[instance].temperature = temperature;
|
|
|
|
sensors[instance].last_update_ms = hal.scheduler->millis();
|
2011-12-11 00:30:24 -04:00
|
|
|
}
|