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

#![warn(missing_docs)]
//! Library to represent an atmospheric sounding.

#![recursion_limit = "1024"]
extern crate chrono;
use chrono::NaiveDateTime;

#[macro_use]
extern crate error_chain;

pub mod error;

/// All the variables stored in the sounding.
///
/// Missing values are -9999.0. The upper air profile variables are stored in parallel vectors.
/// If a profile lacks a certain variable, e.g. cloud fraction, that whole vector has length 0
/// instead of being full of missing values.
pub struct Sounding {
    // Station info section
    /// station number, USAF number, eg 727730
    pub num: i32,
    /// Valid time of sounding
    pub valid_time: NaiveDateTime,
    /// Difference in model initialization time and `valid_time` in hours.
    pub lead_time: i32,
    /// Latitude of grid point used to make sounding.
    pub lat: f32,
    /// Longitude of grid point used to make sounding.
    pub lon: f32,
    /// Elevation of grid point in meters, this is in model terrain, not necessarily the same as
    /// the real world.
    pub elevation: f32,

    // Sounding Indexes
    /// Showalter index
    pub show: f32,
    /// Lifted index
    pub li: f32,
    /// Severe Weather Threat Index
    pub swet: f32,
    /// K-index
    pub kinx: f32,
    /// Lifting Condensation Level, or LCL (hPa), pressure vertical coordinate.
    pub lclp: f32,
    /// Precipitable Water (mm)
    pub pwat: f32,
    /// Total-Totals
    pub totl: f32,
    /// Convective Available Potential Energy, or CAPE. (J/kg)
    pub cape: f32,
    /// Temperature at LCL (K)
    pub lclt: f32,
    /// Convective Inhibitive Energy, or CIN (J/kg)
    pub cins: f32,
    /// Equilibrium Level (hPa), pressure vertical coordinate
    pub eqlv: f32,
    /// Level of Free Convection (hPa), pressure vertical coordinate
    pub lfc: f32,
    /// Bulk Richardson Number
    pub brch: f32,

    // Upper air profile
    /// Pressure (hPa) profile
    pub pressure: Vec<f32>,
    /// Temperature (c) profile
    pub temperature: Vec<f32>,
    /// Wet-bulb (c) profile
    pub wet_bulb: Vec<f32>,
    /// Dew Point (C) profile
    pub dew_point: Vec<f32>,
    /// Equivalent Potential Temperature (K) profile
    pub theta_e: Vec<f32>,
    /// Wind direction (degrees) profile
    pub direction: Vec<f32>,
    /// Wind speed (knots) profile
    pub speed: Vec<f32>,
    /// Vertical velocity (Pa/sec), pressure vertical coordinate
    pub omega: Vec<f32>,
    /// Geopotential Height (m) profile
    pub height: Vec<f32>,
    /// Cloud coverage fraction in percent
    pub cloud_fraction: Vec<f32>,

    // Surface data
    /// Surface pressure reduce to mean sea level (hPa)
    pub mslp: f32,
    /// Surface pressure (hPa)
    pub station_pres: f32,
    /// Low cloud fraction
    pub low_cloud: f32,
    /// Mid cloud fraction
    pub mid_cloud: f32,
    /// Hi cloud fraction
    pub hi_cloud: f32,
    /// U - wind speed (m/s) (West -> East is positive)
    pub uwind: f32,
    /// V - wind speed (m/s) (South -> North is positive)
    pub vwind: f32,
}