PySDM_examples.deJong_Azimi.cloudy_data
1import numpy as np 2 3Z = np.linspace(0.0, 3000.0, 20) / 1000 4NC = np.array( 5 [ 6 [ 7 0.002325942800455995, 8 0.0028045022356462572, 9 0.003452211545804375, 10 0.004361877640877965, 11 0.005702479217508372, 12 0.007812475748456827, 13 0.011468941392117338, 14 0.01888454890111773, 15 0.039528401551653096, 16 0.2356787814538211, 17 9.90222080639631, 18 9.944680372872213, 19 9.948230520394494, 20 9.92828984604994, 21 9.49434950135968, 22 0.060141075852780945, 23 0.023917441369691872, 24 0.013590312401779991, 25 0.00893428827291265, 26 0.006376138306461998, 27 ], 28 [ 29 0.005322014757656553, 30 0.010480022319206528, 31 0.025889114642834388, 32 0.07027637984750808, 33 0.1846578336060549, 34 0.4341420453704166, 35 0.8709339382008977, 36 1.4526980887520684, 37 2.016594898467415, 38 2.4045238186923696, 39 2.579754970394283, 40 2.578563190147033, 41 2.3791690313423275, 42 1.7682308790889112, 43 0.7721032508385711, 44 0.01991858759231357, 45 0.010155146035355484, 46 0.0057644342569305234, 47 0.002836371204245902, 48 0.0008570715559907928, 49 ], 50 [ 51 0.08662893464816847, 52 0.15520311679346102, 53 0.2522813401209703, 54 0.371195396866207, 55 0.49992082338075977, 56 0.6298742026425381, 57 0.7609509414360253, 58 0.8962365323250485, 59 1.0294941603786836, 60 1.1334184385711124, 61 1.1568860052088945, 62 1.0412336653823127, 63 0.7660356607147669, 64 0.4068380577304884, 65 0.12130475619620494, 66 0.007647406381002058, 67 0.0038766052508024296, 68 0.0018366503445636489, 69 0.0006757763437823614, 70 0.00014251192022876204, 71 ], 72 [ 73 0.19820963209310652, 74 0.2528661807124882, 75 0.3158570884138644, 76 0.38364913993931354, 77 0.4493832884815257, 78 0.5027490610084232, 79 0.5311511760149372, 80 0.5227660441149257, 81 0.4715289221976368, 82 0.3822917550688269, 83 0.27222944741946864, 84 0.16541480707209608, 85 0.08242031539019055, 86 0.03156306801519871, 87 0.00795647799253679, 88 0.0011153510331640794, 89 0.0004884291395013965, 90 0.00018766025678891934, 91 5.6861581517651256e-5, 92 1.0878733924819186e-5, 93 ], 94 ] 95) 96QC = np.array( 97 [ 98 [ 99 0.00023259428004559948, 100 0.0002804502235646257, 101 0.0003452211545804375, 102 0.0004361877640877965, 103 0.0005702479217508372, 104 0.0007812475748456827, 105 0.0011468941392117338, 106 0.0018884548901117728, 107 0.003952840155165309, 108 0.02356787814538211, 109 0.9902220806396309, 110 0.9944680372872212, 111 0.9948230520394494, 112 0.992828984604994, 113 0.9494349501359681, 114 0.0060141075852780945, 115 0.002391744136969187, 116 0.0013590312401779991, 117 0.0008934288272912649, 118 0.0006376138306461998, 119 ], 120 [ 121 0.0007446008324123797, 122 0.0017062385212967938, 123 0.0045521266894452, 124 0.01232248611337005, 125 0.030897004687650706, 126 0.06802160484300256, 127 0.1270286320846756, 128 0.1975947205304386, 129 0.25774910322755645, 130 0.29241568674609913, 131 0.30261930788409075, 132 0.29164719358476193, 133 0.2508915891345886, 134 0.16612800717742593, 135 0.061328835766217636, 136 0.0016661182121968014, 137 0.0008661265295733964, 138 0.000476358290888631, 139 0.00021705957692457706, 140 5.7921213943590406e-5, 141 ], 142 [ 143 0.01656031383026699, 144 0.028077798198055417, 145 0.042926091821587964, 146 0.059297486773231455, 147 0.07508578878498474, 148 0.08928990048440437, 149 0.1023228797711026, 150 0.11475197203548103, 151 0.12539200432700545, 152 0.13008943253916275, 153 0.12258521422776973, 154 0.09864887346844652, 155 0.06227244046416078, 156 0.027129943467786286, 157 0.006445011737676357, 158 0.000519163514737142, 159 0.0002532092239246288, 160 0.00010820782161964733, 161 3.4068535970606276e-5, 162 5.79220821416827e-6, 163 ], 164 [ 165 0.030757168729973833, 166 0.0373735991829509, 167 0.04434591786357304, 168 0.05090713814014857, 169 0.05590771651108847, 170 0.05797636144416201, 171 0.05590303149164768, 172 0.04921078599960091, 173 0.038700625460403475, 174 0.026521223162384176, 175 0.015391250719866753, 176 0.0073065234193505635, 177 0.00271332219262275, 178 0.0007486848626273913, 179 0.00015588378375704442, 180 3.626367328902605e-5, 181 1.308012805783289e-5, 182 3.771780816369582e-6, 183 7.63771805417599e-7, 184 8.110312094613359e-8, 185 ], 186 ] 187) 188NR = np.array( 189 [ 190 [ 191 2.3259428004559948e-10, 192 2.804502235646257e-10, 193 3.4522115458043745e-10, 194 4.361877640877965e-10, 195 5.702479217508372e-10, 196 7.812475748456827e-10, 197 1.1468941392117337e-9, 198 1.888454890111773e-9, 199 3.952840155165309e-9, 200 2.356787814538211e-8, 201 9.90222080639631e-7, 202 9.944680372872212e-7, 203 9.948230520394495e-7, 204 9.92828984604994e-7, 205 9.494349501359681e-7, 206 6.014107585278094e-9, 207 2.391744136969187e-9, 208 1.359031240177999e-9, 209 8.934288272912649e-10, 210 6.376138306461997e-10, 211 ], 212 [ 213 0.0003533420532791116, 214 0.0012826071194146278, 215 0.00424553213537039, 216 0.012673002046450213, 217 0.03362575699559721, 218 0.07780646668937528, 219 0.15331141994790556, 220 0.2511957333816394, 221 0.33814936627199305, 222 0.3780523095075203, 223 0.3598088317322889, 224 0.2931296678544589, 225 0.19341647440335039, 226 0.08651190873395816, 227 0.01810784482882837, 228 4.1494189890466366e-6, 229 1.094349113436064e-6, 230 3.546676197032568e-7, 231 9.14315623350356e-8, 232 1.2064015228861402e-8, 233 ], 234 [ 235 0.04286183602944602, 236 0.06852467872215309, 237 0.09903215887493476, 238 0.12997070440033048, 239 0.15687232059974845, 240 0.17740978878087613, 241 0.19118547711395129, 242 0.19691911172103171, 243 0.18984923336031656, 244 0.16372127417779328, 245 0.11889273491240387, 246 0.06789317377731022, 247 0.02793297053418389, 248 0.007207574524451838, 249 0.0008713383759906845, 250 7.702005312220904e-7, 251 2.014372171094809e-7, 252 4.725224070296279e-8, 253 7.766049118140948e-9, 254 6.427033141822792e-10, 255 ], 256 [ 257 0.061874092050410494, 258 0.07128188742881585, 259 0.08111795686596689, 260 0.08980451732153531, 261 0.09449050033043897, 262 0.09134074347469225, 263 0.07752987289413114, 264 0.05487127635724207, 265 0.030981982591521077, 266 0.013660394132669304, 267 0.004693803464682259, 268 0.001247946986806364, 269 0.0002471714862400846, 270 3.290890468105509e-5, 271 2.207068943645234e-6, 272 1.0190141478651551e-8, 273 2.047865946736227e-9, 274 3.3636644705214316e-10, 275 3.774223029783484e-11, 276 1.9828529963408933e-12, 277 ], 278 ] 279) 280QR = np.array( 281 [ 282 [ 283 2.3259428004559948e-10, 284 2.804502235646257e-10, 285 3.4522115458043745e-10, 286 4.361877640877965e-10, 287 5.702479217508372e-10, 288 7.812475748456827e-10, 289 1.1468941392117337e-9, 290 1.888454890111773e-9, 291 3.952840155165309e-9, 292 2.356787814538211e-8, 293 9.90222080639631e-7, 294 9.944680372872212e-7, 295 9.948230520394495e-7, 296 9.92828984604994e-7, 297 9.494349501359681e-7, 298 6.014107585278094e-9, 299 2.391744136969187e-9, 300 1.359031240177999e-9, 301 8.934288272912649e-10, 302 6.376138306461997e-10, 303 ], 304 [ 305 0.00035200151617089407, 306 0.0013111967057274174, 307 0.0044786571531358664, 308 0.013892547736570221, 309 0.03863471127263126, 310 0.09460221788307477, 311 0.19892355317632127, 312 0.34815123733343706, 313 0.491830001583502, 314 0.5511897143312499, 315 0.4932820778423064, 316 0.3591619186873109, 317 0.20663782645801065, 318 0.0795659654091024, 319 0.014340018974944115, 320 2.6031739369300225e-6, 321 6.868825051011121e-7, 322 2.2197596523612783e-7, 323 5.69233748241551e-8, 324 7.456832323236176e-9, 325 ], 326 [ 327 0.11482311944110231, 328 0.1953165800079347, 329 0.2975360397035871, 330 0.40314885321904576, 331 0.4844836094841661, 332 0.5176110996740826, 333 0.49495233107947906, 334 0.4267437023807917, 335 0.3323077368114287, 336 0.23095573846997497, 337 0.1386092144185074, 338 0.06764570083287864, 339 0.024471513975755962, 340 0.005687897829495739, 341 0.0006318064836266766, 342 4.7346216121926047e-7, 343 1.2338769749465383e-7, 344 2.8730188420905227e-8, 345 4.6732350569938025e-9, 346 3.823035518651779e-10, 347 ], 348 [ 349 0.2051516890449052, 350 0.17270749883053185, 351 0.14171673193376189, 352 0.11262150188264747, 353 0.08580121859435079, 354 0.06174295102301346, 355 0.04112156435049741, 356 0.02469744659130083, 357 0.012972304357154157, 358 0.005761888080833301, 359 0.002087794572011011, 360 0.0005919792707091766, 361 0.00012376637910333285, 362 1.7087512894278817e-5, 363 1.1819761014876744e-6, 364 5.503422402098125e-9, 365 1.0561195797504775e-9, 366 1.6226072373027708e-10, 367 1.7126710513565556e-11, 368 9.041044132623166e-13, 369 ], 370 ] 371)
Z =
array([0. , 0.15789474, 0.31578947, 0.47368421, 0.63157895,
0.78947368, 0.94736842, 1.10526316, 1.26315789, 1.42105263,
1.57894737, 1.73684211, 1.89473684, 2.05263158, 2.21052632,
2.36842105, 2.52631579, 2.68421053, 2.84210526, 3. ])
NC =
array([[2.32594280e-03, 2.80450224e-03, 3.45221155e-03, 4.36187764e-03,
5.70247922e-03, 7.81247575e-03, 1.14689414e-02, 1.88845489e-02,
3.95284016e-02, 2.35678781e-01, 9.90222081e+00, 9.94468037e+00,
9.94823052e+00, 9.92828985e+00, 9.49434950e+00, 6.01410759e-02,
2.39174414e-02, 1.35903124e-02, 8.93428827e-03, 6.37613831e-03],
[5.32201476e-03, 1.04800223e-02, 2.58891146e-02, 7.02763798e-02,
1.84657834e-01, 4.34142045e-01, 8.70933938e-01, 1.45269809e+00,
2.01659490e+00, 2.40452382e+00, 2.57975497e+00, 2.57856319e+00,
2.37916903e+00, 1.76823088e+00, 7.72103251e-01, 1.99185876e-02,
1.01551460e-02, 5.76443426e-03, 2.83637120e-03, 8.57071556e-04],
[8.66289346e-02, 1.55203117e-01, 2.52281340e-01, 3.71195397e-01,
4.99920823e-01, 6.29874203e-01, 7.60950941e-01, 8.96236532e-01,
1.02949416e+00, 1.13341844e+00, 1.15688601e+00, 1.04123367e+00,
7.66035661e-01, 4.06838058e-01, 1.21304756e-01, 7.64740638e-03,
3.87660525e-03, 1.83665034e-03, 6.75776344e-04, 1.42511920e-04],
[1.98209632e-01, 2.52866181e-01, 3.15857088e-01, 3.83649140e-01,
4.49383288e-01, 5.02749061e-01, 5.31151176e-01, 5.22766044e-01,
4.71528922e-01, 3.82291755e-01, 2.72229447e-01, 1.65414807e-01,
8.24203154e-02, 3.15630680e-02, 7.95647799e-03, 1.11535103e-03,
4.88429140e-04, 1.87660257e-04, 5.68615815e-05, 1.08787339e-05]])
QC =
array([[2.32594280e-04, 2.80450224e-04, 3.45221155e-04, 4.36187764e-04,
5.70247922e-04, 7.81247575e-04, 1.14689414e-03, 1.88845489e-03,
3.95284016e-03, 2.35678781e-02, 9.90222081e-01, 9.94468037e-01,
9.94823052e-01, 9.92828985e-01, 9.49434950e-01, 6.01410759e-03,
2.39174414e-03, 1.35903124e-03, 8.93428827e-04, 6.37613831e-04],
[7.44600832e-04, 1.70623852e-03, 4.55212669e-03, 1.23224861e-02,
3.08970047e-02, 6.80216048e-02, 1.27028632e-01, 1.97594721e-01,
2.57749103e-01, 2.92415687e-01, 3.02619308e-01, 2.91647194e-01,
2.50891589e-01, 1.66128007e-01, 6.13288358e-02, 1.66611821e-03,
8.66126530e-04, 4.76358291e-04, 2.17059577e-04, 5.79212139e-05],
[1.65603138e-02, 2.80777982e-02, 4.29260918e-02, 5.92974868e-02,
7.50857888e-02, 8.92899005e-02, 1.02322880e-01, 1.14751972e-01,
1.25392004e-01, 1.30089433e-01, 1.22585214e-01, 9.86488735e-02,
6.22724405e-02, 2.71299435e-02, 6.44501174e-03, 5.19163515e-04,
2.53209224e-04, 1.08207822e-04, 3.40685360e-05, 5.79220821e-06],
[3.07571687e-02, 3.73735992e-02, 4.43459179e-02, 5.09071381e-02,
5.59077165e-02, 5.79763614e-02, 5.59030315e-02, 4.92107860e-02,
3.87006255e-02, 2.65212232e-02, 1.53912507e-02, 7.30652342e-03,
2.71332219e-03, 7.48684863e-04, 1.55883784e-04, 3.62636733e-05,
1.30801281e-05, 3.77178082e-06, 7.63771805e-07, 8.11031209e-08]])
NR =
array([[2.32594280e-10, 2.80450224e-10, 3.45221155e-10, 4.36187764e-10,
5.70247922e-10, 7.81247575e-10, 1.14689414e-09, 1.88845489e-09,
3.95284016e-09, 2.35678781e-08, 9.90222081e-07, 9.94468037e-07,
9.94823052e-07, 9.92828985e-07, 9.49434950e-07, 6.01410759e-09,
2.39174414e-09, 1.35903124e-09, 8.93428827e-10, 6.37613831e-10],
[3.53342053e-04, 1.28260712e-03, 4.24553214e-03, 1.26730020e-02,
3.36257570e-02, 7.78064667e-02, 1.53311420e-01, 2.51195733e-01,
3.38149366e-01, 3.78052310e-01, 3.59808832e-01, 2.93129668e-01,
1.93416474e-01, 8.65119087e-02, 1.81078448e-02, 4.14941899e-06,
1.09434911e-06, 3.54667620e-07, 9.14315623e-08, 1.20640152e-08],
[4.28618360e-02, 6.85246787e-02, 9.90321589e-02, 1.29970704e-01,
1.56872321e-01, 1.77409789e-01, 1.91185477e-01, 1.96919112e-01,
1.89849233e-01, 1.63721274e-01, 1.18892735e-01, 6.78931738e-02,
2.79329705e-02, 7.20757452e-03, 8.71338376e-04, 7.70200531e-07,
2.01437217e-07, 4.72522407e-08, 7.76604912e-09, 6.42703314e-10],
[6.18740921e-02, 7.12818874e-02, 8.11179569e-02, 8.98045173e-02,
9.44905003e-02, 9.13407435e-02, 7.75298729e-02, 5.48712764e-02,
3.09819826e-02, 1.36603941e-02, 4.69380346e-03, 1.24794699e-03,
2.47171486e-04, 3.29089047e-05, 2.20706894e-06, 1.01901415e-08,
2.04786595e-09, 3.36366447e-10, 3.77422303e-11, 1.98285300e-12]])
QR =
array([[2.32594280e-10, 2.80450224e-10, 3.45221155e-10, 4.36187764e-10,
5.70247922e-10, 7.81247575e-10, 1.14689414e-09, 1.88845489e-09,
3.95284016e-09, 2.35678781e-08, 9.90222081e-07, 9.94468037e-07,
9.94823052e-07, 9.92828985e-07, 9.49434950e-07, 6.01410759e-09,
2.39174414e-09, 1.35903124e-09, 8.93428827e-10, 6.37613831e-10],
[3.52001516e-04, 1.31119671e-03, 4.47865715e-03, 1.38925477e-02,
3.86347113e-02, 9.46022179e-02, 1.98923553e-01, 3.48151237e-01,
4.91830002e-01, 5.51189714e-01, 4.93282078e-01, 3.59161919e-01,
2.06637826e-01, 7.95659654e-02, 1.43400190e-02, 2.60317394e-06,
6.86882505e-07, 2.21975965e-07, 5.69233748e-08, 7.45683232e-09],
[1.14823119e-01, 1.95316580e-01, 2.97536040e-01, 4.03148853e-01,
4.84483609e-01, 5.17611100e-01, 4.94952331e-01, 4.26743702e-01,
3.32307737e-01, 2.30955738e-01, 1.38609214e-01, 6.76457008e-02,
2.44715140e-02, 5.68789783e-03, 6.31806484e-04, 4.73462161e-07,
1.23387697e-07, 2.87301884e-08, 4.67323506e-09, 3.82303552e-10],
[2.05151689e-01, 1.72707499e-01, 1.41716732e-01, 1.12621502e-01,
8.58012186e-02, 6.17429510e-02, 4.11215644e-02, 2.46974466e-02,
1.29723044e-02, 5.76188808e-03, 2.08779457e-03, 5.91979271e-04,
1.23766379e-04, 1.70875129e-05, 1.18197610e-06, 5.50342240e-09,
1.05611958e-09, 1.62260724e-10, 1.71267105e-11, 9.04104413e-13]])