Note
Go to the end to download the full example code.
CHIANTI IDL Comparison: Continuum#
Compare the free-free and free-bound calculations to that in the CHIANTI IDL routines.
import hissw
import matplotlib.colors
import matplotlib.pyplot as plt
import numpy as np
from astropy.visualization import quantity_support
import fiasco
from fiasco.tests.idl.helpers import read_idl_test_output
quantity_support()
<astropy.visualization.units.quantity_support.<locals>.MplQuantityConverter object at 0x7f51a30dbf90>
Define a function for comparing and plotting the outputs from fiasco and IDL. Note that we have precomputed the IDL outputs.
def plot_idl_comparison(wavelength, temperature, result_fiasco, result_idl):
wave_mesh, temp_mesh = np.meshgrid(wavelength, temperature)
temp_mesh = temp_mesh.to_value('K')
wave_mesh = wave_mesh.to_value('Angstrom')
result_fiasco = result_fiasco.to_value('erg cm3 s-1 Angstrom-1')
result_idl = result_idl.to_value('erg cm3 s-1 Angstrom-1')
difference = result_fiasco - result_idl
difference_norm = difference / result_idl
fig = plt.figure(figsize=(10,10))
axes = fig.subplot_mosaic(
"""
AABB
.CC.
""",
sharex=True,
sharey=True,
)
norm = matplotlib.colors.LogNorm(vmin=1e-30,vmax=1e-26)
# IDL result
im = axes['A'].pcolormesh(temp_mesh, wave_mesh, result_idl, norm=norm)
axes['A'].set_title('IDL')
axes['A'].set_xscale('log')
axes['A'].set_xlabel('Temperature [K]')
axes['A'].set_ylabel(r'Wavelength [$\AA$]')
# fiasco result
im = axes['B'].pcolormesh(temp_mesh, wave_mesh, result_fiasco, norm=norm)
axes['B'].set_title('fiasco')
fig.colorbar(im, ax=[axes['A'], axes['B']], orientation='horizontal')
# Normalized difference
im = axes['C'].pcolormesh(temp_mesh, wave_mesh, difference_norm,
norm=matplotlib.colors.SymLogNorm(1e-3, vmin=-1, vmax=1),
cmap='RdBu')
axes['C'].set_title('Normalized Difference')
fig.colorbar(im, ax=axes['C'])
plt.show()
First, let’s compare the outputs for the free-free continuum emission, i.e. that emission produced by thermal bremsstrahlung.
idl_result_freefree = read_idl_test_output('freefree_all_ions', '8.0.7')
ion_kwargs = {'abundance': idl_result_freefree['abundance'], 'ioneq_filename': idl_result_freefree['ioneq']}
all_ions = [fiasco.Ion(ion_name, idl_result_freefree['temperature'], **ion_kwargs) for ion_name in fiasco.list_ions()]
all_ions = fiasco.IonCollection(*all_ions)
free_free = all_ions.free_free(idl_result_freefree['wavelength'])
plot_idl_comparison(idl_result_freefree['wavelength'],
idl_result_freefree['temperature'],
free_free,
idl_result_freefree['free_free'])
# This is just for printing the code used to produce the IDL result
env = hissw.Environment(ssw_home='', idl_home='')
template = env.env.from_string(idl_result_freefree['idl_script'])
print('IDL code to produce free-free result:')
print(template.render(**idl_result_freefree))
/home/docs/checkouts/readthedocs.org/user_builds/fiasco/conda/latest/lib/python3.11/site-packages/astropy/units/format/utils.py:220: UnitsWarning: The unit 'Angstrom' has been deprecated in the VOUnit standard. Suggested: 0.1nm.
warnings.warn(message, UnitsWarning)
/home/docs/checkouts/readthedocs.org/user_builds/fiasco/conda/latest/lib/python3.11/site-packages/astropy/units/format/utils.py:220: UnitsWarning: The unit 'erg' has been deprecated in the VOUnit standard. Suggested: cm**2.g.s**-2.
warnings.warn(message, UnitsWarning)
/home/docs/checkouts/readthedocs.org/user_builds/fiasco/conda/latest/lib/python3.11/site-packages/asdf/_asdf.py:359: AsdfWarning: File 'file:///home/docs/checkouts/readthedocs.org/user_builds/fiasco/conda/latest/lib/python3.11/site-packages/fiasco/tests/idl/data/freefree_all_ions_v8.0.7.asdf' was created with extension URI 'asdf://astropy.org/core/extensions/core-1.5.0' (from package asdf-astropy==0.5.0), which is not currently installed
warnings.warn(msg, AsdfWarning)
IDL code to produce free-free result:
; set common block
common elements, abund, abund_ref, ioneq, ioneq_logt, ioneq_ref
; read abundance and ionization equilibrium
abund_file = FILEPATH('sun_coronal_1992_feldman_ext.abund', ROOT_DIR=!xuvtop, SUBDIR='abundance')
ioneq_file = FILEPATH('chianti.ioneq', ROOT_DIR=!xuvtop, SUBDIR='ioneq')
read_abund, abund_file, abund, abund_ref
read_ioneq, ioneq_file, ioneq_logt, ioneq, ioneq_ref
; set temperature and wavelength
temperature = [100000d / 1d,3855226053388321d / 34359738368d,4325634777407421d / 34359738368d,1213360511705375d / 8589934592d,5445651543418815d / 34359738368d,3055060763706939d / 17179869184d,6855669111786069d / 34359738368d,1923046815002979d / 8589934592d,8630776059679839d / 34359738368d,1210486251738545d / 4294967296d,1358187913129459d / 4294967296d,3047823805882249d / 8589934592d,3419714555660727d / 8589934592d,3836982839896381d / 8589934592d,8610331110407387d / 17179869184d,4830475201763701d / 8589934592d,2709941159713435d / 4294967296d,6081207982541875d / 8589934592d,6823227580860407d / 8589934592d,1913946815906879d / 2147483648d,1000000d / 1d,4819032566735401d / 4294967296d,1351760867939819d / 1073741824d,6066802558526875d / 4294967296d,6807064429273519d / 4294967296d,1909412977316837d / 1073741824d,4284793194866293d / 2147483648d,600952129688431d / 268435456d,1348558759324975d / 536870912d,6052431258692725d / 2147483648d,106108430713239d / 33554432d,3809779757352811d / 1073741824d,4274643194575909d / 1073741824d,1199057137467619d / 268435456d,5381456944004617d / 1073741824d,3019047001102313d / 536870912d,1693713224820897d / 268435456d,7601509978177343d / 1073741824d,8529034476075509d / 1073741824d,4784867039767197d / 536870912d,10000000d / 1d,94121729819051d / 8388608d,844850542462387d / 67108864d,3791751599079297d / 268435456d,8508830536591899d / 536870912d,1193383110823023d / 67108864d,5355991493582877d / 268435456d,6009521296884309d / 268435456d,6742793796624875d / 268435456d,3782769536682953d / 134217728d,8488674457059119d / 268435456d,297639043543189d / 8388608d,2671651996609943d / 67108864d,5995285687338095d / 134217728d,6726821180005771d / 134217728d,3773808751377891d / 67108864d,4234283062052251d / 67108864d,593867967045105d / 8388608d,5330646547547193d / 67108864d,5981083799708997d / 67108864d,100000000d / 1d]
wavelength = [25d / 1d,26d / 1d,27d / 1d,28d / 1d,29d / 1d,30d / 1d,31d / 1d,32d / 1d,33d / 1d,34d / 1d,35d / 1d,36d / 1d,37d / 1d,38d / 1d,39d / 1d,40d / 1d,41d / 1d,42d / 1d,43d / 1d,44d / 1d,45d / 1d,46d / 1d,47d / 1d,48d / 1d,49d / 1d,50d / 1d,51d / 1d,52d / 1d,53d / 1d,54d / 1d,55d / 1d,56d / 1d,57d / 1d,58d / 1d,59d / 1d,60d / 1d,61d / 1d,62d / 1d,63d / 1d,64d / 1d,65d / 1d,66d / 1d,67d / 1d,68d / 1d,69d / 1d,70d / 1d,71d / 1d,72d / 1d,73d / 1d,74d / 1d,75d / 1d,76d / 1d,77d / 1d,78d / 1d,79d / 1d,80d / 1d,81d / 1d,82d / 1d,83d / 1d,84d / 1d,85d / 1d,86d / 1d,87d / 1d,88d / 1d,89d / 1d,90d / 1d,91d / 1d,92d / 1d,93d / 1d,94d / 1d,95d / 1d,96d / 1d,97d / 1d,98d / 1d,99d / 1d,100d / 1d,101d / 1d,102d / 1d,103d / 1d,104d / 1d,105d / 1d,106d / 1d,107d / 1d,108d / 1d,109d / 1d,110d / 1d,111d / 1d,112d / 1d,113d / 1d,114d / 1d,115d / 1d,116d / 1d,117d / 1d,118d / 1d,119d / 1d,120d / 1d,121d / 1d,122d / 1d,123d / 1d,124d / 1d,125d / 1d,126d / 1d,127d / 1d,128d / 1d,129d / 1d,130d / 1d,131d / 1d,132d / 1d,133d / 1d,134d / 1d,135d / 1d,136d / 1d,137d / 1d,138d / 1d,139d / 1d,140d / 1d,141d / 1d,142d / 1d,143d / 1d,144d / 1d,145d / 1d,146d / 1d,147d / 1d,148d / 1d,149d / 1d,150d / 1d,151d / 1d,152d / 1d,153d / 1d,154d / 1d,155d / 1d,156d / 1d,157d / 1d,158d / 1d,159d / 1d,160d / 1d,161d / 1d,162d / 1d,163d / 1d,164d / 1d,165d / 1d,166d / 1d,167d / 1d,168d / 1d,169d / 1d,170d / 1d,171d / 1d,172d / 1d,173d / 1d,174d / 1d,175d / 1d,176d / 1d,177d / 1d,178d / 1d,179d / 1d,180d / 1d,181d / 1d,182d / 1d,183d / 1d,184d / 1d,185d / 1d,186d / 1d,187d / 1d,188d / 1d,189d / 1d,190d / 1d,191d / 1d,192d / 1d,193d / 1d,194d / 1d,195d / 1d,196d / 1d,197d / 1d,198d / 1d,199d / 1d,200d / 1d,201d / 1d,202d / 1d,203d / 1d,204d / 1d,205d / 1d,206d / 1d,207d / 1d,208d / 1d,209d / 1d,210d / 1d,211d / 1d,212d / 1d,213d / 1d,214d / 1d,215d / 1d,216d / 1d,217d / 1d,218d / 1d,219d / 1d,220d / 1d,221d / 1d,222d / 1d,223d / 1d,224d / 1d,225d / 1d,226d / 1d,227d / 1d,228d / 1d,229d / 1d,230d / 1d,231d / 1d,232d / 1d,233d / 1d,234d / 1d,235d / 1d,236d / 1d,237d / 1d,238d / 1d,239d / 1d,240d / 1d,241d / 1d,242d / 1d,243d / 1d,244d / 1d,245d / 1d,246d / 1d,247d / 1d,248d / 1d,249d / 1d,250d / 1d,251d / 1d,252d / 1d,253d / 1d,254d / 1d,255d / 1d,256d / 1d,257d / 1d,258d / 1d,259d / 1d,260d / 1d,261d / 1d,262d / 1d,263d / 1d,264d / 1d,265d / 1d,266d / 1d,267d / 1d,268d / 1d,269d / 1d,270d / 1d,271d / 1d,272d / 1d,273d / 1d,274d / 1d,275d / 1d,276d / 1d,277d / 1d,278d / 1d,279d / 1d,280d / 1d,281d / 1d,282d / 1d,283d / 1d,284d / 1d,285d / 1d,286d / 1d,287d / 1d,288d / 1d,289d / 1d,290d / 1d,291d / 1d,292d / 1d,293d / 1d,294d / 1d,295d / 1d,296d / 1d,297d / 1d,298d / 1d,299d / 1d,300d / 1d,301d / 1d,302d / 1d,303d / 1d,304d / 1d,305d / 1d,306d / 1d,307d / 1d,308d / 1d,309d / 1d,310d / 1d,311d / 1d,312d / 1d,313d / 1d,314d / 1d,315d / 1d,316d / 1d,317d / 1d,318d / 1d,319d / 1d,320d / 1d,321d / 1d,322d / 1d,323d / 1d,324d / 1d,325d / 1d,326d / 1d,327d / 1d,328d / 1d,329d / 1d,330d / 1d,331d / 1d,332d / 1d,333d / 1d,334d / 1d,335d / 1d,336d / 1d,337d / 1d,338d / 1d,339d / 1d,340d / 1d,341d / 1d,342d / 1d,343d / 1d,344d / 1d,345d / 1d,346d / 1d,347d / 1d,348d / 1d,349d / 1d,350d / 1d,351d / 1d,352d / 1d,353d / 1d,354d / 1d,355d / 1d,356d / 1d,357d / 1d,358d / 1d,359d / 1d,360d / 1d,361d / 1d,362d / 1d,363d / 1d,364d / 1d,365d / 1d,366d / 1d,367d / 1d,368d / 1d,369d / 1d,370d / 1d,371d / 1d,372d / 1d,373d / 1d,374d / 1d,375d / 1d,376d / 1d,377d / 1d,378d / 1d,379d / 1d,380d / 1d,381d / 1d,382d / 1d,383d / 1d,384d / 1d,385d / 1d,386d / 1d,387d / 1d,388d / 1d,389d / 1d,390d / 1d,391d / 1d,392d / 1d,393d / 1d,394d / 1d,395d / 1d,396d / 1d,397d / 1d,398d / 1d,399d / 1d,400d / 1d,401d / 1d,402d / 1d,403d / 1d,404d / 1d,405d / 1d,406d / 1d,407d / 1d,408d / 1d,409d / 1d,410d / 1d,411d / 1d,412d / 1d,413d / 1d]
; calculate free-free
freefree, temperature, wavelength, free_free, /no_setup
Next, let’s compare the outputs for the free-bound continuum emission.
idl_result_freebound = read_idl_test_output('freebound_all_ions', '8.0.7')
ion_kwargs = {'abundance': idl_result_freebound['abundance'], 'ioneq_filename': idl_result_freebound['ioneq']}
all_ions = [fiasco.Ion(ion_name, idl_result_freebound['temperature'], **ion_kwargs) for ion_name in fiasco.list_ions()]
all_ions = fiasco.IonCollection(*all_ions)
free_bound = all_ions.free_bound(idl_result_freebound['wavelength'])
plot_idl_comparison(idl_result_freebound['wavelength'],
idl_result_freebound['temperature'],
free_bound,
idl_result_freebound['free_bound'])
# This is just for printing the code used to produce the IDL result
template = env.env.from_string(idl_result_freebound['idl_script'])
print('IDL code to produce free-free result:')
print(template.render(**idl_result_freebound))
/home/docs/checkouts/readthedocs.org/user_builds/fiasco/conda/latest/lib/python3.11/site-packages/asdf/_asdf.py:359: AsdfWarning: File 'file:///home/docs/checkouts/readthedocs.org/user_builds/fiasco/conda/latest/lib/python3.11/site-packages/fiasco/tests/idl/data/freebound_all_ions_v8.0.7.asdf' was created with extension URI 'asdf://astropy.org/core/extensions/core-1.5.0' (from package asdf-astropy==0.5.0), which is not currently installed
warnings.warn(msg, AsdfWarning)
WARNING: H 2 not included in free-bound emission. _fblvl dataset missing for H 2. [fiasco.collections]
WARNING: He 3 not included in free-bound emission. _fblvl dataset missing for He 3. [fiasco.collections]
WARNING: Li 1 not included in free-bound emission. _fblvl dataset missing for Li 1. [fiasco.collections]
WARNING: Li 2 not included in free-bound emission. _fblvl dataset missing for Li 2. [fiasco.collections]
WARNING: Li 3 not included in free-bound emission. _fblvl dataset missing for Li 3. [fiasco.collections]
WARNING: Li 4 not included in free-bound emission. _fblvl dataset missing for Li 4. [fiasco.collections]
WARNING: Be 1 not included in free-bound emission. _fblvl dataset missing for Be 1. [fiasco.collections]
WARNING: Be 2 not included in free-bound emission. _fblvl dataset missing for Be 2. [fiasco.collections]
WARNING: Be 3 not included in free-bound emission. _fblvl dataset missing for Be 3. [fiasco.collections]
WARNING: Be 4 not included in free-bound emission. _fblvl dataset missing for Be 4. [fiasco.collections]
WARNING: Be 5 not included in free-bound emission. _fblvl dataset missing for Be 5. [fiasco.collections]
WARNING: B 1 not included in free-bound emission. _fblvl dataset missing for B 1. [fiasco.collections]
WARNING: B 2 not included in free-bound emission. _fblvl dataset missing for B 2. [fiasco.collections]
WARNING: B 3 not included in free-bound emission. _fblvl dataset missing for B 3. [fiasco.collections]
WARNING: B 4 not included in free-bound emission. _fblvl dataset missing for B 4. [fiasco.collections]
WARNING: B 5 not included in free-bound emission. _fblvl dataset missing for B 5. [fiasco.collections]
WARNING: B 6 not included in free-bound emission. _fblvl dataset missing for B 6. [fiasco.collections]
WARNING: C 7 not included in free-bound emission. _fblvl dataset missing for C 7. [fiasco.collections]
WARNING: N 8 not included in free-bound emission. _fblvl dataset missing for N 8. [fiasco.collections]
WARNING: O 9 not included in free-bound emission. _fblvl dataset missing for O 9. [fiasco.collections]
WARNING: F 1 not included in free-bound emission. _fblvl dataset missing for F 1. [fiasco.collections]
WARNING: F 2 not included in free-bound emission. _fblvl dataset missing for F 2. [fiasco.collections]
WARNING: F 3 not included in free-bound emission. _fblvl dataset missing for F 3. [fiasco.collections]
WARNING: F 4 not included in free-bound emission. _fblvl dataset missing for F 4. [fiasco.collections]
WARNING: F 5 not included in free-bound emission. _fblvl dataset missing for F 5. [fiasco.collections]
WARNING: F 6 not included in free-bound emission. _fblvl dataset missing for F 6. [fiasco.collections]
WARNING: F 7 not included in free-bound emission. _fblvl dataset missing for F 7. [fiasco.collections]
WARNING: F 8 not included in free-bound emission. _fblvl dataset missing for F 8. [fiasco.collections]
WARNING: F 9 not included in free-bound emission. _fblvl dataset missing for F 9. [fiasco.collections]
WARNING: F 10 not included in free-bound emission. _fblvl dataset missing for F 10. [fiasco.collections]
WARNING: Ne 11 not included in free-bound emission. _fblvl dataset missing for Ne 11. [fiasco.collections]
WARNING: Na 1 not included in free-bound emission. _fblvl dataset missing for Na 1. [fiasco.collections]
WARNING: Na 2 not included in free-bound emission. _fblvl dataset missing for Na 2. [fiasco.collections]
WARNING: Na 3 not included in free-bound emission. _fblvl dataset missing for Na 3. [fiasco.collections]
WARNING: Na 4 not included in free-bound emission. _fblvl dataset missing for Na 4. [fiasco.collections]
WARNING: Na 5 not included in free-bound emission. _fblvl dataset missing for Na 5. [fiasco.collections]
WARNING: Na 6 not included in free-bound emission. _fblvl dataset missing for Na 6. [fiasco.collections]
WARNING: Na 7 not included in free-bound emission. _fblvl dataset missing for Na 7. [fiasco.collections]
WARNING: Na 8 not included in free-bound emission. _fblvl dataset missing for Na 8. [fiasco.collections]
WARNING: Na 9 not included in free-bound emission. _fblvl dataset missing for Na 9. [fiasco.collections]
WARNING: Na 10 not included in free-bound emission. _fblvl dataset missing for Na 10. [fiasco.collections]
WARNING: Na 11 not included in free-bound emission. _fblvl dataset missing for Na 11. [fiasco.collections]
WARNING: Na 12 not included in free-bound emission. _fblvl dataset missing for Na 12. [fiasco.collections]
WARNING: Mg 13 not included in free-bound emission. _fblvl dataset missing for Mg 13. [fiasco.collections]
WARNING: Al 2 not included in free-bound emission. _fblvl dataset missing for Al 2. [fiasco.collections]
WARNING: Al 14 not included in free-bound emission. _fblvl dataset missing for Al 14. [fiasco.collections]
WARNING: Si 15 not included in free-bound emission. _fblvl dataset missing for Si 15. [fiasco.collections]
WARNING: P 1 not included in free-bound emission. _fblvl dataset missing for P 1. [fiasco.collections]
WARNING: P 2 not included in free-bound emission. _fblvl dataset missing for P 2. [fiasco.collections]
WARNING: P 3 not included in free-bound emission. _fblvl dataset missing for P 3. [fiasco.collections]
WARNING: P 4 not included in free-bound emission. _fblvl dataset missing for P 4. [fiasco.collections]
WARNING: P 5 not included in free-bound emission. _fblvl dataset missing for P 5. [fiasco.collections]
WARNING: P 6 not included in free-bound emission. _fblvl dataset missing for P 6. [fiasco.collections]
WARNING: P 7 not included in free-bound emission. _fblvl dataset missing for P 7. [fiasco.collections]
WARNING: P 8 not included in free-bound emission. _fblvl dataset missing for P 8. [fiasco.collections]
WARNING: P 9 not included in free-bound emission. _fblvl dataset missing for P 9. [fiasco.collections]
WARNING: P 10 not included in free-bound emission. _fblvl dataset missing for P 10. [fiasco.collections]
WARNING: P 11 not included in free-bound emission. _fblvl dataset missing for P 11. [fiasco.collections]
WARNING: P 12 not included in free-bound emission. _fblvl dataset missing for P 12. [fiasco.collections]
WARNING: P 13 not included in free-bound emission. _fblvl dataset missing for P 13. [fiasco.collections]
WARNING: P 14 not included in free-bound emission. _fblvl dataset missing for P 14. [fiasco.collections]
WARNING: P 15 not included in free-bound emission. _fblvl dataset missing for P 15. [fiasco.collections]
WARNING: P 16 not included in free-bound emission. _fblvl dataset missing for P 16. [fiasco.collections]
WARNING: S 17 not included in free-bound emission. _fblvl dataset missing for S 17. [fiasco.collections]
WARNING: Cl 1 not included in free-bound emission. _fblvl dataset missing for Cl 1. [fiasco.collections]
WARNING: Cl 2 not included in free-bound emission. _fblvl dataset missing for Cl 2. [fiasco.collections]
WARNING: Cl 3 not included in free-bound emission. _fblvl dataset missing for Cl 3. [fiasco.collections]
WARNING: Cl 4 not included in free-bound emission. _fblvl dataset missing for Cl 4. [fiasco.collections]
WARNING: Cl 5 not included in free-bound emission. _fblvl dataset missing for Cl 5. [fiasco.collections]
WARNING: Cl 6 not included in free-bound emission. _fblvl dataset missing for Cl 6. [fiasco.collections]
WARNING: Cl 7 not included in free-bound emission. _fblvl dataset missing for Cl 7. [fiasco.collections]
WARNING: Cl 8 not included in free-bound emission. _fblvl dataset missing for Cl 8. [fiasco.collections]
WARNING: Cl 9 not included in free-bound emission. _fblvl dataset missing for Cl 9. [fiasco.collections]
WARNING: Cl 10 not included in free-bound emission. _fblvl dataset missing for Cl 10. [fiasco.collections]
WARNING: Cl 11 not included in free-bound emission. _fblvl dataset missing for Cl 11. [fiasco.collections]
WARNING: Cl 12 not included in free-bound emission. _fblvl dataset missing for Cl 12. [fiasco.collections]
WARNING: Cl 13 not included in free-bound emission. _fblvl dataset missing for Cl 13. [fiasco.collections]
WARNING: Cl 14 not included in free-bound emission. _fblvl dataset missing for Cl 14. [fiasco.collections]
WARNING: Cl 15 not included in free-bound emission. _fblvl dataset missing for Cl 15. [fiasco.collections]
WARNING: Cl 16 not included in free-bound emission. _fblvl dataset missing for Cl 16. [fiasco.collections]
WARNING: Cl 17 not included in free-bound emission. _fblvl dataset missing for Cl 17. [fiasco.collections]
WARNING: Cl 18 not included in free-bound emission. _fblvl dataset missing for Cl 18. [fiasco.collections]
WARNING: Ar 2 not included in free-bound emission. _fblvl dataset missing for Ar 2. [fiasco.collections]
WARNING: Ar 3 not included in free-bound emission. _fblvl dataset missing for Ar 3. [fiasco.collections]
WARNING: Ar 5 not included in free-bound emission. _fblvl dataset missing for Ar 5. [fiasco.collections]
WARNING: Ar 6 not included in free-bound emission. _fblvl dataset missing for Ar 6. [fiasco.collections]
WARNING: Ar 19 not included in free-bound emission. _fblvl dataset missing for Ar 19. [fiasco.collections]
WARNING: K 1 not included in free-bound emission. _fblvl dataset missing for K 1. [fiasco.collections]
WARNING: K 2 not included in free-bound emission. _fblvl dataset missing for K 2. [fiasco.collections]
WARNING: K 3 not included in free-bound emission. _fblvl dataset missing for K 3. [fiasco.collections]
WARNING: K 4 not included in free-bound emission. _fblvl dataset missing for K 4. [fiasco.collections]
WARNING: K 5 not included in free-bound emission. _fblvl dataset missing for K 5. [fiasco.collections]
WARNING: K 6 not included in free-bound emission. _fblvl dataset missing for K 6. [fiasco.collections]
WARNING: K 7 not included in free-bound emission. _fblvl dataset missing for K 7. [fiasco.collections]
WARNING: K 8 not included in free-bound emission. _fblvl dataset missing for K 8. [fiasco.collections]
WARNING: K 9 not included in free-bound emission. _fblvl dataset missing for K 9. [fiasco.collections]
WARNING: K 10 not included in free-bound emission. _fblvl dataset missing for K 10. [fiasco.collections]
WARNING: K 11 not included in free-bound emission. _fblvl dataset missing for K 11. [fiasco.collections]
WARNING: K 12 not included in free-bound emission. _fblvl dataset missing for K 12. [fiasco.collections]
WARNING: K 13 not included in free-bound emission. _fblvl dataset missing for K 13. [fiasco.collections]
WARNING: K 14 not included in free-bound emission. _fblvl dataset missing for K 14. [fiasco.collections]
WARNING: K 15 not included in free-bound emission. _fblvl dataset missing for K 15. [fiasco.collections]
WARNING: K 16 not included in free-bound emission. _fblvl dataset missing for K 16. [fiasco.collections]
WARNING: K 17 not included in free-bound emission. _fblvl dataset missing for K 17. [fiasco.collections]
WARNING: K 18 not included in free-bound emission. _fblvl dataset missing for K 18. [fiasco.collections]
WARNING: K 19 not included in free-bound emission. _fblvl dataset missing for K 19. [fiasco.collections]
WARNING: K 20 not included in free-bound emission. _fblvl dataset missing for K 20. [fiasco.collections]
WARNING: Ca 1 not included in free-bound emission. _fblvl dataset missing for Ca 1. [fiasco.collections]
WARNING: Ca 2 not included in free-bound emission. _fblvl dataset missing for Ca 2. [fiasco.collections]
WARNING: Ca 3 not included in free-bound emission. _fblvl dataset missing for Ca 3. [fiasco.collections]
WARNING: Ca 4 not included in free-bound emission. _fblvl dataset missing for Ca 4. [fiasco.collections]
WARNING: Ca 5 not included in free-bound emission. _fblvl dataset missing for Ca 5. [fiasco.collections]
WARNING: Ca 6 not included in free-bound emission. _fblvl dataset missing for Ca 6. [fiasco.collections]
WARNING: Ca 7 not included in free-bound emission. _fblvl dataset missing for Ca 7. [fiasco.collections]
WARNING: Ca 8 not included in free-bound emission. _fblvl dataset missing for Ca 8. [fiasco.collections]
WARNING: Ca 21 not included in free-bound emission. _fblvl dataset missing for Ca 21. [fiasco.collections]
WARNING: Sc 1 not included in free-bound emission. _fblvl dataset missing for Sc 1. [fiasco.collections]
WARNING: Sc 2 not included in free-bound emission. _fblvl dataset missing for Sc 2. [fiasco.collections]
WARNING: Sc 3 not included in free-bound emission. _fblvl dataset missing for Sc 3. [fiasco.collections]
WARNING: Sc 4 not included in free-bound emission. _fblvl dataset missing for Sc 4. [fiasco.collections]
WARNING: Sc 5 not included in free-bound emission. _fblvl dataset missing for Sc 5. [fiasco.collections]
WARNING: Sc 6 not included in free-bound emission. _fblvl dataset missing for Sc 6. [fiasco.collections]
WARNING: Sc 7 not included in free-bound emission. _fblvl dataset missing for Sc 7. [fiasco.collections]
WARNING: Sc 8 not included in free-bound emission. _fblvl dataset missing for Sc 8. [fiasco.collections]
WARNING: Sc 9 not included in free-bound emission. _fblvl dataset missing for Sc 9. [fiasco.collections]
WARNING: Sc 10 not included in free-bound emission. _fblvl dataset missing for Sc 10. [fiasco.collections]
WARNING: Sc 11 not included in free-bound emission. _fblvl dataset missing for Sc 11. [fiasco.collections]
WARNING: Sc 12 not included in free-bound emission. _fblvl dataset missing for Sc 12. [fiasco.collections]
WARNING: Sc 13 not included in free-bound emission. _fblvl dataset missing for Sc 13. [fiasco.collections]
WARNING: Sc 14 not included in free-bound emission. _fblvl dataset missing for Sc 14. [fiasco.collections]
WARNING: Sc 15 not included in free-bound emission. _fblvl dataset missing for Sc 15. [fiasco.collections]
WARNING: Sc 16 not included in free-bound emission. _fblvl dataset missing for Sc 16. [fiasco.collections]
WARNING: Sc 17 not included in free-bound emission. _fblvl dataset missing for Sc 17. [fiasco.collections]
WARNING: Sc 18 not included in free-bound emission. _fblvl dataset missing for Sc 18. [fiasco.collections]
WARNING: Sc 19 not included in free-bound emission. _fblvl dataset missing for Sc 19. [fiasco.collections]
WARNING: Sc 20 not included in free-bound emission. _fblvl dataset missing for Sc 20. [fiasco.collections]
WARNING: Sc 21 not included in free-bound emission. _fblvl dataset missing for Sc 21. [fiasco.collections]
WARNING: Sc 22 not included in free-bound emission. _fblvl dataset missing for Sc 22. [fiasco.collections]
WARNING: Ti 1 not included in free-bound emission. _fblvl dataset missing for Ti 1. [fiasco.collections]
WARNING: Ti 2 not included in free-bound emission. _fblvl dataset missing for Ti 2. [fiasco.collections]
WARNING: Ti 3 not included in free-bound emission. _fblvl dataset missing for Ti 3. [fiasco.collections]
WARNING: Ti 4 not included in free-bound emission. _fblvl dataset missing for Ti 4. [fiasco.collections]
WARNING: Ti 5 not included in free-bound emission. _fblvl dataset missing for Ti 5. [fiasco.collections]
WARNING: Ti 6 not included in free-bound emission. _fblvl dataset missing for Ti 6. [fiasco.collections]
WARNING: Ti 7 not included in free-bound emission. _fblvl dataset missing for Ti 7. [fiasco.collections]
WARNING: Ti 8 not included in free-bound emission. _fblvl dataset missing for Ti 8. [fiasco.collections]
WARNING: Ti 9 not included in free-bound emission. _fblvl dataset missing for Ti 9. [fiasco.collections]
WARNING: Ti 10 not included in free-bound emission. _fblvl dataset missing for Ti 10. [fiasco.collections]
WARNING: Ti 11 not included in free-bound emission. _fblvl dataset missing for Ti 11. [fiasco.collections]
WARNING: Ti 12 not included in free-bound emission. _fblvl dataset missing for Ti 12. [fiasco.collections]
WARNING: Ti 13 not included in free-bound emission. _fblvl dataset missing for Ti 13. [fiasco.collections]
WARNING: Ti 14 not included in free-bound emission. _fblvl dataset missing for Ti 14. [fiasco.collections]
WARNING: Ti 15 not included in free-bound emission. _fblvl dataset missing for Ti 15. [fiasco.collections]
WARNING: Ti 16 not included in free-bound emission. _fblvl dataset missing for Ti 16. [fiasco.collections]
WARNING: Ti 17 not included in free-bound emission. _fblvl dataset missing for Ti 17. [fiasco.collections]
WARNING: Ti 18 not included in free-bound emission. _fblvl dataset missing for Ti 18. [fiasco.collections]
WARNING: Ti 19 not included in free-bound emission. _fblvl dataset missing for Ti 19. [fiasco.collections]
WARNING: Ti 20 not included in free-bound emission. _fblvl dataset missing for Ti 20. [fiasco.collections]
WARNING: Ti 21 not included in free-bound emission. _fblvl dataset missing for Ti 21. [fiasco.collections]
WARNING: Ti 22 not included in free-bound emission. _fblvl dataset missing for Ti 22. [fiasco.collections]
WARNING: Ti 23 not included in free-bound emission. _fblvl dataset missing for Ti 23. [fiasco.collections]
WARNING: V 1 not included in free-bound emission. _fblvl dataset missing for V 1. [fiasco.collections]
WARNING: V 2 not included in free-bound emission. _fblvl dataset missing for V 2. [fiasco.collections]
WARNING: V 3 not included in free-bound emission. _fblvl dataset missing for V 3. [fiasco.collections]
WARNING: V 4 not included in free-bound emission. _fblvl dataset missing for V 4. [fiasco.collections]
WARNING: V 5 not included in free-bound emission. _fblvl dataset missing for V 5. [fiasco.collections]
WARNING: V 6 not included in free-bound emission. _fblvl dataset missing for V 6. [fiasco.collections]
WARNING: V 7 not included in free-bound emission. _fblvl dataset missing for V 7. [fiasco.collections]
WARNING: V 8 not included in free-bound emission. _fblvl dataset missing for V 8. [fiasco.collections]
WARNING: V 9 not included in free-bound emission. _fblvl dataset missing for V 9. [fiasco.collections]
WARNING: V 10 not included in free-bound emission. _fblvl dataset missing for V 10. [fiasco.collections]
WARNING: V 11 not included in free-bound emission. _fblvl dataset missing for V 11. [fiasco.collections]
WARNING: V 12 not included in free-bound emission. _fblvl dataset missing for V 12. [fiasco.collections]
WARNING: V 13 not included in free-bound emission. _fblvl dataset missing for V 13. [fiasco.collections]
WARNING: V 14 not included in free-bound emission. _fblvl dataset missing for V 14. [fiasco.collections]
WARNING: V 15 not included in free-bound emission. _fblvl dataset missing for V 15. [fiasco.collections]
WARNING: V 16 not included in free-bound emission. _fblvl dataset missing for V 16. [fiasco.collections]
WARNING: V 17 not included in free-bound emission. _fblvl dataset missing for V 17. [fiasco.collections]
WARNING: V 18 not included in free-bound emission. _fblvl dataset missing for V 18. [fiasco.collections]
WARNING: V 19 not included in free-bound emission. _fblvl dataset missing for V 19. [fiasco.collections]
WARNING: V 20 not included in free-bound emission. _fblvl dataset missing for V 20. [fiasco.collections]
WARNING: V 21 not included in free-bound emission. _fblvl dataset missing for V 21. [fiasco.collections]
WARNING: V 22 not included in free-bound emission. _fblvl dataset missing for V 22. [fiasco.collections]
WARNING: V 23 not included in free-bound emission. _fblvl dataset missing for V 23. [fiasco.collections]
WARNING: V 24 not included in free-bound emission. _fblvl dataset missing for V 24. [fiasco.collections]
WARNING: Cr 1 not included in free-bound emission. _fblvl dataset missing for Cr 1. [fiasco.collections]
WARNING: Cr 2 not included in free-bound emission. _fblvl dataset missing for Cr 2. [fiasco.collections]
WARNING: Cr 3 not included in free-bound emission. _fblvl dataset missing for Cr 3. [fiasco.collections]
WARNING: Cr 4 not included in free-bound emission. _fblvl dataset missing for Cr 4. [fiasco.collections]
WARNING: Cr 5 not included in free-bound emission. _fblvl dataset missing for Cr 5. [fiasco.collections]
WARNING: Cr 6 not included in free-bound emission. _fblvl dataset missing for Cr 6. [fiasco.collections]
WARNING: Cr 7 not included in free-bound emission. _fblvl dataset missing for Cr 7. [fiasco.collections]
WARNING: Cr 8 not included in free-bound emission. _fblvl dataset missing for Cr 8. [fiasco.collections]
WARNING: Cr 9 not included in free-bound emission. _fblvl dataset missing for Cr 9. [fiasco.collections]
WARNING: Cr 10 not included in free-bound emission. _fblvl dataset missing for Cr 10. [fiasco.collections]
WARNING: Cr 11 not included in free-bound emission. _fblvl dataset missing for Cr 11. [fiasco.collections]
WARNING: Cr 12 not included in free-bound emission. _fblvl dataset missing for Cr 12. [fiasco.collections]
WARNING: Cr 13 not included in free-bound emission. _fblvl dataset missing for Cr 13. [fiasco.collections]
WARNING: Cr 14 not included in free-bound emission. _fblvl dataset missing for Cr 14. [fiasco.collections]
WARNING: Cr 15 not included in free-bound emission. _fblvl dataset missing for Cr 15. [fiasco.collections]
WARNING: Cr 16 not included in free-bound emission. _fblvl dataset missing for Cr 16. [fiasco.collections]
WARNING: Cr 17 not included in free-bound emission. _fblvl dataset missing for Cr 17. [fiasco.collections]
WARNING: Cr 18 not included in free-bound emission. _fblvl dataset missing for Cr 18. [fiasco.collections]
WARNING: Cr 19 not included in free-bound emission. _fblvl dataset missing for Cr 19. [fiasco.collections]
WARNING: Cr 20 not included in free-bound emission. _fblvl dataset missing for Cr 20. [fiasco.collections]
WARNING: Cr 21 not included in free-bound emission. _fblvl dataset missing for Cr 21. [fiasco.collections]
WARNING: Cr 22 not included in free-bound emission. _fblvl dataset missing for Cr 22. [fiasco.collections]
WARNING: Cr 23 not included in free-bound emission. _fblvl dataset missing for Cr 23. [fiasco.collections]
WARNING: Cr 24 not included in free-bound emission. _fblvl dataset missing for Cr 24. [fiasco.collections]
WARNING: Cr 25 not included in free-bound emission. _fblvl dataset missing for Cr 25. [fiasco.collections]
WARNING: Mn 1 not included in free-bound emission. _fblvl dataset missing for Mn 1. [fiasco.collections]
WARNING: Mn 2 not included in free-bound emission. _fblvl dataset missing for Mn 2. [fiasco.collections]
WARNING: Mn 3 not included in free-bound emission. _fblvl dataset missing for Mn 3. [fiasco.collections]
WARNING: Mn 4 not included in free-bound emission. _fblvl dataset missing for Mn 4. [fiasco.collections]
WARNING: Mn 5 not included in free-bound emission. _fblvl dataset missing for Mn 5. [fiasco.collections]
WARNING: Mn 6 not included in free-bound emission. _fblvl dataset missing for Mn 6. [fiasco.collections]
WARNING: Mn 7 not included in free-bound emission. _fblvl dataset missing for Mn 7. [fiasco.collections]
WARNING: Mn 8 not included in free-bound emission. _fblvl dataset missing for Mn 8. [fiasco.collections]
WARNING: Mn 9 not included in free-bound emission. _fblvl dataset missing for Mn 9. [fiasco.collections]
WARNING: Mn 10 not included in free-bound emission. _fblvl dataset missing for Mn 10. [fiasco.collections]
WARNING: Mn 11 not included in free-bound emission. _fblvl dataset missing for Mn 11. [fiasco.collections]
WARNING: Mn 12 not included in free-bound emission. _fblvl dataset missing for Mn 12. [fiasco.collections]
WARNING: Mn 13 not included in free-bound emission. _fblvl dataset missing for Mn 13. [fiasco.collections]
WARNING: Mn 14 not included in free-bound emission. _fblvl dataset missing for Mn 14. [fiasco.collections]
WARNING: Mn 15 not included in free-bound emission. _fblvl dataset missing for Mn 15. [fiasco.collections]
WARNING: Mn 16 not included in free-bound emission. _fblvl dataset missing for Mn 16. [fiasco.collections]
WARNING: Mn 17 not included in free-bound emission. _fblvl dataset missing for Mn 17. [fiasco.collections]
WARNING: Mn 18 not included in free-bound emission. _fblvl dataset missing for Mn 18. [fiasco.collections]
WARNING: Mn 19 not included in free-bound emission. _fblvl dataset missing for Mn 19. [fiasco.collections]
WARNING: Mn 20 not included in free-bound emission. _fblvl dataset missing for Mn 20. [fiasco.collections]
WARNING: Mn 21 not included in free-bound emission. _fblvl dataset missing for Mn 21. [fiasco.collections]
WARNING: Mn 22 not included in free-bound emission. _fblvl dataset missing for Mn 22. [fiasco.collections]
WARNING: Mn 23 not included in free-bound emission. _fblvl dataset missing for Mn 23. [fiasco.collections]
WARNING: Mn 24 not included in free-bound emission. _fblvl dataset missing for Mn 24. [fiasco.collections]
WARNING: Mn 25 not included in free-bound emission. _fblvl dataset missing for Mn 25. [fiasco.collections]
WARNING: Mn 26 not included in free-bound emission. _fblvl dataset missing for Mn 26. [fiasco.collections]
WARNING: Fe 1 not included in free-bound emission. _fblvl dataset missing for Fe 1. [fiasco.collections]
WARNING: Fe 2 not included in free-bound emission. _fblvl dataset missing for Fe 2. [fiasco.collections]
WARNING: Fe 3 not included in free-bound emission. _fblvl dataset missing for Fe 3. [fiasco.collections]
WARNING: Fe 27 not included in free-bound emission. _fblvl dataset missing for Fe 27. [fiasco.collections]
WARNING: Co 1 not included in free-bound emission. _fblvl dataset missing for Co 1. [fiasco.collections]
WARNING: Co 2 not included in free-bound emission. _fblvl dataset missing for Co 2. [fiasco.collections]
WARNING: Co 3 not included in free-bound emission. _fblvl dataset missing for Co 3. [fiasco.collections]
WARNING: Co 4 not included in free-bound emission. _fblvl dataset missing for Co 4. [fiasco.collections]
WARNING: Co 5 not included in free-bound emission. _fblvl dataset missing for Co 5. [fiasco.collections]
WARNING: Co 6 not included in free-bound emission. _fblvl dataset missing for Co 6. [fiasco.collections]
WARNING: Co 7 not included in free-bound emission. _fblvl dataset missing for Co 7. [fiasco.collections]
WARNING: Co 8 not included in free-bound emission. _fblvl dataset missing for Co 8. [fiasco.collections]
WARNING: Co 9 not included in free-bound emission. _fblvl dataset missing for Co 9. [fiasco.collections]
WARNING: Co 10 not included in free-bound emission. _fblvl dataset missing for Co 10. [fiasco.collections]
WARNING: Co 11 not included in free-bound emission. _fblvl dataset missing for Co 11. [fiasco.collections]
WARNING: Co 12 not included in free-bound emission. _fblvl dataset missing for Co 12. [fiasco.collections]
WARNING: Co 13 not included in free-bound emission. _fblvl dataset missing for Co 13. [fiasco.collections]
WARNING: Co 14 not included in free-bound emission. _fblvl dataset missing for Co 14. [fiasco.collections]
WARNING: Co 15 not included in free-bound emission. _fblvl dataset missing for Co 15. [fiasco.collections]
WARNING: Co 16 not included in free-bound emission. _fblvl dataset missing for Co 16. [fiasco.collections]
WARNING: Co 17 not included in free-bound emission. _fblvl dataset missing for Co 17. [fiasco.collections]
WARNING: Co 18 not included in free-bound emission. _fblvl dataset missing for Co 18. [fiasco.collections]
WARNING: Co 19 not included in free-bound emission. _fblvl dataset missing for Co 19. [fiasco.collections]
WARNING: Co 20 not included in free-bound emission. _fblvl dataset missing for Co 20. [fiasco.collections]
WARNING: Co 21 not included in free-bound emission. _fblvl dataset missing for Co 21. [fiasco.collections]
WARNING: Co 22 not included in free-bound emission. _fblvl dataset missing for Co 22. [fiasco.collections]
WARNING: Co 23 not included in free-bound emission. _fblvl dataset missing for Co 23. [fiasco.collections]
WARNING: Co 24 not included in free-bound emission. _fblvl dataset missing for Co 24. [fiasco.collections]
WARNING: Co 25 not included in free-bound emission. _fblvl dataset missing for Co 25. [fiasco.collections]
WARNING: Co 26 not included in free-bound emission. _fblvl dataset missing for Co 26. [fiasco.collections]
WARNING: Co 27 not included in free-bound emission. _fblvl dataset missing for Co 27. [fiasco.collections]
WARNING: Co 28 not included in free-bound emission. _fblvl dataset missing for Co 28. [fiasco.collections]
WARNING: Ni 1 not included in free-bound emission. _fblvl dataset missing for Ni 1. [fiasco.collections]
WARNING: Ni 2 not included in free-bound emission. _fblvl dataset missing for Ni 2. [fiasco.collections]
WARNING: Ni 3 not included in free-bound emission. _fblvl dataset missing for Ni 3. [fiasco.collections]
WARNING: Ni 4 not included in free-bound emission. _fblvl dataset missing for Ni 4. [fiasco.collections]
WARNING: Ni 5 not included in free-bound emission. _fblvl dataset missing for Ni 5. [fiasco.collections]
WARNING: Ni 6 not included in free-bound emission. _fblvl dataset missing for Ni 6. [fiasco.collections]
WARNING: Ni 7 not included in free-bound emission. _fblvl dataset missing for Ni 7. [fiasco.collections]
WARNING: Ni 8 not included in free-bound emission. _fblvl dataset missing for Ni 8. [fiasco.collections]
WARNING: Ni 9 not included in free-bound emission. _fblvl dataset missing for Ni 9. [fiasco.collections]
WARNING: Ni 10 not included in free-bound emission. _fblvl dataset missing for Ni 10. [fiasco.collections]
WARNING: Ni 14 not included in free-bound emission. _fblvl dataset missing for Ni 14. [fiasco.collections]
WARNING: Ni 29 not included in free-bound emission. _fblvl dataset missing for Ni 29. [fiasco.collections]
WARNING: Cu 1 not included in free-bound emission. _fblvl dataset missing for Cu 1. [fiasco.collections]
WARNING: Cu 2 not included in free-bound emission. _fblvl dataset missing for Cu 2. [fiasco.collections]
WARNING: Cu 3 not included in free-bound emission. _fblvl dataset missing for Cu 3. [fiasco.collections]
WARNING: Cu 4 not included in free-bound emission. _fblvl dataset missing for Cu 4. [fiasco.collections]
WARNING: Cu 5 not included in free-bound emission. _fblvl dataset missing for Cu 5. [fiasco.collections]
WARNING: Cu 6 not included in free-bound emission. _fblvl dataset missing for Cu 6. [fiasco.collections]
WARNING: Cu 7 not included in free-bound emission. _fblvl dataset missing for Cu 7. [fiasco.collections]
WARNING: Cu 8 not included in free-bound emission. _fblvl dataset missing for Cu 8. [fiasco.collections]
WARNING: Cu 9 not included in free-bound emission. _fblvl dataset missing for Cu 9. [fiasco.collections]
WARNING: Cu 10 not included in free-bound emission. _fblvl dataset missing for Cu 10. [fiasco.collections]
WARNING: Cu 11 not included in free-bound emission. _fblvl dataset missing for Cu 11. [fiasco.collections]
WARNING: Cu 12 not included in free-bound emission. _fblvl dataset missing for Cu 12. [fiasco.collections]
WARNING: Cu 13 not included in free-bound emission. _fblvl dataset missing for Cu 13. [fiasco.collections]
WARNING: Cu 14 not included in free-bound emission. _fblvl dataset missing for Cu 14. [fiasco.collections]
WARNING: Cu 15 not included in free-bound emission. _fblvl dataset missing for Cu 15. [fiasco.collections]
WARNING: Cu 16 not included in free-bound emission. _fblvl dataset missing for Cu 16. [fiasco.collections]
WARNING: Cu 17 not included in free-bound emission. _fblvl dataset missing for Cu 17. [fiasco.collections]
WARNING: Cu 18 not included in free-bound emission. _fblvl dataset missing for Cu 18. [fiasco.collections]
WARNING: Cu 19 not included in free-bound emission. _fblvl dataset missing for Cu 19. [fiasco.collections]
WARNING: Cu 20 not included in free-bound emission. _fblvl dataset missing for Cu 20. [fiasco.collections]
WARNING: Cu 21 not included in free-bound emission. _fblvl dataset missing for Cu 21. [fiasco.collections]
WARNING: Cu 22 not included in free-bound emission. _fblvl dataset missing for Cu 22. [fiasco.collections]
WARNING: Cu 23 not included in free-bound emission. _fblvl dataset missing for Cu 23. [fiasco.collections]
WARNING: Cu 24 not included in free-bound emission. _fblvl dataset missing for Cu 24. [fiasco.collections]
WARNING: Cu 25 not included in free-bound emission. _fblvl dataset missing for Cu 25. [fiasco.collections]
WARNING: Cu 26 not included in free-bound emission. _fblvl dataset missing for Cu 26. [fiasco.collections]
WARNING: Cu 27 not included in free-bound emission. _fblvl dataset missing for Cu 27. [fiasco.collections]
WARNING: Cu 28 not included in free-bound emission. _fblvl dataset missing for Cu 28. [fiasco.collections]
WARNING: Cu 29 not included in free-bound emission. _fblvl dataset missing for Cu 29. [fiasco.collections]
WARNING: Cu 30 not included in free-bound emission. _fblvl dataset missing for Cu 30. [fiasco.collections]
WARNING: Zn 1 not included in free-bound emission. _fblvl dataset missing for Zn 1. [fiasco.collections]
WARNING: Zn 2 not included in free-bound emission. _fblvl dataset missing for Zn 2. [fiasco.collections]
WARNING: Zn 3 not included in free-bound emission. _fblvl dataset missing for Zn 3. [fiasco.collections]
WARNING: Zn 4 not included in free-bound emission. _fblvl dataset missing for Zn 4. [fiasco.collections]
WARNING: Zn 5 not included in free-bound emission. _fblvl dataset missing for Zn 5. [fiasco.collections]
WARNING: Zn 6 not included in free-bound emission. _fblvl dataset missing for Zn 6. [fiasco.collections]
WARNING: Zn 7 not included in free-bound emission. _fblvl dataset missing for Zn 7. [fiasco.collections]
WARNING: Zn 8 not included in free-bound emission. _fblvl dataset missing for Zn 8. [fiasco.collections]
WARNING: Zn 9 not included in free-bound emission. _fblvl dataset missing for Zn 9. [fiasco.collections]
WARNING: Zn 10 not included in free-bound emission. _fblvl dataset missing for Zn 10. [fiasco.collections]
WARNING: Zn 11 not included in free-bound emission. _fblvl dataset missing for Zn 11. [fiasco.collections]
WARNING: Zn 12 not included in free-bound emission. _fblvl dataset missing for Zn 12. [fiasco.collections]
WARNING: Zn 13 not included in free-bound emission. _fblvl dataset missing for Zn 13. [fiasco.collections]
WARNING: Zn 14 not included in free-bound emission. _fblvl dataset missing for Zn 14. [fiasco.collections]
WARNING: Zn 15 not included in free-bound emission. _fblvl dataset missing for Zn 15. [fiasco.collections]
WARNING: Zn 16 not included in free-bound emission. _fblvl dataset missing for Zn 16. [fiasco.collections]
WARNING: Zn 17 not included in free-bound emission. _fblvl dataset missing for Zn 17. [fiasco.collections]
WARNING: Zn 18 not included in free-bound emission. _fblvl dataset missing for Zn 18. [fiasco.collections]
WARNING: Zn 19 not included in free-bound emission. _fblvl dataset missing for Zn 19. [fiasco.collections]
WARNING: Zn 20 not included in free-bound emission. _fblvl dataset missing for Zn 20. [fiasco.collections]
WARNING: Zn 21 not included in free-bound emission. _fblvl dataset missing for Zn 21. [fiasco.collections]
WARNING: Zn 22 not included in free-bound emission. _fblvl dataset missing for Zn 22. [fiasco.collections]
WARNING: Zn 23 not included in free-bound emission. _fblvl dataset missing for Zn 23. [fiasco.collections]
WARNING: Zn 24 not included in free-bound emission. _fblvl dataset missing for Zn 24. [fiasco.collections]
WARNING: Zn 25 not included in free-bound emission. _fblvl dataset missing for Zn 25. [fiasco.collections]
WARNING: Zn 26 not included in free-bound emission. _fblvl dataset missing for Zn 26. [fiasco.collections]
WARNING: Zn 27 not included in free-bound emission. _fblvl dataset missing for Zn 27. [fiasco.collections]
WARNING: Zn 28 not included in free-bound emission. _fblvl dataset missing for Zn 28. [fiasco.collections]
WARNING: Zn 29 not included in free-bound emission. _fblvl dataset missing for Zn 29. [fiasco.collections]
WARNING: Zn 30 not included in free-bound emission. _fblvl dataset missing for Zn 30. [fiasco.collections]
WARNING: Zn 31 not included in free-bound emission. _fblvl dataset missing for Zn 31. [fiasco.collections]
IDL code to produce free-free result:
; set common block
common elements, abund, abund_ref, ioneq, ioneq_logt, ioneq_ref
; read abundance and ionization equilibrium
abund_file = FILEPATH('sun_coronal_1992_feldman_ext.abund', ROOT_DIR=!xuvtop, SUBDIR='abundance')
ioneq_file = FILEPATH('chianti.ioneq', ROOT_DIR=!xuvtop, SUBDIR='ioneq')
read_abund, abund_file, abund, abund_ref
read_ioneq, ioneq_file, ioneq_logt, ioneq, ioneq_ref
; set temperature and wavelength
temperature = [100000d / 1d,3855226053388321d / 34359738368d,4325634777407421d / 34359738368d,1213360511705375d / 8589934592d,5445651543418815d / 34359738368d,3055060763706939d / 17179869184d,6855669111786069d / 34359738368d,1923046815002979d / 8589934592d,8630776059679839d / 34359738368d,1210486251738545d / 4294967296d,1358187913129459d / 4294967296d,3047823805882249d / 8589934592d,3419714555660727d / 8589934592d,3836982839896381d / 8589934592d,8610331110407387d / 17179869184d,4830475201763701d / 8589934592d,2709941159713435d / 4294967296d,6081207982541875d / 8589934592d,6823227580860407d / 8589934592d,1913946815906879d / 2147483648d,1000000d / 1d,4819032566735401d / 4294967296d,1351760867939819d / 1073741824d,6066802558526875d / 4294967296d,6807064429273519d / 4294967296d,1909412977316837d / 1073741824d,4284793194866293d / 2147483648d,600952129688431d / 268435456d,1348558759324975d / 536870912d,6052431258692725d / 2147483648d,106108430713239d / 33554432d,3809779757352811d / 1073741824d,4274643194575909d / 1073741824d,1199057137467619d / 268435456d,5381456944004617d / 1073741824d,3019047001102313d / 536870912d,1693713224820897d / 268435456d,7601509978177343d / 1073741824d,8529034476075509d / 1073741824d,4784867039767197d / 536870912d,10000000d / 1d,94121729819051d / 8388608d,844850542462387d / 67108864d,3791751599079297d / 268435456d,8508830536591899d / 536870912d,1193383110823023d / 67108864d,5355991493582877d / 268435456d,6009521296884309d / 268435456d,6742793796624875d / 268435456d,3782769536682953d / 134217728d,8488674457059119d / 268435456d,297639043543189d / 8388608d,2671651996609943d / 67108864d,5995285687338095d / 134217728d,6726821180005771d / 134217728d,3773808751377891d / 67108864d,4234283062052251d / 67108864d,593867967045105d / 8388608d,5330646547547193d / 67108864d,5981083799708997d / 67108864d,100000000d / 1d]
wavelength = [25d / 1d,26d / 1d,27d / 1d,28d / 1d,29d / 1d,30d / 1d,31d / 1d,32d / 1d,33d / 1d,34d / 1d,35d / 1d,36d / 1d,37d / 1d,38d / 1d,39d / 1d,40d / 1d,41d / 1d,42d / 1d,43d / 1d,44d / 1d,45d / 1d,46d / 1d,47d / 1d,48d / 1d,49d / 1d,50d / 1d,51d / 1d,52d / 1d,53d / 1d,54d / 1d,55d / 1d,56d / 1d,57d / 1d,58d / 1d,59d / 1d,60d / 1d,61d / 1d,62d / 1d,63d / 1d,64d / 1d,65d / 1d,66d / 1d,67d / 1d,68d / 1d,69d / 1d,70d / 1d,71d / 1d,72d / 1d,73d / 1d,74d / 1d,75d / 1d,76d / 1d,77d / 1d,78d / 1d,79d / 1d,80d / 1d,81d / 1d,82d / 1d,83d / 1d,84d / 1d,85d / 1d,86d / 1d,87d / 1d,88d / 1d,89d / 1d,90d / 1d,91d / 1d,92d / 1d,93d / 1d,94d / 1d,95d / 1d,96d / 1d,97d / 1d,98d / 1d,99d / 1d,100d / 1d,101d / 1d,102d / 1d,103d / 1d,104d / 1d,105d / 1d,106d / 1d,107d / 1d,108d / 1d,109d / 1d,110d / 1d,111d / 1d,112d / 1d,113d / 1d,114d / 1d,115d / 1d,116d / 1d,117d / 1d,118d / 1d,119d / 1d,120d / 1d,121d / 1d,122d / 1d,123d / 1d,124d / 1d,125d / 1d,126d / 1d,127d / 1d,128d / 1d,129d / 1d,130d / 1d,131d / 1d,132d / 1d,133d / 1d,134d / 1d,135d / 1d,136d / 1d,137d / 1d,138d / 1d,139d / 1d,140d / 1d,141d / 1d,142d / 1d,143d / 1d,144d / 1d,145d / 1d,146d / 1d,147d / 1d,148d / 1d,149d / 1d,150d / 1d,151d / 1d,152d / 1d,153d / 1d,154d / 1d,155d / 1d,156d / 1d,157d / 1d,158d / 1d,159d / 1d,160d / 1d,161d / 1d,162d / 1d,163d / 1d,164d / 1d,165d / 1d,166d / 1d,167d / 1d,168d / 1d,169d / 1d,170d / 1d,171d / 1d,172d / 1d,173d / 1d,174d / 1d,175d / 1d,176d / 1d,177d / 1d,178d / 1d,179d / 1d,180d / 1d,181d / 1d,182d / 1d,183d / 1d,184d / 1d,185d / 1d,186d / 1d,187d / 1d,188d / 1d,189d / 1d,190d / 1d,191d / 1d,192d / 1d,193d / 1d,194d / 1d,195d / 1d,196d / 1d,197d / 1d,198d / 1d,199d / 1d,200d / 1d,201d / 1d,202d / 1d,203d / 1d,204d / 1d,205d / 1d,206d / 1d,207d / 1d,208d / 1d,209d / 1d,210d / 1d,211d / 1d,212d / 1d,213d / 1d,214d / 1d,215d / 1d,216d / 1d,217d / 1d,218d / 1d,219d / 1d,220d / 1d,221d / 1d,222d / 1d,223d / 1d,224d / 1d,225d / 1d,226d / 1d,227d / 1d,228d / 1d,229d / 1d,230d / 1d,231d / 1d,232d / 1d,233d / 1d,234d / 1d,235d / 1d,236d / 1d,237d / 1d,238d / 1d,239d / 1d,240d / 1d,241d / 1d,242d / 1d,243d / 1d,244d / 1d,245d / 1d,246d / 1d,247d / 1d,248d / 1d,249d / 1d,250d / 1d,251d / 1d,252d / 1d,253d / 1d,254d / 1d,255d / 1d,256d / 1d,257d / 1d,258d / 1d,259d / 1d,260d / 1d,261d / 1d,262d / 1d,263d / 1d,264d / 1d,265d / 1d,266d / 1d,267d / 1d,268d / 1d,269d / 1d,270d / 1d,271d / 1d,272d / 1d,273d / 1d,274d / 1d,275d / 1d,276d / 1d,277d / 1d,278d / 1d,279d / 1d,280d / 1d,281d / 1d,282d / 1d,283d / 1d,284d / 1d,285d / 1d,286d / 1d,287d / 1d,288d / 1d,289d / 1d,290d / 1d,291d / 1d,292d / 1d,293d / 1d,294d / 1d,295d / 1d,296d / 1d,297d / 1d,298d / 1d,299d / 1d,300d / 1d,301d / 1d,302d / 1d,303d / 1d,304d / 1d,305d / 1d,306d / 1d,307d / 1d,308d / 1d,309d / 1d,310d / 1d,311d / 1d,312d / 1d,313d / 1d,314d / 1d,315d / 1d,316d / 1d,317d / 1d,318d / 1d,319d / 1d,320d / 1d,321d / 1d,322d / 1d,323d / 1d,324d / 1d,325d / 1d,326d / 1d,327d / 1d,328d / 1d,329d / 1d,330d / 1d,331d / 1d,332d / 1d,333d / 1d,334d / 1d,335d / 1d,336d / 1d,337d / 1d,338d / 1d,339d / 1d,340d / 1d,341d / 1d,342d / 1d,343d / 1d,344d / 1d,345d / 1d,346d / 1d,347d / 1d,348d / 1d,349d / 1d,350d / 1d,351d / 1d,352d / 1d,353d / 1d,354d / 1d,355d / 1d,356d / 1d,357d / 1d,358d / 1d,359d / 1d,360d / 1d,361d / 1d,362d / 1d,363d / 1d,364d / 1d,365d / 1d,366d / 1d,367d / 1d,368d / 1d,369d / 1d,370d / 1d,371d / 1d,372d / 1d,373d / 1d,374d / 1d,375d / 1d,376d / 1d,377d / 1d,378d / 1d,379d / 1d,380d / 1d,381d / 1d,382d / 1d,383d / 1d,384d / 1d,385d / 1d,386d / 1d,387d / 1d,388d / 1d,389d / 1d,390d / 1d,391d / 1d,392d / 1d,393d / 1d,394d / 1d,395d / 1d,396d / 1d,397d / 1d,398d / 1d,399d / 1d,400d / 1d,401d / 1d,402d / 1d,403d / 1d,404d / 1d,405d / 1d,406d / 1d,407d / 1d,408d / 1d,409d / 1d,410d / 1d,411d / 1d,412d / 1d,413d / 1d]
; calculate free-bound
freebound, temperature, wavelength, free_bound, /no_setup
Total running time of the script: (1 minutes 3.491 seconds)