Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/102/435 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 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 ! Calculation mode (SST or HFR)
character*8 EP ! Even parity configuration
character*8 OP ! 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_ ! Configuration
character*9 Param_ ! Parameter
integer*4 Fit ! (cm-1) ? Fitted value
real*4 Fit_HFR ! ? Fitted-HFR ratio
character*1 n_Fit_HFR ! [ 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 ! Configuration
character*9 Param__1 ! Parameter
integer*4 Fit_1 ! (cm-1) ? Fitted value
real*4 Fit_HFR_1 ! ? Fitted-HFR ratio
character*1 n_Fit_HFR_1 ! [ 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 ! Configuration
character*5 Level ! Level designation
real*8 E_obs ! (cm-1) ? Observed excitation energy
real*4 DeltaE_obs ! (cm-1) ? Observed fine-structure splitting
real*4 DeltaE_cal ! (cm-1) ? Calculed present fine-structure splitting
real*8 TEC ! (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 ! Configuration
character*9 Level_1 ! Level designation
real*4 E_calc ! (cm-1) HFR energy level
real*4 Delta ! (cm-1) ? E(obs) - E(calc)
character*30 Comp ! (%) 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 ! Configuration
character*9 Level_2 ! Level designation
real*4 E_calc_1 ! (cm-1) HFR energy level
real*4 Delta_1 ! (cm-1) ? E(obs) - E(calc)
character*30 Comp_1 ! (%) 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 !
real*8 Lambda ! (0.1nm) Wavelength calculated with the levels
* of Martin et al. (1990)
real*4 HFR ! Calculated oscillator strength (log gf)
* obtained with HFR method
character*1 n_HFR ! [ *] Cancellation factor (CF) *, as defined
* by Cowan (1981) (1)
real*4 SST ! Calculated oscillator strength (log gf)
* obtained with SST method
character*1 n_SST ! [ +] When +, the length and velocity results
* differ by more than a factor of two
real*4 lg1 ! Oscillator strength (log gf) from other work
character*2 r_lg1 ! Reference to result lg1 (2)
real*4 lg2 ! ? Oscillator strength (log gf) from other work
character*4 r_lg2 ! Reference to result lg2 (2)
real*4 lg3 ! ? Oscillator strength (log gf) from other work
character*2 r_lg3 ! Reference to result lg3 (2)
real*4 lg4 ! ? Oscillator strength (log gf) from other work
character*2 r_lg4 ! Reference to result lg4 (2)
real*4 lg5 ! ? Oscillator strength (log gf) from other work
character*1 r_lg5 ! Reference to result lg5 (2)
real*4 lg6 ! ? Oscillator strength (log gf) from other work
character*1 r_lg6 ! 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 !
real*8 Lambda_1 ! (0.1nm) Wavelength calculated with the levels
* of Martin et al. (1990)
real*4 HFR_1 ! Calculated oscillator strength (log gf)
* obtained with HFR method
character*1 n_HFR_1 ! [ *] Cancellation factor (CF) *, as defined
* by Cowan (1981) (1)
real*4 SST_1 ! Calculated oscillator strength (log gf)
* obtained with SST method
character*1 n_SST_1 ! [ +] When +, the length and velocity results
* differ by more than a factor of two
real*4 lg1_1 ! Oscillator strength (log gf) from other work
character*2 r_lg1_1 ! Reference to result lg1 (2)
real*4 lg2_1 ! ? Oscillator strength (log gf) from other work
character*4 r_lg2_1 ! Reference to result lg2 (2)
real*4 lg3_1 ! ? Oscillator strength (log gf) from other work
character*2 r_lg3_1 ! Reference to result lg3 (2)
real*4 lg4_1 ! ? Oscillator strength (log gf) from other work
character*2 r_lg4_1 ! Reference to result lg4 (2)
real*4 lg5_1 ! ? Oscillator strength (log gf) from other work
character*1 r_lg5_1 ! Reference to result lg5 (2)
real*4 lg6_1 ! ? Oscillator strength (log gf) from other work
character*1 r_lg6_1 ! 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 ! (0.1nm) Laboratory wavelength
character*20 Transition_2 !
real*4 log_gf ! HFR oscillator strength
real*4 Wlambda ! (0.1pm) Equivalent width
real*4 Einf ! (eV)
real*4 A ! Abundance in the logarithmic scale
* calculated with the Holweger & Muller
* (1974) solar model
character*1 u_A ! 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,EP,OP
c ..............Just test output...........
write(6,1)Mode,EP,OP
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_,Param_,Fit,Fit_HFR,n_Fit_HFR
if(ar__1(22:26) .EQ. '') Fit = iNULL__
if(ar__1(28:31) .EQ. '') Fit_HFR = rNULL__
c ..............Just test output...........
write(6,2)Conf_,Param_,Fit,Fit_HFR,n_Fit_HFR
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,Param__1,Fit_1,Fit_HFR_1,n_Fit_HFR_1
if(ar__2(22:26) .EQ. '') Fit_1 = iNULL__
if(ar__2(28:31) .EQ. '') Fit_HFR_1 = rNULL__
c ..............Just test output...........
write(6,3)Conf__1,Param__1,Fit_1,Fit_HFR_1,n_Fit_HFR_1
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,Level,E_obs,DeltaE_obs,DeltaE_cal,TEC
if(ar__3(21:27) .EQ. '') E_obs = rNULL__
if(ar__3(32:36) .EQ. '') DeltaE_obs = rNULL__
if(ar__3(42:47) .EQ. '') DeltaE_cal = rNULL__
if(ar__3(51:58) .EQ. '') TEC = rNULL__
c ..............Just test output...........
write(6,4)Config,Level,E_obs,DeltaE_obs,DeltaE_cal,TEC
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,Level_1,E_calc,Delta,Comp
if(ar__4(32:36) .EQ. '') Delta = rNULL__
c ..............Just test output...........
write(6,5)Conf,Level_1,E_calc,Delta,Comp
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,Level_2,E_calc_1,Delta_1,Comp_1
if(ar__5(32:36) .EQ. '') Delta_1 = rNULL__
c ..............Just test output...........
write(6,6)Conf_1,Level_2,E_calc_1,Delta_1,Comp_1
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,Lambda,HFR,n_HFR,SST,n_SST,lg1,r_lg1,lg2,r_lg2,lg3,
+ r_lg3,lg4,r_lg4,lg5,r_lg5,lg6,r_lg6
if(ar__6(51:55) .EQ. '') lg2 = rNULL__
if(ar__6(61:65) .EQ. '') lg3 = rNULL__
if(ar__6(70:74) .EQ. '') lg4 = rNULL__
if(ar__6(81:85) .EQ. '') lg5 = rNULL__
if(ar__6(88:92) .EQ. '') lg6 = rNULL__
c ..............Just test output...........
write(6,7)
+ Transition,Lambda,HFR,n_HFR,SST,n_SST,lg1,r_lg1,lg2,r_lg2,lg3,
+ r_lg3,lg4,r_lg4,lg5,r_lg5,lg6,r_lg6
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,Lambda_1,HFR_1,n_HFR_1,SST_1,n_SST_1,lg1_1,
+ r_lg1_1,lg2_1,r_lg2_1,lg3_1,r_lg3_1,lg4_1,r_lg4_1,lg5_1,
+ r_lg5_1,lg6_1,r_lg6_1
if(ar__7(51:55) .EQ. '') lg2_1 = rNULL__
if(ar__7(61:65) .EQ. '') lg3_1 = rNULL__
if(ar__7(70:74) .EQ. '') lg4_1 = rNULL__
if(ar__7(81:85) .EQ. '') lg5_1 = rNULL__
if(ar__7(88:92) .EQ. '') lg6_1 = rNULL__
c ..............Just test output...........
write(6,8)
+ Transition_1,Lambda_1,HFR_1,n_HFR_1,SST_1,n_SST_1,lg1_1,
+ r_lg1_1,lg2_1,r_lg2_1,lg3_1,r_lg3_1,lg4_1,r_lg4_1,lg5_1,
+ r_lg5_1,lg6_1,r_lg6_1
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,Transition_2,log_gf,Wlambda,Einf,A,u_A
c ..............Just test output...........
write(6,9)Lambda_lab,Transition_2,log_gf,Wlambda,Einf,A,u_A
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end