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])