iasinat-lib 0.1.4

Reader for IASI NAT L1C and L2 files
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

use super::*;

// EUM/OPS-EPS/MAN/04/0033 v3E p.60
#[derive(Debug)]
pub struct GiadrL2 {
    pub contents:GiadrL2Contents,
    pub error_data:GiadrL2ErrorData,
    // pub brescia:GiadrL2Brescia,
    // pub forli:GiadrL2Forli,
}

// EUM/OPS-EPS/MAN/04/0033 v3E p.60
// Scaling factor applied
#[derive(Debug)]
pub struct GiadrL2Contents {
    /// Pressure levels on which retrieved temperature profiles
    /// are given (Pa)
    pub pressure_levels_temp:Vec<f64>,

    /// Pressure levels on which retrieved humidity profiles
    /// are given (Pa)
    pub pressure_levels_humidity:Vec<f64>,

    /// Pressure levels on which retrieved ozone profiles
    /// are given (Pa)
    pub pressure_levels_ozone:Vec<f64>,

    /// Wavelengths for surface emissivity (micron)
    pub surface_emissivity_wavelengths:Vec<f64>,
}

// EUM/OPS-EPS/MAN/04/0033 v3E p.60
#[derive(Debug)]
pub struct GiadrL2ErrorData {
    pub num_temperature_pcs:u8,
    pub num_water_vapour_pcs:u8,
    pub num_ozone_pcs:u8,
}

impl GiadrL2 {
    pub fn read_bin<R:Read+Seek>(rd:&mut NatReader<R>,rec:&Grh)->Result<Self> {
	let contents = GiadrL2Contents::read_bin(rd,rec)?;
	let error_data = GiadrL2ErrorData::read_bin(rd,rec)?;
	Ok(Self {
	    contents,
	    error_data
	})
    }
}

impl GiadrL2Contents {
    pub fn read_bin<R:Read+Seek>(rd:&mut NatReader<R>,rec:&Grh)->Result<Self> {
	rec.seek_to_record(rd,20)?;
	let pressure_levels_temp = read_vec_and_scale(rd,1e2)?;
	rec.seek_to_record(rd,425)?;
	let pressure_levels_humidity = read_vec_and_scale(rd,1e2)?;
	rec.seek_to_record(rd,830)?;
	let pressure_levels_ozone = read_vec_and_scale(rd,1e2)?;
	rec.seek_to_record(rd,1235)?;
	let surface_emissivity_wavelengths = read_vec_and_scale(rd,1e4)?;
	Ok(Self {
	    pressure_levels_temp,
	    pressure_levels_humidity,
	    pressure_levels_ozone,
	    surface_emissivity_wavelengths
	})
    }
}

impl GiadrL2ErrorData {
    pub fn read_bin<R:Read+Seek>(rd:&mut NatReader<R>,rec:&Grh)->Result<Self> {
	rec.seek_to_record(rd,1284)?;
	let num_temperature_pcs = u8::read_bin(rd)?;
	let num_water_vapour_pcs = u8::read_bin(rd)?;
	let num_ozone_pcs = u8::read_bin(rd)?;
	Ok(Self {
	    num_temperature_pcs,
	    num_water_vapour_pcs,
	    num_ozone_pcs
	})
    }

    fn ntri(n:u8)->usize {
	let n = n as usize;
	(n*(n + 1))/2
    }

    pub fn nerrt(&self)->usize { Self::ntri(self.num_temperature_pcs) }

    pub fn nerrw(&self)->usize { Self::ntri(self.num_water_vapour_pcs) }

    pub fn nerro(&self)->usize { Self::ntri(self.num_ozone_pcs) }
}