Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/526/A115 into FORTRAN code for reading data files line by line.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-16
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/A+A/526/A115 Electron-impact excitation of H-like ions (Malespin+, 2011)
*================================================================================
*Electron-impact excitation of H-like Cr, Mn, Fe, Co and Ni for applications in
*modeling X-ray astrophysical sources.
* Malespin C., Ballance C.P., Pindzola M.S., Witthoeft M.C., Kallman T.R.,
* Loch S.D.
* <Astron. Astrophys. 526, A115 (2011)>
* =2011A&A...526A.115M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cr_legy.dat' ! Cr^23+^ configurations, level labels and energies
integer*4 nr__
parameter (nr__=25) ! Number of records
character*51 ar__ ! Full-size record
integer*4 index ! Level index for bound level in H-like ion
character*2 Conf ! Configuration description (1s, 2s etc)
character*1 v2S_1 ! 2S+1 multiplicity
character*1 L ! L value
real*4 J ! J value
real*8 E ! (cm-1) Level energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cr_negy.dat' ! n-resolved labels and energies for Cr^23+^
integer*4 nr__1
parameter (nr__1=5) ! Number of records
character*51 ar__1 ! Full-size record
integer*4 index_1 ! Level index for bound level in H-like ion
integer*4 n ! [1/5] Bundled-n description
real*4 W ! ((statistical weight of n-shell)-1)/2
real*8 E_1 ! (cm-1) Bundled-n energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cr_lcol.dat' ! Cr^23+^ effective collision strengths
integer*4 nr__2
parameter (nr__2=300) ! Number of records
character*107 ar__2 ! Full-size record
integer*4 i ! Level index for upper level in transition
integer*4 j_1 ! Level index for lower lever in transition
real*4 A ! (s-1) A-value for transitions
real*4 ECS7_29E5 ! Effective collision strength at 7.29E+5K
real*4 ECS1_46E6 ! Effective collision strength at 1.46E+6K
real*4 ECS3_46E6 ! Effective collision strength at 3.46E+6K
real*4 ECS7_29E6 ! Effective collision strength at 7.29E+6K
real*4 ECS1_46E7 ! Effective collision strength at 1.46E+7K
real*4 ECS3_46E7 ! Effective collision strength at 3.46E+7K
real*4 ECS7_29E7 ! Effective collision strength at 7.29E+7K
real*4 ECS1_46E8 ! Effective collision strength at 1.46E+8K
real*4 ECS3_64E8 ! Effective collision strength at 3.64E+8K
real*4 Inf ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cr_ncol.dat' ! Cr^23+^ n-resolved effective collision strengths
integer*4 nr__3
parameter (nr__3=10) ! Number of records
character*107 ar__3 ! Full-size record
integer*4 i_1 ! Level index for upper level in transition
integer*4 j_2 ! Level index for lower lever in transition
real*4 A_1 ! (s-1) A-value for transitions
real*4 ECS7_29E5_1 ! Effective collision strength at 7.29E+5K
real*4 ECS1_46E6_1 ! Effective collision strength at 1.46E+6K
real*4 ECS3_46E6_1 ! Effective collision strength at 3.46E+6K
real*4 ECS7_29E6_1 ! Effective collision strength at 7.29E+6K
real*4 ECS1_46E7_1 ! Effective collision strength at 1.46E+7K
real*4 ECS3_46E7_1 ! Effective collision strength at 3.46E+7K
real*4 ECS7_29E7_1 ! Effective collision strength at 7.29E+7K
real*4 ECS1_46E8_1 ! Effective collision strength at 1.46E+8K
real*4 ECS3_64E8_1 ! Effective collision strength at 3.64E+8K
real*4 Inf_1 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mn_legy.dat' ! Mn^24+^ configurations, level labels and energies
integer*4 nr__4
parameter (nr__4=25) ! Number of records
character*51 ar__4 ! Full-size record
integer*4 index_2 ! Level index for bound level in H-like ion
character*2 Conf_1 ! Configuration description (1s, 2s etc)
character*1 v2S_1_1 ! 2S+1 multiplicity
character*1 L_1 ! L value
real*4 J_3 ! J value
real*8 E_2 ! (cm-1) Level energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mn_negy.dat' ! n-resolved labels and energies for Mn^24+^
integer*4 nr__5
parameter (nr__5=5) ! Number of records
character*51 ar__5 ! Full-size record
integer*4 index_3 ! Level index for bound level in H-like ion
integer*4 n_1 ! [1/5] Bundled-n description
real*4 W_1 ! ((statistical weight of n-shell)-1)/2
real*8 E_3 ! (cm-1) Bundled-n energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mn_lcol.dat' ! Mn^24+^ effective collision strengths
integer*4 nr__6
parameter (nr__6=300) ! Number of records
character*125 ar__6 ! Full-size record
integer*4 i_2 ! Level index for upper level in transition
integer*4 j_4 ! Level index for lower lever in transition
real*4 A_2 ! (s-1) A-value for transitions
real*4 ECS1_35E5 ! Effective collision strength at 1.35E+5K
real*4 ECS3_38E5 ! Effective collision strength at 3.38E+5K
real*4 ECS6_67E5 ! Effective collision strength at 6.67E+5K
real*4 ECS1_35E6 ! Effective collision strength at 1.35E+6K
real*4 ECS3_38E6 ! Effective collision strength at 3.38E+6K
real*4 ECS6_67E6 ! Effective collision strength at 6.67E+6K
real*4 ECS1_35E7 ! Effective collision strength at 1.35E+7K
real*4 ECS3_38E7 ! Effective collision strength at 3.38E+7K
real*4 ECS6_67E7 ! Effective collision strength at 6.67E+7K
real*4 ECS1_35E8 ! Effective collision strength at 1.35E+8K
real*4 ECS3_38E8 ! Effective collision strength at 3.38E+8K
real*4 Inf_2 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mn_ncol.dat' ! Mn^24+^ n-resolved effective collision strengths
integer*4 nr__7
parameter (nr__7=10) ! Number of records
character*125 ar__7 ! Full-size record
integer*4 i_3 ! Level index for upper level in transition
integer*4 j_5 ! Level index for lower lever in transition
real*4 A_3 ! (s-1) A-value for transitions
real*4 ECS1_35E5_1 ! Effective collision strength at 1.35E+5K
real*4 ECS3_38E5_1 ! Effective collision strength at 3.38E+5K
real*4 ECS6_67E5_1 ! Effective collision strength at 6.67E+5K
real*4 ECS1_35E6_1 ! Effective collision strength at 1.35E+6K
real*4 ECS3_38E6_1 ! Effective collision strength at 3.38E+6K
real*4 ECS6_67E6_1 ! Effective collision strength at 6.67E+6K
real*4 ECS1_35E7_1 ! Effective collision strength at 1.35E+7K
real*4 ECS3_38E7_1 ! Effective collision strength at 3.38E+7K
real*4 ECS6_67E7_1 ! Effective collision strength at 6.67E+7K
real*4 ECS1_35E8_1 ! Effective collision strength at 1.35E+8K
real*4 ECS3_38E8_1 ! Effective collision strength at 3.38E+8K
real*4 Inf_3 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fe_legy.dat' ! Fe^25+^ configurations, level labels and energies
integer*4 nr__8
parameter (nr__8=25) ! Number of records
character*51 ar__8 ! Full-size record
integer*4 index_4 ! Level index for bound level in H-like ion
character*2 Conf_2 ! Configuration description (1s, 2s etc)
character*1 v2S_1_2 ! 2S+1 multiplicity
character*1 L_2 ! L value
real*4 J_6 ! J value
real*8 E_4 ! (cm-1) Level energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fe_negy.dat' ! n-resolved labels and energies for Fe^25+^
integer*4 nr__9
parameter (nr__9=5) ! Number of records
character*51 ar__9 ! Full-size record
integer*4 index_5 ! Level index for bound level in H-like ion
integer*4 n_2 ! [1/5] Bundled-n description
real*4 W_2 ! ((statistical weight of n-shell)-1)/2
real*8 E_5 ! (cm-1) Bundled-n energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fe_lcol.dat' ! Fe^25+^ effective collision strengths
integer*4 nr__10
parameter (nr__10=300) ! Number of records
character*107 ar__10 ! Full-size record
integer*4 i_4 ! Level index for upper level in transition
integer*4 j_7 ! Level index for lower lever in transition
real*4 A_4 ! (s-1) A-value for transitions
real*4 ECS6_76E5 ! Effective collision strength at 6.76E+5K
real*4 ECS1_35E6_2 ! Effective collision strength at 1.35E+6K
real*4 ECS3_38E6_2 ! Effective collision strength at 3.38E+6K
real*4 ECS6_76E6 ! Effective collision strength at 6.76E+6K
real*4 ECS1_35E7_2 ! Effective collision strength at 1.35E+7K
real*4 ECS3_38E7_2 ! Effective collision strength at 3.38E+7K
real*4 ECS6_76E7 ! Effective collision strength at 6.76E+7K
real*4 ECS1_35E8_2 ! Effective collision strength at 1.35E+8K
real*4 ECS3_38E8_2 ! Effective collision strength at 3.38E+8K
real*4 Inf_4 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fe_ncol.dat' ! Fe^25+^ n-resolved effective collision strengths
integer*4 nr__11
parameter (nr__11=10) ! Number of records
character*107 ar__11 ! Full-size record
integer*4 i_5 ! Level index for upper level in transition
integer*4 j_8 ! Level index for lower lever in transition
real*4 A_5 ! (s-1) A-value for transitions
real*4 ECS6_76E5_1 ! Effective collision strength at 6.76E+5K
real*4 ECS1_35E6_3 ! Effective collision strength at 1.35E+6K
real*4 ECS3_38E6_3 ! Effective collision strength at 3.38E+6K
real*4 ECS6_76E6_1 ! Effective collision strength at 6.76E+6K
real*4 ECS1_35E7_3 ! Effective collision strength at 1.35E+7K
real*4 ECS3_38E7_3 ! Effective collision strength at 3.38E+7K
real*4 ECS6_76E7_1 ! Effective collision strength at 6.76E+7K
real*4 ECS1_35E8_3 ! Effective collision strength at 1.35E+8K
real*4 ECS3_38E8_3 ! Effective collision strength at 3.38E+8K
real*4 Inf_5 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'co_legy.dat' ! Co^26+^ configurations, level labels and energies
integer*4 nr__12
parameter (nr__12=25) ! Number of records
character*51 ar__12 ! Full-size record
integer*4 index_6 ! Level index for bound level in H-like ion
character*2 Conf_3 ! Configuration description (1s, 2s etc)
character*1 v2S_1_3 ! 2S+1 multiplicity
character*1 L_3 ! L value
real*4 J_9 ! J value
real*8 E_6 ! (cm-1) Level energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'co_negy.dat' ! n-resolved labels and energies for Co^26+^
integer*4 nr__13
parameter (nr__13=5) ! Number of records
character*51 ar__13 ! Full-size record
integer*4 index_7 ! Level index for bound level in H-like ion
integer*4 n_3 ! [1/5] Bundled-n description
real*4 W_3 ! ((statistical weight of n-shell)-1)/2
real*8 E_7 ! (cm-1) Bundled-n energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'co_ncol.dat' ! Co^26+^ n-resolved effective collision strengths
integer*4 nr__14
parameter (nr__14=10) ! Number of records
character*107 ar__14 ! Full-size record
integer*4 i_6 ! Level index for upper level in transition
integer*4 j_10 ! Level index for lower lever in transition
real*4 A_6 ! (s-1) A-value for transitions
real*4 ECS7_29E5_2 ! Effective collision strength at 7.29E+5K
real*4 ECS1_46E6_2 ! Effective collision strength at 1.46E+6K
real*4 ECS3_64E6 ! Effective collision strength at 3.64E+6K
real*4 ECS7_29E6_2 ! Effective collision strength at 7.29E+6K
real*4 ECS1_46E7_2 ! Effective collision strength at 1.46E+7K
real*4 ECS3_64E7 ! Effective collision strength at 3.64E+7K
real*4 ECS7_29E7_2 ! Effective collision strength at 7.29E+7K
real*4 ECS1_46E8_2 ! Effective collision strength at 1.46E+8K
real*4 ECS3_64E8_2 ! Effective collision strength at 3.64E+8K
real*4 Inf_6 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'co_lcol.dat' ! Co^26+^ effective collision strengths
integer*4 nr__15
parameter (nr__15=300) ! Number of records
character*107 ar__15 ! Full-size record
integer*4 i_7 ! Level index for upper level in transition
integer*4 j_11 ! Level index for lower lever in transition
real*4 A_7 ! (s-1) A-value for transitions
real*4 ECS7_29E5_3 ! Effective collision strength at 7.29E+5K
real*4 ECS1_46E6_3 ! Effective collision strength at 1.46E+6K
real*4 ECS3_64E6_1 ! Effective collision strength at 3.64E+6K
real*4 ECS7_29E6_3 ! Effective collision strength at 7.29E+6K
real*4 ECS1_46E7_3 ! Effective collision strength at 1.46E+7K
real*4 ECS3_64E7_1 ! Effective collision strength at 3.64E+7K
real*4 ECS7_29E7_3 ! Effective collision strength at 7.29E+7K
real*4 ECS1_46E8_3 ! Effective collision strength at 1.46E+8K
real*4 ECS3_64E8_3 ! Effective collision strength at 3.64E+8K
real*4 Inf_7 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ni_legy.dat' ! Ni^27+^ configurations, level labels and energies
integer*4 nr__16
parameter (nr__16=25) ! Number of records
character*51 ar__16 ! Full-size record
integer*4 index_8 ! Level index for bound level in H-like ion
character*2 Conf_4 ! Configuration description (1s, 2s etc)
character*1 v2S_1_4 ! 2S+1 multiplicity
character*1 L_4 ! L value
real*4 J_12 ! J value
real*8 E_8 ! (cm-1) Level energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ni_negy.dat' ! n-resolved labels and energies for Ni^27+^
integer*4 nr__17
parameter (nr__17=5) ! Number of records
character*51 ar__17 ! Full-size record
integer*4 index_9 ! Level index for bound level in H-like ion
integer*4 n_4 ! [1/5] Bundled-n description
real*4 W_4 ! ((statistical weight of n-shell)-1)/2
real*8 E_9 ! (cm-1) Bundled-n energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ni_lcol.dat' ! Ni^27+^ effective collision strengths
integer*4 nr__18
parameter (nr__18=300) ! Number of records
character*107 ar__18 ! Full-size record
integer*4 i_8 ! Level index for upper level in transition
integer*4 j_13 ! Level index for lower lever in transition
real*4 A_8 ! (s-1) A-value for transitions
real*4 ECS7_29E5_4 ! Effective collision strength at 7.29E+5K
real*4 ECS1_46E6_4 ! Effective collision strength at 1.46E+6K
real*4 ECS3_64E6_2 ! Effective collision strength at 3.64E+6K
real*4 ECS7_29E6_4 ! Effective collision strength at 7.29E+6K
real*4 ECS1_46E7_4 ! Effective collision strength at 1.46E+7K
real*4 ECS3_64E7_2 ! Effective collision strength at 3.64E+7K
real*4 ECS7_29E7_4 ! Effective collision strength at 7.29E+7K
real*4 ECS1_46E8_4 ! Effective collision strength at 1.46E+8K
real*4 ECS3_64E8_4 ! Effective collision strength at 3.64E+8K
real*4 Inf_8 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ni_ncol.dat' ! Ni^27+^ n-resolved effective collision strengths
integer*4 nr__19
parameter (nr__19=10) ! Number of records
character*107 ar__19 ! Full-size record
integer*4 i_9 ! Level index for upper level in transition
integer*4 j_14 ! Level index for lower lever in transition
real*4 A_9 ! (s-1) A-value for transitions
real*4 ECS7_29E5_5 ! Effective collision strength at 7.29E+5K
real*4 ECS1_46E6_5 ! Effective collision strength at 1.46E+6K
real*4 ECS3_64E6_3 ! Effective collision strength at 3.64E+6K
real*4 ECS7_29E6_5 ! Effective collision strength at 7.29E+6K
real*4 ECS1_46E7_5 ! Effective collision strength at 1.46E+7K
real*4 ECS3_64E7_3 ! Effective collision strength at 3.64E+7K
real*4 ECS7_29E7_5 ! Effective collision strength at 7.29E+7K
real*4 ECS1_46E8_5 ! Effective collision strength at 1.46E+8K
real*4 ECS3_64E8_5 ! Effective collision strength at 3.64E+8K
real*4 Inf_9 ! ? Infinite energy limit point
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pec_cr.dat' ! Cr^23+^ photon emissivity coefficients
integer*4 nr__20
parameter (nr__20=26) ! Number of records
character*48 ar__20 ! Full-size record
real*4 Te ! (eV) Electron temperature
real*4 PEC2p3_2_1s ! (ph.cm+3/s) Photon emissivity coefficient for 2p_3/2_-1s
real*4 PEC2p1_2_1s ! (ph.cm+3/s) Photon emissivity coefficient for 2p_1/2_-1s
real*4 PEC3p3_2_1s ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
real*4 PEC3p1_2_1s ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pec_mn.dat' ! Mn^24+^ photon emissivity coefficients
integer*4 nr__21
parameter (nr__21=26) ! Number of records
character*48 ar__21 ! Full-size record
real*4 Te_1 ! (eV) Electron temperature
real*4 PEC2p3_2_1s_1 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_3/2_-1s
real*4 PEC2p1_2_1s_1 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_1/2_-1s
real*4 PEC3p3_2_1s_1 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
real*4 PEC3p1_2_1s_1 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pec_fe.dat' ! Fe^25+^ photon emissivity coefficients
integer*4 nr__22
parameter (nr__22=26) ! Number of records
character*48 ar__22 ! Full-size record
real*4 Te_2 ! (eV) Electron temperature
real*4 PEC2p3_2_1s_2 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_3/2_-1s
real*4 PEC2p1_2_1s_2 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_1/2_-1s
real*4 PEC3p3_2_1s_2 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
real*4 PEC3p1_2_1s_2 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pec_co.dat' ! Co^26+^ photon emissivity coefficients
integer*4 nr__23
parameter (nr__23=26) ! Number of records
character*48 ar__23 ! Full-size record
real*4 Te_3 ! (eV) Electron temperature
real*4 PEC2p3_2_1s_3 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_3/2_-1s
real*4 PEC2p1_2_1s_3 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_1/2_-1s
real*4 PEC3p3_2_1s_3 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
real*4 PEC3p1_2_1s_3 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pec_ni.dat' ! Ni^27+^ photon emissivity coefficients
integer*4 nr__24
parameter (nr__24=26) ! Number of records
character*48 ar__24 ! Full-size record
real*4 Te_4 ! (eV) Electron temperature
real*4 PEC2p3_2_1s_4 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_3/2_-1s
real*4 PEC2p1_2_1s_4 ! (ph.cm+3/s) Photon emissivity coefficient for 2p_1/2_-1s
real*4 PEC3p3_2_1s_4 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
real*4 PEC3p1_2_1s_4 ! (ph.cm+3/s) Photon emissivity coefficient for 3p_3/2_-1s
C=============================================================================
C Loading file 'cr_legy.dat' ! Cr^23+^ configurations, level labels and energies
C Format for file interpretation
1 format(3X,I2,12X,A2,6X,A1,1X,A1,1X,F4.1,8X,F10.0)
C Effective file loading
open(unit=1,status='old',file=
+'cr_legy.dat')
write(6,*) '....Loading file: cr_legy.dat'
do i__=1,25
read(1,'(A51)')ar__
read(ar__,1)index,Conf,v2S_1,L,J,E
c ..............Just test output...........
write(6,1)index,Conf,v2S_1,L,J,E
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cr_negy.dat' ! n-resolved labels and energies for Cr^23+^
C Format for file interpretation
2 format(3X,I2,3X,I1,20X,F4.1,8X,F10.1)
C Effective file loading
open(unit=1,status='old',file=
+'cr_negy.dat')
write(6,*) '....Loading file: cr_negy.dat'
do i__=1,5
read(1,'(A51)')ar__1
read(ar__1,2)index_1,n,W,E_1
c ..............Just test output...........
write(6,2)index_1,n,W,E_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cr_lcol.dat' ! Cr^23+^ effective collision strengths
C Format for file interpretation
3 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'cr_lcol.dat')
write(6,*) '....Loading file: cr_lcol.dat'
do i__=1,300
read(1,'(A107)')ar__2
read(ar__2,3)
+ i,j_1,A,ECS7_29E5,ECS1_46E6,ECS3_46E6,ECS7_29E6,ECS1_46E7,
+ ECS3_46E7,ECS7_29E7,ECS1_46E8,ECS3_64E8,Inf
if(ar__2(100:107) .EQ. '') Inf = rNULL__
c ..............Just test output...........
write(6,3)
+ i,j_1,A,ECS7_29E5,ECS1_46E6,ECS3_46E6,ECS7_29E6,ECS1_46E7,
+ ECS3_46E7,ECS7_29E7,ECS1_46E8,ECS3_64E8,Inf
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cr_ncol.dat' ! Cr^23+^ n-resolved effective collision strengths
C Format for file interpretation
4 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'cr_ncol.dat')
write(6,*) '....Loading file: cr_ncol.dat'
do i__=1,10
read(1,'(A107)')ar__3
read(ar__3,4)
+ i_1,j_2,A_1,ECS7_29E5_1,ECS1_46E6_1,ECS3_46E6_1,ECS7_29E6_1,
+ ECS1_46E7_1,ECS3_46E7_1,ECS7_29E7_1,ECS1_46E8_1,ECS3_64E8_1,
+ Inf_1
if(ar__3(100:107) .EQ. '') Inf_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ i_1,j_2,A_1,ECS7_29E5_1,ECS1_46E6_1,ECS3_46E6_1,ECS7_29E6_1,
+ ECS1_46E7_1,ECS3_46E7_1,ECS7_29E7_1,ECS1_46E8_1,ECS3_64E8_1,
+ Inf_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'mn_legy.dat' ! Mn^24+^ configurations, level labels and energies
C Format for file interpretation
5 format(3X,I2,12X,A2,6X,A1,1X,A1,1X,F4.1,8X,F10.0)
C Effective file loading
open(unit=1,status='old',file=
+'mn_legy.dat')
write(6,*) '....Loading file: mn_legy.dat'
do i__=1,25
read(1,'(A51)')ar__4
read(ar__4,5)index_2,Conf_1,v2S_1_1,L_1,J_3,E_2
c ..............Just test output...........
write(6,5)index_2,Conf_1,v2S_1_1,L_1,J_3,E_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'mn_negy.dat' ! n-resolved labels and energies for Mn^24+^
C Format for file interpretation
6 format(3X,I2,3X,I1,20X,F4.1,8X,F10.1)
C Effective file loading
open(unit=1,status='old',file=
+'mn_negy.dat')
write(6,*) '....Loading file: mn_negy.dat'
do i__=1,5
read(1,'(A51)')ar__5
read(ar__5,6)index_3,n_1,W_1,E_3
c ..............Just test output...........
write(6,6)index_3,n_1,W_1,E_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'mn_lcol.dat' ! Mn^24+^ effective collision strengths
C Format for file interpretation
7 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'mn_lcol.dat')
write(6,*) '....Loading file: mn_lcol.dat'
do i__=1,300
read(1,'(A125)')ar__6
read(ar__6,7)
+ i_2,j_4,A_2,ECS1_35E5,ECS3_38E5,ECS6_67E5,ECS1_35E6,ECS3_38E6,
+ ECS6_67E6,ECS1_35E7,ECS3_38E7,ECS6_67E7,ECS1_35E8,ECS3_38E8,
+ Inf_2
if(ar__6(118:125) .EQ. '') Inf_2 = rNULL__
c ..............Just test output...........
write(6,7)
+ i_2,j_4,A_2,ECS1_35E5,ECS3_38E5,ECS6_67E5,ECS1_35E6,ECS3_38E6,
+ ECS6_67E6,ECS1_35E7,ECS3_38E7,ECS6_67E7,ECS1_35E8,ECS3_38E8,
+ Inf_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'mn_ncol.dat' ! Mn^24+^ n-resolved effective collision strengths
C Format for file interpretation
8 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'mn_ncol.dat')
write(6,*) '....Loading file: mn_ncol.dat'
do i__=1,10
read(1,'(A125)')ar__7
read(ar__7,8)
+ i_3,j_5,A_3,ECS1_35E5_1,ECS3_38E5_1,ECS6_67E5_1,ECS1_35E6_1,
+ ECS3_38E6_1,ECS6_67E6_1,ECS1_35E7_1,ECS3_38E7_1,ECS6_67E7_1,
+ ECS1_35E8_1,ECS3_38E8_1,Inf_3
if(ar__7(118:125) .EQ. '') Inf_3 = rNULL__
c ..............Just test output...........
write(6,8)
+ i_3,j_5,A_3,ECS1_35E5_1,ECS3_38E5_1,ECS6_67E5_1,ECS1_35E6_1,
+ ECS3_38E6_1,ECS6_67E6_1,ECS1_35E7_1,ECS3_38E7_1,ECS6_67E7_1,
+ ECS1_35E8_1,ECS3_38E8_1,Inf_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fe_legy.dat' ! Fe^25+^ configurations, level labels and energies
C Format for file interpretation
9 format(3X,I2,12X,A2,6X,A1,1X,A1,1X,F4.1,8X,F10.0)
C Effective file loading
open(unit=1,status='old',file=
+'fe_legy.dat')
write(6,*) '....Loading file: fe_legy.dat'
do i__=1,25
read(1,'(A51)')ar__8
read(ar__8,9)index_4,Conf_2,v2S_1_2,L_2,J_6,E_4
c ..............Just test output...........
write(6,9)index_4,Conf_2,v2S_1_2,L_2,J_6,E_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fe_negy.dat' ! n-resolved labels and energies for Fe^25+^
C Format for file interpretation
10 format(3X,I2,3X,I1,20X,F4.1,8X,F10.1)
C Effective file loading
open(unit=1,status='old',file=
+'fe_negy.dat')
write(6,*) '....Loading file: fe_negy.dat'
do i__=1,5
read(1,'(A51)')ar__9
read(ar__9,10)index_5,n_2,W_2,E_5
c ..............Just test output...........
write(6,10)index_5,n_2,W_2,E_5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fe_lcol.dat' ! Fe^25+^ effective collision strengths
C Format for file interpretation
11 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'fe_lcol.dat')
write(6,*) '....Loading file: fe_lcol.dat'
do i__=1,300
read(1,'(A107)')ar__10
read(ar__10,11)
+ i_4,j_7,A_4,ECS6_76E5,ECS1_35E6_2,ECS3_38E6_2,ECS6_76E6,
+ ECS1_35E7_2,ECS3_38E7_2,ECS6_76E7,ECS1_35E8_2,ECS3_38E8_2,
+ Inf_4
if(ar__10(100:107) .EQ. '') Inf_4 = rNULL__
c ..............Just test output...........
write(6,11)
+ i_4,j_7,A_4,ECS6_76E5,ECS1_35E6_2,ECS3_38E6_2,ECS6_76E6,
+ ECS1_35E7_2,ECS3_38E7_2,ECS6_76E7,ECS1_35E8_2,ECS3_38E8_2,
+ Inf_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fe_ncol.dat' ! Fe^25+^ n-resolved effective collision strengths
C Format for file interpretation
12 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'fe_ncol.dat')
write(6,*) '....Loading file: fe_ncol.dat'
do i__=1,10
read(1,'(A107)')ar__11
read(ar__11,12)
+ i_5,j_8,A_5,ECS6_76E5_1,ECS1_35E6_3,ECS3_38E6_3,ECS6_76E6_1,
+ ECS1_35E7_3,ECS3_38E7_3,ECS6_76E7_1,ECS1_35E8_3,ECS3_38E8_3,
+ Inf_5
if(ar__11(100:107) .EQ. '') Inf_5 = rNULL__
c ..............Just test output...........
write(6,12)
+ i_5,j_8,A_5,ECS6_76E5_1,ECS1_35E6_3,ECS3_38E6_3,ECS6_76E6_1,
+ ECS1_35E7_3,ECS3_38E7_3,ECS6_76E7_1,ECS1_35E8_3,ECS3_38E8_3,
+ Inf_5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'co_legy.dat' ! Co^26+^ configurations, level labels and energies
C Format for file interpretation
13 format(3X,I2,12X,A2,6X,A1,1X,A1,1X,F4.1,8X,F10.0)
C Effective file loading
open(unit=1,status='old',file=
+'co_legy.dat')
write(6,*) '....Loading file: co_legy.dat'
do i__=1,25
read(1,'(A51)')ar__12
read(ar__12,13)index_6,Conf_3,v2S_1_3,L_3,J_9,E_6
c ..............Just test output...........
write(6,13)index_6,Conf_3,v2S_1_3,L_3,J_9,E_6
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'co_negy.dat' ! n-resolved labels and energies for Co^26+^
C Format for file interpretation
14 format(3X,I2,3X,I1,20X,F4.1,8X,F10.1)
C Effective file loading
open(unit=1,status='old',file=
+'co_negy.dat')
write(6,*) '....Loading file: co_negy.dat'
do i__=1,5
read(1,'(A51)')ar__13
read(ar__13,14)index_7,n_3,W_3,E_7
c ..............Just test output...........
write(6,14)index_7,n_3,W_3,E_7
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'co_ncol.dat' ! Co^26+^ n-resolved effective collision strengths
C Format for file interpretation
15 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'co_ncol.dat')
write(6,*) '....Loading file: co_ncol.dat'
do i__=1,10
read(1,'(A107)')ar__14
read(ar__14,15)
+ i_6,j_10,A_6,ECS7_29E5_2,ECS1_46E6_2,ECS3_64E6,ECS7_29E6_2,
+ ECS1_46E7_2,ECS3_64E7,ECS7_29E7_2,ECS1_46E8_2,ECS3_64E8_2,
+ Inf_6
if(ar__14(100:107) .EQ. '') Inf_6 = rNULL__
c ..............Just test output...........
write(6,15)
+ i_6,j_10,A_6,ECS7_29E5_2,ECS1_46E6_2,ECS3_64E6,ECS7_29E6_2,
+ ECS1_46E7_2,ECS3_64E7,ECS7_29E7_2,ECS1_46E8_2,ECS3_64E8_2,
+ Inf_6
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'co_lcol.dat' ! Co^26+^ effective collision strengths
C Format for file interpretation
16 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'co_lcol.dat')
write(6,*) '....Loading file: co_lcol.dat'
do i__=1,300
read(1,'(A107)')ar__15
read(ar__15,16)
+ i_7,j_11,A_7,ECS7_29E5_3,ECS1_46E6_3,ECS3_64E6_1,ECS7_29E6_3,
+ ECS1_46E7_3,ECS3_64E7_1,ECS7_29E7_3,ECS1_46E8_3,ECS3_64E8_3,
+ Inf_7
if(ar__15(100:107) .EQ. '') Inf_7 = rNULL__
c ..............Just test output...........
write(6,16)
+ i_7,j_11,A_7,ECS7_29E5_3,ECS1_46E6_3,ECS3_64E6_1,ECS7_29E6_3,
+ ECS1_46E7_3,ECS3_64E7_1,ECS7_29E7_3,ECS1_46E8_3,ECS3_64E8_3,
+ Inf_7
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ni_legy.dat' ! Ni^27+^ configurations, level labels and energies
C Format for file interpretation
17 format(3X,I2,12X,A2,6X,A1,1X,A1,1X,F4.1,8X,F10.0)
C Effective file loading
open(unit=1,status='old',file=
+'ni_legy.dat')
write(6,*) '....Loading file: ni_legy.dat'
do i__=1,25
read(1,'(A51)')ar__16
read(ar__16,17)index_8,Conf_4,v2S_1_4,L_4,J_12,E_8
c ..............Just test output...........
write(6,17)index_8,Conf_4,v2S_1_4,L_4,J_12,E_8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ni_negy.dat' ! n-resolved labels and energies for Ni^27+^
C Format for file interpretation
18 format(3X,I2,3X,I1,20X,F4.1,8X,F10.1)
C Effective file loading
open(unit=1,status='old',file=
+'ni_negy.dat')
write(6,*) '....Loading file: ni_negy.dat'
do i__=1,5
read(1,'(A51)')ar__17
read(ar__17,18)index_9,n_4,W_4,E_9
c ..............Just test output...........
write(6,18)index_9,n_4,W_4,E_9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ni_lcol.dat' ! Ni^27+^ effective collision strengths
C Format for file interpretation
19 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'ni_lcol.dat')
write(6,*) '....Loading file: ni_lcol.dat'
do i__=1,300
read(1,'(A107)')ar__18
read(ar__18,19)
+ i_8,j_13,A_8,ECS7_29E5_4,ECS1_46E6_4,ECS3_64E6_2,ECS7_29E6_4,
+ ECS1_46E7_4,ECS3_64E7_2,ECS7_29E7_4,ECS1_46E8_4,ECS3_64E8_4,
+ Inf_8
if(ar__18(100:107) .EQ. '') Inf_8 = rNULL__
c ..............Just test output...........
write(6,19)
+ i_8,j_13,A_8,ECS7_29E5_4,ECS1_46E6_4,ECS3_64E6_2,ECS7_29E6_4,
+ ECS1_46E7_4,ECS3_64E7_2,ECS7_29E7_4,ECS1_46E8_4,ECS3_64E8_4,
+ Inf_8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ni_ncol.dat' ! Ni^27+^ n-resolved effective collision strengths
C Format for file interpretation
20 format(
+ I4,I4,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'ni_ncol.dat')
write(6,*) '....Loading file: ni_ncol.dat'
do i__=1,10
read(1,'(A107)')ar__19
read(ar__19,20)
+ i_9,j_14,A_9,ECS7_29E5_5,ECS1_46E6_5,ECS3_64E6_3,ECS7_29E6_5,
+ ECS1_46E7_5,ECS3_64E7_3,ECS7_29E7_5,ECS1_46E8_5,ECS3_64E8_5,
+ Inf_9
if(ar__19(100:107) .EQ. '') Inf_9 = rNULL__
c ..............Just test output...........
write(6,20)
+ i_9,j_14,A_9,ECS7_29E5_5,ECS1_46E6_5,ECS3_64E6_3,ECS7_29E6_5,
+ ECS1_46E7_5,ECS3_64E7_3,ECS7_29E7_5,ECS1_46E8_5,ECS3_64E8_5,
+ Inf_9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pec_cr.dat' ! Cr^23+^ photon emissivity coefficients
C Format for file interpretation
21 format(E8.3,2X,E8.3,2X,E8.3,2X,E8.3,2X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'pec_cr.dat')
write(6,*) '....Loading file: pec_cr.dat'
do i__=1,26
read(1,'(A48)')ar__20
read(ar__20,21)
+ Te,PEC2p3_2_1s,PEC2p1_2_1s,PEC3p3_2_1s,PEC3p1_2_1s
c ..............Just test output...........
write(6,21)
+ Te,PEC2p3_2_1s,PEC2p1_2_1s,PEC3p3_2_1s,PEC3p1_2_1s
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pec_mn.dat' ! Mn^24+^ photon emissivity coefficients
C Format for file interpretation
22 format(E8.3,2X,E8.3,2X,E8.3,2X,E8.3,2X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'pec_mn.dat')
write(6,*) '....Loading file: pec_mn.dat'
do i__=1,26
read(1,'(A48)')ar__21
read(ar__21,22)
+ Te_1,PEC2p3_2_1s_1,PEC2p1_2_1s_1,PEC3p3_2_1s_1,PEC3p1_2_1s_1
c ..............Just test output...........
write(6,22)
+ Te_1,PEC2p3_2_1s_1,PEC2p1_2_1s_1,PEC3p3_2_1s_1,PEC3p1_2_1s_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pec_fe.dat' ! Fe^25+^ photon emissivity coefficients
C Format for file interpretation
23 format(E8.3,2X,E8.3,2X,E8.3,2X,E8.3,2X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'pec_fe.dat')
write(6,*) '....Loading file: pec_fe.dat'
do i__=1,26
read(1,'(A48)')ar__22
read(ar__22,23)
+ Te_2,PEC2p3_2_1s_2,PEC2p1_2_1s_2,PEC3p3_2_1s_2,PEC3p1_2_1s_2
c ..............Just test output...........
write(6,23)
+ Te_2,PEC2p3_2_1s_2,PEC2p1_2_1s_2,PEC3p3_2_1s_2,PEC3p1_2_1s_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pec_co.dat' ! Co^26+^ photon emissivity coefficients
C Format for file interpretation
24 format(E8.3,2X,E8.3,2X,E8.3,2X,E8.3,2X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'pec_co.dat')
write(6,*) '....Loading file: pec_co.dat'
do i__=1,26
read(1,'(A48)')ar__23
read(ar__23,24)
+ Te_3,PEC2p3_2_1s_3,PEC2p1_2_1s_3,PEC3p3_2_1s_3,PEC3p1_2_1s_3
c ..............Just test output...........
write(6,24)
+ Te_3,PEC2p3_2_1s_3,PEC2p1_2_1s_3,PEC3p3_2_1s_3,PEC3p1_2_1s_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pec_ni.dat' ! Ni^27+^ photon emissivity coefficients
C Format for file interpretation
25 format(E8.3,2X,E8.3,2X,E8.3,2X,E8.3,2X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'pec_ni.dat')
write(6,*) '....Loading file: pec_ni.dat'
do i__=1,26
read(1,'(A48)')ar__24
read(ar__24,25)
+ Te_4,PEC2p3_2_1s_4,PEC2p1_2_1s_4,PEC3p3_2_1s_4,PEC3p1_2_1s_4
c ..............Just test output...........
write(6,25)
+ Te_4,PEC2p3_2_1s_4,PEC2p1_2_1s_4,PEC3p3_2_1s_4,PEC3p1_2_1s_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end