earths 0.0.4

High-fidelity Earth simulation engine — orbit, atmosphere, geology, hydrology, biosphere, terrain, lighting, rendering, satellites, and temporal systems with full scientific coupling
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

use sciforge::hub::domain::geology::seismology::{
    epicenter_distance, gutenberg_richter, moment_magnitude, omori_aftershock, p_wave_velocity,
    peak_ground_acceleration, s_wave_velocity, seismic_energy, seismic_moment, snell_seismic,
    travel_time,
};
pub struct Earthquake {
    pub latitude: f64,
    pub longitude: f64,
    pub depthkm: f64,
    pub seismicmomentnm: f64,
}
impl Earthquake {
    pub fn frommoment(lat: f64, lon: f64, depthkm: f64, momentnm: f64) -> Self {
        Self {
            latitude: lat,
            longitude: lon,
            depthkm,
            seismicmomentnm: momentnm,
        }
    }
    pub fn frommagnitude(lat: f64, lon: f64, depthkm: f64, mw: f64) -> Self {
        let moment = seismic_moment(mw);
        Self::frommoment(lat, lon, depthkm, moment)
    }
    pub fn momentmagnitude(&self) -> f64 {
        moment_magnitude(self.seismicmomentnm)
    }
    pub fn energyjoules(&self) -> f64 {
        seismic_energy(self.momentmagnitude())
    }
    pub fn pwavearrivaltime(&self, distancekm: f64, vp: f64) -> f64 {
        travel_time(distancekm * 1000.0, vp)
    }
    pub fn swavearrivaltime(&self, distancekm: f64, vs: f64) -> f64 {
        travel_time(distancekm * 1000.0, vs)
    }
    pub fn pgaatdistance(&self, distancekm: f64) -> f64 {
        peak_ground_acceleration(1.0, 1.0, self.momentmagnitude(), distancekm)
    }
}
pub fn aftershockrate(mainshockmagnitude: f64, timedays: f64) -> f64 {
    let a = mainshockmagnitude - 1.2;
    let c = 0.05;
    let p = 1.1;
    omori_aftershock(a.powf(10.0), c, p, timedays)
}
pub fn annualfrequency(magnitude: f64) -> f64 {
    let a = 8.0;
    let b = 1.0;
    gutenberg_richter(a, b, magnitude)
}
pub fn crustpwavevelocity(density: f64, bulkmodulus: f64, shearmodulus: f64) -> f64 {
    p_wave_velocity(bulkmodulus, shearmodulus, density)
}
pub fn crustswavevelocity(density: f64, shearmodulus: f64) -> f64 {
    s_wave_velocity(shearmodulus, density)
}
pub fn distancefromsptime(tp: f64, ts: f64, vp: f64, vs: f64) -> f64 {
    epicenter_distance(vp, vs, ts - tp)
}
pub fn refractionangle(v1: f64, v2: f64, angle1rad: f64) -> f64 {
    snell_seismic(v1, angle1rad, v2)
}