FORTRAN Generation
(/./ftp/cats/J/ApJ/847/31)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/847/31 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-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/847/31    Protoplanetary disk data in Cha I and Lupus    (Mulders+, 2017)
*================================================================================
*Constraints from dust mass and mass accretion rate measurements on angular
*momentum transport in protoplanetary disks.
*    Mulders G.D., Pascucci I., Manara C.F., Testi L., Herczeg G.J., Henning T.,
*    Mohanty S., Lodato G.
*   <Astrophys. J., 847, 31 (2017)>
*   =2017ApJ...847...31M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Stellar and disk properties for Chamaeleon I and Lupus

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

      character*23  ID         (nr__) ! Identifier, 2MASS J or other star name
      character*5   SpT        (nr__) ! MK spectral type
      integer*4     Teff       (nr__) ! (K) [0/5110]? Stellar effective temperature
      real*4        logL       (nr__) ! ([Lsun]) [-2.7/0.8]? log of stellar luminosity
      real*4        logR       (nr__) ! ([Rsun]) [-0.9/0.6]? log of stellar radius
      real*4        logM       (nr__) ! ([Msun]) [-1.6/0.3]? log of stellar mass
      real*4        e_logM     (nr__) ! ([Msun]) [0.04/0.4]? Uncertainty in logM
      real*4        logLacc    (nr__) ! ([Lsun]) [-5.3/0.7]? Log of accretion luminosity
      real*4        logMdust   (nr__) ! ([Mgeo]) [-0.9/2.7]? Log of dust mass
      real*4        e_logMdust (nr__) ! ([Mgeo]) [0.0003/0.2]? Uncertainty in logMdust
      character*5   Det        (nr__) ! [True /False] Source is detected in dust
*                                     continuum emission (>3{sigma})
      real*8        logMacc    (nr__) ! ([Msun/yr]) [-12.3/-6.6]? log of mass accretion rate
      real*4        e_logMacc  (nr__) ! ([Msun/yr]) [0.2/0.5] Uncertainty in logMacc
      character*5   Det2       (nr__) ! [True /False] Accretion luminosity
*                                     > expected chromospheric emission
      character*6   Reg        (nr__) ! Region: Chamaeleon I (ChaI) or Lupus
      character*9   Exc        (nr__) ! Exclude? (1)
*Note (1): Reason why sources are not included in analysis as follows:
*    underlum  = sub-luminous, lying below the ZAMS on HR diagram or
*                 known edge-on disk (9 occurrences);
*    binary    = apparent binary (1 occurrence);
*    faint     = source not detected with VLT/X-shooter (4 occurrences);
*    noMdust   = no dust mass measurement (14 occurrences);
*    notinslit = source not in VLT/X-shooter slit (1 occurrence);

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

C  Loading file 'table1.dat'	! Stellar and disk properties for Chamaeleon I and Lupus

C  Format for file interpretation

    1 format(
     +  A23,1X,A5,1X,I4,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F6.2,1X,
     +  F5.2,1X,F6.4,1X,A5,1X,F7.2,1X,F4.2,1X,A5,1X,A6,1X,A9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,174
        read(1,'(A119)')ar__
        read(ar__,1)
     +  ID(i__),SpT(i__),Teff(i__),logL(i__),logR(i__),logM(i__),
     +  e_logM(i__),logLacc(i__),logMdust(i__),e_logMdust(i__),
     +  Det(i__),logMacc(i__),e_logMacc(i__),Det2(i__),Reg(i__),
     +  Exc(i__)
        if(ar__(31:34) .EQ. '') Teff(i__) = iNULL__
        if(ar__(36:40) .EQ. '') logL(i__) = rNULL__
        if(ar__(42:46) .EQ. '') logR(i__) = rNULL__
        if(ar__(48:52) .EQ. '') logM(i__) = rNULL__
        if(ar__(54:57) .EQ. '') e_logM(i__) = rNULL__
        if(ar__(59:64) .EQ. '') logLacc(i__) = rNULL__
        if(ar__(66:70) .EQ. '') logMdust(i__) = rNULL__
        if(ar__(72:77) .EQ. '') e_logMdust(i__) = rNULL__
        if(ar__(85:91) .EQ. '') logMacc(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),SpT(i__),Teff(i__),logL(i__),logR(i__),logM(i__),
     +  e_logM(i__),logLacc(i__),logMdust(i__),e_logMdust(i__),
     +  Det(i__),logMacc(i__),e_logMacc(i__),Det2(i__),Reg(i__),
     +  Exc(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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