Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/531/A49 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-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+A/531/A49 UBV and spectroscopic variations of V2368 Oph (Harmanec+, 2011)
*================================================================================
*V2368 Ophiuchi: An eclipsing and double-lined spectroscopic binary used as a
*photometric comparison star for U Oph.
* Harmanec P., Bozic H., Mayer P., Eenens P., Broz M., Wolf M., Yang S.,
* Slechta M., Ruzdjak D., Sudar D., Ak H.
* <Astron. Astrophys. 531, A49 (2011)>
* =2011A&A...531A..49H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2a.dat' ! All-sky Johnson V photometry of V2368 Oph
integer*4 nr__
parameter (nr__=158) ! Number of records
character*37 ar__ ! Full-size record
real*8 HJD ! (d) Heliocentric Julian date
real*4 Vmag ! (mag) Johnson V magnitude
real*4 Bmag ! (mag) ? Johnson B magnitude
real*4 Umag ! (mag) ? Johnson U magnitude
real*4 w ! Weight of the data point
integer*4 Obs ! Code of the observatory (1)
*Note (1): Code of the observatory as follows:
* 1 = Hvar observatory
* 30 = San Pedro Martir 0.84-m
* 61 = Hipparcos satellite Hp photometry transformed to Johnson V
* 66 = Tubitak National Observatory
* 93 = ASAS3 automatic telescope
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2b.dat' ! Standard UBV photometry of V2368 Oph
from three stations
integer*4 nr__1
parameter (nr__1=585) ! Number of records
character*37 ar__1 ! Full-size record
real*8 HJD_1 ! (d) Heliocentric Julian date
real*4 Vmag_1 ! (mag) Johnson V magnitude
real*4 Bmag_1 ! (mag) ? Johnson B magnitude
real*4 Umag_1 ! (mag) ? Johnson U magnitude
real*4 w_1 ! Weight of the data point
integer*4 Obs_1 ! Code of the observatory (1)
*Note (1): Code of the observatory as follows:
* 1 = Hvar observatory
* 30 = San Pedro Martir 0.84-m
* 61 = Hipparcos satellite Hp photometry transformed to Johnson V
* 66 = Tubitak National Observatory
* 93 = ASAS3 automatic telescope
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3a.dat' ! H{alpha} radial velocities of the primary
of V2368 Oph measured with the SPEFO program
(Horn et al. 1996A&A...309..521H)
integer*4 nr__2
parameter (nr__2=153) ! Number of records
character*40 ar__2 ! Full-size record
real*8 HJD_2 ! (d) Heliocentric Julian date
real*4 RV1 ! (km/s) Radial velocity of the primary
real*4 RV2 ! (km/s) ? Radial velocity of the secondary (2)
real*4 w_2 ! ? Weight of the data point (2)
character*4 Line ! [MgII Ha SiII] Line (G1) (2)
*Note (2): Only in table3b.dat
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3b.dat' ! SiII, H{alpha}, and MgII radial velocities of both
components of V2368 Oph measured via Gaussian
fits to the line profiles
integer*4 nr__3
parameter (nr__3=135) ! Number of records
character*40 ar__3 ! Full-size record
real*8 HJD_3 ! (d) Heliocentric Julian date
real*4 RV1_1 ! (km/s) Radial velocity of the primary
real*4 RV2_1 ! (km/s) ? Radial velocity of the secondary (2)
real*4 w_3 ! ? Weight of the data point (2)
character*4 Line_1 ! [MgII Ha SiII] Line (G1) (2)
*Note (2): Only in table3b.dat
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rjd.dat' ! List of H{alpha}, SiII and MgII spectra
integer*4 nr__4
parameter (nr__4=217) ! Number of records
character*48 ar__4 ! Full-size record
integer*4 lam_min ! (0.1nm) Lower limit of wavelength range
integer*4 lam_max ! (0.1nm) Upper limit of wavelength range
character*4 Line_2 ! [MgII Ha SiII] Line (G1)
character*12 FileName ! Filename of individual spectra
real*8 HJD_4 ! (d) Heliocentric Julian date
real*4 w_4 ! weights of individual spectra from S/N
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'spectra/*' ! Individual spectra
integer*4 nr__5
parameter (nr__5=188) ! Number of records
character*18 ar__5 ! Full-size record
real*8 lambda ! (0.1nm) Wavelength
real*4 RFlux ! Relative rectified flux
C=============================================================================
C Loading file 'table2a.dat' ! All-sky Johnson V photometry of V2368 Oph
C Format for file interpretation
1 format(F12.4,1X,F5.3,1X,F5.3,1X,F5.3,1X,F3.1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table2a.dat')
write(6,*) '....Loading file: table2a.dat'
do i__=1,158
read(1,'(A37)')ar__
read(ar__,1)HJD,Vmag,Bmag,Umag,w,Obs
if(ar__(20:24) .EQ. '') Bmag = rNULL__
if(ar__(26:30) .EQ. '') Umag = rNULL__
c ..............Just test output...........
write(6,1)HJD,Vmag,Bmag,Umag,w,Obs
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2b.dat' ! Standard UBV photometry of V2368 Oph
* from three stations
C Format for file interpretation
2 format(F12.4,1X,F5.3,1X,F5.3,1X,F5.3,1X,F3.1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table2b.dat')
write(6,*) '....Loading file: table2b.dat'
do i__=1,585
read(1,'(A37)')ar__1
read(ar__1,2)HJD_1,Vmag_1,Bmag_1,Umag_1,w_1,Obs_1
if(ar__1(20:24) .EQ. '') Bmag_1 = rNULL__
if(ar__1(26:30) .EQ. '') Umag_1 = rNULL__
c ..............Just test output...........
write(6,2)HJD_1,Vmag_1,Bmag_1,Umag_1,w_1,Obs_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3a.dat' ! H{alpha} radial velocities of the primary
* of V2368 Oph measured with the SPEFO program
* (Horn et al. 1996A&A...309..521H)
C Format for file interpretation
3 format(F12.4,2X,F6.2,1X,F6.1,2X,F5.3,2X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table3a.dat')
write(6,*) '....Loading file: table3a.dat'
do i__=1,153
read(1,'(A40)')ar__2
read(ar__2,3)HJD_2,RV1,RV2,w_2,Line
if(ar__2(22:27) .EQ. '') RV2 = rNULL__
if(ar__2(30:34) .EQ. '') w_2 = rNULL__
c ..............Just test output...........
write(6,3)HJD_2,RV1,RV2,w_2,Line
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3b.dat' ! SiII, H{alpha}, and MgII radial velocities of both
* components of V2368 Oph measured via Gaussian
* fits to the line profiles
C Format for file interpretation
4 format(F12.4,2X,F6.2,1X,F6.1,2X,F5.3,2X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table3b.dat')
write(6,*) '....Loading file: table3b.dat'
do i__=1,135
read(1,'(A40)')ar__3
read(ar__3,4)HJD_3,RV1_1,RV2_1,w_3,Line_1
if(ar__3(22:27) .EQ. '') RV2_1 = rNULL__
if(ar__3(30:34) .EQ. '') w_3 = rNULL__
c ..............Just test output...........
write(6,4)HJD_3,RV1_1,RV2_1,w_3,Line_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rjd.dat' ! List of H{alpha}, SiII and MgII spectra
C Format for file interpretation
5 format(I4,1X,I4,1X,A4,1X,A12,1X,F12.4,3X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'rjd.dat')
write(6,*) '....Loading file: rjd.dat'
do i__=1,217
read(1,'(A48)')ar__4
read(ar__4,5)lam_min,lam_max,Line_2,FileName,HJD_4,w_4
c ..............Just test output...........
write(6,5)lam_min,lam_max,Line_2,FileName,HJD_4,w_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'spectra/*' ! Individual spectra
C Format for file interpretation
6 format(F10.4,2X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'spectra/*')
write(6,*) '....Loading file: spectra/*'
do i__=1,188
read(1,'(A18)')ar__5
read(ar__5,6)lambda,RFlux
c ..............Just test output...........
write(6,6)lambda,RFlux
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end