Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/405/231 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-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__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/A+A/405/231 Algol-type eclipsing binaries differential photometry (Kim+, 2003)
*================================================================================
*Search for A-F spectral type pulsating components in
*Algol-type eclipsing binary systems.
* Kim S.-L., Lee J.W., Kwon S.-G., Youn J.-H., Mkrtichian D.E., Kim C.
* <Astron. Astrophys. 405, 231 (2003)>
* =2003A&A...405..231K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'stars.dat' ! Star positions and magnitudes (from Simbad)
integer*4 nr__
parameter (nr__=28) ! Number of records
character*76 ar__ ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000.0
real*8 DEdeg ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*15 Star ! Star name
character*1 n_Star ! [c] c for comparison star
integer*4 RAh ! (h) Right ascension (J2000.0)
integer*4 RAm ! (min) Right ascension (J2000.0)
real*4 RAs ! (s) Right ascension (J2000.0)
character*1 DE_ ! Declination sign (J2000.0)
integer*4 DEd ! (deg) Declination (J2000.0)
integer*4 DEm ! (arcmin) Declination (J2000.0)
integer*4 DEs ! (arcsec) Declination (J2000.0)
real*4 Bmag ! (mag) ? B magnitude
real*4 Vmag ! (mag) ? V magnitude
character*8 SpType ! Spectral type
character*9 Table ! Name of the table containing the data
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'abper.dat' ! Data of AB Per (comp. star HD 275605)
integer*4 nr__1
parameter (nr__1=522) ! Number of records
character*28 ar__1 ! Full-size record
real*8 HJD ! (d) Heliocentric Julian Date
real*4 Airmass ! Airmass at a half of each exposure
character*1 Band ! [BV] Band for the magnitude
real*4 Dmag ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'acuma.dat' ! Data of AC UMa (comp. star GSC 04131-00208)
integer*4 nr__2
parameter (nr__2=105) ! Number of records
character*28 ar__2 ! Full-size record
real*8 HJD_1 ! (d) Heliocentric Julian Date
real*4 Airmass_1 ! Airmass at a half of each exposure
character*1 Band_1 ! [BV] Band for the magnitude
real*4 Dmag_1 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cpand.dat' ! Data of CP And (comp. star GSC 03281-01169)
integer*4 nr__3
parameter (nr__3=74) ! Number of records
character*28 ar__3 ! Full-size record
real*8 HJD_2 ! (d) Heliocentric Julian Date
real*4 Airmass_2 ! Airmass at a half of each exposure
character*1 Band_2 ! [BV] Band for the magnitude
real*4 Dmag_2 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rscep.dat' ! Data of RS Cep (comp. star GSC 04519-01555)
integer*4 nr__4
parameter (nr__4=143) ! Number of records
character*28 ar__4 ! Full-size record
real*8 HJD_3 ! (d) Heliocentric Julian Date
real*4 Airmass_3 ! Airmass at a half of each exposure
character*1 Band_3 ! [BV] Band for the magnitude
real*4 Dmag_3 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rvper.dat' ! Data of RV Per (comp. star GSC 02366-00559)
integer*4 nr__5
parameter (nr__5=117) ! Number of records
character*28 ar__5 ! Full-size record
real*8 HJD_4 ! (d) Heliocentric Julian Date
real*4 Airmass_4 ! Airmass at a half of each exposure
character*1 Band_4 ! [BV] Band for the magnitude
real*4 Dmag_4 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rwcrb.dat' ! Data of RW CrB (comp. star GSC 02039-00004)
integer*4 nr__6
parameter (nr__6=86) ! Number of records
character*28 ar__6 ! Full-size record
real*8 HJD_5 ! (d) Heliocentric Julian Date
real*4 Airmass_5 ! Airmass at a half of each exposure
character*1 Band_5 ! [BV] Band for the magnitude
real*4 Dmag_5 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rwper.dat' ! Data of RW Per (comp. star HD 276248)
integer*4 nr__7
parameter (nr__7=176) ! Number of records
character*28 ar__7 ! Full-size record
real*8 HJD_6 ! (d) Heliocentric Julian Date
real*4 Airmass_6 ! Airmass at a half of each exposure
character*1 Band_6 ! [BV] Band for the magnitude
real*4 Dmag_6 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rxhya.dat' ! Data of RX Hya (comp. star BD-07 2718)
integer*4 nr__8
parameter (nr__8=139) ! Number of records
character*28 ar__8 ! Full-size record
real*8 HJD_7 ! (d) Heliocentric Julian Date
real*4 Airmass_7 ! Airmass at a half of each exposure
character*1 Band_7 ! [BV] Band for the magnitude
real*4 Dmag_7 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'stper.dat' ! Data of ST Per (comp. star GSC 02847-00136)
integer*4 nr__9
parameter (nr__9=206) ! Number of records
character*28 ar__9 ! Full-size record
real*8 HJD_8 ! (d) Heliocentric Julian Date
real*4 Airmass_8 ! Airmass at a half of each exposure
character*1 Band_8 ! [BV] Band for the magnitude
real*4 Dmag_8 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'twdra.dat' ! Data of TW Dra (comp. star HD 139549)
integer*4 nr__10
parameter (nr__10=216) ! Number of records
character*28 ar__10 ! Full-size record
real*8 HJD_9 ! (d) Heliocentric Julian Date
real*4 Airmass_9 ! Airmass at a half of each exposure
character*1 Band_9 ! [BV] Band for the magnitude
real*4 Dmag_9 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'vzleo.dat' ! Data of VZ Leo (comp. star GSC 01403-01558)
integer*4 nr__11
parameter (nr__11=64) ! Number of records
character*28 ar__11 ! Full-size record
real*8 HJD_10 ! (d) Heliocentric Julian Date
real*4 Airmass_10 ! Airmass at a half of each exposure
character*1 Band_10 ! [BV] Band for the magnitude
real*4 Dmag_10 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'xtri.dat' ! Data of X Tri (comp. star GSC 01763-02015)
integer*4 nr__12
parameter (nr__12=312) ! Number of records
character*28 ar__12 ! Full-size record
real*8 HJD_11 ! (d) Heliocentric Julian Date
real*4 Airmass_11 ! Airmass at a half of each exposure
character*1 Band_11 ! [BV] Band for the magnitude
real*4 Dmag_11 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'xzuma.dat' ! Data of XZ UMa (comp. star GSC 03429-00621)
integer*4 nr__13
parameter (nr__13=147) ! Number of records
character*28 ar__13 ! Full-size record
real*8 HJD_12 ! (d) Heliocentric Julian Date
real*4 Airmass_12 ! Airmass at a half of each exposure
character*1 Band_12 ! [BV] Band for the magnitude
real*4 Dmag_12 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'zdra.dat' ! Data of Z Dra (comp. star GSC 04396-01170)
integer*4 nr__14
parameter (nr__14=55) ! Number of records
character*28 ar__14 ! Full-size record
real*8 HJD_13 ! (d) Heliocentric Julian Date
real*4 Airmass_13 ! Airmass at a half of each exposure
character*1 Band_13 ! [BV] Band for the magnitude
real*4 Dmag_13 ! (mag) B or V differential magnitudes of variable
* star (Var-Comp)
C=============================================================================
C Loading file 'stars.dat' ! Star positions and magnitudes (from Simbad)
C Format for file interpretation
1 format(
+ A15,2X,A1,4X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,3X,F5.2,1X,
+ F5.2,1X,A8,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'stars.dat')
write(6,*) '....Loading file: stars.dat'
do i__=1,28
read(1,'(A76)')ar__
read(ar__,1)
+ Star,n_Star,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Bmag,Vmag,SpType,Table
if(ar__(46:50) .EQ. '') Bmag = rNULL__
if(ar__(52:56) .EQ. '') Vmag = rNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Star,n_Star,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Bmag,Vmag,SpType,Table
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'abper.dat' ! Data of AB Per (comp. star HD 275605)
C Format for file interpretation
2 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'abper.dat')
write(6,*) '....Loading file: abper.dat'
do i__=1,522
read(1,'(A28)')ar__1
read(ar__1,2)HJD,Airmass,Band,Dmag
c ..............Just test output...........
write(6,2)HJD,Airmass,Band,Dmag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'acuma.dat' ! Data of AC UMa (comp. star GSC 04131-00208)
C Format for file interpretation
3 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'acuma.dat')
write(6,*) '....Loading file: acuma.dat'
do i__=1,105
read(1,'(A28)')ar__2
read(ar__2,3)HJD_1,Airmass_1,Band_1,Dmag_1
c ..............Just test output...........
write(6,3)HJD_1,Airmass_1,Band_1,Dmag_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cpand.dat' ! Data of CP And (comp. star GSC 03281-01169)
C Format for file interpretation
4 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'cpand.dat')
write(6,*) '....Loading file: cpand.dat'
do i__=1,74
read(1,'(A28)')ar__3
read(ar__3,4)HJD_2,Airmass_2,Band_2,Dmag_2
c ..............Just test output...........
write(6,4)HJD_2,Airmass_2,Band_2,Dmag_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rscep.dat' ! Data of RS Cep (comp. star GSC 04519-01555)
C Format for file interpretation
5 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'rscep.dat')
write(6,*) '....Loading file: rscep.dat'
do i__=1,143
read(1,'(A28)')ar__4
read(ar__4,5)HJD_3,Airmass_3,Band_3,Dmag_3
c ..............Just test output...........
write(6,5)HJD_3,Airmass_3,Band_3,Dmag_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rvper.dat' ! Data of RV Per (comp. star GSC 02366-00559)
C Format for file interpretation
6 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'rvper.dat')
write(6,*) '....Loading file: rvper.dat'
do i__=1,117
read(1,'(A28)')ar__5
read(ar__5,6)HJD_4,Airmass_4,Band_4,Dmag_4
c ..............Just test output...........
write(6,6)HJD_4,Airmass_4,Band_4,Dmag_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rwcrb.dat' ! Data of RW CrB (comp. star GSC 02039-00004)
C Format for file interpretation
7 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'rwcrb.dat')
write(6,*) '....Loading file: rwcrb.dat'
do i__=1,86
read(1,'(A28)')ar__6
read(ar__6,7)HJD_5,Airmass_5,Band_5,Dmag_5
c ..............Just test output...........
write(6,7)HJD_5,Airmass_5,Band_5,Dmag_5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rwper.dat' ! Data of RW Per (comp. star HD 276248)
C Format for file interpretation
8 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'rwper.dat')
write(6,*) '....Loading file: rwper.dat'
do i__=1,176
read(1,'(A28)')ar__7
read(ar__7,8)HJD_6,Airmass_6,Band_6,Dmag_6
c ..............Just test output...........
write(6,8)HJD_6,Airmass_6,Band_6,Dmag_6
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rxhya.dat' ! Data of RX Hya (comp. star BD-07 2718)
C Format for file interpretation
9 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'rxhya.dat')
write(6,*) '....Loading file: rxhya.dat'
do i__=1,139
read(1,'(A28)')ar__8
read(ar__8,9)HJD_7,Airmass_7,Band_7,Dmag_7
c ..............Just test output...........
write(6,9)HJD_7,Airmass_7,Band_7,Dmag_7
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'stper.dat' ! Data of ST Per (comp. star GSC 02847-00136)
C Format for file interpretation
10 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'stper.dat')
write(6,*) '....Loading file: stper.dat'
do i__=1,206
read(1,'(A28)')ar__9
read(ar__9,10)HJD_8,Airmass_8,Band_8,Dmag_8
c ..............Just test output...........
write(6,10)HJD_8,Airmass_8,Band_8,Dmag_8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'twdra.dat' ! Data of TW Dra (comp. star HD 139549)
C Format for file interpretation
11 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'twdra.dat')
write(6,*) '....Loading file: twdra.dat'
do i__=1,216
read(1,'(A28)')ar__10
read(ar__10,11)HJD_9,Airmass_9,Band_9,Dmag_9
c ..............Just test output...........
write(6,11)HJD_9,Airmass_9,Band_9,Dmag_9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'vzleo.dat' ! Data of VZ Leo (comp. star GSC 01403-01558)
C Format for file interpretation
12 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'vzleo.dat')
write(6,*) '....Loading file: vzleo.dat'
do i__=1,64
read(1,'(A28)')ar__11
read(ar__11,12)HJD_10,Airmass_10,Band_10,Dmag_10
c ..............Just test output...........
write(6,12)HJD_10,Airmass_10,Band_10,Dmag_10
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'xtri.dat' ! Data of X Tri (comp. star GSC 01763-02015)
C Format for file interpretation
13 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'xtri.dat')
write(6,*) '....Loading file: xtri.dat'
do i__=1,312
read(1,'(A28)')ar__12
read(ar__12,13)HJD_11,Airmass_11,Band_11,Dmag_11
c ..............Just test output...........
write(6,13)HJD_11,Airmass_11,Band_11,Dmag_11
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'xzuma.dat' ! Data of XZ UMa (comp. star GSC 03429-00621)
C Format for file interpretation
14 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'xzuma.dat')
write(6,*) '....Loading file: xzuma.dat'
do i__=1,147
read(1,'(A28)')ar__13
read(ar__13,14)HJD_12,Airmass_12,Band_12,Dmag_12
c ..............Just test output...........
write(6,14)HJD_12,Airmass_12,Band_12,Dmag_12
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'zdra.dat' ! Data of Z Dra (comp. star GSC 04396-01170)
C Format for file interpretation
15 format(F13.5,1X,F5.3,1X,A1,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'zdra.dat')
write(6,*) '....Loading file: zdra.dat'
do i__=1,55
read(1,'(A28)')ar__14
read(ar__14,15)HJD_13,Airmass_13,Band_13,Dmag_13
c ..............Just test output...........
write(6,15)HJD_13,Airmass_13,Band_13,Dmag_13
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end