PySDM_examples.Szumowski_et_al_1998.sounding

 1import numpy as np
 2
 3from PySDM.physics import constants as const
 4from PySDM.physics import si
 5
 6pressure = (
 7    np.array(
 8        (
 9            1014,
10            1010,
11            1000,
12            990,
13            975,
14            960,
15            950,
16            925,
17            900,
18            875,
19            850,
20            825,
21            800,
22            790,
23            775,
24            765,
25            755,
26            745,
27            730,
28            715,
29            700,
30        )
31    )
32    * si.hPa
33)
34
35temperature = (
36    np.array(
37        (
38            25.2,
39            24.8,
40            23.6,
41            22.5,
42            21.8,
43            20.5,
44            19.9,
45            18.2,
46            16.8,
47            14.8,
48            13.3,
49            11.9,
50            11.0,
51            11.3,
52            10.9,
53            11.2,
54            10.2,
55            11.0,
56            11.2,
57            10.0,
58            8.8,
59        )
60    )
61    * si.kelvin
62    + const.T0
63)
64
65mixing_ratio = (
66    np.array(
67        (
68            14.5,
69            14.5,
70            14.5,
71            14.0,
72            13.7,
73            13.9,
74            13.9,
75            10.3,
76            10.3,
77            10.0,
78            9.9,
79            8.9,
80            7.9,
81            4.0,
82            2.3,
83            1.2,
84            1.2,
85            0.9,
86            0.6,
87            2.0,
88            1.6,
89        )
90    )
91    * si.g
92    / si.kg
93)