FORTRAN Generation
(/./ftp/cats/J/A_A/527/A145)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/527/A145 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-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/A+A/527/A145      Chamaeleon I 870um sources               (Belloche+, 2011)
*================================================================================
*The end of star formation in Chamaeleon I ?
*A LABOCA census of starless and protostellar cores.
*   Belloche A., Schuller F., Parise B., Andre P., Hatchell J., Jorgensen J.K.,
*   Bontemps S., Weiss A., Menten , K.M., Muders D.
*  <Astron. Astrophys. 527, A145 (2011)>
*  =2011A&A...527A.145B
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Sources extracted with Gaussclumps in the 870um
                             continuum sum map of Cha I at scale 5, and
                             possible associations found in the SIMBAD database

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     Seq         ! [1/84] Source sequential number (1)
      integer*4     RAh         ! (h) Right ascension (J2000)
      integer*4     RAm         ! (min) Right ascension (J2000)
      real*4        RAs         ! (s) Right ascension (J2000)
      character*1   DE_         ! Declination sign (J2000)
      integer*4     DEd         ! (deg) Declination (J2000)
      integer*4     DEm         ! (arcmin) Declination (J2000)
      real*4        DEs         ! (arcsec) Declination (J2000)
      real*4        Fp          ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam)
      real*4        Ftot        ! (Jy) Total flux density at 870um
      real*4        amaj        ! (arcsec) Major axis (FWHM)
      real*4        amin        ! (arcsec) Minor axis (FWHM)
      real*4        PA          ! (deg) [-90/90] Position angle (east from north)
      real*4        Size        ! (arcsec) ? Mean source size
*                                    (geometrical mean of amaj and amin)
      character*2   Type        ! Type of source based on the possible
*                                   associations found in the SIMBAD database (2)
      character*23  SIMBAD      ! Source found in the SIMBAD database within the
*                                   FWHM ellipse
      character*3   Stype       ! SIMBAD type of the SIMBAD source
      real*4        Dist        ! (arcsec) ? Distance of the possibly associated SIMBAD
*                                    source to the fitted peak position of the
*                                    LABOCA source
*Note (1): Sequential number 1-84 in table2, and
*          Additional sequential number 1-5 in table5.
*Note (2): Type of source based on the associations found in the SIMBAD
*          database as follows:
*  C  = starless core or Class 0 protostar
*  S  = young stellar object (Class I or more evolved)
*  Sc = candidate association with a YSO
*  R  = residual of a stronger nearby source
*  A  = artefact

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

C  Declarations for 'table5.dat'	! Additional compact sources with SIMBAD association
                             with a formal peak flux density above 3.5{sigma}
                             in the 870um continuum map of Cha I filtered up
                             to scale 5

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1     ! (deg) Right Ascension J2000
      real*8        DEdeg_1     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     Seq_1       ! [1/84] Source sequential number (1)
      integer*4     RAh_1       ! (h) Right ascension (J2000)
      integer*4     RAm_1       ! (min) Right ascension (J2000)
      real*4        RAs_1       ! (s) Right ascension (J2000)
      character*1   DE__1       ! Declination sign (J2000)
      integer*4     DEd_1       ! (deg) Declination (J2000)
      integer*4     DEm_1       ! (arcmin) Declination (J2000)
      real*4        DEs_1       ! (arcsec) Declination (J2000)
      real*4        Fp_1        ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam)
      real*4        Ftot_1      ! (Jy) Total flux density at 870um
      real*4        amaj_1      ! (arcsec) Major axis (FWHM)
      real*4        amin_1      ! (arcsec) Minor axis (FWHM)
      real*4        PA_1        ! (deg) [-90/90] Position angle (east from north)
      real*4        Size_1      ! (arcsec) ? Mean source size
*                                    (geometrical mean of amaj and amin)
      character*2   Type_1      ! Type of source based on the possible
*                                   associations found in the SIMBAD database (2)
      character*23  SIMBAD_1    ! Source found in the SIMBAD database within the
*                                   FWHM ellipse
      character*3   Stype_1     ! SIMBAD type of the SIMBAD source
      real*4        Dist_1      ! (arcsec) ? Distance of the possibly associated SIMBAD
*                                    source to the fitted peak position of the
*                                    LABOCA source
*Note (1): Sequential number 1-84 in table2, and
*          Additional sequential number 1-5 in table5.
*Note (2): Type of source based on the associations found in the SIMBAD
*          database as follows:
*  C  = starless core or Class 0 protostar
*  S  = young stellar object (Class I or more evolved)
*  Sc = candidate association with a YSO
*  R  = residual of a stronger nearby source
*  A  = artefact

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

C  Declarations for 'table6.dat'	! Characteristics of starless (or Class 0) sources
                             extracted with Gaussclumps in the 870um
                             continuum map of Cha I filtered up to scale 5

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

      character*8   Name        ! Name, Cha1-CNN
      integer*4     Seq_2       ! [1/84] Numbering of Gaussclumps sources like
*                                     in Table 2
      real*4        FWHM1       ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHM2       ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHMa1      ! (kAU) Size of the fitted Gaussian (1)
      real*4        FWHMa2      ! (kAU) Size of the fitted Gaussian (1)
      real*4        Ra          ! Aspect ratio (1)
      real*4        Npeak       ! (10+21cm-2) Peak H2 column density (2)
      real*4        AV          ! (mag) Visual extinction derived from 2MASS
      real*4        Mpeak       ! (Msun) Mass in the central beam (HPBW = 21.2") (3)
      real*4        Mtot        ! (Msun) Total mass derived from the Gaussian fit (3)
      real*4        M50_        ! (Msun) Mass computed from the flux measured in an
*                                     aperture of 50" in diameter (3)
      integer*4     CM          ! (%) Mass concentration Mpeak/M50"
      integer*4     e_CM        ! (%) Statistical rms uncertainty on CM
      real*4        npeak_1     ! (10+5cm-3) Beam-averaged free-particle density within
*                                     the central beam (4)
      real*4        n50_        ! (10+5cm-3) Mean free-particle density computed for the
*                                     mass M50" in the aperture of diameter
*                                     50" (4)
      real*4        cn          ! Density contrast npeak/n50"
      real*4        e_cn        ! Statistical rms uncertainty on cn
*Note (1): Sizes smaller than 25.4" were set to 25.4" to compute the deconvolved
*  sizes, in order to account for a fit inaccuracy corresponding to a 5{sigma}
*  detection in peak flux density. As a result, the minimum size that can be
*  measured is about 2100AU. The aspect ratio is the ratio of the deconvolved
*  sizes along the major and minor axes.
*Note (2): Peak H_2_ column density computed assuming a dust opacity of
*   0.01cm^2^/g. The statistical rms uncertainty is 1.1x10^21^cm^-2^.
*Note (3): The statistical rms uncertainties of Mpeak and M50" are 0.006 and
*   0.011M_{sun}_, respectively.
*Note (4): The statistical rms uncertainties of npeak and n50" are
*   5.6x10^4^ and 7.2x10^3^cm^-3^, respectively.

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

C  Declarations for 'table8.dat'	! Characteristics of YSOs extracted with Gaussclumps
                             in the 870um continuum map of Cha I filtered
                             up to scale 5

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

      character*8   Name_1      ! Name, Cha1-SNN
      integer*4     Seq_3       ! [1/87] Numbering of Gaussclumps sources like
*                                       in table2 (table8) or table5 (table9)
      real*4        FWHM1_1     ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHM2_1     ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHMa1_1    ! (kAU) Deconvolved source size (1)
      real*4        FWHMa2_1    ! (kAU) Deconvolved source size (1)
      real*4        SNR         ! Peak signal-to-noise ratio of the fitted
*                                     Gaussian
      real*4        Npeak_2     ! (10+21cm-2) Peak H_2_ column density (2)
      real*4        AV_1        ! (mag) Visual extinction derived from 2MASS
      real*4        Mpeak_1     ! (Msun) Mass in the central beam (HPBW=21.2") (3)
      real*4        Mtot_1      ! (Msun) Total mass derived from the Gaussian fit (3)
      real*4        M50__1      ! (Msun) Mass computed from the flux measured in an
*                                     aperture of 50" in diameter (3)
      integer*4     CM_1        ! (%) Mass concentration Mpeak/M50"
      integer*4     e_CM_1      ! (%) Statistical rms uncertainty on CM
      character*4   Class       ! Infrared class of the nearest YSO listed in
*                                     column SIMBAD of table 2 (table8) or
*                                     table5 (table9)
      integer*4     Ref         ! Reference (4)
*Note (1): The deconvolved size is computed from the fitted size multiplied by
*   (1+1/SNR), with SNR the peak signal-to-noise ratio, and therefore is an
*   upper limit in most cases. A zero value means that the uncertainty was
*   still underestimated.
*Note (2): Peak H_2_ column density computed assuming a dust temperature of 20K
*  and a dust opacity of 0.03cm^2^/g. The statistical rms uncertainty is
*  1.5x10^20^cm^-2^.
*Note (3): The statistical rms uncertainties of Mpeak and M50" are 0.0009 and
*   0.0015M_{sun}_, respectively.
*Note (4): References as follows:
*    2 = Luhman et al. (2008, Cat. J/ApJ/675/1375)
*    3 = Furlan et al. (2009, Cat. J/ApJ/703/1964)
*    4 = Henning et al. (1993A&A...276..129H)
*    5 = Luhman & Muench (2008, Cat. J/ApJ/684/654)

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

C  Declarations for 'table9.dat'	! Characteristics of the additional compact sources
                             listed in Table 5

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

      character*8   Name_2      ! Name, Cha1-SNN
      integer*4     Seq_4       ! [1/87] Numbering of Gaussclumps sources like
*                                       in table2 (table8) or table5 (table9)
      real*4        FWHM1_2     ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHM2_2     ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHMa1_2    ! (kAU) Deconvolved source size (1)
      real*4        FWHMa2_2    ! (kAU) Deconvolved source size (1)
      real*4        SNR_1       ! Peak signal-to-noise ratio of the fitted
*                                     Gaussian
      real*4        Npeak_3     ! (10+21cm-2) Peak H_2_ column density (2)
      real*4        AV_2        ! (mag) Visual extinction derived from 2MASS
      real*4        Mpeak_2     ! (Msun) Mass in the central beam (HPBW=21.2") (3)
      real*4        Mtot_2      ! (Msun) Total mass derived from the Gaussian fit (3)
      real*4        M50__2      ! (Msun) Mass computed from the flux measured in an
*                                     aperture of 50" in diameter (3)
      integer*4     CM_2        ! (%) Mass concentration Mpeak/M50"
      integer*4     e_CM_2      ! (%) Statistical rms uncertainty on CM
      character*4   Class_1     ! Infrared class of the nearest YSO listed in
*                                     column SIMBAD of table 2 (table8) or
*                                     table5 (table9)
      integer*4     Ref_1       ! Reference (4)
*Note (1): The deconvolved size is computed from the fitted size multiplied by
*   (1+1/SNR), with SNR the peak signal-to-noise ratio, and therefore is an
*   upper limit in most cases. A zero value means that the uncertainty was
*   still underestimated.
*Note (2): Peak H_2_ column density computed assuming a dust temperature of 20K
*  and a dust opacity of 0.03cm^2^/g. The statistical rms uncertainty is
*  1.5x10^20^cm^-2^.
*Note (3): The statistical rms uncertainties of Mpeak and M50" are 0.0009 and
*   0.0015M_{sun}_, respectively.
*Note (4): References as follows:
*    2 = Luhman et al. (2008, Cat. J/ApJ/675/1375)
*    3 = Furlan et al. (2009, Cat. J/ApJ/703/1964)
*    4 = Henning et al. (1993A&A...276..129H)
*    5 = Luhman & Muench (2008, Cat. J/ApJ/684/654)

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

C  Loading file 'table2.dat'	! Sources extracted with Gaussclumps in the 870um
*                             continuum sum map of Cha I at scale 5, and
*                             possible associations found in the SIMBAD database

C  Format for file interpretation

    1 format(
     +  I2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.3,1X,F5.3,
     +  1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,A2,1X,A23,1X,A3,1X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,118
        read(1,'(A96)')ar__
        read(ar__,1)
     +  Seq,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Fp,Ftot,amaj,amin,PA,Size,
     +  Type,SIMBAD,Stype,Dist
        if(ar__(57:60) .EQ. '') Size = rNULL__
        if(ar__(93:96) .EQ. '') Dist = rNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  Seq,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Fp,Ftot,amaj,amin,PA,Size,
     +  Type,SIMBAD,Stype,Dist
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Additional compact sources with SIMBAD association
*                             with a formal peak flux density above 3.5{sigma}
*                             in the 870um continuum map of Cha I filtered up
*                             to scale 5

C  Format for file interpretation

    2 format(
     +  I2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.3,1X,F5.3,
     +  1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,A2,1X,A23,1X,A3,1X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,9
        read(1,'(A96)')ar__1
        read(ar__1,2)
     +  Seq_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Fp_1,Ftot_1,
     +  amaj_1,amin_1,PA_1,Size_1,Type_1,SIMBAD_1,Stype_1,Dist_1
        if(ar__1(57:60) .EQ. '') Size_1 = rNULL__
        if(ar__1(93:96) .EQ. '') Dist_1 = rNULL__
        RAdeg_1 = rNULL__
        DEdeg_1 = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_1 from input data
c  (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
        if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
        if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
        if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
        if(DEd_1 .GE. 0) DEdeg_1=DEd_1
        if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
        if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
        if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c    ..............Just test output...........
        write(6,2)
     +  Seq_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Fp_1,Ftot_1,
     +  amaj_1,amin_1,PA_1,Size_1,Type_1,SIMBAD_1,Stype_1,Dist_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Characteristics of starless (or Class 0) sources
*                             extracted with Gaussclumps in the 870um
*                             continuum map of Cha I filtered up to scale 5

C  Format for file interpretation

    3 format(
     +  A8,1X,I2,1X,F5.1,1X,F4.1,1X,F4.1,1X,F3.1,1X,F3.1,1X,F4.1,1X,
     +  F4.1,1X,F5.3,1X,F5.3,F5.3,1X,I2,1X,I2,1X,F4.1,1X,F4.1,1X,F3.1,
     +  1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,60
        read(1,'(A86)')ar__2
        read(ar__2,3)
     +  Name,Seq_2,FWHM1,FWHM2,FWHMa1,FWHMa2,Ra,Npeak,AV,Mpeak,Mtot,
     +  M50_,CM,e_CM,npeak_1,n50_,cn,e_cn
c    ..............Just test output...........
        write(6,3)
     +  Name,Seq_2,FWHM1,FWHM2,FWHMa1,FWHMa2,Ra,Npeak,AV,Mpeak,Mtot,
     +  M50_,CM,e_CM,npeak_1,n50_,cn,e_cn
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Characteristics of YSOs extracted with Gaussclumps
*                             in the 870um continuum map of Cha I filtered
*                             up to scale 5

C  Format for file interpretation

    4 format(
     +  A8,1X,I2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.1,1X,F5.2,1X,
     +  F4.1,1X,F6.4,1X,F6.4,1X,F6.4,1X,I3,1X,I3,1X,A4,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,16
        read(1,'(A83)')ar__3
        read(ar__3,4)
     +  Name_1,Seq_3,FWHM1_1,FWHM2_1,FWHMa1_1,FWHMa2_1,SNR,Npeak_2,
     +  AV_1,Mpeak_1,Mtot_1,M50__1,CM_1,e_CM_1,Class,Ref
c    ..............Just test output...........
        write(6,4)
     +  Name_1,Seq_3,FWHM1_1,FWHM2_1,FWHMa1_1,FWHMa2_1,SNR,Npeak_2,
     +  AV_1,Mpeak_1,Mtot_1,M50__1,CM_1,e_CM_1,Class,Ref
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Characteristics of the additional compact sources
*                             listed in Table 5

C  Format for file interpretation

    5 format(
     +  A8,1X,I2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.1,1X,F5.2,1X,
     +  F4.1,1X,F6.4,1X,F6.4,1X,F6.4,1X,I3,1X,I3,1X,A4,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,5
        read(1,'(A83)')ar__4
        read(ar__4,5)
     +  Name_2,Seq_4,FWHM1_2,FWHM2_2,FWHMa1_2,FWHMa2_2,SNR_1,Npeak_3,
     +  AV_2,Mpeak_2,Mtot_2,M50__2,CM_2,e_CM_2,Class_1,Ref_1
c    ..............Just test output...........
        write(6,5)
     +  Name_2,Seq_4,FWHM1_2,FWHM2_2,FWHMa1_2,FWHMa2_2,SNR_1,Npeak_3,
     +  AV_2,Mpeak_2,Mtot_2,M50__2,CM_2,e_CM_2,Class_1,Ref_1
c    .......End.of.Just test output...........
      end do
      close(1)

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