.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "generated/gallery/idl_comparisons/continuum_comparison.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_generated_gallery_idl_comparisons_continuum_comparison.py: CHIANTI IDL Comparison: Continuum ================================== Compare the free-free and free-bound calculations to that in the CHIANTI IDL routines. .. GENERATED FROM PYTHON SOURCE LINES 7-20 .. code-block:: Python 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 from fiasco.util.setup_db import LATEST_VERSION quantity_support() .. rst-class:: sphx-glr-script-out .. code-block:: none .MplQuantityConverter object at 0x7ce9dfee7410> .. GENERATED FROM PYTHON SOURCE LINES 21-24 Define a function for comparing and plotting the outputs from fiasco and IDL. Note that we have precomputed the IDL outputs. .. GENERATED FROM PYTHON SOURCE LINES 24-62 .. code-block:: Python 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() .. GENERATED FROM PYTHON SOURCE LINES 63-66 First, let's compare the outputs for the free-free continuum emission, i.e. that emission produced by thermal bremsstrahlung. .. GENERATED FROM PYTHON SOURCE LINES 66-81 .. code-block:: Python idl_result_freefree = read_idl_test_output('freefree_all_ions', LATEST_VERSION) ion_kwargs = {'abundance': idl_result_freefree['abundance'], 'ionization_fraction': idl_result_freefree['ionization_fraction']} 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']) print(f"CHIANTI database {idl_result_freefree['database_version']}") print(f"CHIANTI IDL {idl_result_freefree['chianti_idl_version']}") print('IDL code to produce free-free result:') print(idl_result_freefree['idl_script']) .. image-sg:: /generated/gallery/idl_comparisons/images/sphx_glr_continuum_comparison_001.png :alt: IDL, fiasco, Normalized Difference :srcset: /generated/gallery/idl_comparisons/images/sphx_glr_continuum_comparison_001.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none CHIANTI database 11.0.2 CHIANTI IDL 11.0.2 IDL code to produce free-free result: ; set common block ; read abundance and ionization equilibrium abundance_subdirs = ['abundance', 'archive'] abund_file = FILEPATH('sun_coronal_1992_feldman_ext.abund', ROOT_DIR=!xuvtop, SUBDIR=abundance_subdirs) ioneq_file = FILEPATH('chianti.ioneq', ROOT_DIR=!xuvtop, SUBDIR='ioneq') ; 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, abund_file=abund_file, ioneq_file=ioneq_file .. GENERATED FROM PYTHON SOURCE LINES 82-84 Next, let's compare the outputs for the free-bound continuum emission. .. GENERATED FROM PYTHON SOURCE LINES 84-99 .. code-block:: Python idl_result_freebound = read_idl_test_output('freebound_all_ions', LATEST_VERSION) ion_kwargs = {'abundance': idl_result_freebound['abundance'], 'ionization_fraction': idl_result_freebound['ionization_fraction']} 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']) print(f"CHIANTI database {idl_result_freebound['database_version']}") print(f"CHIANTI IDL {idl_result_freebound['chianti_idl_version']}") print('IDL code to produce free-bound result:') print(idl_result_freebound['idl_script']) .. image-sg:: /generated/gallery/idl_comparisons/images/sphx_glr_continuum_comparison_002.png :alt: IDL, fiasco, Normalized Difference :srcset: /generated/gallery/idl_comparisons/images/sphx_glr_continuum_comparison_002.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none 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 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 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 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] CHIANTI database 11.0.2 CHIANTI IDL 11.0.2 IDL code to produce free-bound result: ; read abundance and ionization equilibrium abundance_subdirs = ['abundance', 'archive'] abund_file = FILEPATH('sun_coronal_1992_feldman_ext.abund', ROOT_DIR=!xuvtop, SUBDIR=abundance_subdirs) ioneq_file = FILEPATH('chianti.ioneq', ROOT_DIR=!xuvtop, SUBDIR='ioneq') ; 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, abund_file=abund_file, ioneq_file=ioneq_file .. GENERATED FROM PYTHON SOURCE LINES 100-102 Finally, let's compare the outputs for the two-photon continuum emission. .. GENERATED FROM PYTHON SOURCE LINES 102-118 .. code-block:: Python idl_result_twophoton = read_idl_test_output('twophoton_all_ions', LATEST_VERSION) ion_kwargs = {'abundance': idl_result_twophoton['abundance'], 'ionization_fraction': idl_result_twophoton['ionization_fraction']} all_ions = [fiasco.Ion(ion_name, idl_result_twophoton['temperature'], **ion_kwargs) for ion_name in fiasco.list_ions()] all_ions = fiasco.IonCollection(*all_ions) two_photon = all_ions.two_photon(idl_result_twophoton['wavelength'], idl_result_twophoton['density'], use_two_ion_model=False).squeeze() plot_idl_comparison(idl_result_twophoton['wavelength'], idl_result_twophoton['temperature'], two_photon, idl_result_twophoton['two_photon_continuum']) print(f"CHIANTI database {idl_result_twophoton['database_version']}") print(f"CHIANTI IDL {idl_result_twophoton['chianti_idl_version']}") print('IDL code to produce two-photon result:') print(idl_result_twophoton['idl_script']) .. image-sg:: /generated/gallery/idl_comparisons/images/sphx_glr_continuum_comparison_003.png :alt: IDL, fiasco, Normalized Difference :srcset: /generated/gallery/idl_comparisons/images/sphx_glr_continuum_comparison_003.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none WARNING: Li 2 not included in two-photon emission. _elvlc dataset missing for Li 2. [fiasco.collections] WARNING: Li 3 not included in two-photon emission. _elvlc dataset missing for Li 3. [fiasco.collections] WARNING: Be 3 not included in two-photon emission. _elvlc dataset missing for Be 3. [fiasco.collections] WARNING: Be 4 not included in two-photon emission. _elvlc dataset missing for Be 4. [fiasco.collections] WARNING: B 4 not included in two-photon emission. _elvlc dataset missing for B 4. [fiasco.collections] WARNING: B 5 not included in two-photon emission. _elvlc dataset missing for B 5. [fiasco.collections] WARNING: F 8 not included in two-photon emission. _elvlc dataset missing for F 8. [fiasco.collections] WARNING: F 9 not included in two-photon emission. _elvlc dataset missing for F 9. [fiasco.collections] WARNING: Sc 20 not included in two-photon emission. _elvlc dataset missing for Sc 20. [fiasco.collections] WARNING: Sc 21 not included in two-photon emission. _elvlc dataset missing for Sc 21. [fiasco.collections] WARNING: Ti 21 not included in two-photon emission. _elvlc dataset missing for Ti 21. [fiasco.collections] WARNING: Ti 22 not included in two-photon emission. _wgfa dataset missing for Ti 22. [fiasco.collections] WARNING: V 22 not included in two-photon emission. _elvlc dataset missing for V 22. [fiasco.collections] WARNING: V 23 not included in two-photon emission. _elvlc dataset missing for V 23. [fiasco.collections] WARNING: Cr 23 not included in two-photon emission. _elvlc dataset missing for Cr 23. [fiasco.collections] WARNING: Cr 24 not included in two-photon emission. _elvlc dataset missing for Cr 24. [fiasco.collections] WARNING: Mn 24 not included in two-photon emission. _elvlc dataset missing for Mn 24. [fiasco.collections] WARNING: Mn 25 not included in two-photon emission. _elvlc dataset missing for Mn 25. [fiasco.collections] WARNING: Co 26 not included in two-photon emission. _elvlc dataset missing for Co 26. [fiasco.collections] WARNING: Co 27 not included in two-photon emission. _elvlc dataset missing for Co 27. [fiasco.collections] WARNING: Cu 28 not included in two-photon emission. _elvlc dataset missing for Cu 28. [fiasco.collections] WARNING: Cu 29 not included in two-photon emission. _elvlc dataset missing for Cu 29. [fiasco.collections] WARNING: Zn 30 not included in two-photon emission. _elvlc dataset missing for Zn 30. [fiasco.collections] /home/docs/checkouts/readthedocs.org/user_builds/fiasco/checkouts/stable/examples/idl_comparisons/continuum_comparison.py:31: RuntimeWarning: divide by zero encountered in divide difference_norm = difference / result_idl /home/docs/checkouts/readthedocs.org/user_builds/fiasco/checkouts/stable/examples/idl_comparisons/continuum_comparison.py:31: RuntimeWarning: invalid value encountered in divide difference_norm = difference / result_idl CHIANTI database 11.0.2 CHIANTI IDL 11.0.2 IDL code to produce two-photon result: ; read abundance and ionization equilibrium abundance_subdirs = ['abundance', 'archive'] abund_file = FILEPATH('sun_coronal_1992_feldman_ext.abund', ROOT_DIR=!xuvtop, SUBDIR=abundance_subdirs) ioneq_file = FILEPATH('chianti.ioneq', ROOT_DIR=!xuvtop, SUBDIR='ioneq') ; Set ioneq_file this way to get around a bug that always causes the GUI picker to pop up ; even when the file is specified. Seems to only happen in v9 and later. defsysv,'!ioneq_file',ioneq_file ; set temperature and wavelength temperature = [10000d / 1d,3084180842710657d / 274877906944d,3460507821925937d / 274877906944d,7765507274914385d / 549755813888d,4356521234735043d / 274877906944d,1222024305482773d / 68719476736d,1371133822357211d / 68719476736d,384609363000595d / 17179869184d,6904620847743843d / 274877906944d,60524312586927d / 2147483648d,8692402644028503d / 274877906944d,2438259044705789d / 68719476736d,5471543289057141d / 137438953472d,1534793135958543d / 34359738368d,6888264888325867d / 137438953472d,7728760322821873d / 137438953472d,8671811711082939d / 137438953472d,2432483193016735d / 34359738368d,5458582064688281d / 68719476736d,6124629810901963d / 68719476736d,858993459199993d / 8589934592d,7710452106776579d / 68719476736d,2162817388703693d / 17179869184d,4853442046821451d / 34359738368d,5445651543418759d / 34359738368d,763765190926727d / 4294967296d,3427834555892999d / 17179869184d,7692187260011837d / 34359738368d,4315388029839867d / 17179869184d,4841945006954121d / 17179869184d,2716375826258885d / 8589934592d,6095647611764423d / 17179869184d,3419714555660685d / 8589934592d,1918491419948163d / 4294967296d,269072847200227d / 536870912d,4830475201763631d / 8589934592d,677485289928349d / 1073741824d,1520301995635447d / 2147483648d,6823227580860295d / 8589934592d,7655787263627391d / 8589934592d,8589934591999859d / 8589934592d,4819032566735323d / 4294967296d,1351760867939797d / 1073741824d,1516700639631691d / 1073741824d,3403532214636697d / 2147483648d,7637651909267207d / 4294967296d,2142396597433107d / 1073741824d,4807617037507359d / 2147483648d,5394235037299789d / 2147483648d,6052431258692601d / 2147483648d,6790939565647157d / 2147483648d,7619559514705467d / 2147483648d,8549286389151643d / 2147483648d,149882142183449d / 33554432d,336341059000281d / 67108864d,3019047001102245d / 536870912d,6774852899283435d / 1073741824d,1900377494544293d / 268435456d,8529034476075299d / 1073741824d,598108379970885d / 67108864d,1342177279999967d / 134217728d] density = 1000000000d / 1d wavelength = [1d / 1d,2d / 1d,3d / 1d,4d / 1d,5d / 1d,6d / 1d,7d / 1d,8d / 1d,9d / 1d,10d / 1d,11d / 1d,12d / 1d,13d / 1d,14d / 1d,15d / 1d,16d / 1d,17d / 1d,18d / 1d,19d / 1d,20d / 1d,21d / 1d,22d / 1d,23d / 1d,24d / 1d,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,414d / 1d,415d / 1d,416d / 1d,417d / 1d,418d / 1d,419d / 1d,420d / 1d,421d / 1d,422d / 1d,423d / 1d,424d / 1d,425d / 1d,426d / 1d,427d / 1d,428d / 1d,429d / 1d,430d / 1d,431d / 1d,432d / 1d,433d / 1d,434d / 1d,435d / 1d,436d / 1d,437d / 1d,438d / 1d,439d / 1d,440d / 1d,441d / 1d,442d / 1d,443d / 1d,444d / 1d,445d / 1d,446d / 1d,447d / 1d,448d / 1d,449d / 1d,450d / 1d,451d / 1d,452d / 1d,453d / 1d,454d / 1d,455d / 1d,456d / 1d,457d / 1d,458d / 1d,459d / 1d,460d / 1d,461d / 1d,462d / 1d,463d / 1d,464d / 1d,465d / 1d,466d / 1d,467d / 1d,468d / 1d,469d / 1d,470d / 1d,471d / 1d,472d / 1d,473d / 1d,474d / 1d,475d / 1d,476d / 1d,477d / 1d,478d / 1d,479d / 1d,480d / 1d,481d / 1d,482d / 1d,483d / 1d,484d / 1d,485d / 1d,486d / 1d,487d / 1d,488d / 1d,489d / 1d,490d / 1d,491d / 1d,492d / 1d,493d / 1d,494d / 1d,495d / 1d,496d / 1d,497d / 1d,498d / 1d,499d / 1d,500d / 1d,501d / 1d,502d / 1d,503d / 1d,504d / 1d,505d / 1d,506d / 1d,507d / 1d,508d / 1d,509d / 1d,510d / 1d,511d / 1d,512d / 1d,513d / 1d,514d / 1d,515d / 1d,516d / 1d,517d / 1d,518d / 1d,519d / 1d,520d / 1d,521d / 1d,522d / 1d,523d / 1d,524d / 1d,525d / 1d,526d / 1d,527d / 1d,528d / 1d,529d / 1d,530d / 1d,531d / 1d,532d / 1d,533d / 1d,534d / 1d,535d / 1d,536d / 1d,537d / 1d,538d / 1d,539d / 1d,540d / 1d,541d / 1d,542d / 1d,543d / 1d,544d / 1d,545d / 1d,546d / 1d,547d / 1d,548d / 1d,549d / 1d,550d / 1d,551d / 1d,552d / 1d,553d / 1d,554d / 1d,555d / 1d,556d / 1d,557d / 1d,558d / 1d,559d / 1d,560d / 1d,561d / 1d,562d / 1d,563d / 1d,564d / 1d,565d / 1d,566d / 1d,567d / 1d,568d / 1d,569d / 1d,570d / 1d,571d / 1d,572d / 1d,573d / 1d,574d / 1d,575d / 1d,576d / 1d,577d / 1d,578d / 1d,579d / 1d,580d / 1d,581d / 1d,582d / 1d,583d / 1d,584d / 1d,585d / 1d,586d / 1d,587d / 1d,588d / 1d,589d / 1d,590d / 1d,591d / 1d,592d / 1d,593d / 1d,594d / 1d,595d / 1d,596d / 1d,597d / 1d,598d / 1d,599d / 1d,600d / 1d,601d / 1d,602d / 1d,603d / 1d,604d / 1d,605d / 1d,606d / 1d,607d / 1d,608d / 1d,609d / 1d,610d / 1d,611d / 1d,612d / 1d,613d / 1d,614d / 1d,615d / 1d,616d / 1d,617d / 1d,618d / 1d,619d / 1d,620d / 1d,621d / 1d,622d / 1d,623d / 1d,624d / 1d,625d / 1d,626d / 1d,627d / 1d,628d / 1d,629d / 1d,630d / 1d,631d / 1d,632d / 1d,633d / 1d,634d / 1d,635d / 1d,636d / 1d,637d / 1d,638d / 1d,639d / 1d,640d / 1d,641d / 1d,642d / 1d,643d / 1d,644d / 1d,645d / 1d,646d / 1d,647d / 1d,648d / 1d,649d / 1d,650d / 1d,651d / 1d,652d / 1d,653d / 1d,654d / 1d,655d / 1d,656d / 1d,657d / 1d,658d / 1d,659d / 1d,660d / 1d,661d / 1d,662d / 1d,663d / 1d,664d / 1d,665d / 1d,666d / 1d,667d / 1d,668d / 1d,669d / 1d,670d / 1d,671d / 1d,672d / 1d,673d / 1d,674d / 1d,675d / 1d,676d / 1d,677d / 1d,678d / 1d,679d / 1d,680d / 1d,681d / 1d,682d / 1d,683d / 1d,684d / 1d,685d / 1d,686d / 1d,687d / 1d,688d / 1d,689d / 1d,690d / 1d,691d / 1d,692d / 1d,693d / 1d,694d / 1d,695d / 1d,696d / 1d,697d / 1d,698d / 1d,699d / 1d,700d / 1d,701d / 1d,702d / 1d,703d / 1d,704d / 1d,705d / 1d,706d / 1d,707d / 1d,708d / 1d,709d / 1d,710d / 1d,711d / 1d,712d / 1d,713d / 1d,714d / 1d,715d / 1d,716d / 1d,717d / 1d,718d / 1d,719d / 1d,720d / 1d,721d / 1d,722d / 1d,723d / 1d,724d / 1d,725d / 1d,726d / 1d,727d / 1d,728d / 1d,729d / 1d,730d / 1d,731d / 1d,732d / 1d,733d / 1d,734d / 1d,735d / 1d,736d / 1d,737d / 1d,738d / 1d,739d / 1d,740d / 1d,741d / 1d,742d / 1d,743d / 1d,744d / 1d,745d / 1d,746d / 1d,747d / 1d,748d / 1d,749d / 1d,750d / 1d,751d / 1d,752d / 1d,753d / 1d,754d / 1d,755d / 1d,756d / 1d,757d / 1d,758d / 1d,759d / 1d,760d / 1d,761d / 1d,762d / 1d,763d / 1d,764d / 1d,765d / 1d,766d / 1d,767d / 1d,768d / 1d,769d / 1d,770d / 1d,771d / 1d,772d / 1d,773d / 1d,774d / 1d,775d / 1d,776d / 1d,777d / 1d,778d / 1d,779d / 1d,780d / 1d,781d / 1d,782d / 1d,783d / 1d,784d / 1d,785d / 1d,786d / 1d,787d / 1d,788d / 1d,789d / 1d,790d / 1d,791d / 1d,792d / 1d,793d / 1d,794d / 1d,795d / 1d,796d / 1d,797d / 1d,798d / 1d,799d / 1d,800d / 1d,801d / 1d,802d / 1d,803d / 1d,804d / 1d,805d / 1d,806d / 1d,807d / 1d,808d / 1d,809d / 1d,810d / 1d,811d / 1d,812d / 1d,813d / 1d,814d / 1d,815d / 1d,816d / 1d,817d / 1d,818d / 1d,819d / 1d,820d / 1d,821d / 1d,822d / 1d,823d / 1d,824d / 1d,825d / 1d,826d / 1d,827d / 1d,828d / 1d,829d / 1d,830d / 1d,831d / 1d,832d / 1d,833d / 1d,834d / 1d,835d / 1d,836d / 1d,837d / 1d,838d / 1d,839d / 1d,840d / 1d,841d / 1d,842d / 1d,843d / 1d,844d / 1d,845d / 1d,846d / 1d,847d / 1d,848d / 1d,849d / 1d,850d / 1d,851d / 1d,852d / 1d,853d / 1d,854d / 1d,855d / 1d,856d / 1d,857d / 1d,858d / 1d,859d / 1d,860d / 1d,861d / 1d,862d / 1d,863d / 1d,864d / 1d,865d / 1d,866d / 1d,867d / 1d,868d / 1d,869d / 1d,870d / 1d,871d / 1d,872d / 1d,873d / 1d,874d / 1d,875d / 1d,876d / 1d,877d / 1d,878d / 1d,879d / 1d,880d / 1d,881d / 1d,882d / 1d,883d / 1d,884d / 1d,885d / 1d,886d / 1d,887d / 1d,888d / 1d,889d / 1d,890d / 1d,891d / 1d,892d / 1d,893d / 1d,894d / 1d,895d / 1d,896d / 1d,897d / 1d,898d / 1d,899d / 1d,900d / 1d,901d / 1d,902d / 1d,903d / 1d,904d / 1d,905d / 1d,906d / 1d,907d / 1d,908d / 1d,909d / 1d,910d / 1d,911d / 1d,912d / 1d,913d / 1d,914d / 1d,915d / 1d,916d / 1d,917d / 1d,918d / 1d,919d / 1d,920d / 1d,921d / 1d,922d / 1d,923d / 1d,924d / 1d,925d / 1d,926d / 1d,927d / 1d,928d / 1d,929d / 1d,930d / 1d,931d / 1d,932d / 1d,933d / 1d,934d / 1d,935d / 1d,936d / 1d,937d / 1d,938d / 1d,939d / 1d,940d / 1d,941d / 1d,942d / 1d,943d / 1d,944d / 1d,945d / 1d,946d / 1d,947d / 1d,948d / 1d,949d / 1d,950d / 1d,951d / 1d,952d / 1d,953d / 1d,954d / 1d,955d / 1d,956d / 1d,957d / 1d,958d / 1d,959d / 1d,960d / 1d,961d / 1d,962d / 1d,963d / 1d,964d / 1d,965d / 1d,966d / 1d,967d / 1d,968d / 1d,969d / 1d,970d / 1d,971d / 1d,972d / 1d,973d / 1d,974d / 1d,975d / 1d,976d / 1d,977d / 1d,978d / 1d,979d / 1d,980d / 1d,981d / 1d,982d / 1d,983d / 1d,984d / 1d,985d / 1d,986d / 1d,987d / 1d,988d / 1d,989d / 1d,990d / 1d,991d / 1d,992d / 1d,993d / 1d,994d / 1d,995d / 1d,996d / 1d,997d / 1d,998d / 1d,999d / 1d,1000d / 1d,1001d / 1d,1002d / 1d,1003d / 1d,1004d / 1d,1005d / 1d,1006d / 1d,1007d / 1d,1008d / 1d,1009d / 1d,1010d / 1d,1011d / 1d,1012d / 1d,1013d / 1d,1014d / 1d,1015d / 1d,1016d / 1d,1017d / 1d,1018d / 1d,1019d / 1d,1020d / 1d,1021d / 1d,1022d / 1d,1023d / 1d,1024d / 1d,1025d / 1d,1026d / 1d,1027d / 1d,1028d / 1d,1029d / 1d,1030d / 1d,1031d / 1d,1032d / 1d,1033d / 1d,1034d / 1d,1035d / 1d,1036d / 1d,1037d / 1d,1038d / 1d,1039d / 1d,1040d / 1d,1041d / 1d,1042d / 1d,1043d / 1d,1044d / 1d,1045d / 1d,1046d / 1d,1047d / 1d,1048d / 1d,1049d / 1d,1050d / 1d,1051d / 1d,1052d / 1d,1053d / 1d,1054d / 1d,1055d / 1d,1056d / 1d,1057d / 1d,1058d / 1d,1059d / 1d,1060d / 1d,1061d / 1d,1062d / 1d,1063d / 1d,1064d / 1d,1065d / 1d,1066d / 1d,1067d / 1d,1068d / 1d,1069d / 1d,1070d / 1d,1071d / 1d,1072d / 1d,1073d / 1d,1074d / 1d,1075d / 1d,1076d / 1d,1077d / 1d,1078d / 1d,1079d / 1d,1080d / 1d,1081d / 1d,1082d / 1d,1083d / 1d,1084d / 1d,1085d / 1d,1086d / 1d,1087d / 1d,1088d / 1d,1089d / 1d,1090d / 1d,1091d / 1d,1092d / 1d,1093d / 1d,1094d / 1d,1095d / 1d,1096d / 1d,1097d / 1d,1098d / 1d,1099d / 1d,1100d / 1d,1101d / 1d,1102d / 1d,1103d / 1d,1104d / 1d,1105d / 1d,1106d / 1d,1107d / 1d,1108d / 1d,1109d / 1d,1110d / 1d,1111d / 1d,1112d / 1d,1113d / 1d,1114d / 1d,1115d / 1d,1116d / 1d,1117d / 1d,1118d / 1d,1119d / 1d,1120d / 1d,1121d / 1d,1122d / 1d,1123d / 1d,1124d / 1d,1125d / 1d,1126d / 1d,1127d / 1d,1128d / 1d,1129d / 1d,1130d / 1d,1131d / 1d,1132d / 1d,1133d / 1d,1134d / 1d,1135d / 1d,1136d / 1d,1137d / 1d,1138d / 1d,1139d / 1d,1140d / 1d,1141d / 1d,1142d / 1d,1143d / 1d,1144d / 1d,1145d / 1d,1146d / 1d,1147d / 1d,1148d / 1d,1149d / 1d,1150d / 1d,1151d / 1d,1152d / 1d,1153d / 1d,1154d / 1d,1155d / 1d,1156d / 1d,1157d / 1d,1158d / 1d,1159d / 1d,1160d / 1d,1161d / 1d,1162d / 1d,1163d / 1d,1164d / 1d,1165d / 1d,1166d / 1d,1167d / 1d,1168d / 1d,1169d / 1d,1170d / 1d,1171d / 1d,1172d / 1d,1173d / 1d,1174d / 1d,1175d / 1d,1176d / 1d,1177d / 1d,1178d / 1d,1179d / 1d,1180d / 1d,1181d / 1d,1182d / 1d,1183d / 1d,1184d / 1d,1185d / 1d,1186d / 1d,1187d / 1d,1188d / 1d,1189d / 1d,1190d / 1d,1191d / 1d,1192d / 1d,1193d / 1d,1194d / 1d,1195d / 1d,1196d / 1d,1197d / 1d,1198d / 1d,1199d / 1d,1200d / 1d,1201d / 1d,1202d / 1d,1203d / 1d,1204d / 1d,1205d / 1d,1206d / 1d,1207d / 1d,1208d / 1d,1209d / 1d,1210d / 1d,1211d / 1d,1212d / 1d,1213d / 1d,1214d / 1d,1215d / 1d,1216d / 1d,1217d / 1d,1218d / 1d,1219d / 1d,1220d / 1d,1221d / 1d,1222d / 1d,1223d / 1d,1224d / 1d,1225d / 1d,1226d / 1d,1227d / 1d,1228d / 1d,1229d / 1d,1230d / 1d,1231d / 1d,1232d / 1d,1233d / 1d,1234d / 1d,1235d / 1d,1236d / 1d,1237d / 1d,1238d / 1d,1239d / 1d,1240d / 1d,1241d / 1d,1242d / 1d,1243d / 1d,1244d / 1d,1245d / 1d,1246d / 1d,1247d / 1d,1248d / 1d,1249d / 1d,1250d / 1d,1251d / 1d,1252d / 1d,1253d / 1d,1254d / 1d,1255d / 1d,1256d / 1d,1257d / 1d,1258d / 1d,1259d / 1d,1260d / 1d,1261d / 1d,1262d / 1d,1263d / 1d,1264d / 1d,1265d / 1d,1266d / 1d,1267d / 1d,1268d / 1d,1269d / 1d,1270d / 1d,1271d / 1d,1272d / 1d,1273d / 1d,1274d / 1d,1275d / 1d,1276d / 1d,1277d / 1d,1278d / 1d,1279d / 1d,1280d / 1d,1281d / 1d,1282d / 1d,1283d / 1d,1284d / 1d,1285d / 1d,1286d / 1d,1287d / 1d,1288d / 1d,1289d / 1d,1290d / 1d,1291d / 1d,1292d / 1d,1293d / 1d,1294d / 1d,1295d / 1d,1296d / 1d,1297d / 1d,1298d / 1d,1299d / 1d,1300d / 1d,1301d / 1d,1302d / 1d,1303d / 1d,1304d / 1d,1305d / 1d,1306d / 1d,1307d / 1d,1308d / 1d,1309d / 1d,1310d / 1d,1311d / 1d,1312d / 1d,1313d / 1d,1314d / 1d,1315d / 1d,1316d / 1d,1317d / 1d,1318d / 1d,1319d / 1d,1320d / 1d,1321d / 1d,1322d / 1d,1323d / 1d,1324d / 1d,1325d / 1d,1326d / 1d,1327d / 1d,1328d / 1d,1329d / 1d,1330d / 1d,1331d / 1d,1332d / 1d,1333d / 1d,1334d / 1d,1335d / 1d,1336d / 1d,1337d / 1d,1338d / 1d,1339d / 1d,1340d / 1d,1341d / 1d,1342d / 1d,1343d / 1d,1344d / 1d,1345d / 1d,1346d / 1d,1347d / 1d,1348d / 1d,1349d / 1d,1350d / 1d,1351d / 1d,1352d / 1d,1353d / 1d,1354d / 1d,1355d / 1d,1356d / 1d,1357d / 1d,1358d / 1d,1359d / 1d,1360d / 1d,1361d / 1d,1362d / 1d,1363d / 1d,1364d / 1d,1365d / 1d,1366d / 1d,1367d / 1d,1368d / 1d,1369d / 1d,1370d / 1d,1371d / 1d,1372d / 1d,1373d / 1d,1374d / 1d,1375d / 1d,1376d / 1d,1377d / 1d,1378d / 1d,1379d / 1d,1380d / 1d,1381d / 1d,1382d / 1d,1383d / 1d,1384d / 1d,1385d / 1d,1386d / 1d,1387d / 1d,1388d / 1d,1389d / 1d,1390d / 1d,1391d / 1d,1392d / 1d,1393d / 1d,1394d / 1d,1395d / 1d,1396d / 1d,1397d / 1d,1398d / 1d,1399d / 1d,1400d / 1d,1401d / 1d,1402d / 1d,1403d / 1d,1404d / 1d,1405d / 1d,1406d / 1d,1407d / 1d,1408d / 1d,1409d / 1d,1410d / 1d,1411d / 1d,1412d / 1d,1413d / 1d,1414d / 1d,1415d / 1d,1416d / 1d,1417d / 1d,1418d / 1d,1419d / 1d,1420d / 1d,1421d / 1d,1422d / 1d,1423d / 1d,1424d / 1d,1425d / 1d,1426d / 1d,1427d / 1d,1428d / 1d,1429d / 1d,1430d / 1d,1431d / 1d,1432d / 1d,1433d / 1d,1434d / 1d,1435d / 1d,1436d / 1d,1437d / 1d,1438d / 1d,1439d / 1d,1440d / 1d,1441d / 1d,1442d / 1d,1443d / 1d,1444d / 1d,1445d / 1d,1446d / 1d,1447d / 1d,1448d / 1d,1449d / 1d,1450d / 1d,1451d / 1d,1452d / 1d,1453d / 1d,1454d / 1d,1455d / 1d,1456d / 1d,1457d / 1d,1458d / 1d,1459d / 1d,1460d / 1d,1461d / 1d,1462d / 1d,1463d / 1d,1464d / 1d,1465d / 1d,1466d / 1d,1467d / 1d,1468d / 1d,1469d / 1d,1470d / 1d,1471d / 1d,1472d / 1d,1473d / 1d,1474d / 1d,1475d / 1d,1476d / 1d,1477d / 1d,1478d / 1d,1479d / 1d,1480d / 1d,1481d / 1d,1482d / 1d,1483d / 1d,1484d / 1d,1485d / 1d,1486d / 1d,1487d / 1d,1488d / 1d,1489d / 1d,1490d / 1d,1491d / 1d,1492d / 1d,1493d / 1d,1494d / 1d,1495d / 1d,1496d / 1d,1497d / 1d,1498d / 1d,1499d / 1d,1500d / 1d,1501d / 1d,1502d / 1d,1503d / 1d,1504d / 1d,1505d / 1d,1506d / 1d,1507d / 1d,1508d / 1d,1509d / 1d,1510d / 1d,1511d / 1d,1512d / 1d,1513d / 1d,1514d / 1d,1515d / 1d,1516d / 1d,1517d / 1d,1518d / 1d,1519d / 1d,1520d / 1d,1521d / 1d,1522d / 1d,1523d / 1d,1524d / 1d,1525d / 1d,1526d / 1d,1527d / 1d,1528d / 1d,1529d / 1d,1530d / 1d,1531d / 1d,1532d / 1d,1533d / 1d,1534d / 1d,1535d / 1d,1536d / 1d,1537d / 1d,1538d / 1d,1539d / 1d,1540d / 1d,1541d / 1d,1542d / 1d,1543d / 1d,1544d / 1d,1545d / 1d,1546d / 1d,1547d / 1d,1548d / 1d,1549d / 1d,1550d / 1d,1551d / 1d,1552d / 1d,1553d / 1d,1554d / 1d,1555d / 1d,1556d / 1d,1557d / 1d,1558d / 1d,1559d / 1d,1560d / 1d,1561d / 1d,1562d / 1d,1563d / 1d,1564d / 1d,1565d / 1d,1566d / 1d,1567d / 1d,1568d / 1d,1569d / 1d,1570d / 1d,1571d / 1d,1572d / 1d,1573d / 1d,1574d / 1d,1575d / 1d,1576d / 1d,1577d / 1d,1578d / 1d,1579d / 1d,1580d / 1d,1581d / 1d,1582d / 1d,1583d / 1d,1584d / 1d,1585d / 1d,1586d / 1d,1587d / 1d,1588d / 1d,1589d / 1d,1590d / 1d,1591d / 1d,1592d / 1d,1593d / 1d,1594d / 1d,1595d / 1d,1596d / 1d,1597d / 1d,1598d / 1d,1599d / 1d,1600d / 1d,1601d / 1d,1602d / 1d,1603d / 1d,1604d / 1d,1605d / 1d,1606d / 1d,1607d / 1d,1608d / 1d,1609d / 1d,1610d / 1d,1611d / 1d,1612d / 1d,1613d / 1d,1614d / 1d,1615d / 1d,1616d / 1d,1617d / 1d,1618d / 1d,1619d / 1d,1620d / 1d,1621d / 1d,1622d / 1d,1623d / 1d,1624d / 1d,1625d / 1d,1626d / 1d,1627d / 1d,1628d / 1d,1629d / 1d,1630d / 1d,1631d / 1d,1632d / 1d,1633d / 1d,1634d / 1d,1635d / 1d,1636d / 1d,1637d / 1d,1638d / 1d,1639d / 1d,1640d / 1d,1641d / 1d,1642d / 1d,1643d / 1d,1644d / 1d,1645d / 1d,1646d / 1d,1647d / 1d,1648d / 1d,1649d / 1d,1650d / 1d,1651d / 1d,1652d / 1d,1653d / 1d,1654d / 1d,1655d / 1d,1656d / 1d,1657d / 1d,1658d / 1d,1659d / 1d,1660d / 1d,1661d / 1d,1662d / 1d,1663d / 1d,1664d / 1d,1665d / 1d,1666d / 1d,1667d / 1d,1668d / 1d,1669d / 1d,1670d / 1d,1671d / 1d,1672d / 1d,1673d / 1d,1674d / 1d,1675d / 1d,1676d / 1d,1677d / 1d,1678d / 1d,1679d / 1d,1680d / 1d,1681d / 1d,1682d / 1d,1683d / 1d,1684d / 1d,1685d / 1d,1686d / 1d,1687d / 1d,1688d / 1d,1689d / 1d,1690d / 1d,1691d / 1d,1692d / 1d,1693d / 1d,1694d / 1d,1695d / 1d,1696d / 1d,1697d / 1d,1698d / 1d,1699d / 1d,1700d / 1d,1701d / 1d,1702d / 1d,1703d / 1d,1704d / 1d,1705d / 1d,1706d / 1d,1707d / 1d,1708d / 1d,1709d / 1d,1710d / 1d,1711d / 1d,1712d / 1d,1713d / 1d,1714d / 1d,1715d / 1d,1716d / 1d,1717d / 1d,1718d / 1d,1719d / 1d,1720d / 1d,1721d / 1d,1722d / 1d,1723d / 1d,1724d / 1d,1725d / 1d,1726d / 1d,1727d / 1d,1728d / 1d,1729d / 1d,1730d / 1d,1731d / 1d,1732d / 1d,1733d / 1d,1734d / 1d,1735d / 1d,1736d / 1d,1737d / 1d,1738d / 1d,1739d / 1d,1740d / 1d,1741d / 1d,1742d / 1d,1743d / 1d,1744d / 1d,1745d / 1d,1746d / 1d,1747d / 1d,1748d / 1d,1749d / 1d,1750d / 1d,1751d / 1d,1752d / 1d,1753d / 1d,1754d / 1d,1755d / 1d,1756d / 1d,1757d / 1d,1758d / 1d,1759d / 1d,1760d / 1d,1761d / 1d,1762d / 1d,1763d / 1d,1764d / 1d,1765d / 1d,1766d / 1d,1767d / 1d,1768d / 1d,1769d / 1d,1770d / 1d,1771d / 1d,1772d / 1d,1773d / 1d,1774d / 1d,1775d / 1d,1776d / 1d,1777d / 1d,1778d / 1d,1779d / 1d,1780d / 1d,1781d / 1d,1782d / 1d,1783d / 1d,1784d / 1d,1785d / 1d,1786d / 1d,1787d / 1d,1788d / 1d,1789d / 1d,1790d / 1d,1791d / 1d,1792d / 1d,1793d / 1d,1794d / 1d,1795d / 1d,1796d / 1d,1797d / 1d,1798d / 1d,1799d / 1d,1800d / 1d,1801d / 1d,1802d / 1d,1803d / 1d,1804d / 1d,1805d / 1d,1806d / 1d,1807d / 1d,1808d / 1d,1809d / 1d,1810d / 1d,1811d / 1d,1812d / 1d,1813d / 1d,1814d / 1d,1815d / 1d,1816d / 1d,1817d / 1d,1818d / 1d,1819d / 1d,1820d / 1d,1821d / 1d,1822d / 1d,1823d / 1d,1824d / 1d,1825d / 1d,1826d / 1d,1827d / 1d,1828d / 1d,1829d / 1d,1830d / 1d,1831d / 1d,1832d / 1d,1833d / 1d,1834d / 1d,1835d / 1d,1836d / 1d,1837d / 1d,1838d / 1d,1839d / 1d,1840d / 1d,1841d / 1d,1842d / 1d,1843d / 1d,1844d / 1d,1845d / 1d,1846d / 1d,1847d / 1d,1848d / 1d,1849d / 1d,1850d / 1d,1851d / 1d,1852d / 1d,1853d / 1d,1854d / 1d,1855d / 1d,1856d / 1d,1857d / 1d,1858d / 1d,1859d / 1d,1860d / 1d,1861d / 1d,1862d / 1d,1863d / 1d,1864d / 1d,1865d / 1d,1866d / 1d,1867d / 1d,1868d / 1d,1869d / 1d,1870d / 1d,1871d / 1d,1872d / 1d,1873d / 1d,1874d / 1d,1875d / 1d,1876d / 1d,1877d / 1d,1878d / 1d,1879d / 1d,1880d / 1d,1881d / 1d,1882d / 1d,1883d / 1d,1884d / 1d,1885d / 1d,1886d / 1d,1887d / 1d,1888d / 1d,1889d / 1d,1890d / 1d,1891d / 1d,1892d / 1d,1893d / 1d,1894d / 1d,1895d / 1d,1896d / 1d,1897d / 1d,1898d / 1d,1899d / 1d,1900d / 1d,1901d / 1d,1902d / 1d,1903d / 1d,1904d / 1d,1905d / 1d,1906d / 1d,1907d / 1d,1908d / 1d,1909d / 1d,1910d / 1d,1911d / 1d,1912d / 1d,1913d / 1d,1914d / 1d,1915d / 1d,1916d / 1d,1917d / 1d,1918d / 1d,1919d / 1d,1920d / 1d,1921d / 1d,1922d / 1d,1923d / 1d,1924d / 1d,1925d / 1d,1926d / 1d,1927d / 1d,1928d / 1d,1929d / 1d,1930d / 1d,1931d / 1d,1932d / 1d,1933d / 1d,1934d / 1d,1935d / 1d,1936d / 1d,1937d / 1d,1938d / 1d,1939d / 1d,1940d / 1d,1941d / 1d,1942d / 1d,1943d / 1d,1944d / 1d,1945d / 1d,1946d / 1d,1947d / 1d,1948d / 1d,1949d / 1d,1950d / 1d,1951d / 1d,1952d / 1d,1953d / 1d,1954d / 1d,1955d / 1d,1956d / 1d,1957d / 1d,1958d / 1d,1959d / 1d,1960d / 1d,1961d / 1d,1962d / 1d,1963d / 1d,1964d / 1d,1965d / 1d,1966d / 1d,1967d / 1d,1968d / 1d,1969d / 1d,1970d / 1d,1971d / 1d,1972d / 1d,1973d / 1d,1974d / 1d,1975d / 1d,1976d / 1d,1977d / 1d,1978d / 1d,1979d / 1d,1980d / 1d,1981d / 1d,1982d / 1d,1983d / 1d,1984d / 1d,1985d / 1d,1986d / 1d,1987d / 1d,1988d / 1d,1989d / 1d,1990d / 1d,1991d / 1d,1992d / 1d,1993d / 1d,1994d / 1d,1995d / 1d,1996d / 1d,1997d / 1d,1998d / 1d,1999d / 1d] ; calculate two-photon two_photon,temperature,wavelength,two_photon_continuum,$ edensity=density,abund_file=abund_file,ioneq_file=ioneq_file defsysv,'!ioneq_file','' .. rst-class:: sphx-glr-timing **Total running time of the script:** (1 minutes 59.189 seconds) .. _sphx_glr_download_generated_gallery_idl_comparisons_continuum_comparison.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: continuum_comparison.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: continuum_comparison.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: continuum_comparison.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_