Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/145/167 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-Jun-08
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/AJ/145/167 Abundances of LX Per (Kang+, 2013)
*================================================================================
*Signs of accretion in the abundance patterns of the components of the
*RS CVn-type eclipsing binary star LX Persei.
* Kang Y.-W., Yushchenko A.V., Hong K., Guinan E.F., Gopka V.F.
* <Astron. J., 145, 167 (2013)>
* =2013AJ....145..167K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Iron lines in the spectrum of LX Per
integer*4 nr__
parameter (nr__=107) ! Number of records
character*53 ar__ ! Full-size record
character*5 Ion (nr__) ! Iron line identifier (Fe I or Fe II)
real*8 lambda (nr__) ! (0.1nm) [4515/8946] Wavelength; in {AA}
real*4 loggf (nr__) ! ([-]) Log of the oscillator strength
character*2 r_loggf (nr__) ! Reference for loggf; in refs.dat file
real*4 Elow (nr__) ! (eV) Lower excitation energy level
real*4 EWh (nr__) ! (0.1pm) ? Equivalent width of hot component (m{AA})
real*4 EWc (nr__) ! (0.1pm) ? Equivalent width of cool component (m{AA})
real*4 logNh (nr__) ! ([-]) ? Hot component abundance (logN(H)=12)
real*4 logNc (nr__) ! ([-]) ? Cool component abundance (logN(H)=12)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Lines of chemical elements (except iron) in the
spectrum of the hotter component of LX Per
integer*4 nr__1
parameter (nr__1=110) ! Number of records
character*89 ar__1 ! Full-size record
character*5 Ion_1 (nr__1) ! Elemental line identification
real*8 lambda_1 (nr__1) ! (0.1nm) [4497/9238] Line wavelength; in {AA}
real*4 loggf_1 (nr__1) ! ([-]) Log of the oscillator strength
character*1 f_loggf (nr__1) ! [h] Flag on loggf (h) (1)
character*2 r_loggf_1 (nr__1) ! Reference for loggf; in refs.dat file
real*4 Elow_1 (nr__1) ! (eV) Excitation energy of lower level
real*4 logN (nr__1) ! ([-]) Derived abundance in LX Per (H=12)
real*4 logNo (nr__1) ! ([-]) ? Derived abundance in the Sun (H=12)
real*4 Abun (nr__1) ! ([Sun]) LX Per abundance relative to the Sun
real*4 Bl (nr__1) ! [0/1] Blending in LX Per (1=clean line)
real*4 Blo (nr__1) ! [0/1]? Blending in the Sun (1=clean line)
real*4 D (nr__1) ! [0/1] LX Per depth in synthetic spectrum
real*4 Do (nr__1) ! [0/1]? Solar depth in synthetic spectrum
real*4 e_g (nr__1) ! ([-]) LX Per error in abundance for a change in logg
* of +0.2dex
real*4 e_T (nr__1) ! ([-]) LX Per error in abundance for a change in
* effective temperature of +100K
*Note (1): Isotopic splitting and hyperfine structure of the lines were taken
* into account for Li, Sc, V, Mn, Co, Cu, and Ba.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Lines of chemical elements (except iron) in the
spectrum of the cooler component of LX Per
integer*4 nr__2
parameter (nr__2=123) ! Number of records
character*89 ar__2 ! Full-size record
character*5 Ion_2 (nr__2) ! Elemental line identification
real*8 lambda_2 (nr__2) ! (0.1nm) [4497/9238] Line wavelength; in {AA}
real*4 loggf_2 (nr__2) ! ([-]) Log of the oscillator strength
character*1 f_loggf_1 (nr__2) ! [h] Flag on loggf (h) (1)
character*2 r_loggf_2 (nr__2) ! Reference for loggf; in refs.dat file
real*4 Elow_2 (nr__2) ! (eV) Excitation energy of lower level
real*4 logN_1 (nr__2) ! ([-]) Derived abundance in LX Per (H=12)
real*4 logNo_1 (nr__2) ! ([-]) ? Derived abundance in the Sun (H=12)
real*4 Abun_1 (nr__2) ! ([Sun]) LX Per abundance relative to the Sun
real*4 Bl_1 (nr__2) ! [0/1] Blending in LX Per (1=clean line)
real*4 Blo_1 (nr__2) ! [0/1]? Blending in the Sun (1=clean line)
real*4 D_1 (nr__2) ! [0/1] LX Per depth in synthetic spectrum
real*4 Do_1 (nr__2) ! [0/1]? Solar depth in synthetic spectrum
real*4 e_g_1 (nr__2) ! ([-]) LX Per error in abundance for a change in logg
* of +0.2dex
real*4 e_T_1 (nr__2) ! ([-]) LX Per error in abundance for a change in
* effective temperature of +100K
*Note (1): Isotopic splitting and hyperfine structure of the lines were taken
* into account for Li, Sc, V, Mn, Co, Cu, and Ba.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Mean abundances of chemical elements in the
atmosphere of LX Per
integer*4 nr__3
parameter (nr__3=53) ! Number of records
character*70 ar__3 ! Full-size record
character*1 C (nr__3) ! [ch] Component (h=hotter, c=cooler)
integer*4 Z (nr__3) ! [3/56]? Atomic number
character*5 Ion_3 (nr__3) ! Name of ion
integer*4 NL (nr__3) ! [1/72] Number of spectral lines used
real*4 Abest (nr__3) ! ([Sun]) Best value mean abundance (LX Per minus Solar)
real*4 e_Abest (nr__3) ! ([Sun]) ? Uncertainty in Ab-best
real*4 Ag_0_2 (nr__3) ! ([Sun]) Mean abundance (LX Per minus Solar) with
* surface gravity shifted by +0.2cm/s^2^
real*4 e_Ag_0_2 (nr__3) ! ([Sun]) ? Uncertainty in Ag+0.2
real*4 A_100K (nr__3) ! ([Sun]) Mean abundance (LX Per minus Solar) with
* effective temperature altered by +100K
real*4 e_A_100K (nr__3) ! ([Sun]) ? Uncertainty in A+100K
real*4 Nbest (nr__3) ! ([-]) Best value of LX Per absolute mean abundance
* log(N/H)+12
real*4 e_Nbest (nr__3) ! ([-]) ? Uncertainty in N-best
real*4 Ng_0_2 (nr__3) ! ([-]) Absolute mean abundance of LX Per with surface
* gravity shifted by +0.2cm/s^2^
real*4 e_Ng_0_2 (nr__3) ! ([-]) ? Uncertainty in N+0.2
real*4 N_100K (nr__3) ! ([-]) Absolute mean abundance (log(N/H)+12) of LX Per
* with effective temperature altered by +100K
real*4 e_N_100K (nr__3) ! ([-]) ? Uncertainty in N+100K
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__4
parameter (nr__4=15) ! Number of records
character*118 ar__4 ! Full-size record
integer*4 Ref (nr__4) ! [1/27] Reference number
character*19 BibCode (nr__4) ! Bibliographic code
character*21 Aut (nr__4) ! Authors's name
character*73 Com (nr__4) ! Comments
C=============================================================================
C Loading file 'table1.dat' ! Iron lines in the spectrum of LX Per
C Format for file interpretation
1 format(
+ A5,1X,F8.3,1X,F6.3,1X,A2,1X,F5.3,1X,F4.1,1X,F5.1,1X,F5.3,1X,
+ F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,107
read(1,'(A53)')ar__
read(ar__,1)
+ Ion(i__),lambda(i__),loggf(i__),r_loggf(i__),Elow(i__),
+ EWh(i__),EWc(i__),logNh(i__),logNc(i__)
if(ar__(32:35) .EQ. '') EWh(i__) = rNULL__
if(ar__(37:41) .EQ. '') EWc(i__) = rNULL__
if(ar__(43:47) .EQ. '') logNh(i__) = rNULL__
if(ar__(49:53) .EQ. '') logNc(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Ion(i__),lambda(i__),loggf(i__),r_loggf(i__),Elow(i__),
+ EWh(i__),EWc(i__),logNh(i__),logNc(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Lines of chemical elements (except iron) in the
* spectrum of the hotter component of LX Per
C Format for file interpretation
2 format(
+ A5,1X,F8.3,1X,F6.3,A1,1X,A2,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,110
read(1,'(A89)')ar__1
read(ar__1,2)
+ Ion_1(i__),lambda_1(i__),loggf_1(i__),f_loggf(i__),
+ r_loggf_1(i__),Elow_1(i__),logN(i__),logNo(i__),Abun(i__),
+ Bl(i__),Blo(i__),D(i__),Do(i__),e_g(i__),e_T(i__)
if(ar__1(40:44) .EQ. '') logNo(i__) = rNULL__
if(ar__1(59:63) .EQ. '') Blo(i__) = rNULL__
if(ar__1(71:75) .EQ. '') Do(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Ion_1(i__),lambda_1(i__),loggf_1(i__),f_loggf(i__),
+ r_loggf_1(i__),Elow_1(i__),logN(i__),logNo(i__),Abun(i__),
+ Bl(i__),Blo(i__),D(i__),Do(i__),e_g(i__),e_T(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Lines of chemical elements (except iron) in the
* spectrum of the cooler component of LX Per
C Format for file interpretation
3 format(
+ A5,1X,F8.3,1X,F6.3,A1,1X,A2,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,123
read(1,'(A89)')ar__2
read(ar__2,3)
+ Ion_2(i__),lambda_2(i__),loggf_2(i__),f_loggf_1(i__),
+ r_loggf_2(i__),Elow_2(i__),logN_1(i__),logNo_1(i__),
+ Abun_1(i__),Bl_1(i__),Blo_1(i__),D_1(i__),Do_1(i__),
+ e_g_1(i__),e_T_1(i__)
if(ar__2(40:44) .EQ. '') logNo_1(i__) = rNULL__
if(ar__2(59:63) .EQ. '') Blo_1(i__) = rNULL__
if(ar__2(71:75) .EQ. '') Do_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Ion_2(i__),lambda_2(i__),loggf_2(i__),f_loggf_1(i__),
+ r_loggf_2(i__),Elow_2(i__),logN_1(i__),logNo_1(i__),
+ Abun_1(i__),Bl_1(i__),Blo_1(i__),D_1(i__),Do_1(i__),
+ e_g_1(i__),e_T_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Mean abundances of chemical elements in the
* atmosphere of LX Per
C Format for file interpretation
4 format(
+ A1,1X,I2,1X,A5,1X,I2,1X,F5.2,1X,F3.2,1X,F5.2,1X,F3.2,1X,F5.2,
+ 1X,F3.2,1X,F4.2,1X,F3.2,1X,F4.2,1X,F3.2,1X,F4.2,1X,F3.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,53
read(1,'(A70)')ar__3
read(ar__3,4)
+ C(i__),Z(i__),Ion_3(i__),NL(i__),Abest(i__),e_Abest(i__),
+ Ag_0_2(i__),e_Ag_0_2(i__),A_100K(i__),e_A_100K(i__),
+ Nbest(i__),e_Nbest(i__),Ng_0_2(i__),e_Ng_0_2(i__),N_100K(i__),
+ e_N_100K(i__)
if(ar__3(3:4) .EQ. '') Z(i__) = iNULL__
if(ar__3(21:23) .EQ. '') e_Abest(i__) = rNULL__
if(ar__3(31:33) .EQ. '') e_Ag_0_2(i__) = rNULL__
if(ar__3(41:43) .EQ. '') e_A_100K(i__) = rNULL__
if(ar__3(50:52) .EQ. '') e_Nbest(i__) = rNULL__
if(ar__3(59:61) .EQ. '') e_Ng_0_2(i__) = rNULL__
if(ar__3(68:70) .EQ. '') e_N_100K(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ C(i__),Z(i__),Ion_3(i__),NL(i__),Abest(i__),e_Abest(i__),
+ Ag_0_2(i__),e_Ag_0_2(i__),A_100K(i__),e_A_100K(i__),
+ Nbest(i__),e_Nbest(i__),Ng_0_2(i__),e_Ng_0_2(i__),N_100K(i__),
+ e_N_100K(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
5 format(I2,1X,A19,1X,A21,1X,A73)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,15
read(1,'(A118)')ar__4
read(ar__4,5)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,5)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end