FORTRAN Generation
(/./ftp/cats/J/ApJ/704/964)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/704/964 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-Mar-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/ApJ/704/964     Spectroscopic study of Of^+^ supergiants    (De Becker+, 2009)
*================================================================================
*Spectroscopic study of the variability of three northern Of^+^ supergiants.
*    De Becker M., Rauw G., Linder N.
*   <Astrophys. J., 704, 964-974 (2009)>
*   =2009ApJ...704..964D
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! Journal of the observations for HD 14947,
                             HD 15570 and HD 16691

      integer*4 nr__
      parameter (nr__=210)	! Number of records
      character*36 ar__   	! Full-size record

      integer*4     HD         (nr__) ! HD star number (14947, 15570, or 16691)
      integer*4     Seq        (nr__) ! Running index number from original table of
*                                   the paper
      real*8        HJD        (nr__) ! (d) ? Heliocentric Julian date (-240000)
      character*10  Date       (nr__) ! ("YYYY/MM/DD") ? Date of observation
      character*4   Dom        (nr__) ! [blue/red ] Spectral domain of observation

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! Stellar parameters of HD 14947, HD 15570, and
                             HD 16691, and estimated minimum and maximum
                             rotation periods

      integer*4 nr__1
      parameter (nr__1=3)	! Number of records
      character*51 ar__1  	! Full-size record

      integer*4     HD_1       (nr__1) ! HD star number (14947, 15570, or 16691)
      character*5   SpType     (nr__1) ! MK spectral type
      integer*4     Mass       (nr__1) ! (Msun) Stellar mass (1)
      real*4        Rad        (nr__1) ! (Rsun) Stellar radius (1)
      real*4        Lum        (nr__1) ! (10+5Lsun) Stellar luminosity (1)
      integer*4     vsini      (nr__1) ! (km/s) Rotational velocity
      character*1   r_vsini    (nr__1) ! [bc] Vrot reference (2)
      integer*4     Vinf       (nr__1) ! (km/s) Terminal velocity V_{infty}_
      character*1   r_Vinf     (nr__1) ! [ed] Velocity reference (2)
      real*4        Mdot       (nr__1) ! (10-5Msun/yr) Mass-loss rate
      character*1   r_Mdot     (nr__1) ! [fd] M reference (2)
      real*4        Pmin       (nr__1) ! (d) Minimum rotation period
      real*4        Pmax       (nr__1) ! (d) Maximum rotation period
*Note (1): From Martins et al. 2005A&A...436.1049M
*Note (2): Reference flag as follows:
*   b = Penny L. 1996ApJ...463..737P;
*   c = Conti & Ebbets 1977ApJ...213..438C;
*   d = Markova et al. 2005A&A...440.1133M;
*   e = Leitherer C. 1988ApJ...326..356L;
*   f = Herrero et al. 2000A&A...354..193H

C=============================================================================

C  Loading file 'table2.dat'	! Journal of the observations for HD 14947,
*                             HD 15570 and HD 16691

C  Format for file interpretation

    1 format(3X,I5,1X,I2,1X,F8.3,1X,A10,1X,A4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,210
        read(1,'(A36)')ar__
        read(ar__,1)HD(i__),Seq(i__),HJD(i__),Date(i__),Dom(i__)
        if(ar__(13:20) .EQ. '') HJD(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)HD(i__),Seq(i__),HJD(i__),Date(i__),Dom(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3.dat'	! Stellar parameters of HD 14947, HD 15570, and
*                             HD 16691, and estimated minimum and maximum
*                             rotation periods

C  Format for file interpretation

    2 format(
     +  3X,I5,1X,A5,1X,I2,1X,F4.1,1X,F3.1,1X,I3,A1,1X,I4,A1,1X,F4.2,
     +  A1,1X,F3.1,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,3
        read(1,'(A51)')ar__1
        read(ar__1,2)
     +  HD_1(i__),SpType(i__),Mass(i__),Rad(i__),Lum(i__),vsini(i__),
     +  r_vsini(i__),Vinf(i__),r_Vinf(i__),Mdot(i__),r_Mdot(i__),
     +  Pmin(i__),Pmax(i__)
c    ..............Just test output...........
        write(6,2)
     +  HD_1(i__),SpType(i__),Mass(i__),Rad(i__),Lum(i__),vsini(i__),
     +  r_vsini(i__),Vinf(i__),r_Vinf(i__),Mdot(i__),r_Mdot(i__),
     +  Pmin(i__),Pmax(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end