FORTRAN Generation
(/./ftp/cats/J/ApJ/636/1098)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/636/1098 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-20
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/636/1098      Debris disks around solar-type stars         (Bryden+, 2006)
*================================================================================
*Frequency of debris disks around solar-type stars: first results from a Spitzer
*MIPS survey.
*    Bryden G., Beichman C.A., Trilling D.E., Rieke G.H., Holmes E.K.,
*    Lawler S.M., Stapelfeldt K.R., Werner M.W., Gautier T.N., Blaylock M.,
*    Gordon K.D., Stansberry J.A., Su K.Y.L.
*   <Astrophys. J., 636, 1098-1113 (2006)>
*   =2006ApJ...636.1098B
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! FGK survey stars

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

      integer*4     HD          ! The HD number (Cat. III/135)
      character*3   n_HD        ! [ab] a = Observed, b = Known planet-bearing star.
      character*8   SpType      ! SIMBAD MK spectral type
      real*4        Vmag        ! (mag) SIMBAD visual magnitude (2)
      real*4        Age         ! (Gyr) ? The Wright et al. (2004, Cat. J/ApJS/152/261)
*                                 age (3)
      real*4        AgeMin      ! (Gyr) ? Minimum age
      real*4        AgeMax      ! (Gyr) ? Maximum age
      character*14  r_Age       ! Reference(s) for Age
      real*4        v_Fe_H_     ! ([Sun]) ? The average [Fe/H] abundance
      real*4        e__Fe_H_    ! ([Sun]) ? Dispersion in Abund
      integer*4     o__Fe_H_    ! Number of estimates in Abund
      character*21  Refs        ! Additional reference(s)
*Note (2): Typically from the Hipparcos satellite.
*Note (3): Or an average of other estimates if Wright data is unavailable.

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

C  Declarations for 'refs.dat'	! References

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

      character*2   Ref         ! Reference code
      character*19  Bibcode     ! Bibcode
      character*24  Aut         ! Author's name
      character*65  Com         ! Comment

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

C  Declarations for 'table2.dat'	! Measured and predicted fluxes at 24 and 70{mu}m

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

      integer*4     HD_1        ! The HD number
      character*1   n_HD_1      ! [de] Note on IR excess (1)
      real*4        MIPS24um    ! (mJy) Measured MIPS 24{mu}m flux density
      real*4        Star24um    ! (mJy) Predicted 24{mu}m flux density
      real*4        Rat24um     ! The MIPS/predicted 24{mu}m flux density ratio
      real*4        MIPS70um    ! (mJy) Measured MIPS 70{mu}m flux density
      real*4        e_MIPS70um  ! (mJy) Uncertainty in MIPS70um
      real*4        Star70um    ! (mJy) Predicted 70{mu}m flux density
      real*4        Rat70um     ! The MIPS/predicted 70{mu}m flux density ratio
      real*4        S_N         ! The 70{mu}m signal-to-noise ratio
      real*4        chi70       ! Significance of excess (Eq. 1)
      real*4        FDust       ! (mJy) ? Dust flux density at 70{mu}m (2)
      character*1   l_LDust     ! Limit flag on LDust
      real*4        LDust       ! Dust luminosity from 70{mu}m emission (3)
*Note (1): Flag on HD, defined as follows:
*      d = star with excess 70{mu}m emission;
*      e = star with excess 24{mu}m emission.
*Note (2): 70{mu}m dust fluxes have been color corrected by 15%, appropriate
*          for ~50K emission.
*Note (3): See Eq. 3. In units of the stellar luminosity

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

C  Loading file 'table1.dat'	! FGK survey stars

C  Format for file interpretation

    1 format(
     +  3X,I6,1X,A3,1X,A8,1X,F4.2,1X,F4.1,1X,F4.2,1X,F4.1,1X,A14,1X,
     +  F5.2,1X,F4.2,1X,I1,1X,A21)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,127
        read(1,'(A92)')ar__
        read(ar__,1)
     +  HD,n_HD,SpType,Vmag,Age,AgeMin,AgeMax,r_Age,v_Fe_H_,e__Fe_H_,
     +  o__Fe_H_,Refs
        if(ar__(29:32) .EQ. '') Age = rNULL__
        if(ar__(34:37) .EQ. '') AgeMin = rNULL__
        if(ar__(39:42) .EQ. '') AgeMax = rNULL__
        if(ar__(59:63) .EQ. '') v_Fe_H_ = rNULL__
        if(ar__(65:68) .EQ. '') e__Fe_H_ = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  HD,n_HD,SpType,Vmag,Age,AgeMin,AgeMax,r_Age,v_Fe_H_,e__Fe_H_,
     +  o__Fe_H_,Refs
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    2 format(A2,1X,A19,1X,A24,1X,A65)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,31
        read(1,'(A113)')ar__1
        read(ar__1,2)Ref,Bibcode,Aut,Com
c    ..............Just test output...........
        write(6,2)Ref,Bibcode,Aut,Com
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Measured and predicted fluxes at 24 and 70{mu}m

C  Format for file interpretation

    3 format(
     +  I6,1X,A1,1X,F6.1,1X,F6.1,1X,F4.2,1X,F5.1,1X,F4.1,1X,F5.1,1X,
     +  F3.1,1X,F4.1,1X,F4.1,1X,F5.1,1X,A1,E6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,69
        read(1,'(A72)')ar__2
        read(ar__2,3)
     +  HD_1,n_HD_1,MIPS24um,Star24um,Rat24um,MIPS70um,e_MIPS70um,
     +  Star70um,Rat70um,S_N,chi70,FDust,l_LDust,LDust
        if(ar__2(60:64) .EQ. '') FDust = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  HD_1,n_HD_1,MIPS24um,Star24um,Rat24um,MIPS70um,e_MIPS70um,
     +  Star70um,Rat70um,S_N,chi70,FDust,l_LDust,LDust
c    .......End.of.Just test output...........
      end do
      close(1)

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