Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/686/A177 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/A+A/686/A177 VELOCE I. Cepheids High-precision RVs (Anderson+, 2024)
*================================================================================
*VELOCE I. High-precision radial velocities of Cepheids
* Anderson R.I., Viviani G., Shetye S.S., Mowlavi N., Eyer L., Palaversa L.,
* Holl B., Blanco-Cuaresma S., Kravchenko K., Pawlak M., Cruz Reyes M.,
* Khan S., Netzel H.E., Loebling L., Papics P.I., Postel A., Roelens M.,
* Spetsieri Z.T., Thoul A., Zak J., Bonvin V., Martin D.V., Millon M.,
* Saesen S., Wyttenbach A., Figueira P., Marmier M., Prins S., Raskin G.,
* van Winckel H.
* <Astron. Astrophys. 686, A177 (2024)>
* =2024A&A...686A.177A (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Literature data zero-points
integer*4 nr__
parameter (nr__=34) ! Number of records
character*102 ar__ ! Full-size record
character*22 Ref (nr__) ! Literature or instrument reference
integer*4 Nstars (nr__) ! [1/68] Number of stars used to determine
* the zero-point
integer*4 Nclusters (nr__) ! [1/219] Number of clusters used to determine
* the zero-point
real*4 ZP (nr__) ! (km/s) Zero-point difference between VELOCE and
* the reference
real*4 e_ZP (nr__) ! (km/s) Uncertainty on the zero-point difference
character*59 BibCodes (nr__) ! Bibcode for the references
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Linear period change rates (dP/dt)
integer*4 nr__1
parameter (nr__1=146) ! Number of records
character*47 ar__1 ! Full-size record
character*9 Star (nr__1) ! Name of the Cepheid
real*8 dP_dt (nr__1) ! (s/yr) Linear period change rate
real*8 e_dP_dt (nr__1) ! (s/yr) Uncertainty in the linear period change rate
character*19 Notes (nr__1) ! Notes on the results
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Bona fide single-mode classical Cepheids
integer*4 nr__2
parameter (nr__2=258) ! Number of records
character*125 ar__2 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__2) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__2) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*19 Name (nr__2) ! Name of the Cepheid
integer*4 RAh (nr__2) ! (h) Right ascension (J2000)
integer*4 RAm (nr__2) ! (min) Right ascension (J2000)
integer*4 RAs (nr__2) ! (s) Right ascension (J2000)
character*1 DE_ (nr__2) ! Declination sign (J2000)
integer*4 DEd (nr__2) ! (deg) Declination (J2000)
integer*4 DEm (nr__2) ! (arcmin) Declination (J2000)
integer*4 DEs (nr__2) ! (arcsec) Declination (J2000)
character*19 GaiaDR3 (nr__2) ! Gaia DR3 source_id
real*4 Gmag (nr__2) ! (mag) ? Gaia G magnitude
real*4 Bp_Rp (nr__2) ! (mag) ? Gaia G_Bp_-G_Rp_ color
character*2 Pmode (nr__2) ! Pulsation mode classification according to
* Gaia DR3
character*1 Binflag (nr__2) ! [FTU] Binary flag
integer*4 NRV (nr__2) ! [1/929] Number of radial velocity measurements
real*4 NFS (nr__2) ! [2/19]? Number of terms in the Fourier series
real*8 P (nr__2) ! (d) ? Pulsation period
real*4 e_P (nr__2) ! (d) ? Uncertainty on the pulsation period
real*8 E (nr__2) ! (d) ? Reference epoch at which phi=0 (BJD-2400000)
real*4 e_E (nr__2) ! (d) ? Uncertainty on the reference epoch
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Fourier parameters of modeled Cepheids
integer*4 nr__3
parameter (nr__3=219) ! Number of records
character*120 ar__3 ! Full-size record
character*19 Name_1 (nr__3) ! Name of the Cepheid
integer*4 Polydeg (nr__3) ! [1/20]? Degree of the polynomial used in
* conjunction with the Fourier series
integer*4 DoF (nr__3) ! [2/901] Degrees of freedom
real*8 Chi2 (nr__3) ! Chi-squared (Chi^2^)
real*8 vgamma (nr__3) ! (km/s) Radial velocity offset
real*4 P2P (nr__3) ! (km/s) Peak-to-peak radial velocity amplitude
real*4 A1 (nr__3) ! (km/s) Amplitude of the first Fourier component
real*4 phi1 (nr__3) ! (rad) Phase of the first Fourier component
real*4 R21 (nr__3) ! Ratio of the second to the first
* Fourier harmonics
real*4 phi21 (nr__3) ! (rad) Phase difference between the
* second and the first Fourier harmonics
real*8 R31 (nr__3) ! ? Ratio of the third to the first
* Fourier harmonics
real*4 phi31 (nr__3) ! (rad) ? Phase difference between the
* third and the first Fourier harmonics
real*8 R41 (nr__3) ! ? Ratio of the fourth to the first
* Fourier harmonics
real*4 phi41 (nr__3) ! (rad) ? Phase difference between the
* fourth and the first Fourier harmonics
real*8 R51 (nr__3) ! ? Ratio of the fifth to the first
* Fourier harmonics
real*4 phi51 (nr__3) ! (rad) ? Phase difference between the
* fifth and the first Fourier harmonics
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Other targets in VELOCE
integer*4 nr__4
parameter (nr__4=164) ! Number of records
character*142 ar__4 ! Full-size record
character*26 Name_2 (nr__4) ! Name of the target
integer*4 NRV_1 (nr__4) ! [0/199] Number of radial velocity
* measurements
character*1 CCFpeak (nr__4) ! [yn] Peak of the cross-correlation function
* flag
character*1 CCFodd (nr__4) ! [yn?] CCF looks significantly different from
* those of classical Cepheids flag
character*9 Ncomp (nr__4) ! Specify if more than one absorption
* component was noticed
character*5 Emission (nr__4) ! Exhibits signs of emission features
character*2 vRV (nr__4) ! [yn? NA] RV variability flag
character*2 vFWHW (nr__4) ! [yn? NA] FWHM variability flag
character*2 vAsym (nr__4) ! [yn? NA] Line asymmetry variability flag
character*2 vLineDepth (nr__4) ! [yn? NA] Line depth variability flag
character*46 CCFnote (nr__4) ! Additional notes on the CCF
integer*8 GaiaDR3_1 (nr__4) ! Gaia DR3 source_id
character*6 Class (nr__4) ! Gaia DR3 classification
character*5 G22 (nr__4) ! Gavras et al. (2023A&A...674A..22G,
* Cat. J/A+A/674/A22) classification
C=============================================================================
C Loading file 'table5.dat' ! Literature data zero-points
C Format for file interpretation
1 format(A22,1X,I2,1X,I3,1X,F6.3,1X,F5.3,1X,A59)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,34
read(1,'(A102)')ar__
read(ar__,1)
+ Ref(i__),Nstars(i__),Nclusters(i__),ZP(i__),e_ZP(i__),
+ BibCodes(i__)
c ..............Just test output...........
write(6,1)
+ Ref(i__),Nstars(i__),Nclusters(i__),ZP(i__),e_ZP(i__),
+ BibCodes(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Linear period change rates (dP/dt)
C Format for file interpretation
2 format(A9,1X,F9.3,1X,F7.3,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,146
read(1,'(A47)')ar__1
read(ar__1,2)Star(i__),dP_dt(i__),e_dP_dt(i__),Notes(i__)
c ..............Just test output...........
write(6,2)Star(i__),dP_dt(i__),e_dP_dt(i__),Notes(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea1.dat' ! Bona fide single-mode classical Cepheids
C Format for file interpretation
3 format(
+ A19,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,A19,1X,F5.1,1X,
+ F4.1,1X,A2,1X,A1,1X,I3,1X,F5.1,1X,F10.6,1X,E7.2,1X,F13.6,1X,
+ E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,258
read(1,'(A125)')ar__2
read(ar__2,3)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),GaiaDR3(i__),Gmag(i__),Bp_Rp(i__),
+ Pmode(i__),Binflag(i__),NRV(i__),NFS(i__),P(i__),e_P(i__),
+ E(i__),e_E(i__)
if(ar__2(60:64) .EQ. '') Gmag(i__) = rNULL__
if(ar__2(66:69) .EQ. '') Bp_Rp(i__) = rNULL__
if(ar__2(80:84) .EQ. '') NFS(i__) = rNULL__
if(ar__2(86:95) .EQ. '') P(i__) = rNULL__
if(ar__2(97:103) .EQ. '') e_P(i__) = rNULL__
if(ar__2(105:117) .EQ. '') E(i__) = rNULL__
if(ar__2(119:125) .EQ. '') e_E(i__) = rNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,3)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),GaiaDR3(i__),Gmag(i__),Bp_Rp(i__),
+ Pmode(i__),Binflag(i__),NRV(i__),NFS(i__),P(i__),e_P(i__),
+ E(i__),e_E(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Fourier parameters of modeled Cepheids
C Format for file interpretation
4 format(
+ A19,1X,I2,1X,I3,1X,F9.1,1X,F8.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
+ F6.4,1X,F4.2,1X,F7.4,2X,F4.2,1X,F7.4,1X,F5.2,1X,F7.4,2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,219
read(1,'(A120)')ar__3
read(ar__3,4)
+ Name_1(i__),Polydeg(i__),DoF(i__),Chi2(i__),vgamma(i__),
+ P2P(i__),A1(i__),phi1(i__),R21(i__),phi21(i__),R31(i__),
+ phi31(i__),R41(i__),phi41(i__),R51(i__),phi51(i__)
if(ar__3(21:22) .EQ. '') Polydeg(i__) = iNULL__
if(ar__3(80:86) .EQ. '') R31(i__) = rNULL__
if(ar__3(89:92) .EQ. '') phi31(i__) = rNULL__
if(ar__3(94:100) .EQ. '') R41(i__) = rNULL__
if(ar__3(102:106) .EQ. '') phi41(i__) = rNULL__
if(ar__3(108:114) .EQ. '') R51(i__) = rNULL__
if(ar__3(117:120) .EQ. '') phi51(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_1(i__),Polydeg(i__),DoF(i__),Chi2(i__),vgamma(i__),
+ P2P(i__),A1(i__),phi1(i__),R21(i__),phi21(i__),R31(i__),
+ phi31(i__),R41(i__),phi41(i__),R51(i__),phi51(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Other targets in VELOCE
C Format for file interpretation
5 format(
+ A26,1X,I3,1X,A1,1X,A1,1X,A9,1X,A5,1X,A2,1X,A2,1X,A2,1X,A2,1X,
+ A46,1X,I19,1X,A6,1X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,164
read(1,'(A142)')ar__4
read(ar__4,5)
+ Name_2(i__),NRV_1(i__),CCFpeak(i__),CCFodd(i__),Ncomp(i__),
+ Emission(i__),vRV(i__),vFWHW(i__),vAsym(i__),vLineDepth(i__),
+ CCFnote(i__),GaiaDR3_1(i__),Class(i__),G22(i__)
c ..............Just test output...........
write(6,5)
+ Name_2(i__),NRV_1(i__),CCFpeak(i__),CCFodd(i__),Ncomp(i__),
+ Emission(i__),vRV(i__),vFWHW(i__),vAsym(i__),vLineDepth(i__),
+ CCFnote(i__),GaiaDR3_1(i__),Class(i__),G22(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end