mirror of
https://github.com/acemod/ACE3.git
synced 2024-08-30 18:23:18 +00:00
35 lines
1.1 KiB
Plaintext
35 lines
1.1 KiB
Plaintext
#include "..\script_component.hpp"
|
|
/*
|
|
* Author: Ruthberg
|
|
* Calculates the air density
|
|
*
|
|
* Arguments:
|
|
* 0: temperature - degrees celsius <NUMBER>
|
|
* 1: pressure - hPa <NUMBER>
|
|
* 2: relativeHumidity - value between 0.0 and 1.0 <NUMBER>
|
|
*
|
|
* Return Value:
|
|
* density of air - kg * m^(-3) <NUMBER>
|
|
*
|
|
* Example:
|
|
* [0, 1020, 0.5] call ace_weather_fnc_calculateAirDensity
|
|
*
|
|
* Public: No
|
|
*/
|
|
|
|
params ["_temperature", "_pressure", "_relativeHumidity"];
|
|
|
|
_pressure = _pressure * 100; // hPa to Pa
|
|
|
|
if (_relativeHumidity > 0) then {
|
|
// Saturation vapor pressure calculated according to: http://wahiduddin.net/calc/density_algorithms.htm
|
|
// 610.78 gives pressure in Pa - https://en.wikipedia.org/wiki/Density_of_air
|
|
private _pSat = 610.78 * 10 ^ ((7.5 * _temperature) / (_temperature + 237.3));
|
|
private _vaporPressure = _relativeHumidity * _pSat;
|
|
private _partialPressure = _pressure - _vaporPressure;
|
|
|
|
(_partialPressure * DRY_AIR_MOLAR_MASS + _vaporPressure * WATER_VAPOR_MOLAR_MASS) / (UNIVERSAL_GAS_CONSTANT * KELVIN(_temperature))
|
|
} else {
|
|
_pressure / (SPECIFIC_GAS_CONSTANT_DRY_AIR * KELVIN(_temperature))
|
|
};
|