Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/102/435 into FORTRAN code for loading all data files into arrays.
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 2026-Aug-18
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+AS/102/435 Transitions of SI in visible and infra-red (Biemont+ 1993)
*================================================================================
*Delta(n) <= 2 allowed transitions of neutral sulphur within the visible
*and infrared spectral ranges
* BIEMONT E., QUINET P., ZEIPPEN C.J.
* <Astron. Astrophys. Suppl. Ser. 102, 435 (1993)>
* =1993A&AS..102..435B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1' ! Configurations selected for the present SST and
HFR calculations in S I
integer*4 nr__
parameter (nr__=26) ! Number of records
character*27 ar__ ! Full-size record
character*3 Mode (nr__) ! Calculation mode (SST or HFR)
character*8 EP (nr__) ! Even parity configuration
character*8 OP (nr__) ! Odd parity configuration
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2' ! HFR fitted parameters (even configurations)
integer*4 nr__1
parameter (nr__1=36) ! Number of records
character*32 ar__1 ! Full-size record
character*10 Conf_ (nr__1) ! Configuration
character*9 Param_ (nr__1) ! Parameter
integer*4 Fit (nr__1) ! (cm-1) ? Fitted value
real*4 Fit_HFR (nr__1) ! ? Fitted-HFR ratio
character*1 n_Fit_HFR (nr__1) ! [ f] When f, fixed value
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3' ! HFR fitted parameters (odd configurations)
integer*4 nr__2
parameter (nr__2=37) ! Number of records
character*32 ar__2 ! Full-size record
character*10 Conf__1 (nr__2) ! Configuration
character*9 Param__1 (nr__2) ! Parameter
integer*4 Fit_1 (nr__2) ! (cm-1) ? Fitted value
real*4 Fit_HFR_1 (nr__2) ! ? Fitted-HFR ratio
character*1 n_Fit_HFR_1(nr__2) ! [ f] When f, fixed value
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4' ! Excitation energies and fine-structure splittings
integer*4 nr__3
parameter (nr__3=133) ! Number of records
character*58 ar__3 ! Full-size record
character*9 Config (nr__3) ! Configuration
character*5 Level (nr__3) ! Level designation
real*8 E_obs (nr__3) ! (cm-1) ? Observed excitation energy
real*4 DeltaE_obs (nr__3) ! (cm-1) ? Observed fine-structure splitting
real*4 DeltaE_cal (nr__3) ! (cm-1) ? Calculed present fine-structure splitting
real*8 TEC (nr__3) ! (cm-1) ? Semi-empirical corrections to calculated
* non-relativistic energies
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5' ! Calculated HFR energy level (even configuration
integer*4 nr__4
parameter (nr__4=33) ! Number of records
character*73 ar__4 ! Full-size record
character*8 Conf (nr__4) ! Configuration
character*9 Level_1 (nr__4) ! Level designation
real*4 E_calc (nr__4) ! (cm-1) HFR energy level
real*4 Delta (nr__4) ! (cm-1) ? E(obs) - E(calc)
character*30 Comp (nr__4) ! (%) LS coupling composition:
* contributions >= 5% only are given
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6' ! Calculated HFR energy level (odd configuration)
integer*4 nr__5
parameter (nr__5=40) ! Number of records
character*73 ar__5 ! Full-size record
character*8 Conf_1 (nr__5) ! Configuration
character*9 Level_2 (nr__5) ! Level designation
real*4 E_calc_1 (nr__5) ! (cm-1) HFR energy level
real*4 Delta_1 (nr__5) ! (cm-1) ? E(obs) - E(calc)
character*30 Comp_1 (nr__5) ! (%) LS coupling composition:
* contributions >= 5% only are given
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7' ! Calculated oscillator strengths
(3p3(4S)nl - 3p3(4S)n'l')
integer*4 nr__6
parameter (nr__6=114) ! Number of records
character*93 ar__6 ! Full-size record
character*15 Transition (nr__6) !
real*8 Lambda (nr__6) ! (0.1nm) Wavelength calculated with the levels
* of Martin et al. (1990)
real*4 HFR (nr__6) ! Calculated oscillator strength (log gf)
* obtained with HFR method
character*1 n_HFR (nr__6) ! [ *] Cancellation factor (CF) *, as defined
* by Cowan (1981) (1)
real*4 SST (nr__6) ! Calculated oscillator strength (log gf)
* obtained with SST method
character*1 n_SST (nr__6) ! [ +] When +, the length and velocity results
* differ by more than a factor of two
real*4 lg1 (nr__6) ! Oscillator strength (log gf) from other work
character*2 r_lg1 (nr__6) ! Reference to result lg1 (2)
real*4 lg2 (nr__6) ! ? Oscillator strength (log gf) from other work
character*4 r_lg2 (nr__6) ! Reference to result lg2 (2)
real*4 lg3 (nr__6) ! ? Oscillator strength (log gf) from other work
character*2 r_lg3 (nr__6) ! Reference to result lg3 (2)
real*4 lg4 (nr__6) ! ? Oscillator strength (log gf) from other work
character*2 r_lg4 (nr__6) ! Reference to result lg4 (2)
real*4 lg5 (nr__6) ! ? Oscillator strength (log gf) from other work
character*1 r_lg5 (nr__6) ! Reference to result lg5 (2)
real*4 lg6 (nr__6) ! ? Oscillator strength (log gf) from other work
character*1 r_lg6 (nr__6) ! Reference to result lg6 (2)
*Note (1):
* A small value indicates strong cancellation effect. The starred transitions
* are those for which CF<=0.02
*Note (2):
* a: Kurucz R.L. & Peytremann E., Smith. Astron. Obs. Spec. Report 362 (1975)
* b: Lambert D.L. & Warner B., MNRAS 138, 181 (1968)
* c: Biemont E. & Grevesse N., ADNDT 12, 217 (1973)
* d: Lambert D.L. & Luck R.E., MNRAS 183, 79 (1978)
* e: Aymar M., NIM 110, 211 (1973)
* f: Aymar M., Physica 66, 364 (1973)
* g: Wiese W. L. et al., (1969)
* h: Foster E.W., Proc. Phys. Soc. London 90, 275 (1967)
* i: Bridges J.M. & Wiese W.L., Phys. Rev. A159, 31 (1967)
* k: Miller M.H., Wilkerson T.D., Roig R.A. & Bengston R.D.,
* Phys. Rev. A9, 2312 (1974)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8' ! Calculated oscillator strengths
(3p3 (2D)nl - 3p3(2D)n'l')
integer*4 nr__7
parameter (nr__7=16) ! Number of records
character*93 ar__7 ! Full-size record
character*15 Transition_1(nr__7) !
real*8 Lambda_1 (nr__7) ! (0.1nm) Wavelength calculated with the levels
* of Martin et al. (1990)
real*4 HFR_1 (nr__7) ! Calculated oscillator strength (log gf)
* obtained with HFR method
character*1 n_HFR_1 (nr__7) ! [ *] Cancellation factor (CF) *, as defined
* by Cowan (1981) (1)
real*4 SST_1 (nr__7) ! Calculated oscillator strength (log gf)
* obtained with SST method
character*1 n_SST_1 (nr__7) ! [ +] When +, the length and velocity results
* differ by more than a factor of two
real*4 lg1_1 (nr__7) ! Oscillator strength (log gf) from other work
character*2 r_lg1_1 (nr__7) ! Reference to result lg1 (2)
real*4 lg2_1 (nr__7) ! ? Oscillator strength (log gf) from other work
character*4 r_lg2_1 (nr__7) ! Reference to result lg2 (2)
real*4 lg3_1 (nr__7) ! ? Oscillator strength (log gf) from other work
character*2 r_lg3_1 (nr__7) ! Reference to result lg3 (2)
real*4 lg4_1 (nr__7) ! ? Oscillator strength (log gf) from other work
character*2 r_lg4_1 (nr__7) ! Reference to result lg4 (2)
real*4 lg5_1 (nr__7) ! ? Oscillator strength (log gf) from other work
character*1 r_lg5_1 (nr__7) ! Reference to result lg5 (2)
real*4 lg6_1 (nr__7) ! ? Oscillator strength (log gf) from other work
character*1 r_lg6_1 (nr__7) ! Reference to result lg6 (2)
*Note (1):
* A small value indicates strong cancellation effect. The starred transitions
* are those for which CF<=0.02
*Note (2):
* a: Kurucz R.L. & Peytremann E., Smith. Astron. Obs. Spec. Report 362 (1975)
* b: Lambert D.L. & Warner B., MNRAS 138, 181 (1968)
* c: Biemont E. & Grevesse N., ADNDT 12, 217 (1973)
* d: Lambert D.L. & Luck R.E., MNRAS 183, 79 (1978)
* e: Aymar M., NIM 110, 211 (1973)
* f: Aymar M., Physica 66, 364 (1973)
* g: Wiese W. L. et al., (1969)
* h: Foster E.W., Proc. Phys. Soc. London 90, 275 (1967)
* i: Bridges J.M. & Wiese W.L., Phys. Rev. A159, 31 (1967)
* k: Miller M.H., Wilkerson T.D., Roig R.A. & Bengston R.D.,
* Phys. Rev. A9, 2312 (1974)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9' ! S I lines and corresponding gf values
integer*4 nr__8
parameter (nr__8=30) ! Number of records
character*63 ar__8 ! Full-size record
real*8 Lambda_lab (nr__8) ! (0.1nm) Laboratory wavelength
character*20 Transition_2(nr__8) !
real*4 log_gf (nr__8) ! HFR oscillator strength
real*4 Wlambda (nr__8) ! (0.1pm) Equivalent width
real*4 Einf (nr__8) ! (eV)
real*4 A (nr__8) ! Abundance in the logarithmic scale
* calculated with the Holweger & Muller
* (1974) solar model
character*1 u_A (nr__8) ! When :, the value was not retained in the
* final sample as explained in the text
C=============================================================================
C Loading file 'table1' ! Configurations selected for the present SST and
* HFR calculations in S I
C Format for file interpretation
1 format(A3,3X,A8,5X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table1')
write(6,*) '....Loading file: table1'
do i__=1,26
read(1,'(A27)')ar__
read(ar__,1)Mode(i__),EP(i__),OP(i__)
c ..............Just test output...........
write(6,1)Mode(i__),EP(i__),OP(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2' ! HFR fitted parameters (even configurations)
C Format for file interpretation
2 format(A10,1X,A9,1X,I5,1X,F4.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2')
write(6,*) '....Loading file: table2'
do i__=1,36
read(1,'(A32)')ar__1
read(ar__1,2)
+ Conf_(i__),Param_(i__),Fit(i__),Fit_HFR(i__),n_Fit_HFR(i__)
if(ar__1(22:26) .EQ. '') Fit(i__) = iNULL__
if(ar__1(28:31) .EQ. '') Fit_HFR(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Conf_(i__),Param_(i__),Fit(i__),Fit_HFR(i__),n_Fit_HFR(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3' ! HFR fitted parameters (odd configurations)
C Format for file interpretation
3 format(A10,1X,A9,1X,I5,1X,F4.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3')
write(6,*) '....Loading file: table3'
do i__=1,37
read(1,'(A32)')ar__2
read(ar__2,3)
+ Conf__1(i__),Param__1(i__),Fit_1(i__),Fit_HFR_1(i__),
+ n_Fit_HFR_1(i__)
if(ar__2(22:26) .EQ. '') Fit_1(i__) = iNULL__
if(ar__2(28:31) .EQ. '') Fit_HFR_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Conf__1(i__),Param__1(i__),Fit_1(i__),Fit_HFR_1(i__),
+ n_Fit_HFR_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4' ! Excitation energies and fine-structure splittings
C Format for file interpretation
4 format(A9,2X,A5,4X,F7.1,4X,F5.1,5X,F6.1,3X,F8.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4')
write(6,*) '....Loading file: table4'
do i__=1,133
read(1,'(A58)')ar__3
read(ar__3,4)
+ Config(i__),Level(i__),E_obs(i__),DeltaE_obs(i__),
+ DeltaE_cal(i__),TEC(i__)
if(ar__3(21:27) .EQ. '') E_obs(i__) = rNULL__
if(ar__3(32:36) .EQ. '') DeltaE_obs(i__) = rNULL__
if(ar__3(42:47) .EQ. '') DeltaE_cal(i__) = rNULL__
if(ar__3(51:58) .EQ. '') TEC(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Config(i__),Level(i__),E_obs(i__),DeltaE_obs(i__),
+ DeltaE_cal(i__),TEC(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5' ! Calculated HFR energy level (even configuration
C Format for file interpretation
5 format(A8,1X,A9,3X,F6.0,4X,F5.0,7X,A30)
C Effective file loading
open(unit=1,status='old',file=
+'table5')
write(6,*) '....Loading file: table5'
do i__=1,33
read(1,'(A73)')ar__4
read(ar__4,5)
+ Conf(i__),Level_1(i__),E_calc(i__),Delta(i__),Comp(i__)
if(ar__4(32:36) .EQ. '') Delta(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Conf(i__),Level_1(i__),E_calc(i__),Delta(i__),Comp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6' ! Calculated HFR energy level (odd configuration)
C Format for file interpretation
6 format(A8,1X,A9,3X,F6.0,4X,F5.0,7X,A30)
C Effective file loading
open(unit=1,status='old',file=
+'table6')
write(6,*) '....Loading file: table6'
do i__=1,40
read(1,'(A73)')ar__5
read(ar__5,6)
+ Conf_1(i__),Level_2(i__),E_calc_1(i__),Delta_1(i__),
+ Comp_1(i__)
if(ar__5(32:36) .EQ. '') Delta_1(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ Conf_1(i__),Level_2(i__),E_calc_1(i__),Delta_1(i__),
+ Comp_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7' ! Calculated oscillator strengths
* (3p3(4S)nl - 3p3(4S)n'l')
C Format for file interpretation
7 format(
+ A15,2X,F9.3,2X,F5.2,A1,1X,F5.2,A1,1X,F5.2,A2,1X,F5.2,A4,1X,
+ F5.2,A2,2X,F5.2,A2,4X,F5.2,A1,1X,F5.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table7')
write(6,*) '....Loading file: table7'
do i__=1,114
read(1,'(A93)')ar__6
read(ar__6,7)
+ Transition(i__),Lambda(i__),HFR(i__),n_HFR(i__),SST(i__),
+ n_SST(i__),lg1(i__),r_lg1(i__),lg2(i__),r_lg2(i__),lg3(i__),
+ r_lg3(i__),lg4(i__),r_lg4(i__),lg5(i__),r_lg5(i__),lg6(i__),
+ r_lg6(i__)
if(ar__6(51:55) .EQ. '') lg2(i__) = rNULL__
if(ar__6(61:65) .EQ. '') lg3(i__) = rNULL__
if(ar__6(70:74) .EQ. '') lg4(i__) = rNULL__
if(ar__6(81:85) .EQ. '') lg5(i__) = rNULL__
if(ar__6(88:92) .EQ. '') lg6(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ Transition(i__),Lambda(i__),HFR(i__),n_HFR(i__),SST(i__),
+ n_SST(i__),lg1(i__),r_lg1(i__),lg2(i__),r_lg2(i__),lg3(i__),
+ r_lg3(i__),lg4(i__),r_lg4(i__),lg5(i__),r_lg5(i__),lg6(i__),
+ r_lg6(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8' ! Calculated oscillator strengths
* (3p3 (2D)nl - 3p3(2D)n'l')
C Format for file interpretation
8 format(
+ A15,2X,F9.3,2X,F5.2,A1,1X,F5.2,A1,1X,F5.2,A2,1X,F5.2,A4,1X,
+ F5.2,A2,2X,F5.2,A2,4X,F5.2,A1,1X,F5.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table8')
write(6,*) '....Loading file: table8'
do i__=1,16
read(1,'(A93)')ar__7
read(ar__7,8)
+ Transition_1(i__),Lambda_1(i__),HFR_1(i__),n_HFR_1(i__),
+ SST_1(i__),n_SST_1(i__),lg1_1(i__),r_lg1_1(i__),lg2_1(i__),
+ r_lg2_1(i__),lg3_1(i__),r_lg3_1(i__),lg4_1(i__),r_lg4_1(i__),
+ lg5_1(i__),r_lg5_1(i__),lg6_1(i__),r_lg6_1(i__)
if(ar__7(51:55) .EQ. '') lg2_1(i__) = rNULL__
if(ar__7(61:65) .EQ. '') lg3_1(i__) = rNULL__
if(ar__7(70:74) .EQ. '') lg4_1(i__) = rNULL__
if(ar__7(81:85) .EQ. '') lg5_1(i__) = rNULL__
if(ar__7(88:92) .EQ. '') lg6_1(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ Transition_1(i__),Lambda_1(i__),HFR_1(i__),n_HFR_1(i__),
+ SST_1(i__),n_SST_1(i__),lg1_1(i__),r_lg1_1(i__),lg2_1(i__),
+ r_lg2_1(i__),lg3_1(i__),r_lg3_1(i__),lg4_1(i__),r_lg4_1(i__),
+ lg5_1(i__),r_lg5_1(i__),lg6_1(i__),r_lg6_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9' ! S I lines and corresponding gf values
C Format for file interpretation
9 format(F9.3,3X,A20,1X,F5.2,3X,F5.1,4X,F4.2,4X,F4.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table9')
write(6,*) '....Loading file: table9'
do i__=1,30
read(1,'(A63)')ar__8
read(ar__8,9)
+ Lambda_lab(i__),Transition_2(i__),log_gf(i__),Wlambda(i__),
+ Einf(i__),A(i__),u_A(i__)
c ..............Just test output...........
write(6,9)
+ Lambda_lab(i__),Transition_2(i__),log_gf(i__),Wlambda(i__),
+ Einf(i__),A(i__),u_A(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end