satkit 0.3.14

Satellite Toolkit
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

use pyo3::prelude::*;
use pyo3::types::{PyFloat, PyTuple};
use pyo3::wrap_pyfunction;

use crate::nrlmsise;
use crate::AstroTime;

use super::PyAstroTime;
use super::PyITRFCoord;

///
/// NRL MSISE-00 Density Model
///
/// Args:
///     coord (itrfcoord): ITRF coordinate at which to compute density & temperature
///     time (satkit.time): Optional time object representing instant at which to compute
///
/// Returns:
///     tuple: (rho, T) where rho is atmosphere mass density in kg / m^3 and T is temperature in Kelvin
#[pyfunction(name = "nrlmsise")]
#[pyo3(signature=(*args))]
fn pynrlmsise(args: &Bound<'_, PyTuple>) -> PyResult<(f64, f64)> {
    //let args = args.as_gil_ref();

    if args.len() == 0 {
        return Err(pyo3::exceptions::PyTypeError::new_err(
            "Invalid number of arguments",
        ));
    }

    let time: Option<AstroTime> = {
        if args
            .get_item(args.len() - 1)?
            .is_instance_of::<PyAstroTime>()
        {
            Some(
                args.get_item(args.len() - 1)?
                    .extract::<PyAstroTime>()
                    .unwrap()
                    .inner,
            )
        } else {
            None
        }
    };
    if args.get_item(0)?.is_instance_of::<PyITRFCoord>() {
        let itrf = args.get_item(0)?.extract::<PyITRFCoord>().unwrap().inner;
        Ok(nrlmsise::nrlmsise(
            itrf.hae() / 1.0e3,
            Some(itrf.latitude_rad()),
            Some(itrf.longitude_rad()),
            time,
            true,
        ))
    } else if args.get_item(0)?.is_instance_of::<PyFloat>() {
        let altitude = args.get_item(0)?.extract::<f64>().unwrap();
        let latitude: Option<f64> = {
            if args.len() > 1 && args.get_item(1)?.is_instance_of::<PyFloat>() {
                Some(args.get_item(1)?.extract::<f64>().unwrap())
            } else {
                None
            }
        };
        let longitude: Option<f64> = {
            if args.len() > 2 && args.get_item(2)?.is_instance_of::<PyFloat>() {
                Some(args.get_item(2)?.extract::<f64>().unwrap())
            } else {
                None
            }
        };
        Ok(nrlmsise::nrlmsise(
            altitude / 1.0e3,
            latitude,
            longitude,
            time,
            true,
        ))
    } else {
        return Err(pyo3::exceptions::PyTypeError::new_err("Invalid arguments"));
    }
}

#[pymodule]
pub fn density(_py: Python, m: &Bound<'_, PyModule>) -> PyResult<()> {
    m.add_function(wrap_pyfunction!(pynrlmsise, m)?).unwrap();
    Ok(())
}