Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/453/309 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 2024-Apr-16
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/453/309 HARPS observations of nine galactic Cepheids (Nardetto+, 2006)
*================================================================================
*High resolution spectroscopy for Cepheids distance determination.
*I. Line asymmetry.
* Nardetto N., Mourard D., Kervella P., Mathias P., Merand A., Bersier D.
* <Astron. Astrophys., 453, 309-319 (2006)>
* =2006A&A...453..309N
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Observed sample of Cepheids (sorted by
increasing period)
integer*4 nr__
parameter (nr__=9) ! Number of records
character*51 ar__ ! Full-size record
character*8 Name (nr__) ! Star name
integer*4 HD (nr__) ! HD number
real*8 Per (nr__) ! (d) Pulsation period (1)
real*8 T0 (nr__) ! (d) Reference Julian date (1)
integer*4 Nsp (nr__) ! Number of spectra
integer*4 Nc (nr__) ! Number of cycles
real*4 Vmag (nr__) ! (mag) Visible magnitude from Berdnikov et al.
* (2000, Cat. <J/A+AS/143/211>
*Note (1): For l Car, the reference Julian date (T0) and the pulsation period (P)
* used to compute the phase are from Szabados (1989ApJ...337L..29S).
* For others stars we used ephemeris from Berdnikov et al.
* (2001JAD.....7....3B).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table345.dat' ! HARPS observations results
integer*4 nr__1
parameter (nr__1=126) ! Number of records
character*101 ar__1 ! Full-size record
character*8 Name_1 (nr__1) ! Star name
real*8 JD (nr__1) ! (d) Average Julian data of observation
real*4 Phase (nr__1) ! Average pulsation phase of observation.
* For ephemeris, see the Paper
integer*4 Cycle (nr__1) ! Pulsating cycle of the star corresponding
* to observation
integer*4 NSpec (nr__1) ! Number of spectra associated to observation (1)
real*4 RVg (nr__1) ! (km/s) Gaussian fit radial velocity
real*4 e_RVg (nr__1) ! (km/s) rms uncertainty on RVg
real*4 RVm (nr__1) ! (km/s) Minimum radial velocity (2)
real*4 e_RVm (nr__1) ! (km/s) rms uncertainty on RVm
real*4 RVc (nr__1) ! (km/s) Radial velocity corresponding to the first
* moment of the spectral line
real*4 e_RVc (nr__1) ! (km/s) rms uncertainty on RVc
real*4 FWHM (nr__1) ! (0.1nm) Full-Width at Half Maximum (2)
real*4 e_FWHM (nr__1) ! (0.1nm) rms uncertainty to FWHM
real*4 D (nr__1) ! Line depth. Uncertainties are of about 10^-4^ (2)
real*4 A (nr__1) ! (%) Asymmetry (2)
real*4 e_A (nr__1) ! (%) Uncertainty corresponding to the asymmetry
integer*4 SNR (nr__1) ! Observational spectral line signal to noise ratio
real*4 Chi2 (nr__1) ! Reduced Chi2 factor corresponding to the
* bi-Gaussian fit
*Note (1): Results corresponding to these spectra are averaged.
*Note (2): Quantity derived from a bi-Gaussian fit. See paper for more detail.
C=============================================================================
C Loading file 'table1.dat' ! Observed sample of Cepheids (sorted by
* increasing period)
C Format for file interpretation
1 format(A8,1X,I6,1X,F10.7,1X,F12.4,1X,I3,1X,I2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,9
read(1,'(A51)')ar__
read(ar__,1)
+ Name(i__),HD(i__),Per(i__),T0(i__),Nsp(i__),Nc(i__),Vmag(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),HD(i__),Per(i__),T0(i__),Nsp(i__),Nc(i__),Vmag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table345.dat' ! HARPS observations results
C Format for file interpretation
2 format(
+ A8,1X,F10.2,1X,F4.2,1X,I2,1X,I1,1X,F6.2,1X,F4.2,1X,F6.2,1X,
+ F4.2,1X,F6.2,1X,F4.2,1X,F5.3,1X,F5.3,1X,F4.2,1X,F5.1,1X,F3.1,
+ 1X,I3,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table345.dat')
write(6,*) '....Loading file: table345.dat'
do i__=1,126
read(1,'(A101)')ar__1
read(ar__1,2)
+ Name_1(i__),JD(i__),Phase(i__),Cycle(i__),NSpec(i__),RVg(i__),
+ e_RVg(i__),RVm(i__),e_RVm(i__),RVc(i__),e_RVc(i__),FWHM(i__),
+ e_FWHM(i__),D(i__),A(i__),e_A(i__),SNR(i__),Chi2(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),JD(i__),Phase(i__),Cycle(i__),NSpec(i__),RVg(i__),
+ e_RVg(i__),RVm(i__),e_RVm(i__),RVc(i__),e_RVc(i__),FWHM(i__),
+ e_FWHM(i__),D(i__),A(i__),e_A(i__),SNR(i__),Chi2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end