cffdrs 0.11.0

Rust implementation of the Canadian Forest Fire Danger Rating System
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

use std::f64::consts;

use super::FbpFuelType;

/// Surface fuel consumption SFC (km/m^2)
///
/// * `ffmc` - Fine fuel moisture code
/// * `bui` - Buildup index
pub fn surface_fuel_consumption(fuel_type: FbpFuelType, ffmc: f64, bui: f64) -> f64 {
    match fuel_type {
        FbpFuelType::C1 => (1.5 * (1. - consts::E.powf(-0.230 * (ffmc - 81.)))).max(0.),
        FbpFuelType::C2 | FbpFuelType::M3 | FbpFuelType::M4 => {
            5.0 * (1. - consts::E.powf(-0.0115 * bui)).powf(1.)
        }
        FbpFuelType::C3 | FbpFuelType::C4 => 5.0 * (1. - consts::E.powf(-0.0164 * bui)).powf(2.24),
        FbpFuelType::C5 | FbpFuelType::C6 => 5.0 * (1. - consts::E.powf(-0.0149 * bui)).powf(2.48),
        FbpFuelType::C7 => {
            let ffc = (2. * (1. - consts::E.powf(-0.104 * (ffmc - 70.)))).max(0.);
            let wfc = 1.5 * (1. - consts::E.powf(-0.0201 * bui));
            ffc + wfc
        }
        FbpFuelType::D1 => 1.5 * (1. - consts::E.powf(-0.0183 * bui)),
        FbpFuelType::M1 | FbpFuelType::M2 => {
            let pc = 50.;
            let ph = 50.;
            (pc / 100. * surface_fuel_consumption(FbpFuelType::C2, ffmc, bui))
                + (ph / 100. * surface_fuel_consumption(FbpFuelType::D1, ffmc, bui))
        }
        FbpFuelType::O1b | FbpFuelType::O1a => {
            // Grass fuel load (standard = 0.3km/m^2)
            
            0.3
        }
        FbpFuelType::S1 => {
            let ffc = 4.0 * (1. - consts::E.powf(-0.025 * bui));
            let wfc = 4.0 * (1. - consts::E.powf(-0.034 * bui));
            ffc + wfc
        }
        FbpFuelType::S2 => {
            let ffc = 10.0 * (1. - consts::E.powf(-0.013 * bui));
            let wfc = 6.0 * (1. - consts::E.powf(-0.060 * bui));
            ffc + wfc
        }
        FbpFuelType::S3 => {
            let ffc = 12.0 * (1. - consts::E.powf(-0.0166 * bui));
            let wfc = 20.0 * (1. - consts::E.powf(-0.0210 * bui));
            ffc + wfc
        }

        FbpFuelType::NonFuel => 0.,
    }
}