jupiters 0.0.3

Jupiter celestial simulation crate for the MilkyWay SolarSystem workspace
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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115

use sciforge::hub::prelude::constants::elements::atomic_mass;
use sciforge::hub::prelude::constants::{G, N_A, R_GAS};
use std::sync::LazyLock;

pub const JUPITERMASS: f64 = 1.8982e27;
pub const JUPITEREQUATORIALRADIUS: f64 = 71492000.0;
pub const JUPITERPOLARRADIUS: f64 = 66854000.0;
pub const JUPITERMEANRADIUS: f64 = 69911000.0;

pub fn mh2() -> f64 {
    2.0 * atomic_mass(1) * 1e-3
}
pub fn mhe() -> f64 {
    atomic_mass(2) * 1e-3
}
pub fn mch4() -> f64 {
    (atomic_mass(6) + 4.0 * atomic_mass(1)) * 1e-3
}
pub fn mnh3() -> f64 {
    (atomic_mass(7) + 3.0 * atomic_mass(1)) * 1e-3
}
pub fn mh2o() -> f64 {
    (2.0 * atomic_mass(1) + atomic_mass(8)) * 1e-3
}
pub fn mh2s() -> f64 {
    (2.0 * atomic_mass(1) + atomic_mass(16)) * 1e-3
}
pub fn matmosphere() -> f64 {
    0.898 * mh2() + 0.102 * mhe() + 0.003 * mch4()
}
pub fn mwaterice() -> f64 {
    mh2o()
}
pub fn mammoniaice() -> f64 {
    mnh3()
}
pub fn mmethaneice() -> f64 {
    mch4()
}
pub fn msio2() -> f64 {
    (atomic_mass(14) + 2.0 * atomic_mass(8)) * 1e-3
}
pub fn mfe() -> f64 {
    atomic_mass(26) * 1e-3
}

pub static SURFACEGRAVITY: LazyLock<f64> =
    LazyLock::new(|| G * JUPITERMASS / (JUPITEREQUATORIALRADIUS * JUPITEREQUATORIALRADIUS));
pub static SCALEHEIGHT: LazyLock<f64> =
    LazyLock::new(|| R_GAS * ONEBARTEMPK / (matmosphere() * *SURFACEGRAVITY));
pub static ONEBARDENSITY: LazyLock<f64> = LazyLock::new(|| {
    let rspecific = R_GAS / matmosphere();
    ONEBARPRESSURE / (rspecific * ONEBARTEMPK)
});
pub static WATERICEDENSITY: LazyLock<f64> = LazyLock::new(|| {
    let a = 4.5181e-10f64;
    let c = 7.3560e-10;
    let vcell = a * a * c * (60.0f64.to_radians().sin());
    4.0 * mh2o() / (N_A * vcell)
});
pub static IRONDENSITY: LazyLock<f64> = LazyLock::new(|| {
    let a = 2.8665e-10f64;
    2.0 * mfe() / (N_A * a.powi(3))
});
pub static SILICATEDENSITY: LazyLock<f64> = LazyLock::new(|| {
    let a = 4.9133e-10f64;
    let c = 5.4053e-10;
    let vcell = a * a * c * (60.0f64.to_radians().sin());
    3.0 * msio2() / (N_A * vcell)
});

pub const ONEBARPRESSURE: f64 = 101325.0;
pub const ONEBARTEMPK: f64 = 165.0;
pub const CELSIUSTOKELVIN: f64 = 273.15;
pub const SECONDSPERYEAR: f64 = 3.15576e7;
pub const OMEGAJUPITER: f64 = 1.7585e-4;
pub const JUPITERFLATTENING: f64 = 0.06487;
pub const J2JUPITER: f64 = 14.736e-3;
pub const BONDALBEDO: f64 = 0.503;
pub const CLOUDTOPTEMPMINK: f64 = 110.0;
pub const CLOUDTOPTEMPMAXK: f64 = 165.0;
pub const KTTNTTOJOULE: f64 = 4.184e12;
pub const MTTNTTOJOULE: f64 = 4.184e15;
pub const STONYASTEROIDDENSITY: f64 = 3000.0;
pub const METALLICHYDROGENDENSITY: f64 = 1000.0;
pub fn ironasteroiddensity() -> f64 {
    *IRONDENSITY
}
pub fn iceasteroiddensity() -> f64 {
    *WATERICEDENSITY
}
pub const H2LOVE: f64 = 1.15;
pub const K2LOVE: f64 = 0.565;
pub const L2LOVE: f64 = 0.16;
pub const IOECCENTRICITY: f64 = 0.0041;
pub const IOINCLINATIONDEG: f64 = 0.05;
pub const CLOUDEXTINCTIONCOEFF: f64 = 0.15;
pub const CLOUDLWCSCALE: f64 = 0.20;
pub const METHANEABSORPTION: f64 = 0.003;
pub const TROPOPAUSEM: f64 = 50000.0;
pub const ALBEDOEMAALPHA: f64 = 0.001;
pub const INTERNALHEATINGWM2: f64 = 5.44;
pub const MAXWINDSMS: f64 = 180.0;

pub mod atmosphere;
pub mod biosphere;
pub mod geodata;
pub mod geology;
pub mod hydrology;
pub mod lighting;
pub mod physics;
pub mod rendering;
pub mod satellites;
pub mod temporal;
pub mod terrain;