[1]:
import astropy.units as u
import astropy.table
from astropy.visualization import quantity_support
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
import fiasco
from synthesizAR.atomic import equilibrium_ionization, non_equilibrium_ionization, effective_temperature
from synthesizAR.interfaces.ebtel import run_ebtel, read_xml
import pydrad.configure
import pydrad.configure.data
import pydrad.parse
import pydrad.visualize
NEI Simulations#
For the relevant ions in our bandpass, simulate wither non-equilibrium ionization is an important effect for us to consider when thinking about whether this lines will be observable.
Impulsive Heating Simulations#
[5]:
ebtel_config = read_xml('/Users/wtbarnes/Documents/codes/ebtelPlusPlus/config/ebtel.example.cfg.xml')
[6]:
ebtel_config
[6]:
{'total_time': 5000.0,
'tau': 1.0,
'tau_max': 1e+300,
'loop_length': 4000000000.0,
'saturation_limit': 1.0,
'force_single_fluid': False,
'use_c1_loss_correction': True,
'use_c1_grav_correction': True,
'use_flux_limiting': False,
'calculate_dem': False,
'save_terms': False,
'use_adaptive_solver': False,
'output_filename': 'ebtel++_results_file.txt',
'adaptive_solver_error': 1e-06,
'adaptive_solver_safety': 0.5,
'c1_cond0': 2.0,
'c1_rad0': 0.6,
'helium_to_hydrogen_ratio': 0.075,
'surface_gravity': 1.0,
'dem': OrderedDict([('use_new_method', True),
('temperature', {'bins': 451, 'log_min': 4, 'log_max': 8.5})]),
'heating': OrderedDict([('background', 3.5e-05),
('partition', 1.0),
('events',
[{'event': {'magnitude': 0.1,
'rise_start': 0.0,
'rise_end': 50.0,
'decay_start': 50.0,
'decay_end': 100.0}},
{'event': {'magnitude': 0.15,
'rise_start': 200.0,
'rise_end': 250.0,
'decay_start': 250.0,
'decay_end': 300.0}},
{'event': {'magnitude': 0.05,
'rise_start': 1000.0,
'rise_end': 1050.0,
'decay_start': 1050.0,
'decay_end': 1100.0}}])])}
[87]:
total_time = 1.3e4 * u.s
loop_length = 40 * u.Mm
duration = 200 * u.s
t_0 = 100 * u.s
wait_time = 1500 * u.s
number_events = 1
heating_rate = 0.075 * u.Unit('erg cm-3 s-1')
heating_rate_bg = 2.75e-5 * u.Unit('erg cm-3 s-1')
start_times = t_0 + (duration + wait_time) * np.arange(number_events)
ebtel_config['total_time'] = total_time.to_value('s')
ebtel_config['loop_length'] = loop_length.to_value('cm')
ebtel_config['use_flux_limiting']= True
ebtel_config['use_adaptive_solver'] = True
ebtel_config['saturation_limit'] = 1
ebtel_config['adaptive_solver_error'] = 1e-10
#config['adaptive_sover_safety'] = 0.1
events = []
for st in start_times:
events.append(
{'event': {'magnitude': heating_rate.to_value('erg cm-3 s-1'),
'rise_start': st.to_value('s'),
'rise_end': (st + duration/2).to_value('s'),
'decay_start': (st + duration/2).to_value('s'),
'decay_end': (st + duration).to_value('s')}}
)
ebtel_config['heating']['events'] = events
ebtel_config['heating']['background'] = heating_rate_bg
[88]:
simulation = run_ebtel(ebtel_config, '/Users/wtbarnes/Documents/codes/ebtelPlusPlus/')
time = simulation['time'] * u.s
heat = simulation['heat'] * u.Unit('erg cm-3 s-1')
electron_temperature = simulation['electron_temperature'] * u.K
density = simulation['density'] * u.Unit('cm-3')
[83]:
simulation.keys()
[83]:
dict_keys(['time', 'electron_temperature', 'ion_temperature', 'density', 'electron_pressure', 'ion_pressure', 'velocity', 'heat'])
[84]:
electron_temperature[0]
[84]:
$1007727 \; \mathrm{K}$
[85]:
density[0] / 1e9
[85]:
$0.1771306 \; \mathrm{\frac{1}{cm^{3}}}$
[86]:
fig = plt.figure(figsize=(10,10))
ax = fig.add_subplot(221)
ax.plot(time, heat)
ax = fig.add_subplot(222)
ax.plot(time, electron_temperature)
ax = fig.add_subplot(223)
ax.plot(time, density)
ax = fig.add_subplot(224)
ax.plot(electron_temperature, density)
ax.set_yscale('log')
ax.set_xscale('log')
Identify Ions#
[13]:
line_list = astropy.table.QTable.read('../data/moxsi-line-list.asdf')
WARNING: UnitsWarning: The unit 'Angstrom' has been deprecated in the VOUnit standard. Suggested: 0.1nm. [astropy.units.format.utils]
[14]:
line_list[line_list.argsort(keys='intensity_scaled (coronal)_flare_ext', reverse=True)][['ion name', 'wavelength']].pprint_all()
ion name wavelength
Angstrom
-------- ------------------
Fe XVII 15.012999534606934
Fe XVII 17.051000595092773
Fe XXI 12.281999588012695
Fe XVII 17.09600067138672
Fe XVII 16.77560043334961
Fe XXIII 11.737000465393066
Fe XVIII 14.203900337219238
O VIII 18.967100143432617
Fe XXII 11.767499923706055
Fe XXIV 11.170900344848633
Fe XVIII 16.07200050354004
Fe XIX 13.524900436401367
Fe XXIV 10.619000434875488
Fe XX 12.82699966430664
Fe XVII 15.26200008392334
Fe XXIII 12.16100025177002
Fe XVIII 14.208800315856934
Mg XII 8.41919994354248
Fe XX 12.845399856567383
Si XIV 6.1803998947143555
Fe XVIII 17.62179946899414
Fe XVIII 16.0049991607666
Si XIII 6.647900104522705
Fe XXIII 10.979999542236328
Fe XVIII 15.622099876403809
Fe XXIV 11.029000282287598
O VIII 18.97249984741211
Fe XIX 13.506400108337402
Fe XVIII 14.373000144958496
Fe XXIV 11.426300048828125
Fe XIX 16.109899520874023
Fe XXI 13.506999969482422
Fe XXIV 10.663000106811523
Fe XIX 13.798999786376953
Mg XI 9.168700218200684
Fe XXIII 11.017999649047852
Fe XVIII 14.536999702453613
Fe XXII 12.742600440979004
Mg XII 8.424599647521973
Fe XVIII 15.82800006866455
Fe XXI 12.395000457763672
Ne X 12.132100105285645
Fe XXI 14.008000373840332
Fe XX 12.812000274658203
Si XIV 6.185800075531006
Fe XIX 15.081000328063965
Fe XIX 13.471400260925293
Fe XVII 12.124099731445312
Si XIII 6.740300178527832
Fe XVIII 14.258000373840332
Fe XIX 13.455699920654297
Fe XIX 13.690500259399414
Fe XVIII 16.165700912475586
Fe XVII 12.263999938964844
Si XII 44.160301208496094
Ni XIX 12.4350004196167
Fe XIX 14.668999671936035
Fe XXI 12.821999549865723
Fe XXV 1.8503999710083008
Fe XVI 63.71099853515625
Fe XVII 13.824999809265137
Si XII 40.9109992980957
Fe XXI 9.475000381469727
Fe XX 12.951000213623047
Fe XXIV 11.260600090026855
Fe XX 14.263400077819824
C VI 33.73419952392578
Fe XXII 11.933600425720215
Fe XVIII 15.86989974975586
Ni XIX 14.039799690246582
Fe XX 12.905500411987305
Fe XVIII 14.550999641418457
Mg XI 9.314299583435059
Fe XX 13.052000045776367
Fe XVI 66.35700225830078
Fe XXII 13.230999946594238
Ne X 12.137499809265137
Fe XX 13.404999732971191
Fe XIX 13.741299629211426
Fe XIX 13.64799976348877
Fe XIX 14.991999626159668
O VIII 16.00550079345703
Fe XIX 16.2721004486084
Fe XXII 11.535799980163574
Fe XVIII 15.869999885559082
Fe XVIII 14.352499961853027
Si XII 44.01860046386719
Ni XIX 13.777700424194336
Fe XXII 12.935999870300293
Fe XXII 12.19260025024414
Ni XIX 14.07409954071045
Fe XXI 11.975000381469727
Fe XXII 11.920700073242188
Fe XIX 10.815999984741211
Fe XIX 13.721199989318848
Fe XIX 14.738100051879883
Si XII 45.69060134887695
Fe XVIII 11.420000076293945
Fe XXIII 11.442000389099121
Fe XVII 15.45300006866455
Fe XX 14.763699531555176
Fe XVIII 11.524999618530273
Fe XVI 66.2490005493164
Fe XVI 54.709999084472656
Fe XXII 11.499300003051758
Fe XXIII 8.814000129699707
Fe XIX 13.429699897766113
Fe XVII 11.25
Fe XXII 8.976699829101562
Si XII 40.95100021362305
Fe XX 13.780599594116211
Fe XXI 13.255000114440918
Fe XX 13.14050006866455
Fe XXIV 7.984099864959717
Fe XVI 62.87099838256836
C VI 33.739601135253906
Fe XVIII 14.418899536132812
Fe XXI 12.49899959564209
Fe XXIV 11.187899589538574
Fe XX 12.581000328063965
Fe XXIII 8.303000450134277
Fe XX 13.1427001953125
Fe XXIV 8.316399574279785
Fe XVIII 14.258000373840332
S XV 5.038700103759766
Fe XX 14.913299560546875
Ni XX 11.831999778747559
Ni XIX 12.656000137329102
Fe XV 69.68199920654297
Fe XVI 50.361000061035156
Fe XVIII 16.02560043334961
Ni XXIII 10.310999870300293
Ne IX 13.44729995727539
Si XIII 6.688199996948242
N VII 24.779199600219727
Ca XVIII 19.789100646972656
Fe XVIII 16.305700302124023
Fe XVIII 14.344099998474121
Fe XXI 12.42199993133545
Mg XII 7.105800151824951
Fe XIX 13.768600463867188
Fe XX 14.3125
Si XIII 5.680699825286865
Fe XX 12.426300048828125
Fe XX 14.402899742126465
Mg X 63.29499816894531
Fe XXII 11.795499801635742
S XIV 32.560001373291016
Fe XVIII 14.418999671936035
Fe XVIII 14.579999923706055
Fe XVII 11.128999710083008
Fe XIX 13.017999649047852
O VII 21.601499557495117
Fe XIX 14.991999626159668
Fe XX 12.581000328063965
Fe XVII 46.30670166015625
Fe XVIII 15.766400337219238
Fe XVIII 11.32610034942627
Fe XIX 13.843700408935547
O VIII 16.00670051574707
Fe XVIII 13.961999893188477
Fe XVI 15.313899993896484
Si XIV 5.216800212860107
Fe XXIII 11.336999893188477
Fe XX 12.426300048828125
Fe XX 13.516599655151367
Fe XXI 12.043999671936035
Fe XXII 11.510100364685059
Fe XXV 1.8681999444961548
Mg X 57.875999450683594
S XIV 30.427000045776367
Fe XVIII 13.396900177001953
S XVI 4.727399826049805
Fe XXI 12.64900016784668
Fe XIX 10.632599830627441
Fe XVI 15.211400032043457
Fe XVIII 11.32610034942627
Fe XX 13.843999862670898
Fe XIX 13.77400016784668
Fe XXII 11.885199546813965
Fe XVIII 15.396699905395508
Fe XXII 12.5068998336792
Fe XVI 54.125999450683594
Fe XXII 12.090700149536133
Si XII 45.520599365234375
Fe XXII 12.214300155639648
Fe XX 14.962200164794922
Fe XVIII 11.524999618530273
Fe XXII 13.21030044555664
Ca XVIII 18.690900802612305
Fe XXI 8.572999954223633
Fe XVIII 10.522899627685547
Ni XXI 11.318599700927734
Fe XIX 16.34000015258789
Fe XXI 13.052000045776367
Fe XXIV 32.47809982299805
Fe XVII 10.770000457763672
Fe XIX 9.847999572753906
Fe XVIII 14.77180004119873
Fe XXI 13.352999687194824
Fe XXIV 8.2322998046875
Fe XXIII 35.27399826049805
Fe XX 12.965999603271484
Fe XIX 13.671600341796875
Fe XX 12.982000350952148
Fe XIX 14.041999816894531
Fe XXI 13.574000358581543
Fe XXIII 8.906000137329102
Ni XX 13.309000015258789
Fe XX 14.45300006866455
Fe XIX 13.554800033569336
Fe XX 14.44849967956543
Al XIII 7.170899868011475
Fe XXI 35.63100051879883
Fe XVIII 14.487099647521973
Fe XIX 14.01729965209961
Mg XI 7.850500106811523
Fe XXII 12.05739974975586
Fe XIX 13.673299789428711
Ni XX 11.841400146484375
Fe XXIV 7.99429988861084
Fe XX 15.515000343322754
Fe XXII 12.388799667358398
Fe XXI 12.204000473022461
Fe XXI 13.069999694824219
Fe XX 14.932000160217285
Ni XXV 9.967000007629395
Fe XIX 13.937999725341797
Fe XXI 9.805999755859375
Fe XXI 13.782999992370605
Fe XXII 11.44159984588623
Fe XVIII 15.449700355529785
Fe XX 13.972299575805664
Fe XIX 12.923999786376953
Mg XI 9.231200218200684
Fe XIX 16.378700256347656
Fe XVIII 10.433799743652344
O VIII 15.175999641418457
Fe XVIII 13.318900108337402
Fe XVII 46.29999923706055
Fe XIX 14.20259952545166
Fe XIX 13.635700225830078
Fe XX 12.986000061035156
Fe XVIII 18.10700035095215
Fe XX 13.267000198364258
Fe XVII 16.335599899291992
Ni XX 14.470000267028809
Mg X 63.152000427246094
Fe XX 15.0600004196167
Fe XVI 50.564998626708984
Fe XXV 1.8595000505447388
Fe XX 13.267000198364258
Ca XVIII 19.642000198364258
Ca XVII 21.1560001373291
Fe XX 13.393600463867188
Fe XVIII 14.469599723815918
Fe XXIII 11.845999717712402
Fe XIX 14.932000160217285
Fe XXI 13.100000381469727
Ni XXIV 9.946900367736816
Fe XIX 14.017900466918945
Fe XXI 13.972999572753906
Fe XXI 11.937999725341797
S XIV 32.41600036621094
Fe XXI 12.567999839782715
Fe XIX 13.253999710083008
Ne X 10.238499641418457
Fe XXI 12.822999954223633
Fe XIX 10.684000015258789
Fe XXIV 8.372699737548828
Fe XVIII 11.310600280761719
Fe XVIII 14.889599800109863
N VII 24.78459930419922
Fe XXI 9.159000396728516
Fe XVIII 14.70580005645752
Fe XIX 12.923999786376953
Fe XXV 1.855299949645996
Fe XVIII 14.14430046081543
Fe XIX 10.654999732971191
Mg XII 7.106900215148926
Al XII 7.757299900054932
Fe XVIII 14.45300006866455
Fe XXIII 11.298999786376953
Fe XXIV 32.37710189819336
Fe XVIII 13.424099922180176
Fe XX 14.830400466918945
Fe XVIII 43.25550079345703
Ne IX 13.699000358581543
Fe XXII 11.975000381469727
Fe XX 13.505000114440918
Mg X 65.84500122070312
Fe XIX 40.39039993286133
Fe XVIII 43.294700622558594
Fe XXIII 11.883999824523926
Ni XX 13.255999565124512
S XV 5.101500034332275
Fe XIX 16.747699737548828
Fe XXI 14.02400016784668
Fe XXI 12.461999893188477
Fe XVIII 18.52589988708496
O VII 22.097700119018555
Fe XVIII 13.35509967803955
Fe XIX 13.461999893188477
S XIV 33.54949951171875
Fe XIX 13.72089958190918
Fe XVII 13.890000343322754
Fe XX 14.805999755859375
Ca XVIII 20.218799591064453
Fe XVII 50.60599899291992
Fe XIX 15.195799827575684
Fe XXII 11.380599975585938
Fe XXII 11.604999542236328
Fe XVIII 13.46399974822998
Fe XVII 10.503999710083008
Si XIV 5.2179999351501465
Na XI 10.023200035095215
Fe XX 12.824199676513672
Fe XX 12.978500366210938
Ni XX 11.961000442504883
Fe XVII 16.003799438476562
Fe XX 13.292400360107422
Ca XIX 3.177299976348877
Fe XIX 13.67199993133545
Fe XX 13.090999603271484
Ni XX 12.112000465393066
Fe XXI 13.178999900817871
Fe XVIII 14.610400199890137
Fe XIX 40.86800003051758
Mg X 57.91999816894531
Fe XVIII 14.968700408935547
Fe XX 10.986599922180176
Fe XX 12.565999984741211
S XIV 30.4689998626709
Fe XIX 13.45199966430664
Si XII 31.012100219726562
S XVI 4.732800006866455
Fe XIX 14.9350004196167
Fe XXIV 7.167600154876709
Fe XXI 12.77299976348877
Ni XX 12.927000045776367
Fe XX 13.253999710083008
Fe XXII 12.19320011138916
Ni XXI 11.302000045776367
Fe XX 13.114700317382812
Fe XVII 10.656999588012695
Fe XXII 12.324700355529785
Fe XX 13.332900047302246
Fe XXII 11.747200012207031
Fe XIX 16.77549934387207
Fe XV 59.404998779296875
Fe XXI 9.541999816894531
Fe XX 13.525199890136719
Fe XX 12.600000381469727
Fe XIX 15.163100242614746
Ca XVIII 18.73189926147461
Fe XXIV 31.968299865722656
Fe XVIII 14.670499801635742
Fe XX 13.545000076293945
Fe XIX 12.923999786376953
Fe XIX 13.492300033569336
Fe XX 13.12339973449707
Fe XIX 13.553999900817871
Fe XVII 57.709999084472656
Fe XXI 13.112000465393066
Fe XXI 9.972999572753906
Fe XXIV 30.74650001525879
Fe XXII 8.715399742126465
Fe XXIII 7.4720001220703125
Fe XXI 13.107999801635742
Fe XVII 16.237899780273438
Si XII 32.973201751708984
Fe XXI 12.774999618530273
Fe XXII 11.950599670410156
Fe XIX 13.401399612426758
Fe XXIII 11.33899974822998
Fe XXV 1.573199987411499
Fe XXI 12.805999755859375
Fe XVIII 10.526100158691406
Fe XXII 12.077799797058105
Fe XIX 14.294599533081055
Fe XX 11.006500244140625
Fe XXII 35.13759994506836
Fe XXII 9.240699768066406
Fe XXIII 12.437000274658203
Fe XIX 13.76200008392334
Fe XXI 12.50100040435791
Fe XIX 10.757499694824219
Fe XXIV 32.81589889526367
Fe XVII 50.2599983215332
Fe XIX 9.690500259399414
Fe XIX 14.280799865722656
Ni XXV 10.288000106811523
Fe XIX 16.32469940185547
Fe XVIII 16.08930015563965
Fe XXIII 7.882999897003174
Fe XX 14.812800407409668
Fe XXII 13.675700187683105
Fe XVIII 14.487099647521973
Fe XXII 11.511099815368652
Si XIII 5.404600143432617
Ni XIX 10.110199928283691
Fe XVII 68.51000213623047
Fe XXII 11.787699699401855
Fe XXIII 34.803001403808594
Fe XX 14.123299598693848
Fe XXI 12.532999992370605
Ca XVII 22.15999984741211
Al XIII 7.176300048828125
Fe XVIII 17.355899810791016
Fe XXI 14.529000282287598
Ni XVIII 51.03799819946289
Fe XIX 14.410599708557129
Fe XVII 13.15939998626709
Fe XVIII 10.3503999710083
Fe XX 15.062600135803223
Si XII 44.1775016784668
Fe XX 13.395299911499023
Fe XXI 12.79699993133545
Fe XIX 13.936200141906738
Fe XIX 10.63230037689209
Fe XXIII 11.718000411987305
Ar XVI 24.993999481201172
Mg XII 6.73769998550415
Fe XIX 13.48859977722168
C VI 28.465200424194336
Fe XX 13.090999603271484
Ni XIX 9.977100372314453
Fe XXII 11.820699691772461
Fe XXII 11.726400375366211
Fe XVIII 17.803300857543945
Fe XXI 14.26200008392334
Fe XXIV 1.8659000396728516
Fe XX 15.843700408935547
Ni XXV 9.329999923706055
Si XIII 36.43619918823242
Fe XXII 33.91559982299805
Fe XXII 11.668999671936035
Ni XXI 11.539299964904785
Fe XIX 15.041600227355957
O VIII 15.17650032043457
Si XI 49.222999572753906
Ca XV 23.023000717163086
Fe XVII 46.400001525878906
Fe XVI 40.153099060058594
Ni XVIII 52.720001220703125
Fe XV 52.9109992980957
Mg X 65.6729965209961
Ca XVI 21.450000762939453
Fe XVI 15.163000106811523
Fe XVI 46.71820068359375
Fe XVI 42.3036994934082
Fe XVI 15.109199523925781
Si XII 31.0226993560791
Si XII 32.88779830932617
Ni XVIII 44.365299224853516
Si XI 43.762001037597656
Si XI 52.29600143432617
Ni XVIII 52.6150016784668
S XIII 35.66699981689453
Fe XVI 46.661399841308594
Ne IX 11.546600341796875
C VI 28.466299057006836
Fe XVI 54.74700164794922
Ni XVII 55.15999984741211
Ni XVIII 50.255001068115234
Fe XVI 36.749000549316406
O VII 18.62700080871582
Fe XV 63.957000732421875
Al XI 52.44599914550781
N VII 20.909500122070312
Si XII 33.31669998168945
Cr XV 18.49690055847168
Fe XVI 39.82709884643555
Ne IX 13.553099632263184
Fe XVI 17.45319938659668
Al XI 48.297000885009766
Fe XVI 15.158499717712402
Mg X 44.04999923706055
Mg X 47.31019973754883
Fe XVI 15.537599563598633
Mg X 63.31100082397461
O VII 21.803600311279297
C VI 26.989599227905273
Fe XVI 41.932098388671875
Cr XV 21.152799606323242
S XIII 37.724998474121094
S XIII 32.242000579833984
Fe XV 46.2869987487793
Fe XVI 66.37699890136719
Fe XIV 58.9630012512207
Ca XIV 24.132999420166016
Si X 50.52399826049805
Cr XV 20.862899780273438
Fe XVI 35.105899810791016
Al XI 52.29899978637695
Cr XV 21.212099075317383
Mg X 44.04999923706055
Fe XVI 36.803001403808594
Al XI 54.38800048828125
Ca XIV 24.085899353027344
Fe XV 43.78099822998047
Fe XV 53.165000915527344
Mg X 47.22850036621094
Si X 50.691001892089844
Ni XVII 57.57899856567383
Ni XVII 47.6619987487793
Fe XVI 40.24489974975586
Si IX 55.4010009765625
Si VIII 61.06999969482422
Fe XVI 35.999000549316406
O VII 17.76799964904785
C V 40.267398834228516
Ca XI 30.447799682617188
Si IX 55.35599899291992
Fe XV 69.98699951171875
Mg X 47.887699127197266
Mg IX 62.750999450683594
S XII 36.39799880981445
Fe XV 49.4900016784668
Ca XIII 26.033000946044922
Si XI 46.4010009765625
Fe XV 56.20000076293945
Ca XI 35.212398529052734
Si VIII 61.01900100708008
Ca XI 35.57569885253906
Ca XI 35.676300048828125
N VI 28.78700065612793
Fe XIV 59.57899856567383
Fe XV 69.94100189208984
C V 41.471500396728516
Fe XV 64.87799835205078
Fe XIII 62.35300064086914
Fe XV 55.79199981689453
Si XI 37.340999603271484
Ca XIII 25.976999282836914
Ca XII 27.972999572753906
Mg IX 67.23899841308594
Si X 57.3650016784668
Fe XIV 67.14099884033203
Fe XIV 69.68499755859375
O VII 17.395999908447266
Ar XIII 29.319000244140625
Si IX 55.27199935913086
Fe XV 66.2300033569336
Si XI 33.51300048828125
Ca XII 32.65700149536133
Si VIII 69.63200378417969
Fe XIV 58.70899963378906
S XI 39.2400016784668
Fe XIV 58.888999938964844
Si X 61.895999908447266
Si IX 55.30500030517578
Ca XII 31.660999298095703
Si X 61.96799850463867
S X 42.54309844970703
N VI 29.534700393676758
Ca XIII 32.803001403808594
Fe XIII 62.974998474121094
Ca XII 27.972999572753906
Si VIII 60.98899841308594
Si XI 43.685001373291016
Ca XIII 26.7189998626709
Fe XIV 59.625999450683594
Ne VIII 67.38200378417969
Ca XII 32.280399322509766
Si XI 46.297000885009766
Si X 52.611000061035156
Si IX 61.64899826049805
Si X 50.702999114990234
Si IX 69.39800262451172
Fe XIV 55.4370002746582
Si VIII 61.91400146484375
S XII 40.59270095825195
Fe XIV 60.65999984741211
Si VIII 69.79000091552734
Si X 57.20800018310547
Fe XV 65.61499786376953
Si XI 37.9900016784668
Fe XIV 60.39699935913086
Si VII 69.66400146484375
Fe XV 44.981998443603516
Fe XII 66.2969970703125
S X 42.494998931884766
Si XI 49.17599868774414
Fe XIII 62.694000244140625
Ca XII 32.49879837036133
Si X 55.09600067138672
Ca XII 37.604000091552734
Ca XII 28.47800064086914
Si XI 50.47200012207031
S XI 39.29999923706055
Si IX 64.4489974975586
Fe XIV 66.76699829101562
Fe XIII 62.71699905395508
Si VII 68.14800262451172
Si XI 54.38199996948242
Ar XII 31.3882999420166
Fe XII 66.0469970703125
Ni XV 59.59000015258789
Fe XV 50.12099838256836
Fe XIV 48.486000061035156
Fe XIV 58.49599838256836
Fe XV 55.63399887084961
Al X 59.106998443603516
Si X 47.48789978027344
Si XI 46.26499938964844
Si IX 69.48999786376953
S VII 60.1619987487793
S XII 37.714500427246094
S VIII 63.88600158691406
Si IX 64.5780029296875
S XII 43.346099853515625
Fe XII 65.90499877929688
S XII 36.563499450683594
Fe XV 56.23699951171875
Fe XIV 58.75199890136719
Si X 52.48699951171875
Fe XIV 44.57979965209961
Si IX 60.512001037597656
Si X 60.11800003051758
Si X 39.442901611328125
S VIII 61.599998474121094
Fe XIV 55.38999938964844
Fe XIV 56.138999938964844
S XII 36.572898864746094
Mg IX 67.13500213623047
Fe XIV 61.185001373291016
Si IX 55.382999420166016
Si IX 66.70800018310547
Fe XIV 61.150001525878906
Ar XII 31.347000122070312
Si XI 49.051998138427734
Si IX 55.11600112915039
Ni XV 63.56999969482422
Fe XIV 58.10499954223633
Fe XIV 67.25399780273438
Fe XIV 55.2760009765625
Fe XIV 58.209999084472656
Mg VIII 59.15299987792969
S VIII 63.30400085449219
Si VIII 69.90499877929688
Fe XIV 60.19900131225586
S XI 39.31999969482422
Mg VIII 69.46700286865234
Si XI 46.314998626708984
Si IX 56.027000427246094
Si IX 55.23400115966797
S X 42.48500061035156
Al X 63.13399887084961
Si X 54.5989990234375
S X 54.0260009765625
Si X 39.551998138427734
Fe XIV 61.06999969482422
Si X 47.54499816894531
C V 34.972801208496094
Si IX 57.11399841308594
Ni XV 65.9000015258789
S IX 47.43299865722656
Mg IX 55.060001373291016
Fe XIV 58.49599838256836
Si IX 62.974998474121094
Al X 51.979000091552734
Si X 54.8849983215332
Ca XI 25.32670021057129
Si IX 43.97600173950195
Si IX 60.62900161743164
S XI 39.323001861572266
Fe XII 63.3849983215332
Mg VII 67.59400177001953
Fe XIV 56.68600082397461
Fe XIV 55.84199905395508
Ca XI 30.867000579833984
Fe XIV 62.51499938964844
Si X 54.56999969482422
Si IX 55.09400177001953
Fe XIV 59.202999114990234
Ar IX 48.73939895629883
S XI 44.92499923706055
C V 40.730201721191406
Si X 53.595001220703125
Fe XIV 69.55400085449219
Fe XIV 58.0369987487793
Si VIII 61.790000915527344
Mg IX 48.34000015258789
Si IX 55.606998443603516
Si IX 58.582000732421875
Si XI 46.40800094604492
S X 54.22800064086914
Ar IX 49.34510040283203
Al IX 61.069000244140625
Fe XIV 65.11900329589844
Si IX 57.27799987792969
Fe XIV 58.4379997253418
Si XI 35.446998596191406
Si IX 61.84400177001953
Si IX 61.50199890136719
Si VIII 63.731998443603516
Fe XIV 68.85199737548828
Fe XIV 58.91699981689453
Fe XIII 64.25800323486328
Fe XIV 56.9630012512207
S IX 47.249000549316406
Si IX 58.66400146484375
Ar IX 41.47600173950195
Fe XIII 64.35199737548828
Si IX 69.58499908447266
Fe XII 68.84300231933594
S VIII 63.80099868774414
Si VII 69.38500213623047
Ar XI 34.33000183105469
Mg IX 67.08999633789062
Ar IX 49.18550109863281
Fe XIV 56.542999267578125
S XI 47.71500015258789
S IX 49.11899948120117
Si IX 52.8380012512207
Si VII 69.60199737548828
Fe IX 68.50299835205078
Fe XIII 62.459999084472656
Fe XIII 59.24800109863281
Si X 49.700801849365234
Fe XIV 58.31800079345703
Si X 48.549400329589844
Fe XIV 55.505001068115234
Si IX 61.402000427246094
Mg VIII 59.03799819946289
S VIII 52.75600051879883
S IX 64.6969985961914
Si X 59.94499969482422
Si VII 68.37799835205078
Si XI 54.25
S VIII 54.11800003051758
S IX 46.37300109863281
Mg VIII 69.41500091552734
Si IX 64.65299987792969
Si IX 61.599998474121094
Al VIII 67.46399688720703
Si X 52.154998779296875
Fe XII 63.35499954223633
Mg IX 56.483001708984375
Fe XIII 62.154998779296875
Ni XI 62.72999954223633
Fe XIV 56.10300064086914
Si X 49.70100021362305
Fe XI 66.52799987792969
Si IX 60.709999084472656
Ar X 44.630001068115234
Al IX 60.895999908447266
S X 47.65439987182617
Mg IX 62.597999572753906
Fe XII 67.29100036621094
S X 47.79199981689453
Si VIII 63.90299987792969
Si IX 43.939998626708984
Mg VII 67.55699920654297
Si X 53.573001861572266
Fe XII 67.64099884033203
Fe XII 66.5260009765625
N VI 29.084299087524414
Fe XII 65.83300018310547
S IX 56.935001373291016
Mg IX 67.14099884033203
Si IX 58.73099899291992
S VIII 52.790000915527344
Si VII 69.58000183105469
S IX 54.17499923706055
Si IX 52.810001373291016
Fe XIII 63.191001892089844
Si X 48.43579864501953
Si IX 43.948001861572266
Si X 50.75899887084961
Ar X 43.27199935913086
Fe XII 67.97200012207031
S VIII 61.97800064086914
Si VII 69.86000061035156
Si VIII 61.2859992980957
Si IX 69.49199676513672
Fe XII 67.7020034790039
Si X 48.489601135253906
Si X 50.30500030517578
Si VIII 61.904998779296875
Si IX 61.696998596191406
S X 54.332000732421875
Si X 65.11199951171875
Fe XIII 64.08300018310547
Ar X 37.430999755859375
Ar X 44.07699966430664
Fe XII 66.01399993896484
Fe XIII 61.74300003051758
Fe XIII 59.57699966430664
Si VIII 61.89500045776367
Si X 52.06999969482422
Si IX 66.86000061035156
Al VIII 67.40699768066406
Si IX 59.5890007019043
Si IX 52.50400161743164
Si IX 66.84500122070312
Si VIII 63.2239990234375
Mg IX 67.24600219726562
Mg IX 51.65399932861328
Si IX 66.697998046875
Si IX 43.917999267578125
S X 43.00199890136719
Si X 54.70199966430664
Si IX 69.38400268554688
Si X 54.4640007019043
Ni XII 60.02000045776367
Al X 55.375
Si X 47.64609909057617
Fe XIII 63.632999420166016
Si IX 64.5790023803711
Si X 50.15399932861328
Si X 54.52199935913086
Si IX 66.92900085449219
Si IX 57.77799987792969
Si X 61.89500045776367
Si IX 52.20600128173828
Si X 50.31700134277344
Si X 50.33300018310547
Si X 47.04750061035156
Fe XII 68.98400115966797
Fe XII 66.1969985961914
Si X 56.68000030517578
Si IX 52.547000885009766
Si X 48.3484992980957
S IX 56.08100128173828
Si IX 66.8499984741211
Si VIII 62.85100173950195
Si X 54.66400146484375
Si VIII 64.32499694824219
Si X 50.71900177001953
Si IX 55.71799850463867
Si IX 65.22899627685547
Si IX 61.124000549316406
Si IX 56.518001556396484
Si IX 63.07500076293945
[15]:
hot_line_table = line_list[np.logical_and(line_list['max temperature_flare_ext'] > 3 * u.MK,
line_list['max temperature_flare_ext'] < 12 * u.MK)]
hot_line_table[hot_line_table.argsort(keys='intensity_scaled (coronal)_flare_ext', reverse=True)][
['ion name', 'wavelength', 'intensity_scaled (coronal)_flare_ext']].pprint_all()
ion name wavelength intensity_scaled (coronal)_flare_ext
Angstrom
-------- ------------------ ------------------------------------
Fe XVII 15.012999534606934 1.0
Fe XVII 17.051000595092773 0.7990830377280506
Fe XVII 17.09600067138672 0.7228795890158841
Fe XVII 16.77560043334961 0.5949798048043987
Fe XVIII 14.203900337219238 0.4896022804982188
O VIII 18.967100143432617 0.4122510037782132
Fe XVIII 16.07200050354004 0.33495904420023276
Fe XIX 13.524900436401367 0.3077585024598136
Fe XX 12.82699966430664 0.29560189950656485
Fe XVII 15.26200008392334 0.28855152685279967
Fe XVIII 14.208800315856934 0.27566060993420255
Fe XX 12.845399856567383 0.266082030490089
Fe XVIII 17.62179946899414 0.24004924504959327
Fe XVIII 16.0049991607666 0.2387232091747517
Fe XVIII 15.622099876403809 0.21002123717615723
O VIII 18.97249984741211 0.20605731780412825
Fe XIX 13.506400108337402 0.19494433901315555
Fe XVIII 14.373000144958496 0.19125242262590028
Fe XIX 16.109899520874023 0.18051193174733077
Fe XIX 13.798999786376953 0.15840642466078908
Mg XI 9.168700218200684 0.14946556922191895
Fe XVIII 14.536999702453613 0.14459757683859914
Fe XVIII 15.82800006866455 0.1331959162139273
Ne X 12.132100105285645 0.1296553562706885
Fe XX 12.812000274658203 0.1262192366900194
Fe XIX 15.081000328063965 0.11896888085285239
Fe XIX 13.471400260925293 0.11644694023944287
Fe XVII 12.124099731445312 0.11339930939334654
Si XIII 6.740300178527832 0.1106974765714088
Fe XVIII 14.258000373840332 0.11065064895685729
Fe XIX 13.455699920654297 0.10847066476691014
Fe XIX 13.690500259399414 0.10602674971923219
Fe XVIII 16.165700912475586 0.10542858056985398
Fe XVII 12.263999938964844 0.10193770260448132
Si XII 44.160301208496094 0.10068175237526199
Ni XIX 12.4350004196167 0.09868881845621418
Fe XIX 14.668999671936035 0.09764415209687931
Fe XVI 63.71099853515625 0.08924358504936311
Fe XVII 13.824999809265137 0.0879335250168359
Si XII 40.9109992980957 0.08714054397041306
Fe XX 12.951000213623047 0.08480209256050454
Fe XX 14.263400077819824 0.08114561486203685
C VI 33.73419952392578 0.07721653009600904
Fe XVIII 15.86989974975586 0.07606396144080754
Ni XIX 14.039799690246582 0.07532523461561419
Fe XX 12.905500411987305 0.0717605770153224
Fe XVIII 14.550999641418457 0.07013357311839734
Mg XI 9.314299583435059 0.06983122883637835
Fe XX 13.052000045776367 0.06688102035043217
Fe XVI 66.35700225830078 0.06674867602868391
Ne X 12.137499809265137 0.06484078899052988
Fe XX 13.404999732971191 0.06426023061392712
Fe XIX 13.741299629211426 0.06367427982521132
Fe XIX 13.64799976348877 0.06124638966977607
Fe XIX 14.991999626159668 0.061011027254637806
O VIII 16.00550079345703 0.060592315116662104
Fe XIX 16.2721004486084 0.058833518495426326
Fe XVIII 15.869999885559082 0.05656686928372383
Fe XVIII 14.352499961853027 0.056289421279127956
Si XII 44.01860046386719 0.05546283810863876
Ni XIX 13.777700424194336 0.05501988721911229
Ni XIX 14.07409954071045 0.05398649743721928
Fe XIX 10.815999984741211 0.051822278916743965
Fe XIX 13.721199989318848 0.051577213894886405
Fe XIX 14.738100051879883 0.05061878214002545
Si XII 45.69060134887695 0.050449048875881336
Fe XVIII 11.420000076293945 0.049307385293843826
Fe XVII 15.45300006866455 0.047673209009889396
Fe XX 14.763699531555176 0.047331328992566146
Fe XVIII 11.524999618530273 0.046917444607971884
Fe XVI 66.2490005493164 0.04650029613422574
Fe XVI 54.709999084472656 0.04645653558773096
Fe XIX 13.429699897766113 0.044905277232104954
Fe XVII 11.25 0.04458767268142438
Si XII 40.95100021362305 0.04359226165423845
Fe XX 13.780599594116211 0.04334895194900121
Fe XX 13.14050006866455 0.043260958852414315
Fe XVI 62.87099838256836 0.042536662032874
C VI 33.739601135253906 0.04213364999763524
Fe XVIII 14.418899536132812 0.04198042947994366
Fe XX 12.581000328063965 0.040833399615154434
Fe XX 13.1427001953125 0.040412458766163895
Fe XVIII 14.258000373840332 0.039668601438808446
Fe XX 14.913299560546875 0.037709025845725774
Ni XX 11.831999778747559 0.03767134146412616
Ni XIX 12.656000137329102 0.03761767141978405
Fe XVI 50.361000061035156 0.035792486974997886
Fe XVIII 16.02560043334961 0.03570691710580132
Ne IX 13.44729995727539 0.0352391004263809
Si XIII 6.688199996948242 0.03491540183370167
N VII 24.779199600219727 0.03463389674328
Fe XVIII 16.305700302124023 0.03435892010446177
Fe XVIII 14.344099998474121 0.03400298971487221
Fe XIX 13.768600463867188 0.03384284740727596
Fe XX 14.3125 0.033667228940748986
Fe XX 12.426300048828125 0.03291538683062994
Fe XX 14.402899742126465 0.03275909477437533
Mg X 63.29499816894531 0.032448608979999266
S XIV 32.560001373291016 0.0323343197150074
Fe XVIII 14.418999671936035 0.03225707280734121
Fe XVIII 14.579999923706055 0.03211511363279899
Fe XVII 11.128999710083008 0.032050939369853385
Fe XIX 13.017999649047852 0.032003465524315204
O VII 21.601499557495117 0.031891621394518864
Fe XIX 14.991999626159668 0.03140837553704566
Fe XX 12.581000328063965 0.03108486417906092
Fe XVII 46.30670166015625 0.031045749409376988
Fe XVIII 15.766400337219238 0.030708480459568525
Fe XVIII 11.32610034942627 0.030619932238577594
Fe XIX 13.843700408935547 0.03035222441290314
O VIII 16.00670051574707 0.030269172791130783
Fe XVIII 13.961999893188477 0.03002176793136807
Fe XVI 15.313899993896484 0.029879005057668415
Fe XX 12.426300048828125 0.02882088840195315
Fe XX 13.516599655151367 0.028675548737424113
Mg X 57.875999450683594 0.027449896626626776
S XIV 30.427000045776367 0.027413783952741112
Fe XVIII 13.396900177001953 0.027293621695097636
Fe XIX 10.632599830627441 0.02657790716739449
Fe XVI 15.211400032043457 0.026573232245620203
Fe XVIII 11.32610034942627 0.0265622517090096
Fe XX 13.843999862670898 0.026460020798328942
Fe XIX 13.77400016784668 0.026444041260405282
Fe XVIII 15.396699905395508 0.026014637485102746
Fe XVI 54.125999450683594 0.025707734884033894
Si XII 45.520599365234375 0.025486112098778574
Fe XX 14.962200164794922 0.025304056653718054
Fe XVIII 11.524999618530273 0.025290902675160078
Fe XVIII 10.522899627685547 0.024498275951999146
Fe XIX 16.34000015258789 0.024190689640817362
Fe XVII 10.770000457763672 0.023871780993225625
Fe XIX 9.847999572753906 0.0236616996208335
Fe XVIII 14.77180004119873 0.02352350242871807
Fe XX 12.965999603271484 0.023213139186119364
Fe XIX 13.671600341796875 0.023140874774942928
Fe XX 12.982000350952148 0.023081344715554872
Fe XIX 14.041999816894531 0.022942503494296114
Ni XX 13.309000015258789 0.02277110016065875
Fe XX 14.45300006866455 0.022760943509238747
Fe XIX 13.554800033569336 0.022513426342591823
Fe XX 14.44849967956543 0.022229385739754685
Fe XVIII 14.487099647521973 0.02201529349355861
Fe XIX 14.01729965209961 0.021774421861430516
Mg XI 7.850500106811523 0.02168266963225494
Fe XIX 13.673299789428711 0.021602486104542344
Ni XX 11.841400146484375 0.021599196758897325
Fe XX 15.515000343322754 0.02131992669450924
Fe XX 14.932000160217285 0.021020649251908823
Fe XIX 13.937999725341797 0.020968003348795325
Fe XVIII 15.449700355529785 0.020490699473778674
Fe XX 13.972299575805664 0.020469913771521523
Fe XIX 12.923999786376953 0.02039908969646279
Mg XI 9.231200218200684 0.02038227962299861
Fe XIX 16.378700256347656 0.020333225883396835
Fe XVIII 10.433799743652344 0.02031735248149734
O VIII 15.175999641418457 0.02011130251389851
Fe XVIII 13.318900108337402 0.020046235535070574
Fe XVII 46.29999923706055 0.0199500220518005
Fe XIX 14.20259952545166 0.01984383913036219
Fe XIX 13.635700225830078 0.019724684414911462
Fe XX 12.986000061035156 0.01967054231829778
Fe XVIII 18.10700035095215 0.019649977023697075
Fe XX 13.267000198364258 0.019594448853658523
Fe XVII 16.335599899291992 0.019492720971561784
Ni XX 14.470000267028809 0.01940433867493388
Mg X 63.152000427246094 0.019274089570629514
Fe XX 15.0600004196167 0.01907239522433792
Fe XVI 50.564998626708984 0.018985721181499812
Fe XX 13.267000198364258 0.01887679238257305
Ca XVII 21.1560001373291 0.018860163139847688
Fe XX 13.393600463867188 0.018739722639807863
Fe XVIII 14.469599723815918 0.018621514420576586
Fe XIX 14.932000160217285 0.018333352016878086
Fe XIX 14.017900466918945 0.018027781292741074
S XIV 32.41600036621094 0.017926504389995566
Fe XIX 13.253999710083008 0.01780574644376687
Ne X 10.238499641418457 0.017790483436973637
Fe XIX 10.684000015258789 0.017362191374909505
Fe XVIII 11.310600280761719 0.017337463765566087
Fe XVIII 14.889599800109863 0.017311097052677027
N VII 24.78459930419922 0.017303403764060782
Fe XVIII 14.70580005645752 0.017259867159082124
Fe XIX 12.923999786376953 0.01714123814931826
Fe XVIII 14.14430046081543 0.017135969358642256
Fe XIX 10.654999732971191 0.017022177328555523
Al XII 7.757299900054932 0.016955139064576748
Fe XVIII 14.45300006866455 0.0168909515619719
Fe XVIII 13.424099922180176 0.016718225581462066
Fe XX 14.830400466918945 0.016606851217206512
Fe XVIII 43.25550079345703 0.016588307257011183
Ne IX 13.699000358581543 0.016576163534173863
Fe XX 13.505000114440918 0.016552768941263514
Mg X 65.84500122070312 0.01653361640890651
Fe XIX 40.39039993286133 0.016507420508522684
Fe XVIII 43.294700622558594 0.016458807391251582
Ni XX 13.255999565124512 0.016379154293431664
Fe XIX 16.747699737548828 0.01632216998208309
Fe XVIII 18.52589988708496 0.016177357034703313
Fe XVIII 13.35509967803955 0.016094563131033832
Fe XIX 13.461999893188477 0.01605772539442215
S XIV 33.54949951171875 0.016002416230505932
Fe XIX 13.72089958190918 0.015942939108536554
Fe XVII 13.890000343322754 0.015790673092412895
Fe XX 14.805999755859375 0.015590428894662528
Fe XVII 50.60599899291992 0.015582969832999055
Fe XIX 15.195799827575684 0.01548531005655829
Fe XVIII 13.46399974822998 0.014810692728535734
Fe XVII 10.503999710083008 0.014800927582606023
Na XI 10.023200035095215 0.014646703909370198
Fe XX 12.824199676513672 0.014572755931364126
Fe XX 12.978500366210938 0.014533136512329407
Ni XX 11.961000442504883 0.014496092502525343
Fe XVII 16.003799438476562 0.014454884437429179
Fe XX 13.292400360107422 0.014422815457249688
Fe XIX 13.67199993133545 0.01413883202646433
Fe XX 13.090999603271484 0.014091182308016977
Ni XX 12.112000465393066 0.01406360135179621
Fe XVIII 14.610400199890137 0.013956690627896175
Fe XIX 40.86800003051758 0.013910066103314193
Mg X 57.91999816894531 0.013898470247605491
Fe XVIII 14.968700408935547 0.013890713878241794
Fe XX 10.986599922180176 0.013866913734661177
Fe XX 12.565999984741211 0.013858117170158423
S XIV 30.4689998626709 0.013756881576649062
Fe XIX 13.45199966430664 0.013727131530682293
Si XII 31.012100219726562 0.013700636911835576
Fe XIX 14.9350004196167 0.013626028922017664
Ni XX 12.927000045776367 0.013534148832382817
Fe XX 13.253999710083008 0.01349938360321063
Fe XX 13.114700317382812 0.013425533331838636
Fe XVII 10.656999588012695 0.013422460324727601
Fe XX 13.332900047302246 0.013373400903587155
Fe XIX 16.77549934387207 0.013336566510857331
Fe XV 59.404998779296875 0.013263691211093111
Fe XX 13.525199890136719 0.01322561130515175
Fe XX 12.600000381469727 0.013136563064711903
Fe XIX 15.163100242614746 0.013132418422522935
Fe XVIII 14.670499801635742 0.0126369949866854
Fe XX 13.545000076293945 0.012629149636667721
Fe XIX 12.923999786376953 0.012586011808904824
Fe XIX 13.492300033569336 0.012579582363086223
Fe XX 13.12339973449707 0.012577571018220504
Fe XIX 13.553999900817871 0.012563150268335002
Fe XVII 57.709999084472656 0.012539741112450575
Fe XVII 16.237899780273438 0.012330859884277467
Si XII 32.973201751708984 0.01229396730169684
Fe XIX 13.401399612426758 0.012229486454368122
Fe XVIII 10.526100158691406 0.012159272571218514
Fe XIX 14.294599533081055 0.012089911954156606
Fe XX 11.006500244140625 0.012061830487045374
Fe XIX 13.76200008392334 0.011838587682209883
Fe XIX 10.757499694824219 0.011826168868421645
Fe XVII 50.2599983215332 0.011775698705007139
Fe XIX 9.690500259399414 0.01167499546238898
Fe XIX 14.280799865722656 0.011641402960270945
Fe XIX 16.32469940185547 0.011581178170832079
Fe XVIII 16.08930015563965 0.01156749773169528
Fe XX 14.812800407409668 0.011491591693494344
Fe XVIII 14.487099647521973 0.011383214210557274
Ni XIX 10.110199928283691 0.011174044643049382
Fe XVII 68.51000213623047 0.011167641612188216
Fe XX 14.123299598693848 0.011135510592403978
Ca XVII 22.15999984741211 0.011107746159628005
Fe XVIII 17.355899810791016 0.011088773722254157
Ni XVIII 51.03799819946289 0.011038947987024336
Fe XIX 14.410599708557129 0.011018262169210218
Fe XVII 13.15939998626709 0.011007895213708428
Fe XVIII 10.3503999710083 0.010972075068789947
Fe XX 15.062600135803223 0.010960352067074498
Si XII 44.1775016784668 0.010949129247308863
Fe XX 13.395299911499023 0.010938657682671542
Fe XIX 13.936200141906738 0.010815965581656805
Fe XIX 10.63230037689209 0.010812061571201877
Ar XVI 24.993999481201172 0.010749924685804872
Fe XIX 13.48859977722168 0.01069100374181761
C VI 28.465200424194336 0.010689497736800134
Fe XX 13.090999603271484 0.010632412206779684
Ni XIX 9.977100372314453 0.010574439065346662
Fe XVIII 17.803300857543945 0.010471216210668786
Fe XX 15.843700408935547 0.010362115875318184
Fe XIX 15.041600227355957 0.010229387362530212
O VIII 15.17650032043457 0.010048752125656746
Ca XV 23.023000717163086 0.008795925542824056
Fe XVII 46.400001525878906 0.008780780699897775
Fe XVI 40.153099060058594 0.008658515701709894
Ni XVIII 52.720001220703125 0.008438365574888263
Fe XV 52.9109992980957 0.008335206662811163
Mg X 65.6729965209961 0.008326110989277167
Ca XVI 21.450000762939453 0.008183895305343327
Fe XVI 15.163000106811523 0.00816026706457009
Fe XVI 46.71820068359375 0.007995786513907123
Fe XVI 42.3036994934082 0.007540705907685791
Fe XVI 15.109199523925781 0.006980012572200949
Si XII 31.0226993560791 0.006965868419351094
Si XII 32.88779830932617 0.006678871128268564
Ni XVIII 44.365299224853516 0.006320532011585248
Si XI 43.762001037597656 0.006307501965482692
Ni XVIII 52.6150016784668 0.005877772045205512
S XIII 35.66699981689453 0.005637143585472517
Fe XVI 46.661399841308594 0.005602199489966244
Ne IX 11.546600341796875 0.005372643071946554
C VI 28.466299057006836 0.00534384290119645
Fe XVI 54.74700164794922 0.0052392664368922865
Ni XVII 55.15999984741211 0.005238340743766359
Ni XVIII 50.255001068115234 0.005227165401955449
Fe XVI 36.749000549316406 0.005092138238787243
O VII 18.62700080871582 0.005067018980418721
Fe XV 63.957000732421875 0.005059506577822573
Al XI 52.44599914550781 0.005037632910020915
N VII 20.909500122070312 0.004904355032598671
Si XII 33.31669998168945 0.004889309043894498
Cr XV 18.49690055847168 0.004830273246736862
Fe XVI 39.82709884643555 0.004756818115364387
Ne IX 13.553099632263184 0.004626113384010363
Fe XVI 17.45319938659668 0.0043832757643035665
Al XI 48.297000885009766 0.004370497053388531
Fe XVI 15.158499717712402 0.004188417470210847
Mg X 44.04999923706055 0.004121595859765846
Mg X 47.31019973754883 0.004114217896106279
Fe XVI 15.537599563598633 0.004013097893797865
Mg X 63.31100082397461 0.0038203998970354927
C VI 26.989599227905273 0.0035869758505102174
Fe XVI 41.932098388671875 0.0035016236936337916
Cr XV 21.152799606323242 0.0034660087118058573
S XIII 37.724998474121094 0.003432289015229992
S XIII 32.242000579833984 0.0033990524843805223
Fe XV 46.2869987487793 0.0033930831940416923
Fe XVI 66.37699890136719 0.0033157598742392233
Ca XIV 24.132999420166016 0.0032081574532214084
Cr XV 20.862899780273438 0.002961884024453743
Fe XVI 35.105899810791016 0.002936676325049824
Al XI 52.29899978637695 0.002793611099014601
Cr XV 21.212099075317383 0.0027876907943488586
Mg X 44.04999923706055 0.0027268281470441565
Fe XVI 36.803001403808594 0.0025985253004621608
Al XI 54.38800048828125 0.002549125965681085
Ca XIV 24.085899353027344 0.0025432723435804376
Fe XV 43.78099822998047 0.0024154521255133287
Fe XV 53.165000915527344 0.00238864573127441
Mg X 47.22850036621094 0.0023675430118262064
Ni XVII 57.57899856567383 0.0022769191630498003
Ni XVII 47.6619987487793 0.0022235562470589504
Fe XVI 40.24489974975586 0.0020876760983862264
Fe XVI 35.999000549316406 0.001874071528721757
O VII 17.76799964904785 0.0018407626585840568
Mg X 47.887699127197266 0.0016183232071091502
Fe XV 49.4900016784668 0.001556803100273336
Ca XIII 26.033000946044922 0.0015122358779478963
Fe XV 64.87799835205078 0.0011257402135655225
Si XI 37.340999603271484 0.0010144114221865878
Ca XIII 25.976999282836914 0.0010045646969176264
Ca XII 27.972999572753906 0.0010006147032915848
O VII 17.395999908447266 0.0009625862770155171
Ar XIII 29.319000244140625 0.000957939239680539
Si XI 33.51300048828125 0.0009394490121149916
Ca XII 32.65700149536133 0.0009312172353319899
Ca XII 31.660999298095703 0.0007911445825152357
Ca XIII 32.803001403808594 0.0007390908477881756
Ca XII 27.972999572753906 0.0007055862211240275
Ca XIII 26.7189998626709 0.0006683623190949546
Ne VIII 67.38200378417969 0.0006360175420088604
Ca XII 32.280399322509766 0.0006308778896189857
Fe XV 44.981998443603516 0.0005291633614871386
Ca XII 32.49879837036133 0.0005051186655920418
Ca XII 37.604000091552734 0.0004974576019647164
Ca XII 28.47800064086914 0.0004880981904797062
Ar XII 31.3882999420166 0.00042363582396702244
Ar XII 31.347000122070312 0.00030667805071884265
[16]:
elements = np.unique(hot_line_table['element'])
Equilibrium Ionization#
[17]:
temperature_grid = 10**np.arange(4,8,0.001) * u.K
[18]:
for el_name in elements:
el = fiasco.Element(el_name, temperature_grid)
fig = plt.figure(figsize=(15,5))
ax = fig.add_subplot()
for i in el:
ioneq = el.equilibrium_ionization[:, i.charge_state]
l, = ax.plot(el.temperature, ioneq)
ax.text(el.temperature[ioneq.argmax()].value, ioneq.max().value, i.ionization_stage_roman, color=l.get_color())
ax.set_title(el.element_name)
ax.set_xscale('log')
[454]:
fe = fiasco.Element('iron', temperature_grid)
[455]:
istate = 17
cur_ion = fe[istate-1]
with quantity_support():
plt.plot(fe.temperature, cur_ion.ionization_rate, label=f'ionization, {cur_ion.ionization_stage_roman}')
plt.plot(fe.temperature, cur_ion.recombination_rate, label=f'recombination, {cur_ion.ionization_stage_roman}')
plt.plot(fe.temperature, fe[istate-2].ionization_rate, label=f'ionization, {fe[istate-2].ionization_stage_roman}')
plt.plot(fe.temperature, fe[istate].recombination_rate, label=f'recombination, {fe[istate].ionization_stage_roman}')
plt.xscale('log')
plt.yscale('log')
plt.ylim(1e-12,1e-9)
plt.title(cur_ion.ion_name_roman)
plt.legend()
[455]:
<matplotlib.legend.Legend at 0x2bf81c100>
[456]:
with quantity_support():
plt.plot(fe[15].temperature, fe[15].direct_ionization_rate, label='DI')
plt.plot(fe[15].temperature, fe[15].excitation_autoionization_rate, label='EA')
plt.plot(fe[15].temperature, fe[15].direct_ionization_rate+fe[15].excitation_autoionization_rate, label='total')
plt.xscale('log')
plt.yscale('log')
plt.ylim(1e-12,1e-9)
plt.legend()
[456]:
<matplotlib.legend.Legend at 0x29b6a9bb0>
Simulate NEI#
[30]:
for el_name in elements:
el = fiasco.Element(el_name, temperature_grid,)
nei = non_equilibrium_ionization(el, time, electron_temperature, density, check_solution=True)
ieq = equilibrium_ionization(el, electron_temperature)
ionization_stages = np.unique(hot_line_table[hot_line_table['element'] == el.atomic_symbol]['ionization stage'])
fig = plt.figure(figsize=(15,5))
ax = fig.add_subplot()
ax2 = ax.twinx()
ax3 = ax.twinx()
ax3.spines.right.set_position(("axes", 1.05))
for ion in el:
if ion.ionization_stage in ionization_stages:
l, = ax.plot(time, ieq[:, ion.charge_state], ls='--')
ax.plot(time, nei[:, ion.charge_state], ls='-', color=l.get_color(),
label=f'{ion.ion_name_roman} ({ion.formation_temperature.to("MK"):.2f})')
ax2.plot(time, electron_temperature.to('MK'), color='k')
ax3.plot(time, density/1e9, color='k', ls='-.')
ax.legend(loc='upper right', bbox_to_anchor=(1,1.3), ncol=3, frameon=False)
ax.set_xlim(90,750)
ax2.set_xlim(ax.get_xlim())
ax.set_title(el.element_name)
Take the difference between the NEI and IEQ results
[32]:
for el_name in elements:
el = fiasco.Element(el_name, temperature_grid,)
nei = non_equilibrium_ionization(el, time, electron_temperature, density, check_solution=True)
ieq = equilibrium_ionization(el, electron_temperature)
ionization_stages = np.unique(hot_line_table[hot_line_table['element'] == el.atomic_symbol]['ionization stage'])
fig = plt.figure(figsize=(15,5))
ax = fig.add_subplot()
ax2 = ax.twinx()
ax3 = ax.twinx()
ax3.spines.right.set_position(("axes", 1.05))
for ion in el:
if ion.ionization_stage in ionization_stages:
ax.plot(time, ieq[:, ion.charge_state] - nei[:, ion.charge_state], ls='-',
label=f'{ion.ion_name_roman} ({ion.formation_temperature.to("MK"):.2f})')
ax2.plot(time, electron_temperature.to('MK'), color='k')
ax3.plot(time, density/1e9, color='k', ls='-.')
ax.legend(loc='upper right', bbox_to_anchor=(1,1.3), ncol=3, frameon=False)
ax.set_xlim(90,750)
ax.set_ylim(-1, 1)
ax2.set_xlim(ax.get_xlim())
ax.set_title(el.element_name)
ax.axhline(y=0, ls=':', color='k')
Effective temperature#
[474]:
fig = plt.figure(figsize=(10,10))
ax1 = fig.add_subplot(221)
ax1.plot(time, heat, color='k')
ax2 = fig.add_subplot(222)
ax2.plot(time, electron_temperature, color='k')
ax3 = fig.add_subplot(223)
ax3.plot(time, density, color='k')
ax4 = fig.add_subplot(224)
ax4.plot(electron_temperature, density, color='k')
ax4.set_yscale('log')
ax4.set_xscale('log')
for el_name in elements:
el = fiasco.Element(el_name, temperature_grid)
T_eff = effective_temperature(el, time, electron_temperature, density)
ax2.plot(time, T_eff, label=el.atomic_symbol)
ax4.plot(T_eff, density)
ax2.legend()
[474]:
<matplotlib.legend.Legend at 0x2ba763220>
HYDRAD NEI Simulation#
In the EBTEL simulations above, the density can rise unrealistically quickly (due to the lack of spatial dimension) such that the time taken to reach ionization equilibrium (and thus the impact of non-equilibrium ionization) may be underestimated. As such, we’ll run a more sophisticated 1D simulation and compute the ion fractions for the same ions as above out of equilibrium to assess more accurately whether this is something we need to worry about in our observations.
[657]:
elements
[657]:
<Column name='element' dtype='str2' description='' length=14>
| Al |
| Ar |
| C |
| Ca |
| Cr |
| Fe |
| Mg |
| N |
| Na |
| Ne |
| Ni |
| O |
| S |
| Si |
[695]:
hydrad_config = pydrad.configure.data.get_defaults()
WARNING: UnitsWarning: The unit 'erg' has been deprecated in the VOUnit standard. Suggested: cm**2.g.s**-2. [astropy.units.format.utils]
[696]:
hydrad_config['general']['loop_length'] = 2*(loop_length + hydrad_config['general']['footpoint_height'])
hydrad_config['general']['total_time'] = total_time
hydrad_config['general']['heat_flux_timestep_limit'] = 1e-5 * u.s
hydrad_config['general']['write_file_equation_terms'] = False
hydrad_config['general']['write_file_hydrogen_level_populations'] = False
hydrad_config['general']['write_file_ion_populations'] = True
hydrad_config['general']['write_file_physical'] = True
hydrad_config['general']['write_file_timescales'] = False
hydrad_config['initial_conditions']['footpoint_density'] = 2.5e9 * u.cm**(-3)
hydrad_config['initial_conditions']['footpoint_temperature'] = 1e4 * u.K
hydrad_config['initial_conditions']['heating_scale_height'] = 1e300 * u.cm
hydrad_config['radiation']['decouple_ionization_state_solver'] = True
hydrad_config['radiation']['elements_nonequilibrium'] = ['Ca', 'Fe', 'Mg', 'Ne', 'O', 'Si']
hydrad_config['radiation']['rates_dataset'] = 'chianti_v10'
hydrad_config['radiation']['use_power_law_radiative_losses'] = True
hydrad_config['grid']['adapt_every_n_time_steps'] = 10
hydrad_config['grid']['initial_refinement_level'] = 6
hydrad_config['grid']['maximum_refinement_level'] = 6
hydrad_config['grid']['refine_on_hydrogen_energy'] = False
hydrad_config['solver']['cutoff_temperature_fraction'] = 0.2
hydrad_config['solver']['cutoff_ion_fraction'] = 1e-6
hydrad_config['heating']['background']['use_initial_conditions'] = True
hydrad_config['heating']['electron_heating'] = 1.0
events = []
for event in ebtel_config['heating']['events']:
event = event['event']
events.append({
'location': hydrad_config['general']['loop_length'] / 2,
'scale_height': 1e300 * u.cm,
'time_start': event['rise_start'] * u.s,
'total_duration': (event['decay_end'] - event['rise_start']) * u.s,
'rise_duration': (event['rise_end'] - event['rise_start']) * u.s,
'decay_duration': (event['decay_end'] - event['decay_start']) * u.s,
'rate': event['magnitude'] * u.Unit('erg cm-3 s-1'),
})
hydrad_config['heating']['events'] = events
[697]:
hydrad_config
[697]:
{'asdf_library': {'author': 'Space Telescope Science Institute',
'homepage': 'http://github.com/spacetelescope/asdf',
'name': 'asdf',
'version': '2.7.1'},
'history': {'extensions': [{'extension_class': 'astropy.io.misc.asdf.extension.AstropyAsdfExtension',
'software': {'name': 'astropy', 'version': '4.0.1.post1'}},
{'extension_class': 'asdf.extension.BuiltinExtension',
'software': {'name': 'asdf', 'version': '2.7.1'}}]},
'general': {'footpoint_height': <Quantity 5.e+08 cm>,
'force_single_fluid': False,
'heat_flux_limiting_coefficient': 0.167,
'heat_flux_timestep_limit': <Quantity 1.e-05 s>,
'logging_frequency': 1000,
'loop_inclination': <Quantity 0. deg>,
'loop_length': <Quantity 90. Mm>,
'minimum_collisional_coupling_timescale': <Quantity 0.01 s>,
'output_interval': <Quantity 1. s>,
'total_time': <Quantity 13000. s>,
'use_kinetic_model': False,
'write_file_equation_terms': False,
'write_file_hydrogen_level_populations': False,
'write_file_ion_populations': True,
'write_file_physical': True,
'write_file_timescales': False},
'grid': {'adapt': True,
'adapt_every_n_time_steps': 10,
'enforce_conservation': True,
'initial_refinement_level': 6,
'linear_restriction': True,
'maximum_cell_width': <Quantity 0.5 Mm>,
'maximum_fractional_difference': 0.2,
'maximum_refinement_level': 6,
'maximum_variation': 0.1,
'minimum_delta_s': <Quantity 1. cm>,
'minimum_fractional_difference': 0.1,
'refine_on_density': True,
'refine_on_electron_energy': True,
'refine_on_hydrogen_energy': False},
'heating': {'alfven_wave': False,
'background': {'use_initial_conditions': True},
'beam': False,
'electron_heating': 1.0,
'events': [{'location': <Quantity 45. Mm>,
'scale_height': <Quantity 1.e+300 cm>,
'time_start': <Quantity 100. s>,
'total_duration': <Quantity 200. s>,
'rise_duration': <Quantity 100. s>,
'decay_duration': <Quantity 100. s>,
'rate': <Quantity 0.01 erg / (cm3 s)>}]},
'initial_conditions': {'footpoint_density': <Quantity 2.5e+09 1 / cm3>,
'footpoint_temperature': <Quantity 10000. K>,
'heating_location': <Quantity 45. Mm>,
'heating_range_fine_tuning': 10000.0,
'heating_range_lower_bound': <Quantity 1.e-08 erg / (cm3 s)>,
'heating_range_step_size': 0.001,
'heating_range_upper_bound': <Quantity 100. erg / (cm3 s)>,
'heating_scale_height': <Quantity 1.e+300 cm>,
'isothermal': False,
'use_poly_fit_gravity': False,
'use_poly_fit_magnetic_field': False},
'radiation': {'abundance_dataset': 'asplund',
'decouple_ionization_state_solver': True,
'density_dependent_rates': False,
'elements_equilibrium': [],
'elements_nonequilibrium': ['Ca', 'Fe', 'Mg', 'Ne', 'O', 'Si'],
'emissivity_dataset': 'chianti_v7',
'nlte_chromosphere': False,
'optically_thick_radiation': False,
'ranges_dataset': 'ranges',
'rates_dataset': 'chianti_v10',
'use_power_law_radiative_losses': True},
'solver': {'cutoff_ion_fraction': 1e-06,
'epsilon': 0.01,
'epsilon_d': 0.1,
'epsilon_r': 1.8649415311920072,
'maximum_optically_thin_density': <Quantity 1.e+12 1 / cm3>,
'minimum_radiation_temperature': <Quantity 20000. K>,
'minimum_temperature': <Quantity 10000. K>,
'safety_advection': 1.0,
'safety_atomic': 1.0,
'safety_conduction': 1.0,
'safety_radiation': 0.1,
'safety_viscosity': 1.0,
'timestep_increase_limit': 0.05,
'zero_over_temperature_interval': <Quantity 500. K>,
'cutoff_temperature_fraction': 0.2}}
[698]:
conf = pydrad.configure.Configure(hydrad_config)
[700]:
conf.setup_simulation('hydrad_simulations/hydrad-moxsi-nei-higher-cutoff',
'/Users/wtbarnes/Documents/codes/HYDRAD/',
overwrite=True)
WARNING: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError. [astropy.units.quantity]
INFO:
Calculating initial hydrostatic conditions...
Peak heating range = 1.3062e-06 -> 1.3092e-06 erg cm^-3 s^-1
Optimum peak heating rate = 1.3063e-06 erg cm^-3 s^-1
Writing initial conditions file...
Done!
[pydrad.configure.util]
WARNING: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError. [astropy.units.quantity]
WARNING: In file included from ../source/main.cpp:14:
In file included from ../source/mesh.h:13:
../source/eqns.h:38:31: warning: default member initializer for non-static data member is a C++11 extension [-Wc++11-extensions]
double Tc, Te_max, old_Tc = 0.0;
^
1 warning generated.
In file included from ../source/mesh.cpp:17:
In file included from ../source/mesh.h:13:
../source/eqns.h:38:31: warning: default member initializer for non-static data member is a C++11 extension [-Wc++11-extensions]
double Tc, Te_max, old_Tc = 0.0;
^
1 warning generated.
In file included from ../source/eqns.cpp:19:
../source/eqns.h:38:31: warning: default member initializer for non-static data member is a C++11 extension [-Wc++11-extensions]
double Tc, Te_max, old_Tc = 0.0;
^
1 warning generated.
[pydrad.configure.util]
WARNING: UnitsWarning: The unit 'erg' has been deprecated in the VOUnit standard. Suggested: cm**2.g.s**-2. [astropy.units.format.utils]
[723]:
strand = pydrad.parse.Strand('hydrad_simulations/hydrad-moxsi-nei-higher-cutoff/')
[624]:
strand.initial_conditions.peek(limits={'temperature': [0,2], 'density': [1e7, 1e10]})
Compare the initial conditions of HYDRAD run…
[590]:
strand.initial_conditions.spatial_average(
'electron_temperature',
bounds=strand.config['general']['loop_length']*(0.25, 0.75),).to('MK')
[590]:
$0.38303183 \; \mathrm{MK}$
[591]:
strand.initial_conditions.spatial_average(
'electron_density',
bounds=strand.config['general']['loop_length']/2*(0.25, 1.25),) / 1e9
[591]:
$0.019967621 \; \mathrm{\frac{1}{cm^{3}}}$
To that of EBTEL
[592]:
density[0] / 1e9
[592]:
$0.01798226 \; \mathrm{\frac{1}{cm^{3}}}$
[593]:
electron_temperature[0].to('MK')
[593]:
$0.4106181 \; \mathrm{MK}$
[724]:
pydrad.visualize.plot_time_distance(
strand[90:-1],
['electron_temperature', 'electron_density','population_fraction_iron_18', 'population_fraction_iron_17'],
0.5*u.Mm,
time_unit='s',
space_unit='Mm',
units={'electron_temperature': 'MK'},
norm={'population_fraction_iron_18': matplotlib.colors.LogNorm(vmin=1e-6, vmax=1),
'population_fraction_iron_17': matplotlib.colors.LogNorm(vmin=1e-6, vmax=1),
'electron_density': matplotlib.colors.LogNorm(vmin=1e7, vmax=1e11)},
)
[724]:
(<Figure size 1000x1000 with 8 Axes>,
array([<AxesSubplot:ylabel='$s$ [Mm]'>, <AxesSubplot:ylabel='$s$ [Mm]'>,
<AxesSubplot:ylabel='$s$ [Mm]'>,
<AxesSubplot:xlabel='$t$ [s]', ylabel='$s$ [Mm]'>], dtype=object))
[740]:
i_time = -1
plt.plot(strand[i_time].grid_centers, strand[i_time].population_fraction_iron_18, label='NEI')
plt.plot(strand[i_time].grid_centers, equilibrium_ionization(fiasco.Element('Fe', temperature_grid),strand[i_time].electron_temperature)[:,17], label="EQ")
plt.legend()
[740]:
<matplotlib.legend.Legend at 0x3c0465eb0>
[ ]:
pydrad.visualize.plot_time_mesh()
[668]:
strand[:500].peek_time_distance()
[691]:
strand[100:500:2].peek_time_distance(cmap={'velocity':'RdBu'},
norm={'electron_density': matplotlib.colors.LogNorm(vmin=1e8, vmax=1e11),
'velocity': matplotlib.colors.Normalize(vmin=-100, vmax=100),
'electron_temperature': matplotlib.colors.Normalize(vmin=1, vmax=15)},
units={'velocity': 'km / s',
'electron_temperature': "MK"})
[717]:
from sunpy.coordinates.utils import solar_angle_equivalency
from sunpy.coordinates import get_earth
(150 * u.Mm).to('arcsec', equivalencies=solar_angle_equivalency(get_earth()))
INFO: Apparent body location accounts for 500.01 seconds of light travel time [sunpy.coordinates.ephemeris]
[717]:
$206.40326 \; \mathrm{{}^{\prime\prime}}$
[722]:
foo = 'foo \n'
print(foo)
print(foo.strip())
foo
foo
[755]:
strand = pydrad.parse.Strand('hydrad_simulations/hydrad-moxsi-no-nei/')[90:5000]
[757]:
plt.figure(figsize=(10,5))
plt.subplot(121)
plt.plot(strand.time, strand.spatial_average('electron_temperature', bounds=strand.loop_length/2*(0.25, 1.25)))
plt.plot(time, electron_temperature)
plt.xlim(0,5000)
plt.subplot(122)
plt.plot(strand.time, strand.spatial_average('electron_density', bounds=strand.loop_length/2*(0.25, 1.25)))
plt.plot(time, density)
plt.xlim(0,5000)
[757]:
(0.0, 5000.0)
Sandbox#
[3]:
from collections import OrderedDict
[4]:
total_time = 1.3e4 * u.s
loop_length = 40 * u.Mm
heating_background = 2.75e-5 * u.Unit('erg cm-3 s-1')
# setup config
ebtel_config = {
'total_time': total_time.to_value('s'),
'tau': 1.0,
'tau_max': 1e+300,
'loop_length': loop_length.to_value('cm'),
'saturation_limit': 1.0,
'force_single_fluid': False,
'use_c1_loss_correction': True,
'use_c1_grav_correction': True,
'use_flux_limiting': True,
'calculate_dem': False,
'dem': OrderedDict({
'use_new_method': True,
'temperature': {
'bins': 451,
'log_min': 4,
'log_max': 8.5
}
}),
'save_terms': False,
'use_adaptive_solver': True,
'adaptive_solver_error': 1e-10,
'adaptive_solver_safety': 0.5,
'c1_cond0': 6.0,
'c1_rad0': 0.6,
'helium_to_hydrogen_ratio': 0.075,
'surface_gravity': 1.0,
'heating': OrderedDict({
'background': heating_background.to_value('erg cm-3 s-1'),
'partition': 1.0,
}),
}
# setup heating
duration = 200 * u.s
t_0 = 100 * u.s
wait_time = 1500 * u.s
number_events = 1
heating_rate = 0.05 * u.Unit('erg cm-3 s-1')
start_times = t_0 + (duration + wait_time) * np.arange(number_events)
events = []
for st in start_times:
events.append(
{'event': {'magnitude': heating_rate.to_value('erg cm-3 s-1'),
'rise_start': st.to_value('s'),
'rise_end': (st + duration/2).to_value('s'),
'decay_start': (st + duration/2).to_value('s'),
'decay_end': (st + duration).to_value('s')}}
)
ebtel_config['heating']['events'] = events
[5]:
ebtel_config
[5]:
{'total_time': 13000.0,
'tau': 1.0,
'tau_max': 1e+300,
'loop_length': 4000000000.0,
'saturation_limit': 1.0,
'force_single_fluid': False,
'use_c1_loss_correction': True,
'use_c1_grav_correction': True,
'use_flux_limiting': True,
'calculate_dem': False,
'dem': OrderedDict([('use_new_method', True),
('temperature', {'bins': 451, 'log_min': 4, 'log_max': 8.5})]),
'save_terms': False,
'use_adaptive_solver': True,
'adaptive_solver_error': 1e-10,
'adaptive_solver_safety': 0.5,
'c1_cond0': 6.0,
'c1_rad0': 0.6,
'helium_to_hydrogen_ratio': 0.075,
'surface_gravity': 1.0,
'heating': OrderedDict([('background', 2.75e-05),
('partition', 1.0),
('events',
[{'event': {'magnitude': 0.05,
'rise_start': 100.0,
'rise_end': 200.0,
'decay_start': 200.0,
'decay_end': 300.0}}])])}
[6]:
run_ebtel(ebtel_config, '/Users/wtbarnes/Documents/codes/ebtelPlusPlus/')
[6]:
{'time': array([0.00000000e+00, 2.00000000e+00, 6.50000000e+00, ...,
1.25817733e+04, 1.27264025e+04, 1.28897902e+04]),
'electron_temperature': array([1007727., 1007750., 1007848., ..., 1024449., 1024248., 1024022.]),
'ion_temperature': array([1007727., 1007702., 1007596., ..., 1020788., 1020599., 1020385.]),
'density': array([1.771306e+08, 1.771304e+08, 1.771296e+08, ..., 1.598711e+08,
1.601044e+08, 1.603605e+08]),
'electron_pressure': array([0.02463291, 0.02463344, 0.02463572, ..., 0.0226016 , 0.02263016,
0.02266135]),
'ion_pressure': array([0.02302641, 0.02302582, 0.02302329, ..., 0.02105208, 0.02107891,
0.0211082 ]),
'velocity': array([-2068.353, -2037.979, -1907.785, ..., 20524.63 , 19950.89 ,
19326.38 ]),
'heat': array([2.75e-05, 2.75e-05, 2.75e-05, ..., 2.75e-05, 2.75e-05, 2.75e-05])}
[ ]: