PySDM_examples.deJong_Azimi.cloudy_data_0d
1import numpy as np 2 3MOM_data = { 4 "Golovin": { 5 "t": np.array( 6 [ 7 0.0, 8 10.0, 9 20.0, 10 30.0, 11 40.0, 12 50.0, 13 60.0, 14 70.0, 15 80.0, 16 90.0, 17 100.0, 18 110.0, 19 120.0, 20 ] 21 ), 22 "aM0": np.array( 23 [ 24 100.0, 25 60.41666666666108, 26 36.50173611110569, 27 22.05313223379227, 28 13.32376739124922, 29 8.049776132212969, 30 4.863406413211965, 31 2.938308041315548, 32 1.775227774961471, 33 1.072533447372555, 34 0.6479889577875856, 35 0.39149332866333514, 36 0.23652721940076313, 37 ] 38 ), 39 "aM1": np.array( 40 [ 41 10.0, 42 10.0, 43 10.0, 44 10.000000000000002, 45 10.000000000000002, 46 9.999999999999991, 47 10.000000000000007, 48 10.000000000000002, 49 9.999999999999956, 50 9.999999999999966, 51 9.999999999999872, 52 10.000000000000062, 53 10.000000000000108, 54 ] 55 ), 56 "aM2": np.array( 57 [ 58 2.0, 59 5.333333333333704, 60 14.222222222223586, 61 37.92592592592996, 62 101.13580246914695, 63 269.69547325105935, 64 719.1879286694924, 65 1917.8344764519795, 66 5114.225270538599, 67 13637.934054769523, 68 36367.82414605215, 69 96980.86438947004, 70 258615.63837192673, 71 ] 72 ), 73 "aDists": np.array( 74 [ 75 [100.0, 0.1, 1.0], 76 [4.863406413211965, 69.86262088575636, 0.029431646782264554], 77 [0.2365272194007631, 25819.285405029976, 0.0016374749145528168], 78 ] 79 ), 80 "bM0": np.array( 81 [ 82 [ 83 100.0, 84 56.396636657301194, 85 31.895287187247146, 86 18.294957761981376, 87 10.660817516471703, 88 6.290591039424345, 89 3.7443020392268367, 90 2.2413636351719695, 91 1.3464994349989194, 92 0.8107024270698977, 93 0.48877185047542177, 94 0.2949239965078066, 95 0.17804598720980272, 96 ], 97 [ 98 1.0e-6, 99 4.0200190988219875, 100 4.606433069139865, 101 3.758158525366545, 102 2.6629362775448184, 103 1.7591745237414111, 104 1.1190966122081454, 105 0.696938913153295, 106 0.42872454995378795, 107 0.2618284523381802, 108 0.15921539069060084, 109 0.09656819650467219, 110 0.058480487065190184, 111 ], 112 ] 113 ), 114 "bM1": np.array( 115 [ 116 [ 117 10.0, 118 6.354890848885312, 119 3.815306159255158, 120 2.2491773102504062, 121 1.3276630984776567, 122 0.7884211778943261, 123 0.4708461282568907, 124 0.28236149853196507, 125 0.16980173598586254, 126 0.10229476631106786, 127 0.06169480173887945, 128 0.03723418494509775, 129 0.022481089449094008, 130 ], 131 [ 132 1.0e-5, 133 3.6451213684882426, 134 6.184706455055038, 135 7.75083495382188, 136 8.672348657723711, 137 9.211590117984043, 138 9.529164800771182, 139 9.717649154060018, 140 9.83020871507534, 141 9.89771554128402, 142 9.938315405635903, 143 9.962775953501628, 144 9.977529002206483, 145 ], 146 ] 147 ), 148 "bM2": np.array( 149 [ 150 [ 151 2.0, 152 1.3279917682405473, 153 0.815368316209652, 154 0.484406234346711, 155 0.2866707891616959, 156 0.1703612676331259, 157 0.10175389387661611, 158 0.06101894293570185, 159 0.03669236327334466, 160 0.022103697711640027, 161 0.013330436832490668, 162 0.00804502307966626, 163 0.0048573143944901, 164 ], 165 [ 166 0.0002, 167 4.005880364582517, 168 13.40830335654937, 169 37.445419449158464, 170 100.85962493652575, 171 269.5533461366049, 172 719.1621396161328, 173 1917.977828184842, 174 5114.73836146586, 175 13639.390826596753, 176 36371.78858209591, 177 96991.55465118046, 178 258644.4048115109, 179 ], 180 ] 181 ), 182 "bDists1": np.array( 183 [ 184 [100.0, 3.7443020392268367, 0.17804598720980272], 185 [0.1, 0.09035855289429393, 0.08979659938892945], 186 [1.0, 1.3916783213625432, 1.4061293503242225], 187 ] 188 ).T, 189 "bDists2": np.array( 190 [ 191 [1.0e-6, 1.1190998946220907, 0.05848079412663001], 192 [1.0, 66.94723995610593, 25749.457224681548], 193 [3.0, 0.12718999335996053, 0.00662584618581055], 194 ] 195 ).T, 196 }, 197 "Geometric": { 198 "t": np.array( 199 [ 200 0.0, 201 20.0, 202 40.0, 203 60.0, 204 80.0, 205 100.0, 206 120.0, 207 140.0, 208 160.0, 209 180.0, 210 200.0, 211 220.0, 212 240.0, 213 ] 214 ), 215 "aM0": np.array( 216 [ 217 100.0, 218 83.09411530706723, 219 69.0668818692553, 220 57.42833556073367, 221 47.771892319320436, 222 39.760132461015466, 223 33.11299776856557, 224 27.59799298461195, 225 23.022057367028655, 226 19.22483391072762, 227 16.073111976893788, 228 13.456251155748634, 229 11.282412772916492, 230 ] 231 ), 232 "aM1": np.array( 233 [ 234 10.0, 235 9.999999999999998, 236 9.999999999999995, 237 9.999999999999998, 238 10.000000000000002, 239 10.000000000000004, 240 9.999999999999996, 241 9.999999999999996, 242 9.999999999999998, 243 9.999999999999998, 244 9.999999999999998, 245 9.999999999999988, 246 9.99999999999999, 247 ] 248 ), 249 "aM2": np.array( 250 [ 251 2.0, 252 2.8852151953194194, 253 4.151238620053756, 254 5.950500062169887, 255 8.485178624869318, 256 12.0129886133155, 257 16.84418831764611, 258 23.32258967740333, 259 31.78481092193408, 260 42.500028311068455, 261 55.60685101272762, 262 71.07585744751505, 263 88.7218493043981, 264 ] 265 ), 266 "aDists": np.array( 267 [ 268 [100.0, 0.1, 1.0], 269 [33.12307697837833, 1.386085645680255, 0.21781072474918567], 270 [11.288994413248544, 8.001687588099609, 0.11070395858333783], 271 ] 272 ), 273 "bM0": np.array( 274 [ 275 [ 276 100.0, 277 83.08735661284119, 278 69.025406058178, 279 57.28710188741288, 280 47.43565678877952, 281 39.167760941853565, 282 32.289941767938586, 283 26.61554482816344, 284 21.95171830328991, 285 18.122227928945424, 286 14.977334881119972, 287 12.392963329500928, 288 10.267170524593165, 289 ], 290 [ 291 1.0e-6, 292 0.006754965534216085, 293 0.0414689469896406, 294 0.14122472581368004, 295 0.33622210669071656, 296 0.5923365303626719, 297 0.822966528248585, 298 0.9822712716512826, 299 1.0700559173008375, 300 1.1022183815198043, 301 1.0953066097618704, 302 1.062767038894064, 303 1.014700508663851, 304 ], 305 ] 306 ), 307 "bM1": np.array( 308 [ 309 [ 310 10.0, 311 9.998176966812272, 312 9.983636638206804, 313 9.914575608610326, 314 9.648425409616147, 315 9.022885248658554, 316 8.145992984487792, 317 7.201239010259968, 318 6.289684097629613, 319 5.449920021025792, 320 4.694685159768059, 321 4.0259638571565235, 322 3.4404908668086382, 323 ], 324 [ 325 3.0e-6, 326 0.001826033187727366, 327 0.016366361793195474, 328 0.08542739138967238, 329 0.3515775903838529, 330 0.9771177513414488, 331 1.8540100155122126, 332 2.798763989740037, 333 3.7103189023703886, 334 4.55008297897421, 335 5.3053178402319405, 336 5.974039142843473, 337 6.559512133191355, 338 ], 339 ] 340 ), 341 "bM2": np.array( 342 [ 343 [ 344 2.0, 345 2.884419791825751, 346 4.140031572580394, 347 5.841267749589132, 348 7.602242548420331, 349 8.443150132466357, 350 8.291305172948194, 351 7.649776302203622, 352 6.858287226599269, 353 6.056020626291585, 354 5.296417497830004, 355 4.600656776850216, 356 3.976052070235482, 357 ], 358 [ 359 1.2e-5, 360 0.0008129200697375103, 361 0.011232418249945829, 362 0.10926457152459564, 363 0.8829072854441761, 364 3.5698485088870786, 365 8.556214924636981, 366 15.693719076620901, 367 25.002698714778422, 368 36.65139187807145, 369 50.772315702453426, 370 67.35212811578029, 371 86.19438593789097, 372 ], 373 ] 374 ), 375 "bDists1": np.array( 376 [ 377 [100.0, 0.1, 1.0], 378 [10.267170524593165, 0.8205678828026997, 0.4083712124798814], 379 ] 380 ), 381 "bDists2": np.array( 382 [ 383 [1.0e-6, 1.0, 3.0], 384 [1.014700508663851, 6.675884373858597, 0.968333274304494], 385 ] 386 ), 387 }, 388} 389 390for kernel, data in MOM_data.items(): 391 for prefix in ("a", "b"): 392 data[f"{prefix}Moments"] = np.array([data[f"{prefix}M{i}"] for i in range(3)])
MOM_data =
{'Golovin': {'t': array([ 0., 10., 20., 30., 40., 50., 60., 70., 80., 90., 100.,
110., 120.]), 'aM0': array([100. , 60.41666667, 36.50173611, 22.05313223,
13.32376739, 8.04977613, 4.86340641, 2.93830804,
1.77522777, 1.07253345, 0.64798896, 0.39149333,
0.23652722]), 'aM1': array([10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10.]), 'aM2': array([2.00000000e+00, 5.33333333e+00, 1.42222222e+01, 3.79259259e+01,
1.01135802e+02, 2.69695473e+02, 7.19187929e+02, 1.91783448e+03,
5.11422527e+03, 1.36379341e+04, 3.63678241e+04, 9.69808644e+04,
2.58615638e+05]), 'aDists': array([[1.00000000e+02, 1.00000000e-01, 1.00000000e+00],
[4.86340641e+00, 6.98626209e+01, 2.94316468e-02],
[2.36527219e-01, 2.58192854e+04, 1.63747491e-03]]), 'bM0': array([[1.00000000e+02, 5.63966367e+01, 3.18952872e+01, 1.82949578e+01,
1.06608175e+01, 6.29059104e+00, 3.74430204e+00, 2.24136364e+00,
1.34649943e+00, 8.10702427e-01, 4.88771850e-01, 2.94923997e-01,
1.78045987e-01],
[1.00000000e-06, 4.02001910e+00, 4.60643307e+00, 3.75815853e+00,
2.66293628e+00, 1.75917452e+00, 1.11909661e+00, 6.96938913e-01,
4.28724550e-01, 2.61828452e-01, 1.59215391e-01, 9.65681965e-02,
5.84804871e-02]]), 'bM1': array([[1.00000000e+01, 6.35489085e+00, 3.81530616e+00, 2.24917731e+00,
1.32766310e+00, 7.88421178e-01, 4.70846128e-01, 2.82361499e-01,
1.69801736e-01, 1.02294766e-01, 6.16948017e-02, 3.72341849e-02,
2.24810894e-02],
[1.00000000e-05, 3.64512137e+00, 6.18470646e+00, 7.75083495e+00,
8.67234866e+00, 9.21159012e+00, 9.52916480e+00, 9.71764915e+00,
9.83020872e+00, 9.89771554e+00, 9.93831541e+00, 9.96277595e+00,
9.97752900e+00]]), 'bM2': array([[2.00000000e+00, 1.32799177e+00, 8.15368316e-01, 4.84406234e-01,
2.86670789e-01, 1.70361268e-01, 1.01753894e-01, 6.10189429e-02,
3.66923633e-02, 2.21036977e-02, 1.33304368e-02, 8.04502308e-03,
4.85731439e-03],
[2.00000000e-04, 4.00588036e+00, 1.34083034e+01, 3.74454194e+01,
1.00859625e+02, 2.69553346e+02, 7.19162140e+02, 1.91797783e+03,
5.11473836e+03, 1.36393908e+04, 3.63717886e+04, 9.69915547e+04,
2.58644405e+05]]), 'bDists1': array([[1.00000000e+02, 1.00000000e-01, 1.00000000e+00],
[3.74430204e+00, 9.03585529e-02, 1.39167832e+00],
[1.78045987e-01, 8.97965994e-02, 1.40612935e+00]]), 'bDists2': array([[1.00000000e-06, 1.00000000e+00, 3.00000000e+00],
[1.11909989e+00, 6.69472400e+01, 1.27189993e-01],
[5.84807941e-02, 2.57494572e+04, 6.62584619e-03]]), 'aMoments': array([[1.00000000e+02, 6.04166667e+01, 3.65017361e+01, 2.20531322e+01,
1.33237674e+01, 8.04977613e+00, 4.86340641e+00, 2.93830804e+00,
1.77522777e+00, 1.07253345e+00, 6.47988958e-01, 3.91493329e-01,
2.36527219e-01],
[1.00000000e+01, 1.00000000e+01, 1.00000000e+01, 1.00000000e+01,
1.00000000e+01, 1.00000000e+01, 1.00000000e+01, 1.00000000e+01,
1.00000000e+01, 1.00000000e+01, 1.00000000e+01, 1.00000000e+01,
1.00000000e+01],
[2.00000000e+00, 5.33333333e+00, 1.42222222e+01, 3.79259259e+01,
1.01135802e+02, 2.69695473e+02, 7.19187929e+02, 1.91783448e+03,
5.11422527e+03, 1.36379341e+04, 3.63678241e+04, 9.69808644e+04,
2.58615638e+05]]), 'bMoments': array([[[1.00000000e+02, 5.63966367e+01, 3.18952872e+01, 1.82949578e+01,
1.06608175e+01, 6.29059104e+00, 3.74430204e+00, 2.24136364e+00,
1.34649943e+00, 8.10702427e-01, 4.88771850e-01, 2.94923997e-01,
1.78045987e-01],
[1.00000000e-06, 4.02001910e+00, 4.60643307e+00, 3.75815853e+00,
2.66293628e+00, 1.75917452e+00, 1.11909661e+00, 6.96938913e-01,
4.28724550e-01, 2.61828452e-01, 1.59215391e-01, 9.65681965e-02,
5.84804871e-02]],
[[1.00000000e+01, 6.35489085e+00, 3.81530616e+00, 2.24917731e+00,
1.32766310e+00, 7.88421178e-01, 4.70846128e-01, 2.82361499e-01,
1.69801736e-01, 1.02294766e-01, 6.16948017e-02, 3.72341849e-02,
2.24810894e-02],
[1.00000000e-05, 3.64512137e+00, 6.18470646e+00, 7.75083495e+00,
8.67234866e+00, 9.21159012e+00, 9.52916480e+00, 9.71764915e+00,
9.83020872e+00, 9.89771554e+00, 9.93831541e+00, 9.96277595e+00,
9.97752900e+00]],
[[2.00000000e+00, 1.32799177e+00, 8.15368316e-01, 4.84406234e-01,
2.86670789e-01, 1.70361268e-01, 1.01753894e-01, 6.10189429e-02,
3.66923633e-02, 2.21036977e-02, 1.33304368e-02, 8.04502308e-03,
4.85731439e-03],
[2.00000000e-04, 4.00588036e+00, 1.34083034e+01, 3.74454194e+01,
1.00859625e+02, 2.69553346e+02, 7.19162140e+02, 1.91797783e+03,
5.11473836e+03, 1.36393908e+04, 3.63717886e+04, 9.69915547e+04,
2.58644405e+05]]])}, 'Geometric': {'t': array([ 0., 20., 40., 60., 80., 100., 120., 140., 160., 180., 200.,
220., 240.]), 'aM0': array([100. , 83.09411531, 69.06688187, 57.42833556,
47.77189232, 39.76013246, 33.11299777, 27.59799298,
23.02205737, 19.22483391, 16.07311198, 13.45625116,
11.28241277]), 'aM1': array([10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10.]), 'aM2': array([ 2. , 2.8852152 , 4.15123862, 5.95050006, 8.48517862,
12.01298861, 16.84418832, 23.32258968, 31.78481092, 42.50002831,
55.60685101, 71.07585745, 88.7218493 ]), 'aDists': array([[100. , 0.1 , 1. ],
[ 33.12307698, 1.38608565, 0.21781072],
[ 11.28899441, 8.00168759, 0.11070396]]), 'bM0': array([[1.00000000e+02, 8.30873566e+01, 6.90254061e+01, 5.72871019e+01,
4.74356568e+01, 3.91677609e+01, 3.22899418e+01, 2.66155448e+01,
2.19517183e+01, 1.81222279e+01, 1.49773349e+01, 1.23929633e+01,
1.02671705e+01],
[1.00000000e-06, 6.75496553e-03, 4.14689470e-02, 1.41224726e-01,
3.36222107e-01, 5.92336530e-01, 8.22966528e-01, 9.82271272e-01,
1.07005592e+00, 1.10221838e+00, 1.09530661e+00, 1.06276704e+00,
1.01470051e+00]]), 'bM1': array([[1.00000000e+01, 9.99817697e+00, 9.98363664e+00, 9.91457561e+00,
9.64842541e+00, 9.02288525e+00, 8.14599298e+00, 7.20123901e+00,
6.28968410e+00, 5.44992002e+00, 4.69468516e+00, 4.02596386e+00,
3.44049087e+00],
[3.00000000e-06, 1.82603319e-03, 1.63663618e-02, 8.54273914e-02,
3.51577590e-01, 9.77117751e-01, 1.85401002e+00, 2.79876399e+00,
3.71031890e+00, 4.55008298e+00, 5.30531784e+00, 5.97403914e+00,
6.55951213e+00]]), 'bM2': array([[2.00000000e+00, 2.88441979e+00, 4.14003157e+00, 5.84126775e+00,
7.60224255e+00, 8.44315013e+00, 8.29130517e+00, 7.64977630e+00,
6.85828723e+00, 6.05602063e+00, 5.29641750e+00, 4.60065678e+00,
3.97605207e+00],
[1.20000000e-05, 8.12920070e-04, 1.12324182e-02, 1.09264572e-01,
8.82907285e-01, 3.56984851e+00, 8.55621492e+00, 1.56937191e+01,
2.50026987e+01, 3.66513919e+01, 5.07723157e+01, 6.73521281e+01,
8.61943859e+01]]), 'bDists1': array([[100. , 0.1 , 1. ],
[ 10.26717052, 0.82056788, 0.40837121]]), 'bDists2': array([[1.00000000e-06, 1.00000000e+00, 3.00000000e+00],
[1.01470051e+00, 6.67588437e+00, 9.68333274e-01]]), 'aMoments': array([[100. , 83.09411531, 69.06688187, 57.42833556,
47.77189232, 39.76013246, 33.11299777, 27.59799298,
23.02205737, 19.22483391, 16.07311198, 13.45625116,
11.28241277],
[ 10. , 10. , 10. , 10. ,
10. , 10. , 10. , 10. ,
10. , 10. , 10. , 10. ,
10. ],
[ 2. , 2.8852152 , 4.15123862, 5.95050006,
8.48517862, 12.01298861, 16.84418832, 23.32258968,
31.78481092, 42.50002831, 55.60685101, 71.07585745,
88.7218493 ]]), 'bMoments': array([[[1.00000000e+02, 8.30873566e+01, 6.90254061e+01, 5.72871019e+01,
4.74356568e+01, 3.91677609e+01, 3.22899418e+01, 2.66155448e+01,
2.19517183e+01, 1.81222279e+01, 1.49773349e+01, 1.23929633e+01,
1.02671705e+01],
[1.00000000e-06, 6.75496553e-03, 4.14689470e-02, 1.41224726e-01,
3.36222107e-01, 5.92336530e-01, 8.22966528e-01, 9.82271272e-01,
1.07005592e+00, 1.10221838e+00, 1.09530661e+00, 1.06276704e+00,
1.01470051e+00]],
[[1.00000000e+01, 9.99817697e+00, 9.98363664e+00, 9.91457561e+00,
9.64842541e+00, 9.02288525e+00, 8.14599298e+00, 7.20123901e+00,
6.28968410e+00, 5.44992002e+00, 4.69468516e+00, 4.02596386e+00,
3.44049087e+00],
[3.00000000e-06, 1.82603319e-03, 1.63663618e-02, 8.54273914e-02,
3.51577590e-01, 9.77117751e-01, 1.85401002e+00, 2.79876399e+00,
3.71031890e+00, 4.55008298e+00, 5.30531784e+00, 5.97403914e+00,
6.55951213e+00]],
[[2.00000000e+00, 2.88441979e+00, 4.14003157e+00, 5.84126775e+00,
7.60224255e+00, 8.44315013e+00, 8.29130517e+00, 7.64977630e+00,
6.85828723e+00, 6.05602063e+00, 5.29641750e+00, 4.60065678e+00,
3.97605207e+00],
[1.20000000e-05, 8.12920070e-04, 1.12324182e-02, 1.09264572e-01,
8.82907285e-01, 3.56984851e+00, 8.55621492e+00, 1.56937191e+01,
2.50026987e+01, 3.66513919e+01, 5.07723157e+01, 6.73521281e+01,
8.61943859e+01]]])}}