sciforge 0.0.3

A comprehensive scientific computing library in pure Rust with zero dependencies
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93

//! Dispatch handler for atmospheric thermodynamics functions.

use super::super::params::*;
use crate::hub::domain::common::errors::{HubError, HubResult};
use crate::hub::domain::meteorology as meteo;
use crate::hub::engine::experience::runner::RunOutput;

pub(super) fn dispatch(func: &str, p: &Params) -> HubResult<RunOutput> {
    match func {
        "barometric_formula" => Ok(RunOutput::Scalar(meteo::atmosphere::barometric_formula(
            get_f(p, "p0")?,
            get_f(p, "m")?,
            get_f(p, "g")?,
            get_f(p, "h")?,
            get_f(p, "t")?,
        ))),
        "scale_height" => Ok(RunOutput::Scalar(meteo::atmosphere::scale_height(
            get_f(p, "t")?,
            get_f(p, "m")?,
            get_f(p, "g")?,
        ))),
        "lapse_rate_dry" => Ok(RunOutput::Scalar(meteo::atmosphere::lapse_rate_dry())),
        "lapse_rate_moist" => Ok(RunOutput::Scalar(meteo::atmosphere::lapse_rate_moist(
            get_f(p, "t")?,
            get_f(p, "l_v")?,
            get_f(p, "r_s")?,
        ))),
        "potential_temperature" => Ok(RunOutput::Scalar(meteo::atmosphere::potential_temperature(
            get_f(p, "t")?,
            get_f(p, "p")?,
            get_f(p, "p0")?,
        ))),
        "virtual_temperature" => Ok(RunOutput::Scalar(meteo::atmosphere::virtual_temperature(
            get_f(p, "t")?,
            get_f(p, "r")?,
        ))),
        "mixing_ratio" => Ok(RunOutput::Scalar(meteo::atmosphere::mixing_ratio(
            get_f(p, "e")?,
            get_f(p, "p")?,
        ))),
        "saturation_vapor_pressure" => Ok(RunOutput::Scalar(
            meteo::atmosphere::saturation_vapor_pressure(get_f(p, "t_celsius")?),
        )),
        "relative_humidity" => Ok(RunOutput::Scalar(meteo::atmosphere::relative_humidity(
            get_f(p, "e")?,
            get_f(p, "es")?,
        ))),
        "dew_point" => Ok(RunOutput::Scalar(meteo::atmosphere::dew_point(
            get_f(p, "t")?,
            get_f(p, "rh")?,
        ))),
        "density_altitude" => Ok(RunOutput::Scalar(meteo::atmosphere::density_altitude(
            get_f(p, "pressure_altitude")?,
            get_f(p, "temperature_c")?,
        ))),
        "brunt_vaisala_frequency" => Ok(RunOutput::Scalar(
            meteo::atmosphere::brunt_vaisala_frequency(
                get_f(p, "g")?,
                get_f(p, "theta")?,
                get_f(p, "dtheta_dz")?,
            ),
        )),
        "mie_phase" => Ok(RunOutput::Scalar(meteo::atmosphere::mie_phase(
            get_f(p, "cos_theta")?,
            get_f(p, "g")?,
        ))),
        "lifted_condensation_level" => Ok(RunOutput::Scalar(
            meteo::atmosphere::lifted_condensation_level(
                get_f(p, "t_surface")?,
                get_f(p, "dew_point")?,
            ),
        )),
        "dry_adiabatic_temperature" => Ok(RunOutput::Scalar(
            meteo::atmosphere::dry_adiabatic_temperature(
                get_f(p, "t_surface")?,
                get_f(p, "altitude")?,
            ),
        )),
        "convective_available_potential_energy" => Ok(RunOutput::Scalar(
            meteo::atmosphere::convective_available_potential_energy(
                get_v(p, "t_parcel")?,
                get_v(p, "t_env")?,
                get_f(p, "dz")?,
            ),
        )),
        "convective_inhibition" => Ok(RunOutput::Scalar(meteo::atmosphere::convective_inhibition(
            get_v(p, "t_parcel")?,
            get_v(p, "t_env")?,
            get_f(p, "dz")?,
        ))),
        _ => Err(HubError::InvalidInput(format!("unknown function: {func}"))),
    }
}