Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-19
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/MNRAS/378/179 Eclipsing Binary with Eccentric Orbits Catalog (Bulut+, 2007)
*================================================================================
*A New Catalog of Eclipsing Binary Stars with Eccentric Orbits.
* Bulut I., Demircan O.
* <Mon. Not. R. Astron. Soc. 378, 179 (2007)>
* =2007MNRAS.378..179B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Identifications of eclipsing binaries
integer*4 nr__
parameter (nr__=124) ! Number of records
character*103 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension J2000
real*8 Dec ! (deg) Declination J2000
integer*4 Seq ! Serial number of entry
character*13 Name ! Name of star
integer*4 HD ! ? Henry Draper Catalog number
integer*4 HIP ! ? Hipparcos Catalog number
integer*4 RAh ! (h) Right Ascension (J2000)
integer*4 RAm ! (min) Right Ascension (J2000)
real*4 RAs ! (s) Right Ascension (J2000)
character*1 DE_ ! Declination sign (J2000)
integer*4 DEd ! (deg) Declination (J2000)
integer*4 DEm ! (arcmin) Declination (J2000)
real*4 DEs ! (arcsec) Declination (J2000)
real*4 mag ! (mag) ? Magnitude at maximum brightness
real*4 Dmag ! (mag) ? Magnitude changes brightness
character*2 Filt ! Photometric system for magnitudes
character*16 SpType ! MK spectral type
character*10 Type ! Type of the variability (1)
character*5 Ref1 ! Reference, in refs.dat file
*Note (1): In this column contains various abbreviations from GCVS
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Apsidal Motion Parameters
integer*4 nr__1
parameter (nr__1=124) ! Number of records
character*171 ar__1 ! Full-size record
integer*4 Seq_1 ! Serial number of entry
character*13 Name_1 ! Name of star
real*8 T0 ! (d) ? HJD(2400000+) Epoch of primary minimum
character*1 n_T0 ! []] T0 indicates periastron epoch (1)
real*8 e_T0 ! (d) ? rms uncertainty on T0
real*8 Ps ! (d) ? Sidereal period
real*8 e_Ps ! (d) ? rms uncertainty on Ps
real*8 Pa ! (d) ? Anomalistic period (periastre passages)
real*8 e_Pa ! (d) ? rms uncertainty on Pa
character*1 l_e ! [~] Limit flag on e
real*8 e ! ? Orbital eccentricity
real*8 e_e ! ? rms uncertainty on e
real*4 omega ! (deg) ? Position of periastron in epoch of reference
real*4 e_omega ! (deg) ? rms uncertainty on omega
character*6 n_U ! Note on U (2)
character*1 l_U ! [=>] Limit flag on U
real*8 U ! (yr) ? Apsidal motion period
real*8 e_U ! (yr) ? rms uncertainty on U
real*4 logk2 ! ? Mean internal structure constant of system
real*4 e_logk2 ! ? rms uncertainty on logk2
character*10 Ref2 ! Reference numbers for above parameters,
* in refs.dat file
*Note (1): The values in brackets (]) are the time of periastron passage
*Note (2): If a value of the apsidal motion period is not known, the phase
* of secondary minimum (phi_II) is given in U column.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Photometric Parameters
integer*4 nr__2
parameter (nr__2=238) ! Number of records
character*132 ar__2 ! Full-size record
integer*4 Seq_2 ! Serial number of entry
character*13 Name_2 ! Name of star
real*4 L1 ! ? Fractional luminosity of primary component
real*4 e_L1 ! ? rms uncertainty on L1
character*2 n_L1 ! Filters for L1 and L2 (1)
real*4 L2 ! ? Fractional luminosity of secondary component
real*4 e_L2 ! ? rms uncertainty on L2
real*4 r1 ! ? Relative mean radius of primary component
real*4 e_r1 ! ? rms uncertainty on r1
real*4 r2 ! ? Relative mean radius of secondary component
real*4 e_r2 ! ? rms uncertainty on r2
real*4 i ! (deg) ? Orbital inclination
real*4 e_i ! (deg) ? rms uncertainty on i
real*4 e_1 ! ? Orbital eccentricity obtained by the light
* curve solution.
real*4 e_e_1 ! ? rms uncertainty on e
real*4 q ! ? Mass ratio
real*4 e_q ! ? rms uncertainty on q
character*9 Filt_1 ! Wavelength of the light curve (1)
character*6 Ref3 ! Reference for above parameters,
* in refs.dat file
*Note (1): UBVRIJ are Johnson standard filters, ubvy are Stroemgren filters
* and some values given in numbers are in Angstroem unit.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Absolute Parameters
integer*4 nr__3
parameter (nr__3=124) ! Number of records
character*124 ar__3 ! Full-size record
integer*4 Seq_3 ! Serial number of entry
character*13 Name_3 ! Name of star
real*4 M1 ! (solMass) ? Primary mass
real*4 e_M1 ! (solMass) ? rms uncertainty on M1
real*4 M2 ! (solMass) ? Secondary mass
real*4 e_M2 ! (solMass) ? rms uncertainty on M2
real*4 R1_1 ! (solRad) ? Primary radius
real*4 e_R1_1 ! (solRad) ? rms uncertainty on R1
real*4 R2_1 ! (solRad) ? Secondary radius
real*4 e_R2_1 ! (solRad) ? rms uncertainty on R2
integer*4 T1 ! ([K]) ? Primary Effective temperature
integer*4 e_T1 ! ([K]) ? rms uncertainty on T1
integer*4 T2 ! ([K]) ? Secondary Effective temperature
integer*4 e_T2 ! ([K]) ? rms uncertainty on T2
real*4 logg1 ! ([cm/s2]) ? Primary surface gravity
real*4 e_logg1 ! ([cm/s2]) ? rms uncertainty on logg1
real*4 logg2 ! ([cm/s2]) ? Secondary surface gravity
real*4 e_logg2 ! ([cm/s2]) ? rms uncertainty on logg2
character*8 Ref4 ! Reference for above parameters,
* in refs.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Candidates
integer*4 nr__4
parameter (nr__4=150) ! Number of records
character*108 ar__4 ! Full-size record
real*8 RA__1 ! (deg) Right Ascension J2000
real*8 Dec_1 ! (deg) Declination J2000
integer*4 Seq_4 ! Serial number of entry
character*13 Name_4 ! Name of star
integer*4 HIP_1 ! ? Hipparcos Catalogue number
integer*4 RAh_1 ! (h) Right Ascension (J2000)
integer*4 RAm_1 ! (min) Right Ascension (J2000)
real*4 RAs_1 ! (s) Right Ascension (J2000)
character*1 DE__1 ! Declination sign (J2000)
integer*4 DEd_1 ! (deg) Declination (J2000)
integer*4 DEm_1 ! (arcmin) Declination (J2000)
real*4 DEs_1 ! (arcsec) Declination (J2000)
real*8 T0_1 ! (d) ? HJD Epoch of primary minimum
real*8 Ps_1 ! (d) ? Orbital period
real*4 phi2 ! ? Phase of secondary minimum
real*4 mag_1 ! (mag) ? Magnitude at maximum brightness
real*4 Dmag_1 ! (mag) ? Magnitude changes brightness
character*2 Filt_2 ! Photometric system for magnitudes
character*10 SpType_1 ! MK spectral type
character*6 Ref5 ! Reference, in refs.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! Bibliography
integer*4 nr__5
parameter (nr__5=157) ! Number of records
character*222 ar__5 ! Full-size record
integer*4 Ref ! Reference number
character*19 BibCode ! BibCode
character*25 Aut ! Author's name
character*173 Com ! Comments
C=============================================================================
C Loading file 'table1.dat' ! Identifications of eclipsing binaries
C Format for file interpretation
1 format(
+ I3,1X,A13,1X,I6,1X,I6,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,
+ F4.1,1X,F5.2,1X,F4.2,1X,A2,1X,A16,1X,A10,1X,A5)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,124
read(1,'(A103)')ar__
read(ar__,1)
+ Seq,Name,HD,HIP,RAh,RAm,RAs,DE_,DEd,DEm,DEs,mag,Dmag,Filt,
+ SpType,Type,Ref1
if(ar__(19:24) .EQ. '') HD = iNULL__
if(ar__(26:31) .EQ. '') HIP = iNULL__
if(ar__(57:61) .EQ. '') mag = rNULL__
if(ar__(63:66) .EQ. '') Dmag = rNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAs .GE. 0) RA_=RA_+RAs/240.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DEs .GE. 0) Dec=Dec+DEs/3600.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)
+ Seq,Name,HD,HIP,RAh,RAm,RAs,DE_,DEd,DEm,DEs,mag,Dmag,Filt,
+ SpType,Type,Ref1
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Apsidal Motion Parameters
C Format for file interpretation
2 format(
+ I3,1X,A13,1X,F11.5,1X,A1,1X,F7.5,1X,F11.8,1X,F10.8,1X,F11.8,
+ 1X,F10.8,1X,A1,1X,F7.5,1X,F7.5,1X,F6.2,1X,F5.2,1X,A6,1X,A1,1X,
+ F10.3,2X,F9.3,1X,F6.3,1X,F5.3,1X,A10)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,124
read(1,'(A171)')ar__1
read(ar__1,2)
+ Seq_1,Name_1,T0,n_T0,e_T0,Ps,e_Ps,Pa,e_Pa,l_e,e,e_e,omega,
+ e_omega,n_U,l_U,U,e_U,logk2,e_logk2,Ref2
if(ar__1(19:29) .EQ. '') T0 = rNULL__
if(ar__1(33:39) .EQ. '') e_T0 = rNULL__
if(ar__1(41:51) .EQ. '') Ps = rNULL__
if(ar__1(53:62) .EQ. '') e_Ps = rNULL__
if(ar__1(64:74) .EQ. '') Pa = rNULL__
if(ar__1(76:85) .EQ. '') e_Pa = rNULL__
if(ar__1(89:95) .EQ. '') e = rNULL__
if(ar__1(97:103) .EQ. '') e_e = rNULL__
if(ar__1(105:110) .EQ. '') omega = rNULL__
if(ar__1(112:116) .EQ. '') e_omega = rNULL__
if(ar__1(127:136) .EQ. '') U = rNULL__
if(ar__1(139:147) .EQ. '') e_U = rNULL__
if(ar__1(149:154) .EQ. '') logk2 = rNULL__
if(ar__1(156:160) .EQ. '') e_logk2 = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1,Name_1,T0,n_T0,e_T0,Ps,e_Ps,Pa,e_Pa,l_e,e,e_e,omega,
+ e_omega,n_U,l_U,U,e_U,logk2,e_logk2,Ref2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Photometric Parameters
C Format for file interpretation
3 format(
+ I3,1X,A13,1X,F6.4,1X,F6.4,1X,A2,1X,F6.4,1X,F6.4,1X,F6.4,1X,
+ F6.4,1X,F6.4,1X,F6.4,1X,F6.3,1X,F5.3,1X,F6.4,1X,F6.4,1X,F5.3,
+ 1X,F5.3,1X,A9,1X,A6)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,238
read(1,'(A132)')ar__2
read(ar__2,3)
+ Seq_2,Name_2,L1,e_L1,n_L1,L2,e_L2,r1,e_r1,r2,e_r2,i,e_i,e_1,
+ e_e_1,q,e_q,Filt_1,Ref3
if(ar__2(19:24) .EQ. '') L1 = rNULL__
if(ar__2(26:31) .EQ. '') e_L1 = rNULL__
if(ar__2(36:41) .EQ. '') L2 = rNULL__
if(ar__2(43:48) .EQ. '') e_L2 = rNULL__
if(ar__2(50:55) .EQ. '') r1 = rNULL__
if(ar__2(57:62) .EQ. '') e_r1 = rNULL__
if(ar__2(64:69) .EQ. '') r2 = rNULL__
if(ar__2(71:76) .EQ. '') e_r2 = rNULL__
if(ar__2(78:83) .EQ. '') i = rNULL__
if(ar__2(85:89) .EQ. '') e_i = rNULL__
if(ar__2(91:96) .EQ. '') e_1 = rNULL__
if(ar__2(98:103) .EQ. '') e_e_1 = rNULL__
if(ar__2(105:109) .EQ. '') q = rNULL__
if(ar__2(111:115) .EQ. '') e_q = rNULL__
c ..............Just test output...........
write(6,3)
+ Seq_2,Name_2,L1,e_L1,n_L1,L2,e_L2,r1,e_r1,r2,e_r2,i,e_i,e_1,
+ e_e_1,q,e_q,Filt_1,Ref3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Absolute Parameters
C Format for file interpretation
4 format(
+ I3,1X,A13,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,I5,1X,I4,1X,I5,1X,I4,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,A8)
C Effective file loading
open(unit=1,file='table4.dat', status='old')
write(6,*) '....Loading file: table4.dat'
do i__=1,124
read(1,'(A124)')ar__3
read(ar__3,4)
+ Seq_3,Name_3,M1,e_M1,M2,e_M2,R1_1,e_R1_1,R2_1,e_R2_1,T1,e_T1,
+ T2,e_T2,logg1,e_logg1,logg2,e_logg2,Ref4
if(ar__3(19:24) .EQ. '') M1 = rNULL__
if(ar__3(26:30) .EQ. '') e_M1 = rNULL__
if(ar__3(32:37) .EQ. '') M2 = rNULL__
if(ar__3(39:43) .EQ. '') e_M2 = rNULL__
if(ar__3(45:50) .EQ. '') R1_1 = rNULL__
if(ar__3(52:56) .EQ. '') e_R1_1 = rNULL__
if(ar__3(58:63) .EQ. '') R2_1 = rNULL__
if(ar__3(65:69) .EQ. '') e_R2_1 = rNULL__
if(ar__3(71:75) .EQ. '') T1 = iNULL__
if(ar__3(77:80) .EQ. '') e_T1 = iNULL__
if(ar__3(82:86) .EQ. '') T2 = iNULL__
if(ar__3(88:91) .EQ. '') e_T2 = iNULL__
if(ar__3(93:97) .EQ. '') logg1 = rNULL__
if(ar__3(99:103) .EQ. '') e_logg1 = rNULL__
if(ar__3(105:109) .EQ. '') logg2 = rNULL__
if(ar__3(111:115) .EQ. '') e_logg2 = rNULL__
c ..............Just test output...........
write(6,4)
+ Seq_3,Name_3,M1,e_M1,M2,e_M2,R1_1,e_R1_1,R2_1,e_R2_1,T1,e_T1,
+ T2,e_T2,logg1,e_logg1,logg2,e_logg2,Ref4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Candidates
C Format for file interpretation
5 format(
+ I3,1X,A13,1X,I6,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ F9.3,1X,F10.6,1X,F5.3,1X,F5.2,1X,F5.2,1X,A2,1X,A10,1X,A6)
C Effective file loading
open(unit=1,file='table5.dat', status='old')
write(6,*) '....Loading file: table5.dat'
do i__=1,150
read(1,'(A108)')ar__4
read(ar__4,5)
+ Seq_4,Name_4,HIP_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ T0_1,Ps_1,phi2,mag_1,Dmag_1,Filt_2,SpType_1,Ref5
if(ar__4(19:24) .EQ. '') HIP_1 = iNULL__
if(ar__4(50:58) .EQ. '') T0_1 = rNULL__
if(ar__4(60:69) .EQ. '') Ps_1 = rNULL__
if(ar__4(71:75) .EQ. '') phi2 = rNULL__
if(ar__4(77:81) .EQ. '') mag_1 = rNULL__
if(ar__4(83:87) .EQ. '') Dmag_1 = rNULL__
c Derive coordinates RA__1 and Dec_1 from input data
c (RA__1 and Dec_1 are set to rNULL_ when unknown)
RA__1 = RAh_1
if(RAh_1 .GE. 0) RA__1=RAh_1*15.
if(RAm_1 .GE. 0) RA__1=RA__1+RAm_1/4.
if(RAs_1 .GE. 0) RA__1=RA__1+RAs_1/240.
Dec_1 = DEd_1
if(DEm_1 .GE. 0) Dec_1=Dec_1+DEm_1/60.
if(DEs_1 .GE. 0) Dec_1=Dec_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.Dec_1.GE.0) Dec_1=-Dec_1
c ..............Just test output...........
write(6,5)
+ Seq_4,Name_4,HIP_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ T0_1,Ps_1,phi2,mag_1,Dmag_1,Filt_2,SpType_1,Ref5
write(6,'(6H Pos: 2F8.4)') RA__1,Dec_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! Bibliography
C Format for file interpretation
6 format(I3,1X,A19,1X,A25,A173)
C Effective file loading
open(unit=1,file='refs.dat', status='old')
write(6,*) '....Loading file: refs.dat'
do i__=1,157
read(1,'(A222)')ar__5
read(ar__5,6)Ref,BibCode,Aut,Com
c ..............Just test output...........
write(6,6)Ref,BibCode,Aut,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end