import numpy as np import pylab as pl import scipy.integrate as sp import scipy.interpolate as si import sys import math #This is just to plot the ionization rate results altitude=[0.15, 0.45, 0.75, 1.05, 1.35, 1.65, 1.95, 2.25, 2.55, 2.85, 3.15, 3.45, 3.75, 4.05, 4.35, 4.65, 4.95, 5.25, 5.55, 5.85, 6.15, 6.45, 6.75, 7.05, 7.35, 7.65, 7.95, 8.25, 8.55, 8.85, 9.15, 9.45, 9.75, 10.05, 10.35, 10.65, 10.95, 11.25, 11.55, 11.85, 12.15, 12.45, 12.75, 13.05, 13.35, 13.65, 13.95, 14.25, 14.55, 14.85, 15.15, 15.45, 15.75, 16.05, 16.35, 16.65, 16.95, 17.25, 17.55, 17.85, 18.15, 18.45, 18.75, 19.05, 19.35, 19.65, 19.95, 20.25, 20.55, 20.85, 21.15, 21.45, 21.75, 22.05, 22.35, 22.65, 22.95, 23.25, 23.55, 23.85, 24.15, 24.45, 24.75, 25.05, 25.35, 25.65, 25.95, 26.25, 26.55, 26.85, 27.15, 27.45, 27.75, 28.05, 28.35, 28.65, 28.95, 29.25, 29.55, 29.85, 30.15, 30.45, 30.75, 31.05, 31.35, 31.65, 31.95, 32.25, 32.55, 32.85, 33.15, 33.45, 33.75, 34.05, 34.35, 34.65, 34.95, 35.25, 35.55, 35.85, 36.15, 36.45, 36.75, 37.05, 37.35, 37.65, 37.95, 38.25, 38.55, 38.85, 39.15, 39.45, 39.75, 40.05, 40.35, 40.65, 40.95, 41.25, 41.55, 41.85, 42.15, 42.45, 42.75, 43.05, 43.35, 43.65, 43.95, 44.25, 44.55, 44.85, 45.15, 45.45, 45.75, 46.05, 46.35, 46.65, 46.95, 47.25, 47.55, 47.85, 48.15, 48.45, 48.75, 49.05, 49.35, 49.65, 49.95, 50.25, 50.55, 50.85, 51.15, 51.45, 51.75, 52.05, 52.35, 52.65, 52.95, 53.25, 53.55, 53.85, 54.15, 54.45, 54.75, 55.05, 55.35, 55.65, 55.95, 56.25, 56.55, 56.85, 57.15, 57.45, 57.75, 58.05, 58.35, 58.65, 58.95, 59.25, 59.55, 59.85, 60.15, 60.45, 60.75, 61.05, 61.35, 61.65, 61.95, 62.25, 62.55, 62.85, 63.15, 63.45, 63.75, 64.05, 64.35, 64.65, 64.95, 65.25, 65.55, 65.85, 66.15, 66.45, 66.75, 67.05, 67.35, 67.65, 67.95, 68.25, 68.55, 68.85, 69.15, 69.45, 69.75, 70.05, 70.35, 70.65, 70.95, 71.25, 71.55, 71.85, 72.15, 72.45, 72.75, 73.05, 73.35, 73.65, 73.95, 74.25, 74.55, 74.85, 75.15, 75.45, 75.75, 76.05, 76.35, 76.65, 76.95, 77.25, 77.55, 77.85, 78.15, 78.45, 78.75, 79.05, 79.35, 79.65, 79.95, 80.25, 80.55, 80.85, 81.15, 81.45, 81.75, 82.05, 82.35, 82.65, 82.95, 83.25, 83.55, 83.85, 84.15, 84.45, 84.75, 85.05, 85.35, 85.65, 85.95, 86.25, 86.55, 86.85, 87.15, 87.45, 87.75, 88.05, 88.35, 88.65, 88.95, 89.25, 89.55, 89.85, 90.15, 90.45, 90.75, 91.05, 91.35, 91.65, 91.95, 92.25, 92.55, 92.85, 93.15, 93.45, 93.75, 94.05, 94.35, 94.65, 94.95, 95.25, 95.55, 95.85, 96.15, 96.45, 96.75, 97.05, 97.35, 97.65, 97.95, 98.25, 98.55, 98.85, 99.15, 99.45, 99.75, 100.05, 100.35, 100.65, 100.95, 101.25, 101.55, 101.85, 102.15, 102.45, 102.75, 103.05, 103.35, 103.65, 103.95, 104.25, 104.55, 104.85, 105.15, 105.45, 105.75, 106.05, 106.35, 106.65, 106.95, 107.25, 107.55, 107.85, 108.15, 108.45, 108.75, 109.05, 109.35, 109.65, 109.95, 110.25, 110.55, 110.85, 111.15, 111.45, 111.75, 112.05, 112.35, 112.65, 112.95, 113.25, 113.55, 113.85, 114.15, 114.45, 114.75, 115.05, 115.35, 115.65, 115.95, 116.25, 116.55, 116.85, 117.15, 117.45, 117.75, 118.05, 118.35, 118.65, 118.95, 119.25, 119.55, 119.85, 120.15, 120.45, 120.75, 121.05, 121.35, 121.65, 121.95, 122.25, 122.55, 122.85, 123.15, 123.45, 123.75, 124.05, 124.35, 124.65, 124.95, 125.25, 125.55, 125.85, 126.15, 126.45, 126.75, 127.05, 127.35, 127.65, 127.95, 128.25, 128.55, 128.85, 129.15, 129.45, 129.75, 130.05, 130.35, 130.65, 130.95, 131.25, 131.55, 131.85, 132.15, 132.45, 132.75, 133.05, 133.35, 133.65, 133.95, 134.25, 134.55, 134.85, 135.15, 135.45, 135.75, 136.05, 136.35, 136.65, 136.95, 137.25, 137.55, 137.85, 138.15, 138.45, 138.75, 139.05, 139.35, 139.65, 139.95, 140.25, 140.55, 140.85, 141.15, 141.45, 141.75, 142.05, 142.35, 142.65, 142.95, 143.25, 143.55, 143.85, 144.15, 144.45, 144.75, 145.05, 145.35, 145.65, 145.95, 146.25, 146.55, 146.85, 147.15, 147.45, 147.75, 148.05, 148.35, 148.65, 148.95, 149.25, 149.55, 149.85] #Duration of extreme SEP events like Carrington from Dartnell paper# #The energy deposition profiles from that paper are given in terms of the entire duration of the event (~20 hours) #therefore we divide the total energy deposition by the event duration to get an approximation of a time-averaged energy deposition rate for comparison sepduration=72000 #=20 hours in seconds 20*3600 ######solmax solmax_protons=[0.000738870953930608, 0.000765023969177284, 0.000943540027211705, 0.00103454117995736, 0.00116816952059565, 0.00131848958509299, 0.00145442836036734, 0.00155154640340788, 0.00172343983509566, 0.00190991791928832, 0.00214303385939624, 0.00246290286204431, 0.00267822517896217, 0.00290905259028938, 0.00313933668008722, 0.00332582314907538, 0.00329232894404287, 0.00333206409513652, 0.0033938184069831, 0.00343661196442298, 0.00348862062727063, 0.00366118812868924, 0.00379084745972372, 0.00395968864830963, 0.00413412013907627, 0.00434737871927903, 0.00465774303114273, 0.00505368458193123, 0.00538995229541244, 0.00584464521298909, 0.00615322453669069, 0.00650351890093948, 0.00680549748291309, 0.00717265370909232, 0.00765715409348715, 0.00808782645263579, 0.00849009082785298, 0.00905617543335359, 0.00952813780893442, 0.00986573443787794, 0.0105846903209452, 0.0115436205564248, 0.0124300693305263, 0.0132735001389618, 0.0143158350484735, 0.0154288864509121, 0.0164255279461838, 0.0173482281094457, 0.0182520077828451, 0.0190821528447273, 0.019975102862374, 0.0205972619874851, 0.0216153891567157, 0.0238338182294683, 0.0252091376820048, 0.0269814654482321, 0.0292589185837113, 0.030892124458142, 0.0316737359977783, 0.0340769861920228, 0.035915073963338, 0.0377606304476836, 0.0398375316515905, 0.0416531424162555, 0.0430473895477147, 0.0442134463065147, 0.0460917344712146, 0.0477771081919003, 0.0501020265280786, 0.0522210414319188, 0.0545546094282249, 0.0563898892277383, 0.0588566586980159, 0.0612179569878436, 0.0640148725546346, 0.0672371976922123, 0.0703530207872816, 0.0737272730450813, 0.0768144872339965, 0.0810932754620581, 0.0850413309186392, 0.0889644851149316, 0.0930154707493342, 0.0976469893974552, 0.10074770761751, 0.105862385260499, 0.110123166623688, 0.114009075865737, 0.118132320366499, 0.122267361816056, 0.125130098828886, 0.13010521981427, 0.135100124305128, 0.140666420518903, 0.146793959284439, 0.153528677654247, 0.161016605274222, 0.169721161294499, 0.177921749731814, 0.186106595697018, 0.19387481170453, 0.201741306692935, 0.208617834485494, 0.215699178301201, 0.223475830749465, 0.23213445466658, 0.240172624458046, 0.248299862061865, 0.25839105417055, 0.269370205763932, 0.279945517596338, 0.290878962921047, 0.30345813694388, 0.314539049245983, 0.325008265130359, 0.336977586471382, 0.348182043188148, 0.35831823310562, 0.369336536379756, 0.381412850458502, 0.392598593416012, 0.406235669468765, 0.420996207735619, 0.436606932549193, 0.451788921547465, 0.468577466485015, 0.483888042935833, 0.499651554908505, 0.515381026345506, 0.531701581493449, 0.550344094341071, 0.567697345722421, 0.585904178125296, 0.604367129425554, 0.624185204169274, 0.642964474431782, 0.661606621458066, 0.681451352990675, 0.703374567317495, 0.72472284961479, 0.744433644465199, 0.766509728837402, 0.787153495290528, 0.806872109132023, 0.825887553574338, 0.843209676346685, 0.865407327714449, 0.886611441511918, 0.906315842250796, 0.928289429218658, 0.953193574107926, 0.97380329399731, 0.998023956588412, 1.02552257161507, 1.05723595799028, 1.09425493789437, 1.13692765623587, 1.1819864837396, 1.22915645698416, 1.27183441675292, 1.3132610743142, 1.35279697125262, 1.39278733037382, 1.43851565375338, 1.49677943126475, 1.55887618185574, 1.63567552982332, 1.72209391424397, 1.83457745542925, 1.95357099403756, 2.09683859255152, 2.27119767477093, 2.45144091351398, 2.66918283934078, 2.92883452862924, 3.25097236181134, 3.59671774200491, 4.03033762550808, 4.48354016393092, 5.03263373362543, 5.65124599602438, 6.34825743303165, 7.08619614892104, 7.91904610998351, 8.80204559454261, 9.77190534903751, 10.8195145268718, 12.0793647196959, 13.4255350428126, 14.8431142337067, 16.3430636322848, 17.8530032502428, 19.2382551501173, 20.621304896453, 21.9717662126965, 23.2832555435986, 24.6433101540461, 26.0904776995051, 27.5149741960075, 28.8735664068281, 30.1319328552052, 31.3050644256876, 32.2507883336318, 32.9300533826836, 33.4185393383268, 33.727521345397, 33.8351904994875, 33.8918320773321, 33.8373741965093, 33.695904974513, 33.4329076377289, 32.9404457081432, 32.2647983176071, 31.6305056468409, 30.8819281965157, 30.0906925097377, 29.3000663325875, 28.378501290218, 27.3269162278175, 26.3485721705446, 25.2787484859996, 24.3339698973375, 23.4223690923757, 22.4628048297421, 21.4309625508638, 20.4531334726019, 19.414313120136, 18.394363680222, 17.4557943580271, 16.5605903277935, 15.6906430280023, 14.8248984859909, 14.0181587810922, 13.2652628814537, 12.527621368597, 11.8221724406142, 11.1402939931047, 10.5001989700489, 9.84001836574513, 9.23959864714747, 8.67121330191289, 8.12854303244544, 7.61260467558334, 7.15925696354767, 6.71212666222119, 6.28634713127138, 5.90298970164804, 5.54081536391587, 5.2117298295683, 4.88631691960866, 4.57768432898433, 4.28868339866747, 4.01237079852719, 3.74011907336656, 3.50411967370684, 3.29051966547619, 3.08142352567955, 2.87561280020022, 2.68543198234259, 2.50822671663981, 2.32302815435144, 2.16043554574136, 2.02327685906093, 1.88457713441172, 1.75689743401945, 1.64761800319523, 1.53968589878337, 1.43347409225185, 1.34122552832175, 1.25026121047894, 1.16034242973773, 1.07106967310446, 0.990717245813113, 0.919246531083721, 0.856317257145057, 0.804273762323153, 0.755469452255038, 0.707618627351862, 0.656529324442127, 0.60957070756951, 0.561769486758795, 0.521568696399102, 0.484084895374626, 0.452883104620091, 0.418086331019061, 0.389792592410859, 0.360097627410792, 0.333538805094973, 0.308104011818179, 0.286483197543835, 0.265149067690194, 0.244880321780559, 0.226026439724841, 0.209339260469531, 0.192836555130393, 0.177771621195737, 0.164196697848192, 0.151610236127161, 0.139027123327581, 0.127828883174487, 0.118424460939353, 0.109421012602765, 0.101184624723068, 0.0936237591751, 0.0859526604301489, 0.0779297335146196, 0.0710550516267227, 0.0651756253774091, 0.060225325138915, 0.0558498454995947, 0.051358178032016, 0.0476790468834663, 0.0441403951664479, 0.0404557230630913, 0.0374740227172844, 0.0348671979122506, 0.0319419418815747, 0.0291672289782067, 0.0267199264343578, 0.0244047351209679, 0.0224479194300777, 0.0208119996555331, 0.0194843403814661, 0.0179578633299566, 0.0165742661575621, 0.0152529506863094, 0.0139828153278465, 0.0128234103923706, 0.0122562881169922, 0.0112911645244336, 0.0103284393464149, 0.00953533403781125, 0.00875405502863426, 0.00770722001226549, 0.00717622785630395, 0.00683138105020875, 0.00661573141185991, 0.00629664315542521, 0.00595934173899367, 0.0056424312929548, 0.00530448193483031, 0.00486056339414519, 0.00451454299864182, 0.00417193842765477, 0.00382283286508264, 0.00350083281049278, 0.00319518615689888, 0.002877338646414, 0.00258212882387362, 0.00229163091931222, 0.0020643131259124, 0.00189695306109654, 0.00178866760660943, 0.00169699669334488, 0.0016128897500057, 0.00152561512198148, 0.00143820897245462, 0.00135019483836009, 0.00125905641082304, 0.00116325153772513, 0.00107363975731981, 0.00101014411661349, 0.000949938241606227, 0.000860059320732006, 0.00077821251819875, 0.000706031509714641, 0.000619866269833311, 0.000538204380550527, 0.000502735601018009, 0.000474522716174884, 0.000449421151760153, 0.000426318798664625, 0.000403348403650839, 0.000381014980036107, 0.000358494189635663, 0.000336828848672062, 0.000323784891311383, 0.000309830689357953, 0.000288646393067129, 0.000269340662397789, 0.000244296672671823, 0.000209296339230207, 0.00018162741119501, 0.000162538352661669, 0.000145084156932719, 0.000136302008321617, 0.000130458840246201, 0.000123821814700506, 0.000118098004259087, 0.000112070321269493, 0.000105604763070809, 9.89123960206331e-05, 9.22590569691477e-05, 8.52336346994032e-05, 7.83011079293013e-05, 7.12678281814168e-05, 6.46758829226059e-05, 5.8182967324332e-05, 5.27005706600168e-05, 4.89431744809887e-05, 4.58068527398541e-05, 4.25478847191458e-05, 4.01154777609986e-05, 3.78839947479475e-05, 3.47915388860144e-05, 3.20019159996345e-05, 3.00450003857527e-05, 2.8079875188598e-05, 2.6178333986288e-05, 2.44483096814196e-05, 2.26301311310914e-05, 2.10557791611742e-05, 1.95741242370865e-05, 1.79968547925418e-05, 1.62229890094985e-05, 1.46030379072779e-05, 1.32288724768896e-05, 1.21212256044881e-05, 1.10935482342463e-05, 1.05022795923782e-05, 1.00102050614514e-05, 9.49969505805254e-06, 8.97433724794653e-06, 8.65109029279609e-06, 8.24189310669699e-06, 7.8582137578216e-06, 7.3848049011242e-06, 6.88124836414291e-06, 6.29767849432073e-06, 5.79942846784382e-06, 5.2560026236726e-06, 4.78561427662759e-06, 4.34934822854941e-06, 4.01599976841783e-06, 3.71173453109183e-06, 3.48660876174115e-06, 3.37216922046732e-06, 3.23289630990097e-06, 3.06969009537525e-06, 2.95096285207264e-06, 2.77228789720792e-06, 2.58527605041159e-06, 2.40468166589484e-06, 2.24841346034789e-06, 2.05517305592678e-06, 1.92188274468332e-06, 1.77173395440739e-06, 1.63546326895815e-06, 1.50753226670251e-06, 1.39798914211338e-06, 1.29761279133444e-06, 1.23238769996175e-06, 1.17799233871948e-06, 1.13617206494067e-06, 1.07955545528086e-06, 1.02112083464194e-06, 9.67000450045424e-07, 9.23629822104933e-07, 8.78325414259297e-07, 8.33145911389671e-07, 7.78464225409577e-07, 7.19532390634518e-07, 6.61995667003987e-07, 5.95153104422119e-07, 5.4460074601601e-07, 5.05735951094589e-07, 4.72735912946739e-07, 4.37623941435669e-07, 4.10275273038267e-07, 3.91885366404755e-07, 3.76038400254941e-07, 3.56022132088885e-07, 3.39979957215658e-07, 3.28128533743393e-07, 3.13470340616856e-07, 2.96763885659749e-07, 2.83455431004346e-07, 2.7162479837708e-07, 2.53882466525665e-07, 2.36522378458901e-07, 2.18409018159607e-07, 2.01283181162278e-07, 1.85730392467918e-07, 1.74648934690584e-07, 1.6465236539161e-07, 1.56016933541038e-07, 1.49047945444748e-07, 1.41826095450965e-07, 1.34835636405156e-07, 1.2843214888531e-07, 1.24069140096426e-07, 1.19710150838772e-07, 1.15008438250115e-07, 1.09861115202908e-07, 1.04578884706168e-07, 9.78742394077617e-08, 9.13072157925449e-08, 8.51925229177522e-08, 7.81155344209624e-08, 7.17957088056459e-08] solmax_HZE=[0.000214564255932794, 0.0002220707051897, 0.000273814994353869, 0.000300151006539606, 0.000338926989041028, 0.000382588464944393, 0.000422106918021267, 0.000450354020477539, 0.000500314408098908, 0.000554729554730863, 0.000622657225113376, 0.000715857127703108, 0.000778668979787929, 0.000845947401752176, 0.000913153654473259, 0.000967034820562101, 0.000957235660237084, 0.000968624753464856, 0.000986510357341186, 0.00099859302596883, 0.00101376729471048, 0.00106387381648718, 0.00110164694994432, 0.00115062481890885, 0.00120171603163767, 0.00126460690512873, 0.00135536905193233, 0.00147109976200382, 0.0015697457068759, 0.00170333275048959, 0.00179331947399527, 0.00189538159852076, 0.00198363118555431, 0.00209089863015865, 0.00223234669782025, 0.00235848929920634, 0.00247629666410585, 0.00264173843221103, 0.00277929464010793, 0.00287665603530738, 0.00308554757630749, 0.00336609107912791, 0.00362593136591712, 0.00387385989750139, 0.00418011871559169, 0.00450927382351706, 0.00480281271266942, 0.00507407677387648, 0.00533954889350182, 0.00558504236282106, 0.00584663063446837, 0.00602879870739444, 0.00632697501655301, 0.0069791531388927, 0.00738236830235405, 0.00790294468326847, 0.00857234212765228, 0.00905459346827126, 0.00928726555043449, 0.00999916568943308, 0.0105446918594748, 0.0110949105667061, 0.0117120623040772, 0.0122510624610247, 0.0126621303194278, 0.0130084766095763, 0.0135641497733831, 0.0140630701770245, 0.0147500251210206, 0.0153785999111824, 0.0160698225846111, 0.016616261480822, 0.0173475235283072, 0.0180513242394954, 0.0188852017925403, 0.0198477094993357, 0.0207771863229257, 0.0217891293412386, 0.0227107333958847, 0.0239806295436972, 0.0251542046876482, 0.0263241796492651, 0.0275299932048434, 0.0289164500790112, 0.029849528073459, 0.0313811008442041, 0.0326548224261478, 0.0338134566914465, 0.0350375464243369, 0.0362734605027679, 0.0371314408919471, 0.0386227444194149, 0.0401271601945065, 0.0418061081385889, 0.0436621603701749, 0.0457008034927899, 0.0479673712179364, 0.0505986321907098, 0.0530876660908492, 0.0555670314122683, 0.0579189811385097, 0.0603047138069905, 0.0623976656561193, 0.0645415409832168, 0.0668984647751351, 0.0695354957747661, 0.0719849561443071, 0.0744608517354263, 0.0775488005103308, 0.0808967614888685, 0.0841265010819424, 0.0874896128404953, 0.0913441107640725, 0.0947344100153772, 0.0979485678826733, 0.101628741993569, 0.105071599035594, 0.108200954441364, 0.111600136359413, 0.115325595627351, 0.118777458680764, 0.122986448845799, 0.127559713887575, 0.132413180303758, 0.137134729904256, 0.142350882618348, 0.147111759886753, 0.152005065808021, 0.156888441457846, 0.161969104420555, 0.167795657512085, 0.173242435131868, 0.178952598698966, 0.184781374203876, 0.191039097070625, 0.196982691344441, 0.202896604628956, 0.209208545904343, 0.216170627143299, 0.222926719284325, 0.229210153998219, 0.236237693814977, 0.242826315061715, 0.249154157058154, 0.25532385560306, 0.260979678280019, 0.268184322425606, 0.275114844740386, 0.281641267979029, 0.28892381533423, 0.297154312134548, 0.304118767110488, 0.312167170097925, 0.32095384746445, 0.331230780421701, 0.343247110700097, 0.357091593339281, 0.371945982383512, 0.387140055113769, 0.401434274270357, 0.414418523363584, 0.426395118875299, 0.438378197469061, 0.452308936837533, 0.469876030165625, 0.490098266304199, 0.513716296134774, 0.541188677297003, 0.575030817279154, 0.610582398824034, 0.654135921035606, 0.705003785766563, 0.763700229504396, 0.833755791746967, 0.91620394030163, 1.01028761028309, 1.11314929342784, 1.23694201335528, 1.375808099992, 1.54059570193392, 1.73186598414202, 1.95631677360911, 2.2025461155964, 2.47670370516881, 2.77642652148377, 3.10746783665881, 3.47082102537375, 3.88814458697799, 4.33725021625334, 4.81557260838837, 5.32511519359275, 5.84736246678605, 6.35006642839091, 6.85316618226151, 7.35669681162956, 7.84817975555576, 8.35016164045621, 8.86740888168873, 9.38744398260142, 9.87486570887157, 10.3385489713753, 10.7517308348966, 11.1043124845736, 11.3856724193057, 11.6073049788113, 11.7778121428745, 11.8996677557341, 11.9852461454335, 12.0250628817017, 12.0351604082316, 11.9755526727126, 11.8668096232981, 11.7248801152153, 11.5699795014301, 11.3954544686072, 11.2068958700348, 10.9932322963517, 10.7287609163555, 10.4288683342404, 10.110297730885, 9.78797240420205, 9.47268707213908, 9.16123305277203, 8.84640488962329, 8.51521522096598, 8.19112985735204, 7.86568817992104, 7.55433576937178, 7.25347055779616, 6.96665453278117, 6.66762564773189, 6.36802239528925, 6.07725513609327, 5.80437241665307, 5.55105042514715, 5.31801064117671, 5.09600224007057, 4.87372865510356, 4.64898727103917, 4.43395783754379, 4.22989893728414, 4.04141115276635, 3.86584279889282, 3.69463440842486, 3.514123551755, 3.34314700953154, 3.18203316288857, 3.03185129324757, 2.89773524419241, 2.77654126611823, 2.65267625650541, 2.52429579476962, 2.40341276381217, 2.28749262956294, 2.17980423804366, 2.07676487816731, 1.97882027462626, 1.87648698374635, 1.77600033318522, 1.68218009449418, 1.59893341551134, 1.521776457523, 1.45095243066423, 1.38114126695102, 1.30946672753721, 1.23953050281748, 1.17444742709508, 1.11268508441549, 1.05515569016496, 1.00184960448687, 0.947830573030427, 0.892802255123724, 0.84067771175933, 0.791822150268154, 0.746175264735906, 0.703558262789327, 0.663723448869926, 0.62337412411559, 0.583013722555305, 0.545271303593673, 0.511214766018135, 0.479571631089385, 0.451155402327268, 0.424209168804345, 0.396822227307038, 0.370132048874049, 0.3457938156615, 0.322927248419579, 0.301994681470983, 0.283077295562266, 0.264677339545629, 0.246689212903896, 0.229948111973999, 0.214507114113468, 0.200600161953302, 0.187918309865671, 0.175606775668757, 0.163116265575339, 0.151067809676654, 0.139682409988993, 0.129341979019258, 0.12010493483722, 0.111861110319553, 0.104177260603076, 0.0964175526822602, 0.0891177974163523, 0.082271479516295, 0.0760071312458209, 0.0702844089604763, 0.0652722164551272, 0.0604334725900288, 0.055607095139263, 0.0511783494430446, 0.047095790276968, 0.0434090490629183, 0.0400869088087855, 0.0371215555516053, 0.0341422401808228, 0.031326545968967, 0.0287682730001068, 0.0264785997769017, 0.0244303439622615, 0.0226499439798286, 0.0209625463824721, 0.0192841357586932, 0.017711832703974, 0.0162951479508684, 0.0150190253074125, 0.0138813923162482, 0.0128608937507815, 0.0118667969686093, 0.0109129669948362, 0.0100664942142141, 0.00932367383400969, 0.00871127900121423, 0.0082107299994914, 0.0078010740945748, 0.00739515322048558, 0.00700804453256367, 0.00662517068233032, 0.0062383577552663, 0.00584854093686225, 0.00545666023058241, 0.00506299244897686, 0.00466409135640402, 0.00427283502947, 0.00388842108057329, 0.00350400199685103, 0.00314009820515466, 0.00283583278284133, 0.00258093850970941, 0.00238576216236799, 0.00224540293223715, 0.00213336720602003, 0.00202237170822075, 0.00191231836977105, 0.00179713284387105, 0.00168423713638481, 0.00157498265037223, 0.00146383484537263, 0.00135798493111508, 0.00124949813327058, 0.00113982291308973, 0.00102854373794697, 0.000918048201298466, 0.000819708492296997, 0.000739826367387055, 0.00067624148677433, 0.000626625680269098, 0.000592818819950286, 0.000561766584611583, 0.000531545081024479, 0.000500758691598802, 0.000472417120846119, 0.000444065281504103, 0.000414632528972163, 0.000385246091938591, 0.000356229530207621, 0.000326143879957471, 0.000296433079404939, 0.000267625264275524, 0.00023986530982189, 0.000216061845486557, 0.000196114313160602, 0.00018049929191948, 0.000169013338748562, 0.000160924286961568, 0.000152947492776494, 0.00014526347825435, 0.000137425583761948, 0.000129617424640564, 0.000121447548149585, 0.000113337027071918, 0.000104898793882698, 9.64182958597043e-05, 8.79283259310713e-05, 7.99323123843179e-05, 7.17972424997436e-05, 6.45747452907533e-05, 5.88353721916549e-05, 5.42453975535872e-05, 5.02742957354414e-05, 4.74718365348928e-05, 4.51246479020014e-05, 4.24997034421664e-05, 3.99774298273005e-05, 3.77301864384953e-05, 3.55225878657388e-05, 3.33758357580065e-05, 3.12183358835796e-05, 2.90102918234942e-05, 2.68363765121989e-05, 2.46462604866906e-05, 2.24211098330509e-05, 2.02317642346696e-05, 1.82716133527127e-05, 1.6656323974723e-05, 1.53518113283201e-05, 1.43488039824256e-05, 1.36366971952321e-05, 1.29870519482241e-05, 1.23290355183664e-05, 1.16648448516911e-05, 1.10168042955135e-05, 1.03731248017034e-05, 9.74625368665129e-06, 9.10584635719118e-06, 8.48220192408841e-06, 7.86391469139723e-06, 7.23783040435895e-06, 6.61204892582105e-06, 6.0079843405194e-06, 5.48997309802222e-06, 5.0487666146117e-06, 4.69970394640656e-06, 4.42978422432427e-06, 4.23063431952162e-06, 4.02827536171225e-06, 3.82622204564261e-06, 3.6299207932452e-06, 3.4339036566908e-06, 3.24354427208748e-06, 3.0473752949235e-06, 2.85224263184283e-06, 2.65595225823676e-06, 2.46063630322729e-06, 2.26418390591658e-06, 2.07657147114531e-06, 1.90000225998774e-06, 1.75340693339833e-06, 1.62683262857224e-06, 1.52480239513731e-06, 1.44668977665807e-06, 1.38503270031062e-06, 1.32067995313356e-06, 1.26109715101993e-06, 1.19912857127575e-06, 1.13687033844244e-06, 1.07290734637358e-06, 1.00680831681445e-06, 9.42304707657255e-07, 8.77710455636046e-07, 8.14521099306392e-07, 7.53502897705908e-07, 6.93977089388061e-07, 6.42074488033448e-07, 5.99638558120732e-07, 5.62652079110848e-07, 5.32305094376071e-07, 5.09143402385589e-07, 4.85367282720327e-07, 4.61664252330085e-07, 4.39683592104784e-07, 4.18121239732556e-07, 3.96915040011308e-07, 3.76003576065474e-07, 3.5503231124795e-07, 3.34770110503113e-07, 3.14673705230888e-07, 2.93823791035766e-07, 2.73092991270612e-07, 2.52992343489066e-07, 2.35549192320967e-07, 2.1940850010675e-07, 2.06634543574926e-07, 1.96666309610498e-07, 1.88469169393843e-07, 1.79727048119213e-07, 1.72090426772031e-07, 1.64133484801864e-07, 1.56147303295504e-07, 1.49518435159347e-07, 1.42590494090564e-07, 1.35254088144542e-07, 1.27954861316201e-07, 1.20729879965343e-07, 1.12447582072019e-07, 1.02597696706178e-07, 9.35658951335981e-08, 7.92665071492615e-08] solmax_total = [solmax_protons[i] + solmax_HZE[i] for i in range(len(solmax_HZE))] solmax_total_edep=[(solmax_total[i]*33.5*sepduration) for i in range(len(solmax_total)) ] #Scaled up for duration of SEP event #####Solmin solmin_protons=[0.000738870953930608, 0.000765023969177284, 0.000943540027211705, 0.00103454117995736, 0.00116816952059565, 0.00131848958509299, 0.00145442836036734, 0.00155154640340788, 0.00172343983509566, 0.00190991791928832, 0.00214303385939624, 0.00246290286204431, 0.00267822517896217, 0.00290905259028938, 0.00313933668008722, 0.00332582314907538, 0.00329232894404287, 0.00333206409513652, 0.0033938184069831, 0.00343661196442298, 0.00348862062727063, 0.00366118812868924, 0.00379084745972372, 0.00395968864830963, 0.00413412013907627, 0.00434737871927903, 0.00465774303114273, 0.00505368458193123, 0.00538995229541244, 0.00584464521298909, 0.00615322453669069, 0.00650351890093948, 0.00680549748291309, 0.00717265370909232, 0.00765715409348715, 0.00808782645263579, 0.00849009082785298, 0.00905617543335359, 0.00952813780893442, 0.00986573443787794, 0.0105846903209452, 0.0115436205564248, 0.0124300693305263, 0.0132735001389618, 0.0143158350484735, 0.0154288864509121, 0.0164255279461838, 0.0173482281094457, 0.0182520077828451, 0.0190821528447273, 0.019975102862374, 0.0205972619874851, 0.0216153891567157, 0.0238338182294683, 0.0252091376820048, 0.0269814654482321, 0.0292589185837113, 0.030892124458142, 0.0316737359977783, 0.0340769861920228, 0.035915073963338, 0.0377606304476836, 0.0398375316515905, 0.0416531424162555, 0.0430473895477147, 0.0442134463065147, 0.0460917344712146, 0.0477771081919003, 0.0501020265280786, 0.0522210414319188, 0.0545546094282249, 0.0563898892277383, 0.0588566586980159, 0.0612183291609871, 0.0640168798578654, 0.0672418605688722, 0.0703613292117788, 0.0737396853124848, 0.0768309051656451, 0.0811126198413161, 0.0850637983064337, 0.0889904901569319, 0.0930446445453892, 0.0976816537487363, 0.1007873030103, 0.105907928178112, 0.110170989928044, 0.114060201331841, 0.118186869913839, 0.122327669584048, 0.125196563304401, 0.130183456203388, 0.135194238680255, 0.140777373732367, 0.146937240593671, 0.153709668894179, 0.161238565589738, 0.169988111495465, 0.17823870101652, 0.186466945696297, 0.19427495223094, 0.202195872453714, 0.20912349435542, 0.216249126188294, 0.224073924029652, 0.232799042851758, 0.240899839002133, 0.249098486523308, 0.259291269365573, 0.270366129540021, 0.28104592388327, 0.292117764554896, 0.304827592549278, 0.316018524225403, 0.326607039838188, 0.338723228914909, 0.350060473309104, 0.360337888789129, 0.371522084600595, 0.383782008489603, 0.395136188922913, 0.408973272546987, 0.423987813176736, 0.439889444678595, 0.455363961096748, 0.472475095048393, 0.488080243906372, 0.504123924361435, 0.520152918548971, 0.536775672831642, 0.555803658346469, 0.573587796464715, 0.59225051876146, 0.611221201019791, 0.631612530395964, 0.650977092831653, 0.670217855365064, 0.690726687546449, 0.713372621368572, 0.735396520752447, 0.755820018137032, 0.778665857379274, 0.800050013694365, 0.820550923632567, 0.84041817714232, 0.85858475073616, 0.881770554711703, 0.904028224830576, 0.924865186832601, 0.948086185920319, 0.974290433802909, 0.996139278538405, 1.02184947877563, 1.05119685142873, 1.08507545377861, 1.12473198919536, 1.17039763743399, 1.21869712017127, 1.26950971102997, 1.31563711946438, 1.3609460123122, 1.40524781882143, 1.44984837842166, 1.50088895042205, 1.56611111823401, 1.63550649333152, 1.72139162882094, 1.81960321201684, 1.94900683332745, 2.0860508096717, 2.25217561593034, 2.45685861412438, 2.66503302166405, 2.91778197274027, 3.22002968785052, 3.60262289215101, 4.0179328968104, 4.54273213006809, 5.08950817280919, 5.75456938463732, 6.50570356217659, 7.34644321479446, 8.23947731138519, 9.25415600490668, 10.3413364257214, 11.5268726851381, 12.8082709111333, 14.3635742917563, 16.0352725703196, 17.799593268642, 19.6946020124891, 21.5960289601186, 23.3376022864878, 25.1069053755574, 26.8361208050039, 28.5281687998903, 30.3215029795529, 32.2561046116057, 34.175012290869, 36.0322073410051, 37.7647029830572, 39.4425205603601, 40.8329550193581, 41.8786204644497, 42.7032548770274, 43.3236385885393, 43.6567752076515, 43.9547856926043, 44.1277986896507, 44.1788512916315, 44.1092761699131, 43.7170627029989, 43.0619899129875, 42.4579566884785, 41.7122881399662, 40.8658870150692, 40.0371605013809, 39.0400459717674, 37.8522775133792, 36.7738060336672, 35.5276853287635, 34.4703076174892, 33.4405056352563, 32.3065806134407, 31.0476229583642, 29.8839950921532, 28.5896594776625, 27.2921340931806, 26.1161015706986, 24.9711913551375, 23.8388134023276, 22.7013810831054, 21.6479596316546, 20.6587766449046, 19.6688157940895, 18.7132535827201, 17.7805606691363, 16.8810934005667, 15.9357377337329, 15.0801335438805, 14.2669690939224, 13.4695484091884, 12.7249314602738, 12.0572031831496, 11.4004931093872, 10.7621411220382, 10.190098274578, 9.64053105730094, 9.15214919929088, 8.64373442392589, 8.15976134161436, 7.70698424913202, 7.26120540801157, 6.81615547693825, 6.43038685180625, 6.08679778460798, 5.74667065999708, 5.40793331553528, 5.07985922002208, 4.7791289924093, 4.45227675430714, 4.16077523743324, 3.91917056869602, 3.67479898322262, 3.44190320734681, 3.24314962992878, 3.04452485715231, 2.84769364874931, 2.67889279294671, 2.51223347290313, 2.34373746954369, 2.17398507872848, 2.01784014828678, 1.88173083447269, 1.76017194097133, 1.65947325641873, 1.5676657293027, 1.47397588791733, 1.37260726929286, 1.27765613669077, 1.18160441294965, 1.09932871226533, 1.02259659649051, 0.957445419009663, 0.888246300970526, 0.83262702626703, 0.772210655575141, 0.718497983320404, 0.6668739134943, 0.622568904825272, 0.577179756260747, 0.53401991999239, 0.494242357106675, 0.458756318958953, 0.423199121337238, 0.390464521350072, 0.362910446491791, 0.335757111730371, 0.308714793757138, 0.28435303329454, 0.263929494942614, 0.243309043161813, 0.225750426708203, 0.209181946399937, 0.192409302613702, 0.174953988330323, 0.159604723759727, 0.146027928997603, 0.13493790583772, 0.125113868133133, 0.115128441656402, 0.106645913764915, 0.0987164522057611, 0.0904574319750014, 0.0835942069032484, 0.0775600549816853, 0.0712677787576247, 0.0651386532547404, 0.0596440954617204, 0.0545678048690872, 0.0502372865447619, 0.0464977019739079, 0.0433016327015682, 0.0399721364905899, 0.0370405421647394, 0.0341641753900992, 0.0314502950836284, 0.0289506134904131, 0.0275031725613075, 0.0252932024072047, 0.0231431968840082, 0.0213053635535368, 0.0195955276856886, 0.0173621057577608, 0.0160667760293763, 0.0152854054400791, 0.0147833845888734, 0.0140405285304543, 0.0133101375154718, 0.0126022097750195, 0.0118241479500197, 0.010861894559113, 0.0100971477898177, 0.00934061343634902, 0.00859419038298183, 0.00789769460210723, 0.0072014186997908, 0.00649583180301069, 0.00582889528993329, 0.00519775174030715, 0.00467715712836004, 0.00430380977914137, 0.00405447899721075, 0.00384949152335654, 0.00364935313457827, 0.00345101215256895, 0.00325814076862378, 0.00305993854265449, 0.00285638605601121, 0.00264236488782665, 0.0024370865905654, 0.00228653650938084, 0.00217814031189444, 0.001974051552729, 0.00178603755961174, 0.00161920955017852, 0.00142709715439167, 0.00121402710517023, 0.00113231426904221, 0.00106962040878078, 0.00101473980104183, 0.000961368283839604, 0.000909838883916712, 0.000857220943876255, 0.000805694769221443, 0.000754934071269701, 0.000730891551525626, 0.000700867322662193, 0.000652182750391902, 0.000606539938794777, 0.000551663841107625, 0.000468740651796324, 0.000402962391926897, 0.000362324742867159, 0.000326368140781801, 0.000305925033064123, 0.000292724813455532, 0.000278147524162242, 0.000264801980162397, 0.000251116561312644, 0.000236909556455757, 0.00022200064854465, 0.000207368965390498, 0.000191747505847606, 0.000176423446219692, 0.000161119541127311, 0.00014627767943725, 0.000131367210202912, 0.000118651136686966, 0.000108689927206961, 0.000100333982453062, 9.3067142193718e-05, 8.76139207850027e-05, 8.28901407462176e-05, 7.71992682312244e-05, 7.20236743864417e-05, 6.75267280728305e-05, 6.30806692582894e-05, 5.87146381216723e-05, 5.47421367554424e-05, 5.06024060874521e-05, 4.70614920439901e-05, 4.37700705682005e-05, 4.03131704361327e-05, 3.64560422257333e-05, 3.28956093787337e-05, 2.98140482770953e-05, 2.72814547932125e-05, 2.4903086125119e-05, 2.3564552939431e-05, 2.24505466379758e-05, 2.13337181935045e-05, 2.01399700647912e-05, 1.94201426142894e-05, 1.84849841921969e-05, 1.75155807168804e-05, 1.64126721603811e-05, 1.53059526661431e-05, 1.39912618037152e-05, 1.28480757353409e-05, 1.16447671104683e-05, 1.0605278082001e-05, 9.62737921318869e-06, 8.87890041497823e-06, 8.17658435444719e-06, 7.70656691464966e-06, 7.43199045292135e-06, 7.13741584784386e-06, 6.80095409739482e-06, 6.5484008410707e-06, 6.1618710164366e-06, 5.76762755480636e-06, 5.37695518133402e-06, 5.0221980563324e-06, 4.61083114385196e-06, 4.30872766596869e-06, 3.97164636397431e-06, 3.64658221056118e-06, 3.36203159276435e-06, 3.11341062695511e-06, 2.89832984307123e-06, 2.7480336901913e-06, 2.62485558319928e-06, 2.51722473135376e-06, 2.3870865955195e-06, 2.26227261233117e-06, 2.15067529441746e-06, 2.06926707935153e-06, 1.97897329937145e-06, 1.87949477857664e-06, 1.75676016157116e-06, 1.62964371452627e-06, 1.49608623619376e-06, 1.34642104794595e-06, 1.23164706583833e-06, 1.13703024498922e-06, 1.05159267159062e-06, 9.68483358649007e-07, 9.08074872616465e-07, 8.67372627289979e-07, 8.32207167527526e-07, 7.9170239045064e-07, 7.56941634287617e-07, 7.28325987707865e-07, 6.94421590972681e-07, 6.57393623707012e-07, 6.26246539629973e-07, 5.98612739721144e-07, 5.60001964663384e-07, 5.21268916672079e-07, 4.81505143750853e-07, 4.43833762789777e-07, 4.09895846596259e-07, 3.85436156422762e-07, 3.63576174275145e-07, 3.45017438810014e-07, 3.30061797809957e-07, 3.14715462319996e-07, 3.00235590175951e-07, 2.8673859221165e-07, 2.7665224763859e-07, 2.66644914250457e-07, 2.55481634018733e-07, 2.43634458659994e-07, 2.31808243149602e-07, 2.17034026184242e-07, 2.0251775965073e-07, 1.89201095312474e-07, 1.73787203817336e-07, 1.59332972949381e-07] solmin_HZE=[0.000214564255932794, 0.0002220707051897, 0.000273814994353869, 0.000300151006539606, 0.000338926989041028, 0.000382588464944393, 0.000422106918021267, 0.000450354020477539, 0.000500314408098908, 0.000554729554730863, 0.000622657225113376, 0.000715857127703108, 0.000778668979787929, 0.000845947401752176, 0.000913153654473259, 0.000967034820562101, 0.000957235660237084, 0.000968624753464856, 0.000986510357341186, 0.00099859302596883, 0.00101376729471048, 0.00106387381648718, 0.00110164694994432, 0.00115062481890885, 0.00120171603163767, 0.00126460690512873, 0.00135536905193233, 0.00147109976200382, 0.0015697457068759, 0.00170333275048959, 0.00179331947399527, 0.00189538159852076, 0.00198363118555431, 0.00209089863015865, 0.00223234669782025, 0.00235848929920634, 0.00247629666410585, 0.00264173843221103, 0.00277929464010793, 0.00287665603530738, 0.00308554757630749, 0.00336609107912791, 0.00362593136591712, 0.00387385989750139, 0.00418011871559169, 0.00450927382351706, 0.00480281271266942, 0.00507407677387648, 0.00533954889350182, 0.00558504236282106, 0.00584663063446837, 0.00602879870739444, 0.00632697501655301, 0.0069791531388927, 0.00738236830235405, 0.00790294468326847, 0.00857234212765228, 0.00905459346827126, 0.00928726555043449, 0.00999916568943308, 0.0105446918594748, 0.0110949105667061, 0.0117120623040772, 0.0122510624610247, 0.0126621303194278, 0.0130084766095763, 0.0135641497733831, 0.0140630701770245, 0.0147500251210206, 0.0153785999111824, 0.0160698225846111, 0.016616261480822, 0.0173475235283072, 0.0180513815101455, 0.0188855106798456, 0.0198484270308914, 0.0207784648377053, 0.021791039362486, 0.0227132598157029, 0.0239836062903542, 0.025157662008279, 0.026328181350262, 0.0275344825191774, 0.0289217860898329, 0.0298556250399179, 0.0313881153240033, 0.0326621897255238, 0.0338213345634186, 0.0350459532054037, 0.0362827561926899, 0.0371416876720653, 0.0386348090970137, 0.0401416820768685, 0.041823238797461, 0.0436843145575525, 0.0457288170502505, 0.0480017501496143, 0.0506399965299901, 0.0531367919401178, 0.0556228928370107, 0.0579810219031708, 0.0603752009971007, 0.0624760955328979, 0.0646268610352981, 0.0669912713919652, 0.0696386854337891, 0.0720979393196274, 0.0745850155625288, 0.0776888670368998, 0.0810518195154629, 0.0842979167749528, 0.0876827718398196, 0.0915577565025639, 0.0949653244987626, 0.0981982404959759, 0.101901511594913, 0.105365232003343, 0.108516847888778, 0.111942173858989, 0.115696601493006, 0.119175052206619, 0.123415653872957, 0.12802900309868, 0.132928481751152, 0.137696390256706, 0.142963859216492, 0.14777159659881, 0.152709628824408, 0.157640663082535, 0.162769311948002, 0.168657151584575, 0.174172910079832, 0.179955751198685, 0.185866027024302, 0.192215993789953, 0.198253906682679, 0.204264466888564, 0.210683783187876, 0.217762400192331, 0.224627480432708, 0.231026309839595, 0.238178440655628, 0.244887957483572, 0.251345863889339, 0.257656747749017, 0.263454275042928, 0.270825202583938, 0.277932534315691, 0.284651423760285, 0.292148938572119, 0.300606515796536, 0.307795581736833, 0.316099699290545, 0.32516397873817, 0.335783375685102, 0.348218400147042, 0.362542833074625, 0.377940186256751, 0.393704264520378, 0.408594003009793, 0.422123612408801, 0.434679157906757, 0.44727113107965, 0.461937477467131, 0.480440777456119, 0.50185104236025, 0.526928976507665, 0.556186933905194, 0.59228874884709, 0.630541048590304, 0.677550747857123, 0.73244079255957, 0.79602237285, 0.872103167484794, 0.961640300741023, 1.06420044575251, 1.17712530197953, 1.31318525926817, 1.46640502171861, 1.6487535868373, 1.86100775722792, 2.110017834614, 2.3844613430284, 2.69064774589509, 3.0258904983929, 3.39711171765528, 3.80564845591909, 4.27560114002517, 4.78406224675778, 5.32721460319418, 5.90769952112289, 6.50525634949568, 7.08264671128574, 7.66385815323889, 8.24889824780902, 8.82495346414842, 9.41617941470799, 10.0301503428538, 10.6497308591243, 11.2359754922977, 11.7981115520471, 12.307396436486, 12.7528881520937, 13.118523154053, 13.4197342327564, 13.6656128364052, 13.857207562509, 14.00681011617, 14.10663319579, 14.1713081571795, 14.1558797281412, 14.0826372432053, 13.9710593622514, 13.8443160648191, 13.6956239222207, 13.5280359672289, 13.3302416880344, 13.0708450865318, 12.7674092451101, 12.4383023090865, 12.1035867477991, 11.7767666592991, 11.4532875441794, 11.1216498010872, 10.7674396624373, 10.4176432446587, 10.0615391190834, 9.71985035233195, 9.39028716821752, 9.07422341217167, 8.73954619248433, 8.40013323569895, 8.06743149506856, 7.75379159912222, 7.46598831620271, 7.20507623384439, 6.95683292783334, 6.70169782498313, 6.43889902432281, 6.18281527175558, 5.93451659651466, 5.70507840232409, 5.4933571334181, 5.28547854147456, 5.0611384202079, 4.84851959313199, 4.64648180790352, 4.4571595959217, 4.28988022415599, 4.13963072291817, 3.98110308677422, 3.81360150349464, 3.65487547867327, 3.50060196647365, 3.3579724318301, 3.22101111160955, 3.08872669441342, 2.94817655592483, 2.80690550789473, 2.6723245816267, 2.5543820690361, 2.44525709826644, 2.34412881006104, 2.24401628172594, 2.13986336450758, 2.03551415764641, 1.93804341096133, 1.84629741082057, 1.76041538563592, 1.68029970666559, 1.59871319063658, 1.51358745929849, 1.43124952911511, 1.35439915196167, 1.28260392879189, 1.21465438903863, 1.15097741972213, 1.08531956573329, 1.01873421627855, 0.956110498318738, 0.899709207611178, 0.847578762543609, 0.801094645518261, 0.75689206954292, 0.711714924361884, 0.667284800470128, 0.62645450731153, 0.587899599336694, 0.552345311195063, 0.520075959624212, 0.488619457047784, 0.457772283258974, 0.429110843582238, 0.40240159459993, 0.37792011081744, 0.355580275143821, 0.333562242952709, 0.310920260359475, 0.289082428498136, 0.268510293452998, 0.249779251130496, 0.233203686639462, 0.218466871918184, 0.204637479707422, 0.190488677637909, 0.177030526598657, 0.164240175049485, 0.15234151716576, 0.141601154884173, 0.132121365899601, 0.122867207654557, 0.113525812054514, 0.104942870756292, 0.096872471000999, 0.0895571779367192, 0.0829875230646381, 0.0771148788623314, 0.0711999705605684, 0.0656299895496941, 0.0605983961114196, 0.0560677497624837, 0.0520096129537598, 0.0484657266572277, 0.0450472259028317, 0.0415942444695759, 0.0384015791960406, 0.035512037485516, 0.0329067684129079, 0.0305687321635362, 0.0284778631285451, 0.0263962197725322, 0.0243742699152447, 0.022585686522192, 0.0210415688360585, 0.0197740453729641, 0.0187708034969524, 0.0179611045659179, 0.0171426426199955, 0.0163404430766561, 0.015531477253069, 0.0146862042641202, 0.0138331646073546, 0.012971951665373, 0.0121011390767524, 0.0112193264418243, 0.0103473882173873, 0.00946965865607302, 0.00857717959321993, 0.00772832278865701, 0.00701559960815179, 0.00641910040310046, 0.00596419030287563, 0.00563956019803972, 0.00538571491354776, 0.0051265585446647, 0.00486404319921053, 0.00458930349168972, 0.004317571041554, 0.00404676769769398, 0.00377307931735834, 0.00350815446600993, 0.00323490920196816, 0.00295856780962721, 0.00267861112295851, 0.00239696898149799, 0.00214442435998431, 0.00193896363924479, 0.00177538619287212, 0.0016498451747402, 0.00156580401216882, 0.00148896448479232, 0.00141250529188782, 0.00133438063420606, 0.00126031787627835, 0.00118582859734232, 0.00110808113417455, 0.00103078028698015, 0.00095365729734873, 0.000873890691238572, 0.000794512649615267, 0.000717559439971153, 0.000642610469736107, 0.000578504907520296, 0.000525240869844578, 0.000483356193907314, 0.000452295680346081, 0.000430478444725833, 0.000408915795899749, 0.000387674047451901, 0.000366372118766268, 0.000345191252477877, 0.000323575810322741, 0.000302110151837082, 0.00028012830616539, 0.000258125253400382, 0.000236041195362336, 0.000214749055729901, 0.000193150415623937, 0.000173600310310749, 0.000157352755007337, 0.000144315801771664, 0.000133626382956978, 0.000126105089560199, 0.000120084990533732, 0.000113722421542704, 0.000107410311028105, 0.000101436814110885, 9.5491714605726e-05, 8.95548719966522e-05, 8.36601384749298e-05, 7.77400505086735e-05, 7.1853882264682e-05, 6.59678870655649e-05, 6.00529251094327e-05, 5.41966210499324e-05, 4.89729691022986e-05, 4.46393122858375e-05, 4.11247425504994e-05, 3.84128656519503e-05, 3.65076043153501e-05, 3.47677517408919e-05, 3.30203204556327e-05, 3.12716953312064e-05, 2.95510024267247e-05, 2.78415030722707e-05, 2.6150898280602e-05, 2.44353391472847e-05, 2.27503279709102e-05, 2.10725964184173e-05, 1.93894626006477e-05, 1.77090055536432e-05, 1.60811875956513e-05, 1.46843970825689e-05, 1.35013403018494e-05, 1.25564472895374e-05, 1.18342683683306e-05, 1.12978894901595e-05, 1.07565957098946e-05, 1.02194321609319e-05, 9.69016270132361e-06, 9.16309944154498e-06, 8.6516889592805e-06, 8.12546151221063e-06, 7.60480492635986e-06, 7.08804967491378e-06, 6.57020489698572e-06, 6.04714059458919e-06, 5.54272019250938e-06, 5.06821794918582e-06, 4.6685592520356e-06, 4.32825203247063e-06, 4.05458964428214e-06, 3.84876294008039e-06, 3.68107376810744e-06, 3.5098732032176e-06, 3.34677782675e-06, 3.17947692338983e-06, 3.01201633462291e-06, 2.84311991600183e-06, 2.6681630987386e-06, 2.49907666111461e-06, 2.3301001580662e-06, 2.16341448029708e-06, 2.0000778958616e-06, 1.84102274041735e-06, 1.70012357587858e-06, 1.5829685947412e-06, 1.4830611739445e-06, 1.40269415391589e-06, 1.34046900550859e-06, 1.27924548365456e-06, 1.21851063070868e-06, 1.16042890630479e-06, 1.10297424892277e-06, 1.04666591414129e-06, 9.91038977954367e-07, 9.35344102898426e-07, 8.808669210885e-07, 8.27068855209095e-07, 7.72372281661101e-07, 7.17879424961941e-07, 6.64922969591134e-07, 6.18512701051646e-07, 5.7694435460194e-07, 5.43254823677967e-07, 5.16808743025708e-07, 4.95316047822744e-07, 4.73541988914402e-07, 4.53420252005361e-07, 4.32907629849069e-07, 4.12400470899587e-07, 3.94389503867196e-07, 3.75485140621769e-07, 3.5588922649084e-07, 3.36402318495566e-07, 3.17029654870887e-07, 2.95664464632665e-07, 2.72309595119033e-07, 2.50252443605867e-07, 2.20213082952832e-07] solmin_total=[solmin_protons[i] + solmin_HZE[i] for i in range (len(solmin_HZE)) ] solmin_total_edep=[(solmin_total[i]*33.5*sepduration) for i in range(len(solmin_total)) ] #Scaled up for duration of SEP event ######Worstweek worstweek_protons=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2.86028325985846e-07, 6.28955250998595e-07, 9.94036630747278e-07, 1.29744051256935e-06, 1.68353668608819e-06, 2.00216967757635e-06, 2.45081635420052e-06, 3.03633818090247e-06, 3.8611088779347e-06, 5.37605495818172e-06, 7.10776575908786e-06, 1.15219482728308e-05, 1.67770107036141e-05, 2.26332369792461e-05, 2.88226091084889e-05, 3.56252439792457e-05, 4.09040143671155e-05, 4.56631435142096e-05, 5.15137683211964e-05, 5.74767296589002e-05, 6.32645924527257e-05, 6.95445577286617e-05, 7.87143912532061e-05, 8.83926846279208e-05, 9.96942878044734e-05, 0.000114565933299497, 0.000128855012958064, 0.000144404119181741, 0.000166186958029207, 0.000185737191449502, 0.000202485332103291, 0.000221000715751156, 0.000243255860952163, 0.000264179853925988, 0.000287142084575638, 0.000314107151111042, 0.000344467465523661, 0.000373778519763821, 0.000407578352994409, 0.000449907052696428, 0.00049969780316557, 0.000549872166655445, 0.00060596417367566, 0.000656179387827632, 0.000705863114283909, 0.000757741889792345, 0.000809595765380039, 0.000876630801999847, 0.00095376058396684, 0.00103500073160775, 0.0011273049858888, 0.00123223213687669, 0.00133994055072059, 0.00145283616755781, 0.00157667650367556, 0.0017108316241989, 0.0018377755159775, 0.00197356957551301, 0.00211826464013439, 0.00225999412531775, 0.00242071790169196, 0.00259118683966553, 0.00276843924183865, 0.00297057011937411, 0.00318914313301343, 0.00342917642186646, 0.00369679356041438, 0.0039701330942172, 0.00424069233017652, 0.00458219897488312, 0.00504539180327845, 0.0055885863457675, 0.00626295836316171, 0.00702298092533008, 0.00782606861200058, 0.00877937362737303, 0.00967038086462671, 0.0108780642758261, 0.0126889337735713, 0.0141870316950035, 0.015953864782233, 0.0184041617244612, 0.0210339004711877, 0.0241601883411295, 0.0305996049133079, 0.040288986019271, 0.0483607722535677, 0.0608584416846927, 0.0809636314231148, 0.0947961466262335, 0.114303803605448, 0.139718077132657, 0.359724849420804, 0.414448322240086, 0.489509707388423, 0.570770864247108, 0.68885216543942, 0.673448551432495, 0.866504846351324, 1.09785483959559, 1.40293359480959, 1.8049979983166, 2.26315652757884, 2.77328171013824, 3.47527229674754, 4.06161559878637, 5.00580412826207, 6.04047225173087, 7.01784583471431, 7.77043823895699, 9.09838499628067, 10.1276370928348, 11.3420539564938, 12.7521634408826, 14.5257105300626, 17.0469758057671, 19.3709993293026, 21.8860995921107, 25.11462763687, 28.4281285597778, 31.3597387738737, 35.7152117999452, 40.5025459229058, 44.4809744868138, 48.8565419090791, 53.1459155950246, 58.0255933468757, 63.1362910961024, 68.2722839352834, 73.9628678384299, 80.2896859487849, 84.348161557995, 89.3055869886465, 96.392518085848, 105.527062334717, 114.982180332315, 125.007778941344, 133.084333876812, 143.340235049245, 155.989689409746, 166.029868126494, 180.737782149736, 195.350727041393, 209.526580865909, 218.537704481129, 235.498281474206, 249.904364663177, 267.721981423254, 287.255087831004, 308.739992092865, 325.183758299492, 342.895385995233, 372.603577224724, 401.942234399088, 430.472429709417, 460.429062758487, 488.246468474276, 502.705055275842, 540.98493530212, 584.035648543183, 628.10898188256, 674.037144611886, 717.146960977642, 732.213992745017, 755.404212980427, 808.968153760637, 858.062653472111, 913.95904290021, 971.321228943617, 1016.91711384903, 1032.47052602701, 1087.00888694548, 1155.86251477657, 1243.70233165659, 1324.63668950792, 1399.64159528028, 1426.76765681154, 1486.15736197678, 1574.47006056645, 1685.1620673542, 1805.82937899019, 1920.85881971011, 1984.46341298435, 2038.16361728684, 2149.41485772755, 2284.58554588084, 2414.0056160845, 2544.77070214811, 2590.42003670295, 2646.81529131269, 2740.39530468232, 2881.47154944564, 3019.06463852069, 3174.60684805532, 3204.53551673662, 3179.12140562159, 3224.55481676753, 3315.69560880528, 3459.96008586184, 3639.56715344507, 3784.38036550272, 4014.06757796562, 4223.86765582876, 4416.45933866711, 4601.49986106955, 4772.58931630842, 4892.31280282292, 5046.29068357172, 5160.29443504327, 5278.67516550654, 5372.30491701373, 5311.7216573136, 5338.11303888318, 5404.71118260587, 5461.00878024564, 5558.34957742502, 5657.55276331985, 5815.5649095674, 6039.10842924661, 6258.2855825624, 6431.61479810562, 6557.94918087545, 6616.85037194249, 6670.1820360926, 6713.76874786342, 6746.62649104775, 6783.06764623807, 6687.83107437402, 6559.0098965929, 6458.91507131881, 6434.93395785813, 6460.07197676848, 6471.54909135659, 6535.39480245182, 6604.77325500318, 6624.88617902801, 6599.13827792081, 6620.2459273576, 6603.36286302086, 6572.72746284, 6526.58971633283, 6451.42035432463, 6355.77462895393, 6229.84394539044, 6120.52287067405, 6036.26476642367, 5986.67442545103, 5844.20453827081, 5699.31384458828, 5601.1798932951, 5541.04105770448, 5523.0179979764, 5578.03782698649, 5641.38474551446, 5668.71974953537, 5709.36623456288, 5724.10841518771, 5695.36458552089, 5623.5803737494, 5557.91398676593, 5462.96338281165, 5351.72933534192, 5222.1473372674, 5086.03802063677, 4863.24154376974, 4606.20044264136, 4407.44641162378, 4273.71799036943, 4186.88720462098, 4151.49874124896, 4149.28660314636, 4145.53825066454, 4108.28288316475, 4038.58597160255, 3969.52028117338, 3890.07822183009, 3774.3401026788, 3648.09792391649, 3507.50507396449, 3344.5081943806, 3160.09033446685, 2955.6579083016, 2767.81624308528, 2593.24163229304, 2441.78128171869, 2315.31638158784, 2228.67722814594, 2146.56719335316, 2055.42190204842, 1955.87655881021, 1855.28665678531, 1749.86815456624, 1645.00792561445, 1534.95644461691, 1423.32716070205, 1312.07493597192, 1197.64160186527, 1082.18114780716, 970.442336318806, 874.545835113494, 796.825258268512, 736.052029887358, 689.790939499927, 655.881302892353, 622.133742316209, 587.286259402882, 552.247372879598, 517.323558843683, 484.768315280916, 450.645965068815, 418.544661608126, 390.126819105351, 359.937587083279, 329.517400798856, 299.560049014163, 269.940861510125, 243.024005962512, 222.087567257221, 205.466173694718, 192.816787370332, 183.370893515008, 173.57946117821, 164.217803735144, 155.866128429654, 146.903421223343, 137.802665405494, 128.692917942345, 119.567535036267, 109.660225870751, 100.906556854798, 92.2073782061917, 82.909668925573, 74.9472417411824, 68.0931828640379, 62.4987370275471, 58.221106693536, 55.3816513804674, 52.8995305534153, 50.413519119504, 47.7976652720948, 45.1753630754478, 42.5352520277917, 39.8876816945055, 37.3589102865919, 35.0074305557439, 32.6052052799275, 30.1066112055912, 27.5565039151351, 25.0124836079256, 22.5956196297243, 20.7022160975732, 19.2977569600342, 18.2374906331174, 17.4310173478096, 16.6367899124931, 15.8049814460245, 14.9572736001011, 14.154303803695, 13.3975997299471, 12.6128989388517, 11.814849020079, 11.0211461313871, 10.1799350087753, 9.330512943795, 8.49661000183548, 7.73259639498216, 7.11496565014195, 6.59992315822473, 6.18753472135218, 5.86484972225439, 5.60319473306205, 5.3315386819218, 5.07760977698893, 4.82491828274235, 4.58987722221338, 4.33804533277014, 4.08188933922112, 3.82106488668063, 3.56250170776916, 3.30597333307046, 3.06306806085813, 2.81940605143064, 2.58958537647987, 2.38703112284355, 2.22660599124282, 2.10120402086838, 2.00158861414455, 1.91525173231636, 1.81999995335599, 1.73104880672304, 1.63832540096253, 1.55526294179679, 1.47069391454736, 1.39869119949772, 1.31893008859068, 1.23725131126772, 1.14909326956184, 1.05993871269406, 0.976929279867227, 0.906110525412422, 0.846356730623899, 0.800477653877721, 0.771327939913363, 0.735981203347851, 0.700002654722743, 0.668227152496343, 0.634036511704591, 0.595669272115745, 0.567288382127278, 0.54090406367739, 0.51146056768832, 0.48830429643214, 0.468619667673374, 0.437520862176884, 0.404669121118841, 0.37472572692231, 0.325446954544036] worstweek_HZE=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.95652278619485e-08, 8.7001025982178e-08, 1.37501366912238e-07, 1.79470090384329e-07, 2.32877329087915e-07, 2.76952638304317e-07, 3.39012254004943e-07, 4.20005296954855e-07, 5.30187658480321e-07, 7.25903861101142e-07, 9.44106383517508e-07, 1.48756146924103e-06, 2.13972873555782e-06, 2.87282735075114e-06, 3.65453812918419e-06, 4.51990178744558e-06, 5.18004549124766e-06, 5.76182853713362e-06, 6.4717356908496e-06, 7.17462506797624e-06, 7.86458144039386e-06, 8.62560630845938e-06, 9.6515440728023e-06, 1.07346325683779e-05, 1.19740063659834e-05, 1.35496597343728e-05, 1.50548906135237e-05, 1.67026107927153e-05, 1.88831023951018e-05, 2.08923801544487e-05, 2.26080123571398e-05, 2.44426617813434e-05, 2.66437194605552e-05, 2.88280243297517e-05, 3.10848449452789e-05, 3.37408596711159e-05, 3.66987067299441e-05, 3.95721693642682e-05, 4.27491388803651e-05, 4.68231468669218e-05, 5.1513225590564e-05, 5.61033514174518e-05, 6.10200191170808e-05, 6.53923029364957e-05, 6.95609517182836e-05, 7.4060874424345e-05, 7.87233902270331e-05, 8.47306368049866e-05, 9.096859956703e-05, 9.80013377162278e-05, 0.000105292034856585, 0.000113225373292366, 0.000121105904876894, 0.000129443144264905, 0.000138203181511262, 0.000148047386602009, 0.00015751963459105, 0.000167046924896581, 0.000177076659222859, 0.000186350910045693, 0.000195783193502305, 0.000206143389534016, 0.000216013670706126, 0.000228259631154115, 0.000242021026567144, 0.000255840600644028, 0.000270576586382607, 0.0002863738960925, 0.000300860883736193, 0.000316458121736617, 0.000333357280883316, 0.000354242856865403, 0.0003793664162914, 0.000407980883907387, 0.00043985639359735, 0.000473383428081145, 0.000506832037618163, 0.000538262097341357, 0.000569613959038602, 0.000601336021648352, 0.00064053382246316, 0.00069166332917011, 0.000751471299825134, 0.000826970771766271, 0.000918281652081181, 0.00103510735338668, 0.00119197057453536, 0.00137759774985001, 0.00160548172561877, 0.00187888648354598, 0.00222420206400015, 0.00260641591785279, 0.00309861252386903, 0.0037169558104028, 0.00451626740696359, 0.00547056927426999, 0.00658449269257924, 0.00799879563399874, 0.00961679378927091, 0.0114868383946167, 0.0135802531265115, 0.0162544389504115, 0.0194159906595499, 0.0228527570687765, 0.0271625209531828, 0.0322669762910326, 0.0379485673237429, 0.0440234166211507, 0.0512560727094492, 0.058113482815831, 0.0665401147806068, 0.0751442565213302, 0.0851257406371426, 0.0953170874460371, 0.107824854649114, 0.120238005988883, 0.134679281676683, 0.152026618806798, 0.174238098947702, 0.194203574609256, 0.213825039334004, 0.232389649655168, 0.248778258393268, 0.266859926138636, 0.289015358585104, 0.340396586735725, 0.392846623094227, 0.42248057053995, 0.448469020737354, 0.478138218396422, 0.478392933158949, 0.491306026117061, 0.538257133315788, 0.647770833934796, 0.685527506469181, 0.752589594144697, 0.837077879446532, 0.852925180854581, 0.818302534983329, 0.852219216628992, 0.869824954168121, 0.883770468147445, 0.945889463386293, 0.972330920386783, 1.00521480306241, 1.03769379795514, 1.16084597010536, 1.24418566173023, 1.32301766373196, 1.44193377650339, 1.53733671605181, 1.49926494256848, 1.54798107747102, 1.63167140947159, 1.69028227140486, 1.81708590050845, 1.91603302325532, 1.97377445108975, 2.1113402912017, 2.24620454394193, 2.32832799989155, 2.42665334185866, 2.56792797688725, 2.62370258964965, 2.68921428134054, 2.87691662407712, 3.11918760954209, 3.31316115053295, 3.53035623300512, 3.72083295226367, 3.79862779028624, 4.05245233952655, 4.39001099019702, 4.77141898604413, 5.17558596504425, 5.54130679400552, 5.66262724241286, 5.99825697394172, 6.46231167829312, 6.99132388025137, 7.60409871080114, 8.24493977231283, 8.55969490780778, 9.00175233605764, 9.80747478597796, 10.65320890474, 11.5238585466545, 12.4511301701983, 12.9712563289029, 13.64964652762, 14.8042240045303, 16.1528885545131, 17.4542848573647, 18.7511266285984, 19.2265387957732, 19.6228437681608, 20.6523903185826, 22.2243521185281, 24.0504171909243, 26.1493430323682, 28.3265590883362, 31.3954691621173, 34.5757321232571, 37.8183846377531, 40.9120664079399, 43.5148600836059, 47.2839231893588, 51.6242971874537, 56.417332041902, 61.5013141551153, 66.3234081040085, 68.3015873920014, 73.2699665999226, 78.755102017325, 84.5376655989226, 90.9119559865107, 98.0197367687097, 108.742527227451, 122.187172743395, 136.968329040405, 150.693846400354, 164.730035322939, 181.367340616901, 198.294142758507, 215.457185590219, 234.897888405094, 250.987800607855, 263.827243076059, 276.639294669981, 289.510628528034, 301.487983437578, 315.251441487638, 328.347888977877, 349.181134570037, 370.678228432873, 392.692940229545, 414.909787134231, 440.921046045955, 460.789287795055, 478.37560901461, 492.706212912901, 505.421657943608, 519.302673724632, 529.890469761587, 540.584025827399, 552.54703533625, 563.664583723671, 559.24771838256, 559.788047449274, 569.241227359692, 584.492796315767, 606.903627547982, 644.32805461885, 681.967223560781, 711.488333005498, 735.599712071275, 754.048512666661, 758.928578633315, 764.185622648273, 773.959424913358, 784.434348626092, 792.425864866206, 802.250339504824, 800.596537935691, 783.857844926019, 762.346148139794, 746.261064039221, 736.727950827039, 739.887351559631, 754.606562993234, 774.261626913762, 791.302941124503, 803.760464990055, 805.866137240989, 803.362267395928, 797.289907559301, 787.297004957124, 771.085160442277, 752.395612662453, 731.199470941137, 705.424150363019, 674.056845449577, 646.968794980936, 622.966282319927, 602.451933825093, 587.591482699929, 580.397191692827, 570.926020224501, 558.728896767569, 544.227168745353, 527.6183950388, 508.844849267613, 487.847725504039, 464.020362326097, 437.350305532489, 408.253943616226, 377.614212966802, 345.367605751784, 312.908265748315, 285.170731968772, 262.036416954363, 243.566119648975, 230.079257774205, 220.855817362476, 211.1232016905, 201.144633263167, 190.928707522089, 180.42388604504, 169.712536306469, 158.99241001106, 148.296091161634, 137.332305886106, 126.425450285975, 115.44864175696, 104.333895504623, 93.9070054018743, 85.1204947534815, 77.8537177338357, 72.2390770601125, 68.3333173252258, 65.1787979709367, 62.0935853864289, 58.9497208106116, 55.8118391834523, 52.6846784424522, 49.566232552787, 46.3766235937467, 43.2031462051995, 40.0050567846018, 36.7259927935309, 33.4923554010657, 30.2297723669003, 27.3454529459943, 24.9032939273206, 22.9307705361867, 21.4162664898216, 20.3789949090224, 19.4259501423356, 18.4817483786573, 17.5315997955683, 16.5715612104495, 15.6268398020027, 14.690168193653, 13.7536859918071, 12.8150432544233, 11.8639056014956, 10.9199748302025, 9.97394348747561, 9.05004036398756, 8.26699743988741, 7.61003203433448, 7.08732414872706, 6.69021732606113, 6.39730876190375, 6.09961239091604, 5.80294349750999, 5.4955758349578, 5.19374209979069, 4.90425224609204, 4.61083810299339, 4.32095708752987, 4.04697879251546, 3.76293037402991, 3.47154165374398, 3.18467912861886, 2.91404163954917, 2.66955346427647, 2.46894927409881, 2.30883370801866, 2.18912590297651, 2.08990643412916, 1.98985979172, 1.89191035455068, 1.79279656838163, 1.69474768936763, 1.60263101308183, 1.50676933868478, 1.41299549189253, 1.32116375403087, 1.22688269106713, 1.13005249439605, 1.0371588919181, 0.955073931861287, 0.883061721724214, 0.82382857473827, 0.776480258926177, 0.740272206438354, 0.706346754578455, 0.673770037110264, 0.641735826308507, 0.610019586700792, 0.57864844451338, 0.547219348448417, 0.516210522046125, 0.485510225900386, 0.45498048949825, 0.424523402948534, 0.39427014886776, 0.364796726862836, 0.338625785784933, 0.316453518451272, 0.29819047136565, 0.28380801003865, 0.272088642658147, 0.260756063355113, 0.249700505849016, 0.23858077288906, 0.227604333193398, 0.217069267838771, 0.206018049464131, 0.194616188943648, 0.18341398156056, 0.172272435281836, 0.161233938108338, 0.150706722373812, 0.140398155458899, 0.129858039275435] worstweek_total=[worstweek_protons[i] + worstweek_HZE[i] for i in range(len(worstweek_HZE)) ] worstweek_total_edep=[(worstweek_total[i]*33.5*sepduration) for i in range(len(worstweek_total)) ] ######Worst5min worst5min_protons=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.56284401161132e-07, 7.83443192862361e-07, 1.2381981557165e-06, 1.61612600594782e-06, 2.09705754830147e-06, 2.49395517783071e-06, 3.05280127100576e-06, 3.7821426489072e-06, 4.82708648639101e-06, 6.77647876691297e-06, 9.02964820231158e-06, 1.48304664628792e-05, 2.17129520386769e-05, 2.93542858784966e-05, 3.73992032459394e-05, 4.62228593921407e-05, 5.31410990290021e-05, 5.94443133802063e-05, 6.72234243528258e-05, 7.52580068318829e-05, 8.30167156303916e-05, 9.13831079586985e-05, 0.000104112748999231, 0.000117515483293409, 0.00013329771326114, 0.000154312634101251, 0.000174566447492714, 0.000196510514745697, 0.000228016360461838, 0.000256016809358855, 0.000280156471105975, 0.000307199335891593, 0.000339908110427919, 0.000370123982209055, 0.000404077554575379, 0.000443894171816678, 0.000488998920199935, 0.000532265519990325, 0.000583035425454419, 0.000645978394397, 0.000721041394826248, 0.000797951324898159, 0.00088640415874093, 0.000966015551282737, 0.00104613288360911, 0.00112814903268517, 0.00120893717858263, 0.00131406317721179, 0.00144082466129608, 0.00157028947570344, 0.00172450630169582, 0.00190252968542957, 0.00208685988697966, 0.00228161653168968, 0.00249791293926526, 0.00272788898643009, 0.00294545013449869, 0.00318248183195171, 0.00343507144760347, 0.00369118185149715, 0.00400662422148938, 0.00433368518617237, 0.00469238564929907, 0.00508862966727057, 0.00551208428581326, 0.00599299146489751, 0.00653831555312295, 0.00708417008305305, 0.00764595248132247, 0.00840918671182078, 0.0095253354901128, 0.0108278687098823, 0.0124718390925885, 0.0143121705934174, 0.0162094769926974, 0.0185857039510102, 0.0207857929338462, 0.0242350406282876, 0.0300324451002141, 0.0343795139204797, 0.0395739536224244, 0.0470316188529566, 0.0553322062256887, 0.0649958954866029, 0.0913126811258592, 0.131577821314371, 0.160295596564966, 0.209549517116049, 0.295527448540616, 0.347242011155626, 0.42749388997611, 0.539038432821775, 2.22209801529804, 2.47018863278219, 2.81305984052473, 3.19873310443454, 3.8047375448595, 3.18882884940672, 4.32236114041117, 5.66728824339012, 7.51812289976635, 10.0141218436478, 13.0500282556153, 16.3912992481238, 21.1837737700704, 24.8030794623684, 31.2636677361674, 38.120285765692, 44.5322384517724, 49.0183741073642, 58.1936237010754, 64.8628361211227, 73.1927091115101, 82.390479731944, 94.4471737502323, 112.641982863224, 129.317420687647, 147.277641387028, 172.077288866551, 197.785224215735, 220.078089302257, 257.796125171275, 296.203699705678, 327.537091341235, 363.536611978518, 399.056369161986, 438.897962332074, 480.26902133259, 527.353473215133, 575.924324527734, 635.583292638808, 664.935778457156, 705.946076802001, 761.068898672584, 852.857099145668, 939.825830565947, 1031.85137133852, 1108.92885670182, 1215.22882442829, 1323.81432000862, 1401.57323748754, 1540.8812275732, 1666.83798737372, 1774.03978109112, 1850.67007607874, 2001.91361348834, 2113.84145996462, 2268.69415394431, 2451.01126348402, 2643.74878116303, 2798.24202365222, 2964.66800602451, 3291.09589760027, 3581.33031881478, 3856.45098685863, 4141.52728919788, 4404.15171451372, 4504.75618742107, 4924.39968479547, 5349.06784079857, 5788.8829507087, 6248.42621949126, 6678.54815589119, 6785.65626379438, 7028.76256197407, 7602.95210673522, 8152.99521865557, 8772.62079011463, 9402.34980331056, 9915.09899842631, 10120.5986332737, 10775.954657694, 11590.9157258742, 12642.7822312417, 13629.7201522822, 14568.2545755921, 14929.2918129389, 15725.2212963966, 16902.0243044089, 18357.7038042646, 19942.9055901876, 21499.6418515651, 22403.2111775581, 23253.1839875588, 24909.9579757033, 26906.2630679708, 28799.5782713208, 30761.2747521946, 31578.2986615201, 32634.0258508642, 34301.0751069838, 36643.3168985514, 38955.2748095194, 41481.6640980523, 42212.0270232661, 42175.5839763826, 43332.4630145709, 45216.6847062157, 47938.838931886, 51187.7826640475, 53909.3482598789, 58225.7169170669, 62294.2347347545, 66101.3299139315, 69792.5386120283, 73212.219654656, 75970.2151574801, 79456.9200208059, 82328.4239608184, 85271.0244023593, 87738.4596626639, 87315.1485143215, 88648.842336537, 90761.7748378561, 92754.7928684973, 95529.7114853511, 98389.6683296072, 102569.489332622, 108063.172740193, 113442.619449181, 117852.282830584, 121391.190287978, 123762.416842328, 126101.376107274, 128242.140632249, 130128.601266656, 131948.716832794, 131050.830191985, 129563.357319821, 128689.032531167, 129407.423189782, 131166.346982078, 132531.438462788, 135192.044994244, 138033.760333813, 139713.49870428, 140304.691517669, 142055.081416011, 142894.628737572, 143295.077375035, 143371.997693247, 142771.313423255, 141654.432164253, 139718.975312853, 138206.630002451, 137220.177266588, 136945.061817162, 134309.469354951, 131710.862150858, 130278.413230948, 129703.649135424, 130072.070867811, 132226.950181388, 134501.712254176, 135849.210456236, 137556.411806328, 138655.48392505, 138584.26243191, 137456.156152053, 136570.839757497, 134953.765112711, 132901.586544754, 130398.621176561, 127698.109109984, 122700.526943234, 116762.45231348, 112237.364721575, 109334.761407379, 107583.452232195, 107098.48863758, 107444.93168054, 107757.790655807, 107179.620782428, 105757.130559744, 104332.402106316, 102606.919395367, 99896.3340985813, 96885.8045688151, 93441.3111639789, 89383.7630804278, 84740.8574680245, 79536.7243577845, 74723.09961075, 70206.9856652021, 66244.9130756408, 62908.5125725355, 60614.3474425985, 58433.1080885662, 56003.136661259, 53328.5911052139, 50612.1616448993, 47764.7706834518, 44910.9279045405, 41913.0426698711, 38870.8016756012, 35836.381803197, 32710.0250665399, 29567.8145377535, 26519.0266696449, 23899.8308506199, 21778.6406698949, 20119.1299756397, 18853.9651243091, 17925.4749510054, 17003.3125702715, 16049.8649182238, 15091.526009588, 14137.4435031291, 13249.1509511884, 12314.7265456223, 11436.4148327955, 10661.8658629754, 9837.85520801964, 9006.41024416474, 8189.08443256406, 7379.77138430041, 6642.60071965049, 6069.92376639326, 5615.24830645877, 5268.18332801561, 5009.58060088038, 4741.74752567012, 4485.53803214268, 4258.27225225381, 4013.91811093565, 3765.52022288282, 3516.58514894122, 3267.47221144477, 2996.55854856029, 2757.65585776128, 2520.17320464795, 2266.29029987379, 2048.80758081668, 1861.20624408031, 1707.88877433066, 1590.46137993623, 1512.57861006832, 1444.73098136355, 1376.6537278221, 1305.34503474118, 1233.81838414605, 1161.61724675658, 1089.27399881431, 1020.11521357604, 955.585647789701, 889.963623589109, 821.913225466485, 752.437358428212, 683.03920734916, 617.21549964825, 565.682915878091, 527.320922543769, 498.366128011416, 476.368142580471, 454.795573057729, 432.093377695891, 408.934043702238, 386.965693146366, 366.336957820923, 344.949671944457, 323.115311448679, 301.404561806802, 278.381168148572, 255.168263415901, 232.385992545053, 211.534034332682, 194.65460616614, 180.595113973868, 169.324201484345, 160.510943625751, 153.351255375331, 145.927301228263, 139.011882368208, 132.070830729415, 125.578448209875, 118.669738729044, 111.652131060829, 104.483005999505, 97.412465149679, 90.4208218887156, 83.7950551830936, 77.1412092341892, 70.8639179836669, 65.336179427906, 60.9470820486378, 57.5288922225776, 54.8148437947098, 52.4576668942006, 49.8356894604909, 47.3806846304788, 44.8166231667487, 42.5264023571742, 40.2052287803077, 38.2492620901959, 36.0826964365749, 33.8552393895999, 31.4422029999541, 28.9979087864671, 26.7226291447302, 24.7888223196218, 23.1505888959906, 21.8919907319499, 21.0897375491104, 20.1191005500085, 19.1218639524229, 18.2554215733698, 17.3242649605224, 16.2793412221755, 15.5030008919633, 14.7894275578973, 13.9796274349032, 13.3472552165281, 12.8155162463381, 11.9683347666421, 11.0707053467937, 10.2622333282661, 8.90269529053292] worst5min_HZE=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.06615397464453e-09, 2.34439417290092e-09, 3.70521381461293e-09, 4.83613416459082e-09, 6.27528522969585e-09, 7.46297120155137e-09, 9.13527563450568e-09, 1.13177742406391e-08, 1.444510650178e-08, 2.02799936546615e-08, 2.70247120691509e-08, 4.43903943692564e-08, 6.49936963669847e-08, 8.78680318890784e-08, 1.11949790139216e-07, 1.38362365519252e-07, 1.59073232348387e-07, 1.77944495739982e-07, 2.01235234627425e-07, 2.25293410988563e-07, 2.48524577969001e-07, 2.73574215511767e-07, 3.11701454165539e-07, 3.51843458476176e-07, 3.9911519512476e-07, 4.62065197503181e-07, 5.22736972237091e-07, 5.88469012418371e-07, 6.82861853180319e-07, 7.66746246043889e-07, 8.39069995329723e-07, 9.20100506262135e-07, 1.01811628885595e-06, 1.10865147335041e-06, 1.21040439910693e-06, 1.32972631788829e-06, 1.46490353798078e-06, 1.59456241406766e-06, 1.7467288007621e-06, 1.93536403112775e-06, 2.16035495273273e-06, 2.39092338284691e-06, 2.65618066835711e-06, 2.89493828604764e-06, 3.13524980494054e-06, 3.38120597610237e-06, 3.62344774368048e-06, 3.93869153685606e-06, 4.31898404767491e-06, 4.70727771701623e-06, 5.17002904723512e-06, 5.70428012627389e-06, 6.25748826315808e-06, 6.8420900062669e-06, 7.49138942498746e-06, 8.18161237454059e-06, 8.83459853712685e-06, 9.54701942282956e-06, 1.03065766961485e-05, 1.10840677878567e-05, 1.20492133808857e-05, 1.30541224072432e-05, 1.41579575068014e-05, 1.5406319366088e-05, 1.67361442879246e-05, 1.82647512647266e-05, 2.00340716434154e-05, 2.18812466880985e-05, 2.38648150350289e-05, 2.61512462259365e-05, 2.86649816097255e-05, 3.19331139554969e-05, 3.60693425689826e-05, 4.08735551061395e-05, 4.63409360890504e-05, 5.21925457328674e-05, 5.89805091180348e-05, 6.54555634179648e-05, 7.27176565091943e-05, 8.18899396787287e-05, 9.62920986424411e-05, 0.000115713410792072, 0.000140714712313456, 0.000195533357334558, 0.000256880592105724, 0.000342609830115771, 0.00067484585895294, 0.000912324213158317, 0.00130022278070665, 0.00171135399921814, 0.00235830079608711, 0.00283847257226379, 0.00381627651580621, 0.00532880463989984, 0.00756954705396023, 0.0102716330384427, 0.0131767937695765, 0.0174650690248059, 0.0220824116428129, 0.0274122163748642, 0.0334990611965686, 0.0430214728145712, 0.0553695451792497, 0.0676293632388389, 0.08526009836397, 0.106239898415866, 0.130493163657681, 0.155723493854201, 0.190434740402189, 0.220409437245219, 0.267970065444951, 0.313071187231292, 0.370000586200844, 0.426913545402759, 0.504301471228727, 0.577217675261697, 0.674985365131961, 0.836622414713972, 1.10103188331236, 1.28015373540507, 1.46233658964034, 1.63499445715041, 1.75931028541137, 1.89674687067403, 2.1374141849372, 2.98737038747549, 3.84392498215883, 4.21702798222333, 4.50524761429712, 4.84705855780161, 4.57431306951647, 4.53031138584701, 5.2028239188242, 7.1604892745265, 7.6797521655077, 8.82488255497889, 10.3710593103183, 10.4573444792376, 9.43172263771736, 9.76735140689314, 9.71798478064291, 9.64330342839283, 10.5830224643131, 10.8132130782936, 11.1825102341194, 11.4116750905777, 13.4961860054077, 14.7102194909299, 15.8745433932543, 17.9033489264611, 19.5175651744332, 18.1046102582141, 18.1528805032259, 18.8023859555656, 18.9169944634403, 20.5352930567625, 21.7639083336488, 22.4934280496295, 24.4674358544023, 26.186081404846, 26.575187999048, 27.3924048629865, 28.9533807940312, 29.0957258427772, 29.2060399548637, 31.1849296841871, 33.9123554145807, 35.7341210400464, 38.1203295077019, 40.3027389592387, 40.6813086311567, 43.7000613934052, 47.8178655326366, 52.7357447769366, 58.0061768272843, 62.914606141536, 64.2165531094161, 68.6386149821166, 74.5873966082071, 81.6195952665133, 90.2284882334077, 99.4818595650548, 103.984804741336, 110.520947266389, 122.228939739247, 134.270582070738, 146.799082134552, 160.704851785749, 168.548547128401, 179.836709725459, 199.048535305133, 221.159314427325, 242.614956074574, 264.159285122879, 272.768805468288, 281.197129191335, 301.358601212234, 330.86557879364, 364.62463339258, 402.688415801551, 443.144045314063, 503.133645852797, 566.401187832207, 631.30462020658, 693.505156112618, 746.701392036818, 828.27785520957, 924.26554093376, 1030.1033605197, 1141.76201861892, 1247.9831827475, 1294.89373238376, 1415.59042954914, 1553.13309469571, 1699.2672918166, 1858.89421665896, 2035.7758172564, 2315.10990457018, 2661.89406258178, 3043.09907547913, 3399.12977169138, 3757.85120794639, 4146.06743805087, 4527.36379488268, 4904.66998000718, 5330.33759073324, 5678.94703920663, 5937.52991552876, 6184.19609556959, 6426.73636840867, 6645.32292224732, 6904.80970454397, 7152.24927671185, 7565.75766995711, 7996.79298324503, 8441.82694916291, 8892.35753590745, 9427.10006795372, 9834.86852620987, 10195.0864492569, 10485.4783477822, 10740.2991099976, 11023.1192003506, 11235.3934553549, 11447.2742178227, 11687.3030130863, 11908.4180336444, 11798.1765178512, 11797.0897013465, 11981.9395003438, 12284.4392318543, 12740.0156526351, 13513.6440276065, 14288.5251721546, 14896.0963332427, 15394.5606956581, 15775.9352022839, 15870.9820514765, 15974.6448819592, 16176.5660165172, 16397.189845188, 16561.8281817892, 16765.8178961487, 16737.1283422443, 16394.5966403304, 15952.0840028957, 15622.0540852774, 15431.812538801, 15508.1843689942, 15824.7064834491, 16241.1991091579, 16607.5808783237, 16879.7006820222, 16930.5626973919, 16885.9694974543, 16769.0501669355, 16570.0889118994, 16237.8769192099, 15853.9845237933, 15416.9434691481, 14883.2905928809, 14230.8701075395, 13667.0722104652, 13163.0884134409, 12728.6278107768, 12411.5711163046, 12255.3330420628, 12049.3782338972, 11786.1226229681, 11474.750188711, 11119.0967594476, 10718.9549050239, 10274.3918036641, 9771.03079716314, 9208.7990089186, 8596.40934028527, 7952.02952546442, 7273.22503783544, 6590.24361098412, 6006.57115072531, 5519.00319076629, 5129.62319413149, 4845.44189022028, 4651.1781893942, 4446.26594631034, 4236.77852822187, 4022.00902365824, 3800.81295843199, 3575.21535802834, 3349.50097072079, 3124.24570867906, 2893.30516941058, 2663.90480450426, 2432.86987464986, 2198.80726168049, 1979.19253532816, 1794.24506866223, 1641.0527272847, 1522.84672133705, 1440.61571622657, 1374.05273237939, 1309.15796579475, 1243.20210120605, 1177.06689916434, 1111.29503907182, 1045.78983688889, 978.581895515358, 911.601128728967, 844.193560909199, 774.920106605361, 706.549460996403, 637.632828278645, 576.714614213244, 525.205865075434, 483.648980364067, 451.767887073351, 429.924011505056, 409.843850038728, 389.937078125655, 369.884564313386, 349.669851529455, 329.75335540686, 310.010570357905, 290.249920252478, 270.440951624839, 250.365010716315, 230.445161024877, 210.479445208471, 190.996908298878, 174.495750546938, 160.648349990365, 149.640700896452, 141.279635073056, 135.114382799144, 128.834134863775, 122.577665518795, 116.081600178667, 109.705643113539, 103.578188830145, 97.3679888319558, 91.2382914836963, 85.4089429774323, 79.3977628495419, 73.2372012795847, 67.1737295128092, 61.4510154933737, 56.3156342651922, 52.078573805026, 48.7070849887974, 46.1802574931124, 44.0885777442235, 41.9835221315094, 39.9179669647422, 37.8282234872108, 35.7697171654088, 33.8309665137778, 31.8045581625309, 29.8228947531113, 27.8787503568705, 25.8841401862525, 23.8364193892071, 21.8781468774392, 20.1473214730114, 18.6313172456287, 17.3832244899526, 16.3836649595695, 15.6202340299788, 14.9047192134511, 14.2182795440926, 13.5435720575451, 12.8743917064438, 12.213111061365, 11.5507044104887, 10.8953377689524, 10.2467676600456, 9.60234612846831, 8.95886107884115, 8.31962818803531, 7.69685787274963, 7.14355764338937, 6.67519263606825, 6.2897116669284, 5.98617976264011, 5.73938319485629, 5.50066499771681, 5.26765730312834, 5.03303086488684, 4.80152717383138, 4.57930003443306, 4.34652853290066, 4.10625402307811, 3.87032210779145, 3.63564042466405, 3.40295799855149, 3.18064730399215, 2.96298240842642, 2.74061659817492] worst5min_total=[worst5min_protons[i]+ worst5min_HZE[i] for i in range(len(worst5min_HZE)) ] #TOtal ion production from the Dartnell paper #Keep in mind here that this is assumed to be over 20 hours, so must be divided by 20 hours worth of seconds to give a rate/s #Carrington-like event# Carrington_Aug1972=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.93e-269, 5.21e-268, 8.47e-267, 1.52e-266, 2.67e-266, 3.8e-266, 5.05e-266, 5.98e-266, 6.86e-266, 7.23e-266, 7.69e-266, 8.69e-157, 3.11e-156, 5.16e-156, 7.75e-156, 9.94e-156, 7.93e-92, 1.1e-90, 1.89e-90, 7.09e-90, 2.1e-89, 2.3e-89, 7.74e-71, 1.41e-68, 8.87e-38, 3.76e-20, 3.69e-19, 9.02e-19, 1.17e-14, 2.35e-14, 2.88e-13, 4.71e-13, 1.24e-12, 0.790882986, 40.43292993, 40.56532568, 40.83902695, 41.3710513, 42.50937776, 7.428190481, 151.7715222, 206.6370463, 246.7863562, 254.2317375, 350.1649568, 3359.951344, 15009.44967, 17565.14939, 21502.30418, 53806.50637, 69160.45443, 60434.65793, 82304.11084, 109832.0131, 352361.7955, 370991.3923, 416535.2775, 541048.913, 1107841.045, 1192467.661, 2545643.968, 3790582.408, 4481719.01, 7780517.953, 8858188.85, 8865285.186, 9087329.699, 10693339.88, 11313821.92, 11923987.8, 14829945.86, 17451266.82, 22144287.47, 19591656.24, 19629800.74, 18208239.06, 25055172.25, 29406569.86, 35668058.09, 40471398.05, 49995563.36, 52668865.33, 48590490.08, 53835066.23, 55806496.51, 50819442.14, 50663761.17, 53104319.18, 47673085.16, 46971560.83, 53914776.13, 58035887.05, 66149563.36, 72732681.27, 93665413.23, 107494112.3, 115217729.1, 125533603.6, 138098056.7, 135717653, 166513758.7, 202436669.4, 245124999, 287902958.6, 333628505.7, 354677135.4, 400837551.1, 510640409.2, 633245208, 771042205.5, 905142644.4, 1023689573, 1092768371, 1283200661, 1525819686, 1841819721, 2153874443L, 2453260126L, 2607254688L, 2904456352L, 3338211889L, 3850028015L, 4399356945L, 4939031714L, 5287092456L, 5619486386L, 6200747587L, 6880747383L, 7522608395L, 8188328304L, 8502200506L, 8828563067L, 9303928725L, 9953239081L, 10570279799L, 11244630797L, 11436866005L, 11371855691L, 11517748661L, 11796129185L, 12244903376L, 12809594109L, 13236840307L, 13789046850L, 14231850397L, 14607201219L, 14943476503L, 15252178677L, 15303683810L, 15368255036L, 15297214521L, 15224342398L, 15106084117L, 14711947120L, 14409849618L, 14150520100L, 13844729462L, 13598809748L, 13333610039L, 13095376621L, 12947420766L, 12801292225L, 12634588224L, 12381272601L, 11993274508L, 11583370583L, 11170946423L, 10755942374L, 10403907199L, 9913765623L, 9370429564L, 8863499037L, 8438580411L, 8064854731L, 7728480308L, 7412968160L, 7093204358L, 6766614987L, 6432910162L, 6116957815L, 5806900231L, 5530837025L, 5236154561L, 4934334710L, 4640226527L, 4367968759L, 4102764434L, 3869164036L, 3675070737L, 3470576061L, 3251832761L, 3055028612L, 2887515932L, 2753865873L, 2649128061L, 2567632546L, 2486732812L, 2405571928L, 2310059327L, 2214691397L, 2105281947, 1986610715, 1865199571, 1746952447, 1625867904, 1509773468, 1384044037, 1257897199, 1154306364, 1071188307, 1007789819, 962526940, 928958674.8, 893559599.1, 853119062.4, 806692013.8, 763459407.1, 720862090, 675431545.9, 629980634.1, 585159764.2, 539360537.6, 491273394.1, 442369256.6, 400187621.9, 364057513.5, 335204560.5, 312576673.3, 297706973, 283566825.2, 268380636.3, 252920615, 238402463.8, 223402309, 209434336, 195194525.4, 180788427.6, 166489514.1, 152055752.3, 136963573.5, 122526800.7, 110221838.8, 100222308.2, 92344908.14, 86452805.81, 82356312.65, 78242195.16, 73978427.99, 69637771.13, 65220728.33, 61016567.67, 56825600.07, 52852683.15, 49089305.48, 45284287.35, 41478228.13, 37631259.36, 33880928.02, 30574849.12, 27884190.19, 25765258.94, 24239563.3, 23095632.57, 21934189.63, 20796037.45, 19665230.95, 18501195.71, 17324967.38, 16183084.81, 15009902.56, 13789827.94, 12652609.99, 11544961.04, 10351993.49, 9344948.706, 8512434.044, 7838016.518, 7334060.191, 7012215.586, 6695440.594, 6382120.029, 6039094.064, 5703931.271, 5368690.827, 5042143.287, 4724798.031, 4430175.589, 4119884.276, 3796596.586, 3457111.47, 3133152.751, 2835216.568, 2591677.221, 2408978.238, 2274552.735, 2171439.627, 2060691.104, 1953856.285, 1847798.148, 1750309.066, 1653032.53, 1553964.154, 1457514.901, 1360405.022, 1257261.586, 1151116.781, 1047160.615, 950510.487, 875350.935, 812911.6644, 764642.5086, 724353.2555, 692359.4304, 658532.9092, 627201.2404, 596851.8809, 570233.2582, 539289.9526, 507072.9075, 474449.3447, 441715.6072, 409288.492, 378075.3236, 347743.527, 319158.605, 292930.907, 271560.7421, 255249.5411, 242287.3799, 231554.5215, 221923.9988, 213785.3309, 204892.5683, 195969.2766, 186310.7069, 175952.425, 164342.4117, 152743.4399, 141744.7895, 130647.8729, 120554.2683, 111897.1275, 105074.9086, 99291.32055, 96031.04256, 92242.74801, 88883.63037, 85072.3956, 80902.80182, 76029.95406, 72229.55092, 68097.73568, 64589.26637, 61716.93412, 58844.45031, 55006.80183, 51046.5224, 46416.94102, 41212.08205] Carrington_Aug1972_edep=[(Carrington_Aug1972[i]*33.5) for i in range(len(Carrington_Aug1972)) ] Carrington_Aug1972_ionrate=[(Carrington_Aug1972[i]/sepduration) for i in range(len(Carrington_Aug1972)) ] Carrington_Feb1956=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9.83e-10, 2.16e-09, 3.42e-09, 4.46e-09, 5.78e-09, 6.88e-09, 8.42e-09, 1.04e-08, 2.83e-08, 8.68e-08, 1.75e-07, 4.49e-07, 7.55e-07, 1.07e-06, 1.38e-06, 2.25e-06, 4.13e-06, 6.21e-06, 9.16e-06, 1.3e-05, 1.64e-05, 2.22e-05, 6.97e-05, 0.000111558, 0.000157722, 0.000209654, 0.000262948, 0.000296926, 0.000386319, 0.000469136, 0.000593289, 0.000758108, 0.001508155, 0.002289824, 0.003128301, 0.004277321, 0.005755554, 0.006767786, 0.008267174, 0.010159836, 0.017874119, 0.050264398, 0.137798501, 0.221894234, 0.308584209, 0.384191505, 0.447195196, 0.529029669, 0.687613651, 0.833029118, 1.110880567, 1.425140396, 1.735642894, 2.67553732, 3.423632494, 4.099468339, 4.81330344, 5.560697455, 6.032591137, 8.474297647, 32.91028388, 53.73589862, 88.71180938, 122.9651858, 160.67948, 196.2236129, 234.9050291, 260.4467492, 292.4326667, 514.8169738, 1113.953416, 1718.683963, 2518.626364, 3325.262422, 3990.757535, 5898.672308, 7512.742979, 14965.74231, 33839.57717, 40454.37459, 48856.27327, 64791.24109, 92812.95474, 117885.1532, 301088.9382, 584610.6844, 717717.9617, 1040527.08, 1694208.894, 1981305.328, 2483255.402, 3057175.799, 10242719.72, 12156351.5, 14935806.18, 17967040.44, 22824094.16, 24551047.63, 32330796.69, 43512851.92, 58072399.88, 78191560.64, 96706092.47, 119952597.2, 145971112, 169550490.9, 207907409, 254961238.3, 299300982.4, 341979257.7, 402832135.4, 453927936.4, 514428548.8, 592173247.2, 681379600.5, 801251367, 923324012.4, 1047612252, 1196884950, 1349707829, 1493734611, 1662277831, 1868888258, 2043975066, 2247597104L, 2450979817L, 2667210915L, 2905529372L, 3132210708L, 3362490307L, 3597490329L, 3827416585L, 4033790940L, 4311622244L, 4591687581L, 4868150834L, 5162535555L, 5395547272L, 5646648095L, 5980926844L, 6277959458L, 6607507159L, 6981687989L, 7358956936L, 7574826540L, 7922013794L, 8228023284L, 8541645960L, 8831058536L, 9182640596L, 9422201541L, 9642623884L, 9918567658L, 10236319879L, 10508989002L, 10795537672L, 11084387540L, 11264556823L, 11549549622L, 11941274923L, 12286762061L, 12623767097L, 12909474763L, 12948424681L, 12987589919L, 13275613791L, 13435493964L, 13646299823L, 13893889032L, 14015903279L, 13879874198L, 13987960738L, 14210675617L, 14484999615L, 14678805396L, 14800986923L, 14695750251L, 14655836071L, 14725858794L, 14933984809L, 15188336969L, 15357355045L, 15332005650L, 15187956445L, 15209116277L, 15311313714L, 15437502213L, 15487760689L, 15250194827L, 15000781407L, 14796228695L, 14744723009L, 14701995225L, 14769744675L, 14456479120L, 13997714626L, 13637712840L, 13383757698L, 13247406290L, 13219740242L, 13122532654L, 13122383434L, 13069658116L, 12978328721L, 12884925775L, 12794106226L, 12578539729L, 12373679148L, 12078591242L, 11815709788L, 11544093032L, 11122142336L, 10783034237L, 10522651180L, 10221466603L, 9974188065L, 9717629456L, 9515957020L, 9382688450L, 9274555220L, 9141964939L, 8954333117L, 8690153480L, 8412620033L, 8133023819L, 7867145655L, 7642903779L, 7318790138L, 6966216955L, 6645940338L, 6389712662L, 6173446775L, 5976112253L, 5804178392L, 5634232904L, 5448345659L, 5250984782L, 5075646226L, 4894069451L, 4725943861L, 4546713109L, 4356216549L, 4164492346L, 3977333728L, 3798377395L, 3641672929L, 3517776510L, 3365240052L, 3203771091L, 3065683245L, 2952540965L, 2866653019L, 2815677072L, 2779734523L, 2734685067L, 2692504523L, 2637044274L, 2571569882L, 2488746331L, 2400643652L, 2303985948L, 2205991722L, 2102164371, 1997750111, 1870085865, 1735197783, 1625998393, 1542421744, 1481646404, 1443059632, 1418274124, 1391993497, 1356303676, 1310575046, 1266620974, 1221146611, 1166605548, 1110085751, 1051820297, 988625099.1, 919839468, 846827227, 781658441.8, 723326014.3, 674608248.4, 635346674, 608897772.2, 583965603.5, 556708703.4, 527897343.8, 499499030.7, 469876763.9, 441261640, 411530517.6, 381418534.8, 351455961.2, 320866122.1, 289525153.6, 259372842.5, 233563719.6, 212642770.9, 196274670.5, 183889718.1, 174957593.4, 166017771.7, 156803660.7, 147512868.5, 138180246, 129399689.3, 120368829.7, 111834973.3, 104104202.1, 96017355.16, 87902657.66, 79836069.73, 71916021.54, 64803169.43, 59203017.19, 54773075.66, 51461480.6, 48974400.26, 46401219.99, 43929898.43, 41641294.09, 39219248, 36767587.32, 34339468.46, 31888864.64, 29263383.3, 26906646.81, 24575262.9, 22082704.53, 19952916.8, 18142784.65, 16672495.82, 15555820.83, 14817930.71, 14154573.41, 13493673.67, 12786426.57, 12083317.6, 11378236.29, 10671852.48, 9997574.865, 9375293.314, 8728699.479, 8052894.046, 7360154.223, 6676871.772, 6029309.679, 5517904.033, 5139960.928, 4856192.135, 4639309.456, 4420046.389, 4196510.472, 3970828.562, 3758785.403, 3554997.975, 3344411.022, 3133659.621, 2923494.219, 2701295.637, 2475390.215, 2253552.497, 2049322.737, 1885457.949, 1748455.477, 1639730.125, 1553368.694, 1483794.351, 1411597.919, 1343941.756, 1278218.326, 1218638.05, 1152722.921, 1084950.472, 1016145.613, 947117.1845, 878151.8306, 812647.0417, 747610.932, 686459.9409, 631955.2107, 588612.5165, 554773.487, 527839.6207, 504692.65, 480587.9748, 458640.4795, 435707.6642, 414543.6942, 392611.2957, 372571.3623, 350376.6894, 327905.3987, 304527.8923, 280988.7438, 259158.0818, 240339.4707, 224810.2341, 212679.2777, 205165.0076, 196121.1657, 187314.2274, 178900.0352, 169776.586, 159460.3212, 151825.8282, 144286.7665, 136612.3529, 130423.1376, 124858.485, 116533.3269, 107857.4989, 99316.26834, 86907.47226] Carrington_Feb1956_edep=[(Carrington_Feb1956[i]*33.5) for i in range(len(Carrington_Feb1956)) ] Carrington_Feb1956_ionrate=[(Carrington_Feb1956[i]/sepduration) for i in range(len(Carrington_Feb1956)) ] Carrington_Oct1989=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6.55e-11, 1.44e-10, 2.28e-10, 2.97e-10, 3.86e-10, 4.59e-10, 5.62e-10, 6.96e-10, 1.27e-09, 3e-09, 5.53e-09, 1.32e-08, 2.19e-08, 3.09e-08, 3.98e-08, 5.56e-08, 8.31e-08, 1.13e-07, 1.55e-07, 2.09e-07, 2.57e-07, 3.24e-07, 7.22e-07, 1.08e-06, 1.48e-06, 1.96e-06, 2.46e-06, 2.82e-06, 3.67e-06, 4.43e-06, 5.47e-06, 6.83e-06, 1.11e-05, 1.56e-05, 2.05e-05, 2.7e-05, 3.53e-05, 4.12e-05, 4.98e-05, 6.02e-05, 9.31e-05, 0.000202277, 0.000482867, 0.000752891, 0.001033235, 0.00127523, 0.001478251, 0.00174586, 0.00225943, 0.002726926, 0.00361033, 0.004614805, 0.005612159, 0.008328426, 0.010552974, 0.012551841, 0.014660929, 0.016884806, 0.018376964, 0.025353199, 0.093444667, 0.151498358, 0.248878602, 0.344299287, 0.449309235, 0.547704928, 0.654451679, 0.7246436, 0.812427984, 1.548538504, 3.543826063, 5.52875609, 8.141864277, 10.73625164, 13.00102066, 21.40778897, 28.40702599, 64.8029823, 159.0812369, 186.5316877, 221.2115294, 294.5807127, 537.3067035, 886.0618598, 3793.11128, 8367.549846, 9837.818834, 13711.43108, 22647.64397, 25740.05231, 34703.83274, 49521.39986, 615877.0248, 644943.6592, 686122.3211, 740638.8328, 850094.3245, 486783.3171, 760696.2732, 1078877.348, 1544828.088, 2191128.306, 3049019.575, 3984056.293, 5406071.849, 6325263.94, 8230517.927, 10154701.61, 11964213.24, 13041814.86, 15754068.09, 17555405.25, 19858410.77, 22329941.1, 25799728.1, 31517871.54, 36408374.06, 41950974.55, 49476398.75, 57474651.42, 64493230.91, 76130869.7, 89158141.09, 100185340.1, 112023411.5, 123555682.8, 137348073.9, 151576502.2, 166015961.1, 182949674.1, 201764539, 213230168.4, 228157023.4, 249706605.9, 277299898.2, 306767063.3, 337848271.3, 363176719.7, 395091171.1, 435501733.8, 467890677.4, 515346306.2, 562350657.4, 609259529.8, 639277202.3, 695670046.2, 744686138.6, 805326607.4, 871784973.8, 944738799.5, 1000638586, 1061760055, 1162233020, 1262866536, 1362127848, 1466601799, 1563156143, 1615671882, 1746537193, 1894789101, 2046938384, 2205954515L, 2355659477L, 2411610313L, 2492834114L, 2675027644L, 2841719114L, 3031415615L, 3226176534L, 3381333812L, 3435049850L, 3616880836L, 3845315815L, 4135869296L, 4402865608L, 4649779583L, 4739025658L, 4931324989L, 5215100024L, 5570758261L, 5958310448L, 6325835777L, 6526682755L, 6691072777L, 7035664339L, 7455097464L, 7857139203L, 8261082582L, 8394994905L, 8557878831L, 8831351987L, 9254276939L, 9665057109L, 10134322936L, 10211149249L, 10114602016L, 10231169244L, 10487579704L, 10906595321L, 11435398134L, 11857380083L, 12531909846L, 13143858322L, 13703583059L, 14240453505L, 14737105444L, 15076326281L, 15516914487L, 15834871449L, 16167578183L, 16426907974L, 16225656603L, 16287903791L, 16474123045L, 16628879602L, 16908071683L, 17192689039L, 17664298767L, 18338065188L, 19000661941L, 19526415474L, 19911340889L, 20103079955L, 20283498818L, 20436618436L, 20557969418L, 20689981486L, 20421659491L, 20061560119L, 19795363733L, 19766591340L, 19891824498L, 19977185859L, 20247931832L, 20544610953L, 20683768256L, 20676285288L, 20836495459L, 20880793834L, 20876527091L, 20825122601L, 20680856501L, 20473687657L, 20164951979L, 19920772694L, 19758098719L, 19702588434L, 19312950564L, 18934805280L, 18735871140L, 18664775233L, 18733548123L, 19066212785L, 19426555025L, 19660108199L, 19951185339L, 20154729857L, 20184877637L, 20065320528L, 19997094715L, 19830691917L, 19604604872L, 19319224463L, 19013277653L, 18359409473L, 17558902744L, 16967130262L, 16620218479L, 16453733309L, 16474865912L, 16624910541L, 16772885542L, 16779411975L, 16659847843L, 16546556690L, 16382996968L, 16057668882L, 15681267623L, 15222173462L, 14668398096L, 14019010603L, 13271784673L, 12574589657L, 11904233983L, 11295990670L, 10768914411L, 10403976068L, 10052081544L, 9654067237L, 9207139767L, 8750309968L, 8270870472L, 7782648290L, 7268724147L, 6744851124L, 6220416467L, 5678280798L, 5137360221L, 4608363170L, 4153289223L, 3785423833L, 3496994799L, 3276061240L, 3114429153L, 2954921243L, 2789228088L, 2622772496L, 2457074481L, 2302944977L, 2139751095, 1986609813, 1852376172, 1709683737, 1565499043, 1424340729, 1283821334, 1154967858, 1054994328, 975590744.6, 914583031.3, 869503047.5, 823232376.7, 778693308.5, 739478997.2, 697167985.9, 654078342.6, 610751212.8, 567488856.4, 520376905.3, 478910857.4, 437742173.4, 393726217.8, 356000266.4, 323325722.6, 296574971.1, 276008944, 262451432.1, 250635412.9, 238730040.6, 226377864.9, 214000903.8, 201410340, 188883786.8, 176863659.5, 165524732.7, 154102901.9, 142374278.4, 130352660.9, 118346454.1, 107051793.8, 98216906.25, 91551801.15, 86524450.26, 82720591.83, 79003521.85, 75051967.73, 71022838.51, 67206762.26, 63647463.45, 59963960.47, 56177154.86, 52408491.64, 48390189.47, 44355175.43, 40407098.19, 36797935.63, 33874997.34, 31456153.37, 29508082.72, 27973823.1, 26726686.7, 25434615.7, 24246772.32, 23028737.45, 21878282.49, 20665293.46, 19437700.96, 18169669.18, 16938195.08, 15733358.92, 14585697.41, 13432475.57, 12342724.24, 11384216.39, 10613875.11, 10023987.2, 9556034.849, 9149296.439, 8691125.204, 8263116.486, 7810358.607, 7405616.655, 7000194.969, 6661952.772, 6285557.138, 5894104.751, 5472238.532, 5043179.954, 4644727.372, 4311318.402, 4027189.639, 3807192.696, 3667202.6, 3498664.682, 3320286.781, 3171695.555, 3011973.941, 2831241.712, 2694569.798, 2572227.8, 2429586.05, 2320135.408, 2230262.778, 2085947.407, 1931315.377, 1793295.509, 1553664.257] Carrington_Oct1989_edep=[(Carrington_Oct1989[i]*33.5) for i in range(len(Carrington_Oct1989)) ] Carrington_Oct1989_ionrate=[(Carrington_Oct1989[i]/sepduration) for i in range(len(Carrington_Oct1989)) ] ##775 AD event## AD_Aug1972=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.19e-266, 3.22e-265, 5.23e-264, 9.38e-264, 1.65e-263, 2.35e-263, 3.12e-263, 3.69e-263, 4.24e-263, 4.46e-263, 4.75e-263, 5.37e-154, 1.92e-153, 3.19e-153, 4.79e-153, 6.14e-153, 4.9e-89, 6.78e-88, 1.17e-87, 4.38e-87, 1.29e-86, 1.42e-86, 4.78e-68, 8.7e-66, 5.48e-35, 2.32e-17, 2.28e-16, 5.57e-16, 7.24e-12, 1.45e-11, 1.78e-10, 2.91e-10, 7.64e-10, 488.6257919, 24980.3988, 25062.19596, 25231.29493, 25559.99187, 26263.27627, 4589.307797, 93767.95492, 127665.1441, 152470.3159, 157070.2446, 216340.0051, 2075855.613, 9273184.981, 10852155.35, 13284620.59, 33242903.49, 42728927.54, 37337928.74, 50849382.31, 67856756.67, 217697262.8, 229207058.4, 257345123.5, 334272525.9, 684449807.9, 736734087.6, 1572757691, 2341909439L, 2768909609L, 4806983854L, 5472793847L, 5477178127L, 5614362360L, 6606592575L, 6989940725L, 7366915310L, 9162283376L, 10781796065L, 13681252708L, 12104178124L, 12127744681L, 11249471013L, 15479664621L, 18168058656L, 22036550831L, 25004165296L, 30888414795L, 32540042553L, 30020328040L, 33260548428L, 34478543628L, 31397426152L, 31301242846L, 32809076004L, 29453534072L, 29020116128L, 33309795046L, 35855912640L, 40868729428L, 44935931857L, 57868657563L, 66412347529L, 71184176510L, 77557562267L, 85320171851L, 83849503431L, 102876000000.0, 125070000000.0, 151444000000.0, 177873000000.0, 206123000000.0, 219128000000.0, 247647000000.0, 315485000000.0, 391234000000.0, 476368000000.0, 559218000000.0, 632459000000.0, 675138000000.0, 792791000000.0, 942687000000.0, 1137920000000.0, 1330710000000.0, 1515680000000.0, 1610820000000.0, 1794440000000.0, 2062420000000.0, 2378640000000.0, 2718020000000.0, 3051450000000.0, 3266490000000.0, 3471850000000.0, 3830970000000.0, 4251080000000.0, 4647640000000.0, 5058940000000.0, 5252860000000.0, 5454490000000.0, 5748180000000.0, 6149340000000.0, 6530560000000.0, 6947190000000.0, 7065960000000.0, 7025800000000.0, 7115930000000.0, 7287920000000.0, 7565180000000.0, 7914060000000.0, 8178030000000.0, 8519190000000.0, 8792770000000.0, 9024670000000.0, 9232430000000.0, 9423150000000.0, 9454970000000.0, 9494860000000.0, 9450970000000.0, 9405950000000.0, 9332890000000.0, 9089380000000.0, 8902740000000.0, 8742520000000.0, 8553590000000.0, 8401660000000.0, 8237810000000.0, 8090630000000.0, 7999220000000.0, 7908930000000.0, 7805940000000.0, 7649440000000.0, 7409720000000.0, 7156470000000.0, 6901670000000.0, 6645270000000.0, 6427770000000.0, 6124950000000.0, 5789270000000.0, 5476070000000.0, 5213550000000.0, 4982650000000.0, 4774830000000.0, 4579900000000.0, 4382350000000.0, 4180570000000.0, 3974400000000.0, 3779200000000.0, 3587640000000.0, 3417080000000.0, 3235020000000.0, 3048550000000.0, 2866840000000.0, 2698630000000.0, 2534780000000.0, 2390460000000.0, 2270540000000.0, 2144200000000.0, 2009060000000.0, 1887470000000.0, 1783970000000.0, 1701400000000.0, 1636690000000.0, 1586340000000.0, 1536360000000.0, 1486220000000.0, 1427210000000.0, 1368290000000.0, 1300690000000.0, 1227370000000.0, 1152360000000.0, 1079310000000.0, 1004500000000.0, 932773000000.0, 855094000000.0, 777158000000.0, 713157000000.0, 661805000000.0, 622636000000.0, 594671000000.0, 573932000000.0, 552062000000.0, 527077000000.0, 498393000000.0, 471683000000.0, 445365000000.0, 417297000000.0, 389217000000.0, 361525000000.0, 333229000000.0, 303520000000.0, 273306000000.0, 247245000000.0, 224923000000.0, 207097000000.0, 193117000000.0, 183930000000.0, 175194000000.0, 165812000000.0, 156260000000.0, 147291000000.0, 138023000000.0, 129393000000.0, 120596000000.0, 111695000000.0, 102861000000.0, 93943558876L, 84619261906L, 75699889916L, 68097600047L, 61919658823L, 57052818992L, 53412541963L, 50881633790L, 48339836900L, 45705582979L, 43023824824L, 40294873676L, 37697446032L, 35108169365L, 32653609447L, 30328507724L, 27977679562L, 25626208196L, 23249461956L, 20932420558L, 18889848586L, 17227497303L, 15918372593L, 14975762553L, 14269015702L, 13551449406L, 12848272682L, 12149634286L, 11430466401L, 10703765348L, 9998283904L, 9273464786L, 8519674479L, 7817074945L, 7132743819L, 6395700888L, 5773525338L, 5259178535L, 4842507797L, 4531151928L, 4332308891L, 4136597978L, 3943021290L, 3731091919L, 3524020597L, 3316901301L, 3115152682L, 2919089447L, 2737064892L, 2545359745L, 2345625137L, 2135883384, 1935734199, 1751662338, 1601198099, 1488322444, 1405271261, 1341565599, 1273142601, 1207137580, 1141612411, 1081381403, 1021281710, 960074977.7, 900486399.2, 840489671, 776765271.7, 711186558, 646960035, 587247351.9, 540812043.2, 502235618.4, 472413818.1, 447522186.2, 427755661.4, 406856854.7, 387499425.4, 368748889.5, 352303289, 333185799.5, 313281364.3, 293125773, 272902113.3, 252867891.9, 233583674.9, 214843989.8, 197183564.2, 180979486.1, 167776504.2, 157699067.1, 149690744.3, 143059736.3, 137109776.7, 132081519.5, 126587365.2, 121074349.3, 115107061.7, 108707475.7, 101534541, 94368428.12, 87573207.68, 80717276.06, 74481213.79, 69132632.13, 64917707.57, 61344473.16, 59330198.05, 56989702.1, 54914361.57, 52559692.62, 49983621.2, 46973063.21, 44625086.3, 42072355.32, 39904741.88, 38130148.62, 36355461.71, 33984473.77, 31537721.59, 28677459.18, 25461777] AD_Aug1972_edep=[(AD_Aug1972[i]*33.5) for i in range(len(AD_Aug1972)) ] AD_Aug1972_ionrate=[(AD_Aug1972[i]/sepduration) for i in range(len(AD_Aug1972)) ] AD_Feb1956=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2.75e-09, 6.05e-09, 9.57e-09, 1.25e-08, 1.62e-08, 1.93e-08, 2.36e-08, 2.92e-08, 7.93e-08, 2.43e-07, 4.9e-07, 1.26e-06, 2.11e-06, 3e-06, 3.86e-06, 6.3e-06, 1.16e-05, 1.74e-05, 2.56e-05, 3.64e-05, 4.59e-05, 6.21e-05, 0.000195, 0.000312429, 0.000441716, 0.000587159, 0.000736413, 0.000831574, 0.001081926, 0.001313865, 0.001661569, 0.00212316, 0.004223746, 0.006412891, 0.008761136, 0.011979085, 0.016119034, 0.018953895, 0.023153088, 0.028453685, 0.050058344, 0.140770716, 0.385919152, 0.621438069, 0.864222434, 1.075968597, 1.25241704, 1.481603061, 1.925734133, 2.332985397, 3.111137515, 3.991255117, 4.860849922, 7.493122818, 9.588241798, 11.48099095, 13.480161, 15.57331631, 16.89490405, 23.73315917, 92.16870094, 150.4930188, 248.4467244, 344.376897, 449.9997324, 549.5448038, 657.8761658, 729.4084312, 818.9883472, 1441.798918, 3119.743349, 4813.354658, 7053.677234, 9312.746098, 11176.53495, 16519.85033, 21040.22449, 41913.13065, 94771.28426, 113296.7181, 136827.1163, 181454.6646, 259932.412, 330149.7328, 843231.1433, 1637263.524, 2010044.411, 2914105.197, 4744809.66, 5548853.329, 6954617.145, 8561941.398, 28685810.61, 34045137.08, 41829291.35, 50318580.75, 63921269.02, 68757783.3, 90545786.3, 121862304.5, 162637845.3, 218983664.7, 270835552.5, 335939826.6, 408807405.5, 474843928.6, 582266500, 714045682.9, 838223785.7, 957748770.9, 1128173636, 1271272784, 1440711093, 1658443273, 1908275020, 2243988472L, 2585865716L, 2933947964L, 3352001807L, 3779998304L, 4183360410L, 4655383371L, 5234017537L, 5724366502L, 6294631369L, 6864225985L, 7469803849L, 8137239678L, 8772084529L, 9417006691L, 10075148897L, 10719081485L, 11297054507L, 12075150206L, 12859502551L, 13633766880L, 14458222160L, 15110795907L, 15814030092L, 16750212780L, 17582083762L, 18505016655L, 19552949307L, 20609530553L, 21214096015L, 22186430325L, 23043441990L, 23921775179L, 24732305448L, 25716947874L, 26387863420L, 27005179305L, 27777988779L, 28667887192L, 29431525664L, 30234035262L, 31042989604L, 31547572573L, 32345724792L, 33442792575L, 34410365531L, 35354183472L, 36154337749L, 36263421080L, 36373107435L, 37179748490L, 37627509681L, 38217893594L, 38911293136L, 39253006829L, 38872043123L, 39174750813L, 39798487183L, 40566760307L, 41109533719L, 41451715897L, 41156989562L, 41045205715L, 41241311694L, 41824190422L, 42536530317L, 43009883163L, 42938889522L, 42535464620L, 42594724949L, 42880939586L, 43234343710L, 43375097844L, 42709769736L, 42011261299L, 41438389980L, 41294142904L, 41174479265L, 41364218706L, 40486885670L, 39202067597L, 38193844848L, 37482616841L, 37100750426L, 37023268757L, 36751028716L, 36750610810L, 36602948027L, 36347170479L, 36085585735L, 35831236066L, 35227519485L, 34653785947L, 33827361322L, 33091134247L, 32330443041L, 31148725870L, 30199018080L, 29469788026L, 28626289036L, 27933759562L, 27215240274L, 26650435470L, 26277202865L, 25974364415L, 25603031430L, 25077548837L, 24337686062L, 23560424542L, 22777386027L, 22032766348L, 21404753461L, 20497039222L, 19509621053L, 18612652803L, 17895060330L, 17289385034L, 16736729014L, 16255210208L, 15779260049L, 15258663326L, 14705933496L, 14214879484L, 13706354685L, 13235501340L, 12733546825L, 12200041234L, 11663097499L, 11138940166L, 10637754190L, 10198886900L, 9851901986L, 9424707641L, 8972496885L, 8585767392L, 8268900573L, 8028362375L, 7885598890L, 7784938013L, 7658772287L, 7540641250L, 7385319006L, 7201951108L, 6969995067L, 6723254276L, 6452554243L, 6178111137L, 5887331752L, 5594908667L, 5237371562L, 4859603343L, 4553778999L, 4319713830L, 4149506120L, 4041439818L, 3972025403L, 3898423751L, 3798470666L, 3670402843L, 3547304855L, 3419949133L, 3267201166L, 3108911547L, 2945733034L, 2768748256L, 2576106884L, 2371628447L, 2189116430L, 2025750350, 1889311138, 1779354982, 1705282059, 1635456907, 1559121100, 1478431866, 1398899413, 1315939148, 1235799493, 1152534367, 1068202602, 984289272.8, 898619220.6, 810845551.9, 726400842.9, 654119687.9, 595528377.6, 549687795.9, 515002438.1, 489987086.1, 464950178, 439145094.1, 413125256.3, 386988268.5, 362397398.6, 337105529.3, 313205569.6, 291554734.2, 268906671.4, 246180610.2, 223589284.8, 201408359.4, 181488071.2, 165804257.6, 153397741.9, 144123272.6, 137157943.4, 129951482.2, 123030287, 116620810.6, 109837616.5, 102971483.7, 96171282.23, 89308109.26, 81955173.45, 75354885.8, 68825600.68, 61844929.62, 55880235.75, 50810770.85, 46693072.81, 43565707.41, 41499168.75, 39641367.08, 37790448.05, 35809728.38, 33840597.99, 31865943.85, 29887641.91, 27999256.75, 26256492.01, 24445638.18, 22552974.2, 20612883.68, 18699279.53, 16885713.98, 15453468.85, 14394999.55, 13600275.29, 12992872.6, 12378803.39, 11752767.62, 11120721.77, 10526872.66, 9956144.602, 9366373.759, 8776142.361, 8187552.115, 7565261.688, 6932589.863, 6311310.072, 5739343.214, 5280422.694, 4896732.905, 4592236.154, 4350371.913, 4155521.669, 3953327.991, 3763849.812, 3579784.454, 3412923.644, 3228321.413, 3038517.562, 2845822.342, 2652500.989, 2459356.284, 2275903.254, 2093762.809, 1922503.046, 1769856.834, 1648471.074, 1553701.323, 1478270.206, 1413444.688, 1345937.017, 1284470.755, 1220245.001, 1160973.084, 1099549.103, 1043425.168, 981266.659, 918333.4529, 852862.2948, 786938.4414, 725799.3831, 673095.8894, 629604.6341, 595630.6186, 574586.1173, 549257.8895, 524593.1355, 501028.3078, 475477.1317, 446585.3505, 425204.1522, 404090.2194, 382597.2192, 365263.673, 349679.2798, 326363.8016, 302066.2356, 278145.6237, 243393.4891] AD_Feb1956_edep=[(AD_Feb1956[i]*33.5) for i in range(len(AD_Feb1956)) ] AD_Feb1956_ionrate=[(AD_Feb1956[i]/sepduration) for i in range(len(AD_Feb1956)) ] AD_Oct1989=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.28e-09, 2.82e-09, 4.45e-09, 5.81e-09, 7.54e-09, 8.97e-09, 1.1e-08, 1.36e-08, 2.49e-08, 5.87e-08, 1.08e-07, 2.59e-07, 4.28e-07, 6.04e-07, 7.77e-07, 1.09e-06, 1.62e-06, 2.21e-06, 3.03e-06, 4.09e-06, 5.03e-06, 6.32e-06, 1.41e-05, 2.11e-05, 2.9e-05, 3.84e-05, 4.8e-05, 5.51e-05, 7.17e-05, 8.66e-05, 0.000107, 0.000133, 0.000218, 0.000305, 0.000401, 0.000528, 0.00069, 0.000805, 0.000973, 0.00118, 0.00182, 0.003952375, 0.009434932, 0.014711026, 0.020188766, 0.024917211, 0.028884111, 0.034113034, 0.044147872, 0.05328244, 0.070543609, 0.090170435, 0.10965811, 0.162732297, 0.206198581, 0.245255217, 0.286465474, 0.329918656, 0.359074514, 0.495385815, 1.825850981, 2.960184183, 4.862933908, 6.727395061, 8.779224479, 10.70181545, 12.78758092, 14.15908765, 15.87434019, 30.25748437, 69.24416887, 108.0284736, 159.086991, 209.779715, 254.0319005, 418.2949525, 555.0557137, 1266.210184, 3108.348953, 3644.713778, 4322.336431, 5755.924883, 10498.64058, 17313.10057, 74115.04779, 163496.7474, 192224.8938, 267912.8806, 442520.952, 502944.6979, 678091.4221, 967617.5166, 12033856.05, 12601800.1, 13406405.68, 14471624.59, 16610317.19, 9511444.867, 14863534.57, 21080595.92, 30184985.13, 42813291.55, 59575956.21, 77845995.21, 105631299.4, 123591744.2, 160819228.5, 198416587.2, 233773325, 254828994, 307824744.2, 343021757.8, 388021061, 436313234.6, 504110726.1, 615839711.2, 711397105.1, 819696089.6, 966738222.9, 1123019132, 1260157833, 1487550095, 1742094919, 1957559565, 2188868157L, 2414201604L, 2683696392L, 2961711080L, 3243849174L, 3574723450L, 3942354271L, 4166385575L, 4458047087L, 4879112597L, 5418268458L, 5994038634L, 6601346210L, 7096248423L, 7719837060L, 8509434454L, 9142294375L, 10069548003L, 10987983945L, 11904554292L, 12491081042L, 13592962325L, 14550706446L, 15735583691L, 17034139056L, 18459611677L, 19551858920L, 20746134616L, 22709314186L, 24675630836L, 26615135474L, 28656491836L, 30543103989L, 31569229026L, 34126256243L, 37023006816L, 39995909683L, 43102986500L, 46028128851L, 47121373558L, 48708436387L, 52268385255L, 55525433441L, 59231985726L, 63037493580L, 66069170807L, 67118748964L, 70671613946L, 75135092107L, 80812327383L, 86029270132L, 90853816448L, 92597629540L, 96355039495L, 101900000000.0, 108849000000.0, 116422000000.0, 123603000000.0, 127527000000.0, 130739000000.0, 137473000000.0, 145668000000.0, 153524000000.0, 161416000000.0, 164033000000.0, 167216000000.0, 172559000000.0, 180823000000.0, 188849000000.0, 198018000000.0, 199520000000.0, 197633000000.0, 199911000000.0, 204921000000.0, 213108000000.0, 223441000000.0, 231686000000.0, 244866000000.0, 256823000000.0, 267760000000.0, 278250000000.0, 287954000000.0, 294582000000.0, 303191000000.0, 309404000000.0, 315904000000.0, 320972000000.0, 317039000000.0, 318256000000.0, 321894000000.0, 324918000000.0, 330373000000.0, 335935000000.0, 345149000000.0, 358314000000.0, 371261000000.0, 381534000000.0, 389055000000.0, 392802000000.0, 396327000000.0, 399319000000.0, 401690000000.0, 404269000000.0, 399027000000.0, 391990000000.0, 386789000000.0, 386227000000.0, 388674000000.0, 390342000000.0, 395632000000.0, 401429000000.0, 404148000000.0, 404002000000.0, 407132000000.0, 407998000000.0, 407914000000.0, 406910000000.0, 404091000000.0, 400043000000.0, 394011000000.0, 389240000000.0, 386061000000.0, 384976000000.0, 377363000000.0, 369974000000.0, 366087000000.0, 364698000000.0, 366042000000.0, 372542000000.0, 379583000000.0, 384146000000.0, 389834000000.0, 393811000000.0, 394400000000.0, 392064000000.0, 390731000000.0, 387479000000.0, 383062000000.0, 377486000000.0, 371508000000.0, 358732000000.0, 343090000000.0, 331527000000.0, 324749000000.0, 321496000000.0, 321909000000.0, 324840000000.0, 327732000000.0, 327859000000.0, 325523000000.0, 323309000000.0, 320114000000.0, 313757000000.0, 306402000000.0, 297432000000.0, 286611000000.0, 273923000000.0, 259322000000.0, 245700000000.0, 232601000000.0, 220717000000.0, 210418000000.0, 203287000000.0, 196411000000.0, 188635000000.0, 179902000000.0, 170976000000.0, 161608000000.0, 152068000000.0, 142026000000.0, 131790000000.0, 121543000000.0, 110950000000.0, 100381000000.0, 90044565360L, 81152701965L, 73964839833L, 68329114937L, 64012209883L, 60854018893L, 57737333014L, 54499791284L, 51247351987L, 48009715283L, 44998120134L, 41809412623L, 38817126728L, 36194284425L, 33406162511L, 30588882803L, 27830736669L, 25085074630L, 22567357411L, 20613936490L, 19062439597L, 17870386623L, 16989551627L, 16085451346L, 15215185507L, 14448962125L, 13622231140L, 12780286167L, 11933700854L, 11088381176L, 10167842796L, 9357621874L, 8553211257L, 7693166716L, 6956024965L, 6317584592L, 5794891459L, 5393044012L, 5128138617L, 4897260912L, 4664637302L, 4423283431L, 4181445267L, 3935433440L, 3690672341L, 3455806489L, 3234250875L, 3011075366L, 2781905319L, 2547010352L, 2312416498L, 2091725822, 1919097586, 1788865556, 1690634229, 1616309189, 1543679941, 1466469018, 1387742326, 1313178557, 1243632060, 1171658691, 1097666852, 1024029504, 945514365.4, 866672687.4, 789529700.6, 719008897.1, 661896491.2, 614633776.5, 576569681, 546591197.4, 522222923.6, 496976655.7, 473766930.8, 449967283, 427488104.8, 403787049.6, 379800651.6, 355024095.1, 330961852.9, 307420099.7, 284995504, 262462262.6, 241169195.9, 222440545.3, 207388553.2, 195862508.4, 186719009.1, 178771592.2, 169819209.7, 161456184.1, 152609575.3, 144701167.9, 136779479, 130170435.7, 122815897.9, 115167160.4, 106924155.5, 98540616.32, 90755099.37, 84240494.37, 78688794.12, 74390189.94, 71654870.08, 68361743.43, 64876349.6, 61972968.96, 58852107.45, 55320711.49, 52650226.85, 50259739.9, 47472608.34, 45334010.52, 43577954.91, 40758121.85, 37736707.66, 35039884.82, 30357638.41] AD_Oct1989_edep=[(AD_Oct1989[i]*33.5) for i in range(len(AD_Oct1989)) ] AD_Oct1989_ionrate=[(AD_Oct1989[i]/sepduration) for i in range(len(AD_Oct1989)) ] #Photoionization rates from Peter et al 2014 #these ionrates are in 10log peter_photoionrateslog10=[1.01481, 1.05926, 1.08889, 1.1037, 1.13333, 1.14815, 1.17778, 1.19259, 1.20741, 1.23704, 1.25185, 1.28148, 1.31111, 1.34074, 1.38519, 1.41481, 1.45926, 1.5037, 1.53333, 1.57778, 1.60741, 1.66667, 1.6963, 1.8, 1.87407, 1.94815, 2.03704, 2.0963, 2.14074, 2.21481, 2.27407, 2.34815, 2.42222, 2.4963, 2.55556, 2.64444, 2.73333, 2.82222, 2.86667, 2.95556, 3.07407, 3.16296, 3.23704, 3.41481, 3.56296, 3.68148, 3.78519, 3.88889, 4.02222, 4.18519, 4.34815, 4.43704, 4.51111, 4.58519, 4.64444, 4.7037, 4.76296, 4.77778, 4.79259, 4.79259, 4.77778, 4.74815, 4.7037, 4.64444, 4.6, 4.55556, 4.51111, 4.45185, 4.39259, 4.33333, 4.27407, 4.2, 4.12593, 4.06667, 4.02222, 3.99259, 3.96296, 3.94815, 3.91852, 3.9037, 3.87407, 3.84444, 3.82963, 3.8, 3.77037, 3.71111, 3.65185, 3.56296, 3.45926, 3.32593, 3.22222, 3.1037, 3.0, 2.83704, 2.71852, 2.62963, 2.48148, 2.39259, 2.31852, 2.22963, 2.14074, 2.05185, 1.96296, 1.85926, 1.75556, 1.66667, 1.57778, 1.47407, 1.32593, 1.22222, 1.08889, 1.04444, 1.01481] peter_altitude=[247.494, 246.6, 245.261, 244.367, 243.474, 242.581, 241.687, 240.794, 240.347, 239.007, 238.561, 237.667, 236.328, 234.988, 233.648, 231.414, 230.074, 228.734, 227.395, 225.608, 223.821, 222.035, 220.248, 216.228, 213.102, 209.975, 207.295, 204.615, 202.829, 200.596, 197.916, 195.682, 192.556, 189.876, 187.643, 185.409, 182.283, 179.603, 178.263, 176.03, 172.457, 170.67, 168.437, 164.864, 161.737, 159.95, 158.164, 156.377, 154.144, 151.464, 148.784, 146.998, 146.104, 144.764, 143.871, 142.084, 141.191, 139.851, 139.404, 138.065, 135.831, 135.385, 134.491, 133.598, 133.151, 132.258, 131.811, 131.365, 130.918, 130.471, 130.025, 129.578, 128.685, 128.238, 127.792, 127.792, 127.345, 126.452, 126.005, 125.112, 124.665, 123.772, 123.325, 122.878, 122.432, 121.985, 121.538, 120.645, 119.752, 119.305, 118.859, 118.412, 117.52, 117.072, 117.072, 116.625, 116.179, 116.179, 115.732, 115.285, 114.839, 114.839, 114.392, 113.945, 113.499, 113.052, 113.052, 112.605, 111.712, 111.712, 111.266, 110.819, 110.372] peter_photoionrates=[math.pow(10,peter_photoionrateslog10[i]) for i in range (len(peter_photoionrateslog10)) ] #Previous studies# ###Dubach### altitude_dubach=[83.1054, 82.2596, 81.7432, 80.9013, 80.0581, 78.7169, 77.8737, 77.0279, 76.3547, 75.6788, 74.5023, 73.9911, 72.9846, 71.8067, 71.1322, 70.4577, 69.2877, 68.1217, 66.7897, 65.961, 65.131, 64.3932, 63.9768, 62.99, 61.8372, 61.1837, 60.3655, 59.5526, 59.0638, 58.5789, 57.928, 57.119, 56.314, 55.5063] ionrate_dubach=[1.00538, 1.28994, 1.65562, 2.03251, 2.53218, 3.29533, 4.10545, 5.26746, 6.18844, 7.48755, 9.74583, 11.7938, 14.4761, 19.1213, 22.7975, 27.1805, 32.8693, 38.0326, 44.6515, 47.3169, 50.8846, 51.5941, 51.5762, 50.7703, 50.7088, 47.7778, 45.0084, 39.9764, 37.6723, 33.9682, 31.0768, 26.4107, 21.4761, 17.985] ###Boroucki### altitude_boroucki=[79.422, 78.1402, 77.0414, 75.7628, 73.3751, 71.1765, 68.9666, 66.5727, 63.7934, 61.5515, 60.7991, 60.0457, 59.2891, 58.7186, 58.3414, 57.9569, 56.6422, 56.4428, 55.5012, 55.1178, 54.5453, 53.5986, 53.21, 52.0803, 51.1377, 50.1951, 48.3213, 46.0753, 44.5787, 43.0832, 41.9586, 40.0868, 38.03, 36.5314, 34.2864, 32.4136, 30.1656, 28.4819, 26.236, 23.987, 22.4883, 18.7437, 15.9324, 14.0606, 12.7459, 11.0602, 9.56149, 8.25195, 7.12641, 6.1869, 4.87737, 3.56475, 2.62422, 1.31365, 0.376212] ionrate_boroucki=[3.34184, 4.69447, 6.26683, 9.26381, 15.4238, 27.023, 39.2732, 59.0497, 71.1866, 61.0911, 53.3258, 45.7632, 37.3212, 30.4365, 27.9573, 22.7999, 18.594, 14.9084, 12.3666, 10.2582, 8.08629, 6.16124, 4.69447, 3.76397, 3.06962, 2.50336, 2.00716, 1.60932, 1.40476, 1.24721, 1.08868, 0.903066, 0.761937, 0.642862, 0.524273, 0.427559, 0.331356, 0.284364, 0.227999, 0.173721, 0.146572, 0.0991539, 0.0705844, 0.0585501, 0.0477493, 0.0396083, 0.0334184, 0.0296705, 0.0254627, 0.0218516, 0.0194009, 0.016369, 0.0138109, 0.0120554, 0.0107034] ###Upadhyayay### altitude_upad=[0, 2.38775510204082, 4.77551020408163, 7.16326530612245, 9.55102040816327, 11.9387755102041, 14.3265306122449, 16.7142857142857, 19.1020408163265, 21.4897959183673, 23.8775510204082, 26.265306122449, 28.6530612244898, 31.0408163265306, 33.4285714285714, 35.8163265306122, 38.2040816326531, 40.5918367346939, 42.9795918367347, 45.3673469387755, 47.7551020408163, 50.1428571428571, 52.530612244898, 54.9183673469388, 57.3061224489796, 59.6938775510204, 62.0816326530612, 64.4693877551021, 66.8571428571429, 69.2448979591837, 71.6326530612245, 74.0204081632653, 76.4081632653061, 78.7959183673469, 81.1836734693878, 83.5714285714286, 85.9591836734694, 88.3469387755102, 90.734693877551, 93.1224489795918, 95.5102040816327, 97.8979591836735, 100.285714285714, 102.673469387755, 105.061224489796, 107.448979591837, 109.836734693878, 112.224489795918, 114.612244897959, 117] protons_upad=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0582219317770242, 0.983802543159385, 18.2568215396693, 106.29239236701, 272.401525447983, 653.534468507615, 892.359420375615, 942.088916773812, 867.339163197218, 709.141515608641, 536.135823251239, 387.025791266918, 252.575917234365, 176.799586919844, 118.27709497018, 72.7408961306672, 41.370945098076, 24.389216757934, 15.0168829616447, 7.76848057831586, 3.45967693839115, 1.85035188358631, 0.987236819373512, 0.544928468172109, 0.289623963250251, 0.138946999766423, 0.0806145535537336, 0.0401827906173959, 0.0158125327564783, 0] muons_upad=[175.105, 191.494440112899, 206.767730834902, 220.710807618506, 235.42887188844, 250.070503653219, 259.615108465842, 260.66548108998, 258.957825871563, 255.853063830465, 251.620496619538, 246.529425891634, 239.625848653946, 226.835341491189, 210.798114144224, 194.571102082328, 176.414948558318, 157.786272367194, 141.067969969711, 124.855988163736, 106.395766246134, 87.3807374086759, 72.2275604372746, 59.4154571336199, 46.436536594224, 33.244900250847, 22.1226250346854, 13.2237385552601, 5.98554862764091, 3.0003081041862, 1.34787030876326, 0.534245525772378, 0.198167125788632, 0.0690496821860641, 0.0319627703527496, 0.0117945659429756, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] total_upad=[protons_upad[i]+muons_upad[i] for i in range(len(protons_upad)) ] #Set size of output plot fig = pl.figure(figsize=(10.5, 7.7)) pl.rcParams.update({'font.size': 16}) #fig = pl.figure(figsize=(6.5, 7.7)) ####Solmax#### # #pl.plot(solmax_protons,altitude,'b--',label='Protons (Z=1)',linewidth=2) #pl.plot(solmax_HZE,altitude,'r--',label='HZE (Z=2-28)',linewidth=2) #pl.plot(solmax_total,altitude,'k-',label='Total',linewidth=2) #pl.title('Ion production rate - Solar Maximum') ####Solmin##### #pl.plot(solmin_protons,altitude,'b--',label='This study: Protons (Z=1)',linewidth=2) #pl.plot(solmin_HZE,altitude,'r--',label='This study: HZE (Z=2-28)',linewidth=2) #pl.plot(solmin_total,altitude,'k-',label='This study (Z=1-28)',linewidth=2) #pl.title('Ion production rate - Solar Minimum') ### ##pl.plot(total_upad,altitude_upad,'k:',label='Upadhyay & Singh, 1995',linewidth=2) ##pl.plot(ionrate_dubach,altitude_dubach,'r--',label='Dubach et al.,1974',linewidth=2) ##pl.plot(ionrate_boroucki,altitude_boroucki,'b--',label='Boroucki et al.,1982',linewidth=2) ####Worstweek#### # #pl.plot(worstweek_protons,altitude,'b--',label='Protons (Z=1)',linewidth=2) #pl.plot(worstweek_HZE,altitude,'r--',label='HZE (Z=2-28)',linewidth=2) #pl.plot(worstweek_total,altitude,'k-',label='Total',linewidth=2) #pl.title('Ion production rate - October 1989 Worst Week') ####Worst5min#### #pl.plot(worst5min_protons,altitude,'b--',label='Protons (Z=1)',linewidth=2) #pl.plot(worst5min_HZE,altitude,'r--',label='HZE (Z=2-28)',linewidth=2) #pl.plot(worst5min_total,altitude,'k-',label='Total',linewidth=2) #pl.title('Ion production rate - October 1989 Peak') ####All + Peter et al photoionization flxues#### #pl.plot(solmax_total,altitude,'k-',label='GCR:Solar Maximum',linewidth=2) #pl.plot(solmin_total,altitude,'k--',label='GCR:Solar Minimum',linewidth=2) #pl.plot(worstweek_total,altitude,'b-.',label='SEP:October 1989, Worst Week',linewidth=3) #pl.plot(worst5min_total,altitude,'r--',label='SEP:October 1989, Peak',linewidth=2) #pl.plot(peter_photoionrates,peter_altitude,'g:',label='EUV/X-Ray',linewidth=3) ####All + extreme carrington SEP + Peter et al photoionization flxues#### pl.plot(Carrington_Aug1972_ionrate,altitude,color='orange',label='SEP:Carrington (Aug1972 shape)', linewidth=2,linestyle='--') #pl.plot(Carrington_Feb1956_ionrate,altitude,color='blue',label='SEP:Carrington (Feb1956 shape)', linewidth=2,linestyle='--') #pl.plot(Carrington_Oct1989_ionrate,altitude,color='green',label='SEP:Carrington (Oct1989 shape)', linewidth=2,linestyle='--') #pl.plot(AD_Aug1972_ionrate,altitude,color='orange',label='SEP:775 AD (Aug1972 shape)', linewidth=3,linestyle='-') pl.plot(AD_Feb1956_ionrate,altitude,color='blue',label='SEP:775 AD (Feb1956 shape)', linewidth=3,linestyle='-') #pl.plot(AD_Oct1989_ionrate,altitude,color='green',label='SEP:775 AD (Oct1989 shape)', linewidth=3,linestyle='-') pl.plot(worstweek_total,altitude,color='red',label='SEP:October 1989, Worst Week',linewidth=2, linestyle='--') #pl.plot(worst5min_total,altitude,color='red',label='SEP:October 1989, Peak',linewidth=2, linestyle='-') #pl.plot(solmax_total,altitude,'k--',label='GCR:Solar Maximum',linewidth=2) pl.plot(solmin_total,altitude,'k-',label='GCR:Solar Minimum',linewidth=2) pl.plot(peter_photoionrates,peter_altitude,'g:',label='EUV/X-Ray',linewidth=3) ####Extreme carrington_edep_comparison#### #pl.plot(Carrington_Aug1972_edep,altitude,color='orange',label='Carrington Aug1972', linewidth=2,linestyle='--') #pl.plot(Carrington_Feb1956_edep,altitude,color='blue',label='Carrington Feb1956', linewidth=2,linestyle='--') #pl.plot(Carrington_Oct1989_edep,altitude,color='green',label='Carrington Oct1989', linewidth=2,linestyle='--') # #pl.plot(AD_Aug1972_edep,altitude,color='orange',label='775 AD Aug1972', linewidth=3,linestyle='-') #pl.plot(AD_Feb1956_edep,altitude,color='blue',label='775 AD Feb1956', linewidth=3,linestyle='-') #pl.plot(AD_Oct1989_edep,altitude,color='green',label='775 AD Oct1989', linewidth=3,linestyle='-') # # #pl.plot(solmax_total_edep,altitude,color='black',label='GCR Solar max', linewidth=2,linestyle='--') #pl.plot(solmin_total_edep,altitude,color='black',label='GCR Solar min', linewidth=2,linestyle='-') #pl.plot(worstweek_total_edep,altitude,color='gray',label='Worst Week SEP', linewidth=2,linestyle='--') #pl.xlim(1e-3,1e2) pl.title('Ion production rate - Summary') #pl.title('Ion production rate') #pl.title('Total energy deposition over 20 hours') pl.yscale('linear') pl.xscale('log') pl.ylabel("Altitude [km]") pl.xlabel("Ion production rate [#/cm3-sec]") #pl.xlabel("Total energy deposition [eV/cm3]") pl.legend() pl.legend(loc='upper left') #pl.xlim(1e-8,1.2e3) pl.xlim(1e-8,1.5e6) #pl.ylim(0,150) pl.ylim(0,250) pl.minorticks_on() pl.show()