PySDM_examples.Jensen_and_Nugent_2017.table_3

 1import numpy as np
 2
 3from PySDM.physics import si
 4
 5NA = (
 6    np.asarray(
 7        [
 8            111800,
 9            68490,
10            38400,
11            21820,
12            13300,
13            8496,
14            5486,
15            3805,
16            2593,
17            1919,
18            1278,
19            998.4,
20            777.9,
21            519.5,
22            400.5,
23            376.9,
24            265.3,
25            212.4,
26            137.8,
27            121.4,
28            100.9,
29            122.2,
30            50.64,
31            38.3,
32            55.47,
33            21.45,
34            12.95,
35            43.23,
36            26.26,
37            30.5,
38            4.385,
39            4.372,
40            4.465,
41            4.395,
42            4.427,
43            4.411,
44            0,
45            0,
46            0,
47            4.522,
48            0,
49            4.542,
50        ]
51    )
52    / si.m**3
53)
54
55RD = np.linspace(0.8, 9, num=len(NA), endpoint=True) * si.um
NA = array([1.118e+05, 6.849e+04, 3.840e+04, 2.182e+04, 1.330e+04, 8.496e+03, 5.486e+03, 3.805e+03, 2.593e+03, 1.919e+03, 1.278e+03, 9.984e+02, 7.779e+02, 5.195e+02, 4.005e+02, 3.769e+02, 2.653e+02, 2.124e+02, 1.378e+02, 1.214e+02, 1.009e+02, 1.222e+02, 5.064e+01, 3.830e+01, 5.547e+01, 2.145e+01, 1.295e+01, 4.323e+01, 2.626e+01, 3.050e+01, 4.385e+00, 4.372e+00, 4.465e+00, 4.395e+00, 4.427e+00, 4.411e+00, 0.000e+00, 0.000e+00, 0.000e+00, 4.522e+00, 0.000e+00, 4.542e+00])
RD = array([8.0e-07, 1.0e-06, 1.2e-06, 1.4e-06, 1.6e-06, 1.8e-06, 2.0e-06, 2.2e-06, 2.4e-06, 2.6e-06, 2.8e-06, 3.0e-06, 3.2e-06, 3.4e-06, 3.6e-06, 3.8e-06, 4.0e-06, 4.2e-06, 4.4e-06, 4.6e-06, 4.8e-06, 5.0e-06, 5.2e-06, 5.4e-06, 5.6e-06, 5.8e-06, 6.0e-06, 6.2e-06, 6.4e-06, 6.6e-06, 6.8e-06, 7.0e-06, 7.2e-06, 7.4e-06, 7.6e-06, 7.8e-06, 8.0e-06, 8.2e-06, 8.4e-06, 8.6e-06, 8.8e-06, 9.0e-06])