FORTRAN Generation
(/./ftp/cats/J/ApJ/831/125)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/831/125 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-18
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/831/125  ALMA 887{mu}m obs. of ChaI star-forming region  (Pascucci+, 2016)
*================================================================================
*A steeper than linear disk mass-stellar mass scaling relation.
*    Pascucci I., Testi L., Herczeg G.J., Long F., Manara C.F., Hendler N.,
*    Mulders G.D., Krijt S., Ciesla F., Henning T., Mohanty S.,
*    Drabek-Maunder E., Apai D., Szucs L., Sacco G., Olofsson J.
*   <Astrophys. J., 831, 125-125 (2016)>
*   =2016ApJ...831..125P    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Source properties

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

      character*17  v2MASS     (nr__) ! 2MASS identifier (JHHMMSSss+DDMMSSs)
      character*1   f_2MASS    (nr__) ! [*] Flag on 2MASS (1)
      character*10  Name       (nr__) ! Other identifier
      real*4        Multi      (nr__) ! (arcsec) [0/29]?=0 Multiplicity information from lit.
      character*4   r_Multi    (nr__) ! Reference for Multi (2)
      character*5   SpT        (nr__) ! Luhman+ 2007, J/ApJS/173/104;
*                                     2008, J/ApJ/675/1375 MK spectral type
      character*4   SED        (nr__) ! SED type ("FS" (flat spectrum),
*                                     Class "II" (82 instances), or "II/T")
      character*4   ALMA       (nr__) ! ALMA sample type ("Cool" or "Hot")
      character*5   ASpT       (nr__) ! Adopted MK spectral type
      character*4   r_ASpT     (nr__) ! Reference for ASpType (2)
      real*4        logM_      (nr__) ! ([Lsun]) [-1.6/0.3] Log stellar mass
      real*4        b_logM_    (nr__) ! ([Lsun]) [-1.8/0.3]? Lower boundary in logM* (3)
      real*4        B_logM__1  (nr__) ! ([Lsun]) [-1.1/0.6]? Upper boundary in logM* (3)
*Note (1):
*     * =  T30 is the secondary of T31 at a separation of 16.52". ISO 138 is
*          the secondary of ISO 143 at 18.16". T53 is the secondary of T52
*          at 11.18". Hn18 is the secondary of CHXR60 (not included in our
*          ALMA survey) at a separation of 28.28".
*Note (2): Reference as follows:
*   A15 = Anthonioz et al. (2015A&A...574A..41A);
*   D13 = Daemgen et al. (2013A&A...554A..43D);
*  KH07 = Kraus and Hillenbrand 2007, J/ApJ/662/413;
*  La08 = Lafreniere et al. (2008ApJ...683..844L);
*   L07 = Luhman 2007, J/ApJS/173/104;
*   M14 = Manara et al. (2014A&A...568A..18M);
*  M16a = Manara et al. (2016A&A...585A.136M);
*  M16b = Manara et al. (2016A&A...591L...3M);
*   N12 = Nguyen et al. 2012, J/ApJ/745/119.
*Note (3): When we fixed the isochrone there are no uncertainties associated
*          with the estimated stellar mass, see Section 2.1.

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

C  Declarations for 'table3.dat'	! Measured continuum flux densities

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

      character*17  v2MASS_1   (nr__1) ! 2MASS identifier (JHHMMSSss+DDMMSSs)
      character*1   n_2MASS    (nr__1) ! [d*] Flag on 2MASS (4)
      real*8        Fnu        (nr__1) ! (mJy) [-0.7/1364] Integrated flux density at 887um
      real*4        e_Fnu      (nr__1) ! (mJy) [0.1/2] Uncertainty in Fnu
      real*8        DelRA      (nr__1) ! (arcsec) [-0.6/1.1] Offsets from the phase center
*                                   in RA (5)
      real*4        e_DelRA    (nr__1) ! (arcsec) [0/0.2] Uncertainty in DelRA (5)
      real*8        DelDE      (nr__1) ! (arcsec) [-0.4/0.5] Offsets from the phase center
*                                   in DE (5)
      real*4        e_DelDE    (nr__1) ! (arcsec) [0/0.2] Uncertainty in DelDE (5)
      real*4        MajAxis    (nr__1) ! (arcsec) [0/0.8] FWHM of major axis (5)
      real*4        MinAxis    (nr__1) ! (arcsec) [0/0.5] FWHM of minor axis (5)
*Note (4): Flag as follows:
*    d = Sources with rings. Integrated flux density measured on image
*        within the 3{sigma} contour.
*    * = Sources that have additional mm detections in their exposures:
*        J11044258-7741571 (ISO 52) at ~6", coordinates
*        (11:04:40.59;-77:41:56.9); J11082238-7730277 (ISO 143) at ~10",
*        coordinates (11:08:21.11;-77:30:18.9); and J11123092-7644241 (T53)
*        at ~11", coordinates (11:12:27.7;-76:44:22.3). In the first two cases
*        there is no object in the SIMBAD Astronomical Database associated with
*        the mm emission. In the case of T53 we detect the disk from the
*        companion T52. Fluxes from these additional detections are not
*        reported in the table.
*Note (5): Sources with a FWHM reported in the last column of the table are
*          those that were fitted with an elliptical gaussian. Undetected
*          sources have values of 0.0 in all columns following the Fnu columns.
*          For these sources flux densities are measured assuming a point source
*          model and fixed DelRA and DelDE to the median values of the detected
*          sources.

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

C  Declarations for 'figure8.dat'	! The M_dust_-M_*_ relation

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

      integer*4     Loc        (nr__2) ! [1/4] Location code (1)
      character*17  ID         (nr__2) ! Source identifier
      integer*4     Det_       (nr__2) ! Detection code; 1=yes, 0=3{sigma} upper limit
      real*4        logM__1    (nr__2) ! ([Msun]) [-1.8/0.6] Log stellar mass
      real*4        logdust    (nr__2) ! ([Mgeo]) [-0.9/2.5] Log dust mass in Earth units
      real*4        logd20     (nr__2) ! ([Mgeo]) [-1.1/2.7] Log dust mass in Earth units (2)
*Note (1): Location code as follows:
*    1 = Taurus;
*    2 = Lupus;
*    3 = ChamaeleonI;
*    4 = Upper Sco.
*Note (2): Assuming a fixed dust temperature of 20K.

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

C  Loading file 'table1.dat'	! Source properties

C  Format for file interpretation

    1 format(
     +  A17,A1,1X,A10,1X,F6.3,1X,A4,1X,A5,1X,A4,1X,A4,1X,A5,1X,A4,1X,
     +  F5.2,1X,F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,93
        read(1,'(A86)')ar__
        read(ar__,1)
     +  v2MASS(i__),f_2MASS(i__),Name(i__),Multi(i__),r_Multi(i__),
     +  SpT(i__),SED(i__),ALMA(i__),ASpT(i__),r_ASpT(i__),logM_(i__),
     +  b_logM_(i__),B_logM__1(i__)
        if(ar__(76:80) .EQ. '') b_logM_(i__) = rNULL__
        if(ar__(82:86) .EQ. '') B_logM__1(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  v2MASS(i__),f_2MASS(i__),Name(i__),Multi(i__),r_Multi(i__),
     +  SpT(i__),SED(i__),ALMA(i__),ASpT(i__),r_ASpT(i__),logM_(i__),
     +  b_logM_(i__),B_logM__1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Measured continuum flux densities

C  Format for file interpretation

    2 format(
     +  A17,A1,1X,F7.2,1X,F4.2,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,
     +  F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,93
        read(1,'(A71)')ar__1
        read(ar__1,2)
     +  v2MASS_1(i__),n_2MASS(i__),Fnu(i__),e_Fnu(i__),DelRA(i__),
     +  e_DelRA(i__),DelDE(i__),e_DelDE(i__),MajAxis(i__),MinAxis(i__)
c    ..............Just test output...........
        write(6,2)
     +  v2MASS_1(i__),n_2MASS(i__),Fnu(i__),e_Fnu(i__),DelRA(i__),
     +  e_DelRA(i__),DelDE(i__),e_DelDE(i__),MajAxis(i__),MinAxis(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'figure8.dat'	! The M_dust_-M_*_ relation

C  Format for file interpretation

    3 format(I1,1X,A17,1X,I1,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'figure8.dat')
      write(6,*) '....Loading file: figure8.dat'
      do i__=1,446
        read(1,'(A42)')ar__2
        read(ar__2,3)
     +  Loc(i__),ID(i__),Det_(i__),logM__1(i__),logdust(i__),
     +  logd20(i__)
c    ..............Just test output...........
        write(6,3)
     +  Loc(i__),ID(i__),Det_(i__),logM__1(i__),logdust(i__),
     +  logd20(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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