FORTRAN Generation
(/./ftp/cats/J/ApJ/648/580)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/648/580 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-29
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/648/580       Rotational velocities in 19 open clusters     (Huang+, 2006)
*================================================================================
*Stellar rotation in young clusters.
*I. Evolution of projected rotational velocity distributions.
*    Huang W., Gies D.R.
*   <Astrophys. J., 648, 580-590 (2006)>
*   =2006ApJ...648..580H
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! ZAMS model star parameters

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

      integer*4     Index      (nr__) ! Index
      character*6   PSpType    (nr__) ! Pseudospectral SubType spectral type
      real*4        Teff       (nr__) ! (kK) Effective temperature
      real*4        M          (nr__) ! (solMass) Mass
      real*4        R          (nr__) ! (solRad) Radius

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

C  Declarations for 'table3.dat'	! Projected Rotational Velocities of B Stars
                             in 19 Open Clusters

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

      character*11  Cluster    (nr__1) ! Cluster name
      integer*4     WEBDA      (nr__1) ! Webda ID (http://obswww.unige.ch/web) (1)
      integer*4     vsini      (nr__1) ! (km/s) Derived projected rotational velocity
      integer*4     e_vsini    (nr__1) ! (km/s) Measurement error in vsini
      real*4        RV1        (nr__1) ! (km/s) Radial velocity from first night,
*                                 HJD specified in obs.dat file
      real*4        RV2        (nr__1) ! (km/s) ? Radial velocity from second night if
*                                  available, HJD specified in obs.dat file
      real*4        RV3        (nr__1) ! (km/s) ? Radial velocity from third night if available,
*                                  HJD specified in obs.dat file
      real*4        Index_1    (nr__1) ! Index of estimated pseudo-spectral type,
*                                 detailed in table2.dat
      character*95  Notes      (nr__1) ! Additional notes
*Note (1): For stars not in the WEBDA database, we assign a number > 10000
*     and give an identification and coordinates for epoch 2000.0 in the
*     "Notes" column; the other stars are listed using the Webda numbering
*     systems from references detailed in refs.dat file.

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

C  Declarations for 'obs.dat'	! The heliocentric Julian date of observation
                             for each night

      integer*4 nr__2
      parameter (nr__2=19)	! Number of records
      character*49 ar__2  	! Full-size record

      character*11  Cluster_1  (nr__2) ! Cluster name
      real*8        HJD1       (nr__2) ! (d) Heliocentric Julian date of night 1
      real*8        HJD2       (nr__2) ! (d) ? Heliocentric Julian date of night 2
      real*8        HJD3       (nr__2) ! (d) ? Heliocentric Julian date of night 3

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

C  Declarations for 'refs.dat'	! Webda numbering systems references

      integer*4 nr__3
      parameter (nr__3=19)	! Number of records
      character*93 ar__3  	! Full-size record

      character*11  Cluster_2  (nr__3) ! Cluster name
      character*19  BibCode    (nr__3) ! BibCode
      character*18  Aut        (nr__3) ! Author's name
      character*20  Com        (nr__3) ! Comments
      character*21  Acronym    (nr__3) ! Acronym used in Simbad name for stars in this
*                                   cluster with number <10000

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

C  Declarations for 'table4.dat'	! Summary data for 19 open clusters

      integer*4 nr__4
      parameter (nr__4=19)	! Number of records
      character*51 ar__4  	! Full-size record

      character*11  Cluster_3  (nr__4) ! Cluster name
      real*4        logAge     (nr__4) ! ([yr]) Age
      integer*4     N          (nr__4) ! Number of stars
      integer*4     v_vsini_   (nr__4) ! (km/s) Mean rotational velocity from all stars
      integer*4     NB3        (nr__4) ! Number of stars with SpType<=B3
      integer*4     v_vsini_B3 (nr__4) ! (km/s) Mean rotational velocity from stars with
*                                   SpType<=B3
      integer*4     N300       (nr__4) ! Number of stars with vsini>300km/s
      integer*4     N300p      (nr__4) ! (%) Number of stars with vsini>300km/s
*                                   in percentage
      integer*4     Nbin       (nr__4) ! ? Number of binaries
      integer*4     Nbinp      (nr__4) ! (%) ? Number of binaries in percentage

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

C  Loading file 'table2.dat'	! ZAMS model star parameters

C  Format for file interpretation

    1 format(I2,1X,A6,2X,F4.1,1X,F4.1,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,21
        read(1,'(A24)')ar__
        read(ar__,1)Index(i__),PSpType(i__),Teff(i__),M(i__),R(i__)
c    ..............Just test output...........
        write(6,1)Index(i__),PSpType(i__),Teff(i__),M(i__),R(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Projected Rotational Velocities of B Stars
*                             in 19 Open Clusters

C  Format for file interpretation

    2 format(
     +  A11,2X,I5,1X,I3,1X,I2,1X,F6.1,1X,F6.1,1X,F6.1,1X,F4.1,1X,A95)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,496
        read(1,'(A147)')ar__1
        read(ar__1,2)
     +  Cluster(i__),WEBDA(i__),vsini(i__),e_vsini(i__),RV1(i__),
     +  RV2(i__),RV3(i__),Index_1(i__),Notes(i__)
        if(ar__1(34:39) .EQ. '') RV2(i__) = rNULL__
        if(ar__1(41:46) .EQ. '') RV3(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Cluster(i__),WEBDA(i__),vsini(i__),e_vsini(i__),RV1(i__),
     +  RV2(i__),RV3(i__),Index_1(i__),Notes(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'obs.dat'	! The heliocentric Julian date of observation
*                             for each night

C  Format for file interpretation

    3 format(A11,1X,F11.3,2X,F11.3,2X,F11.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'obs.dat')
      write(6,*) '....Loading file: obs.dat'
      do i__=1,19
        read(1,'(A49)')ar__2
        read(ar__2,3)Cluster_1(i__),HJD1(i__),HJD2(i__),HJD3(i__)
        if(ar__2(26:36) .EQ. '') HJD2(i__) = rNULL__
        if(ar__2(39:49) .EQ. '') HJD3(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)Cluster_1(i__),HJD1(i__),HJD2(i__),HJD3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! Webda numbering systems references

C  Format for file interpretation

    4 format(A11,1X,A19,1X,A18,1X,A20,1X,A21)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,19
        read(1,'(A93)')ar__3
        read(ar__3,4)
     +  Cluster_2(i__),BibCode(i__),Aut(i__),Com(i__),Acronym(i__)
c    ..............Just test output...........
        write(6,4)
     +  Cluster_2(i__),BibCode(i__),Aut(i__),Com(i__),Acronym(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Summary data for 19 open clusters

C  Format for file interpretation

    5 format(
     +  A11,2X,F4.2,1X,I2,1X,I3,1X,I2,1X,I3,1X,I1,2X,I2,2X,I2,2X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,19
        read(1,'(A51)')ar__4
        read(ar__4,5)
     +  Cluster_3(i__),logAge(i__),N(i__),v_vsini_(i__),NB3(i__),
     +  v_vsini_B3(i__),N300(i__),N300p(i__),Nbin(i__),Nbinp(i__)
        if(ar__4(40:41) .EQ. '') Nbin(i__) = iNULL__
        if(ar__4(44:45) .EQ. '') Nbinp(i__) = iNULL__
c    ..............Just test output...........
        write(6,5)
     +  Cluster_3(i__),logAge(i__),N(i__),v_vsini_(i__),NB3(i__),
     +  v_vsini_B3(i__),N300(i__),N300p(i__),Nbin(i__),Nbinp(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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