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 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-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      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     Seq        (nr__) ! [1/84] Source sequential number (1)
      integer*4     RAh        (nr__) ! (h) Right ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Right ascension (J2000)
      real*4        RAs        (nr__) ! (s) Right ascension (J2000)
      character*1   DE_        (nr__) ! Declination sign (J2000)
      integer*4     DEd        (nr__) ! (deg) Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Declination (J2000)
      real*4        DEs        (nr__) ! (arcsec) Declination (J2000)
      real*4        Fp         (nr__) ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam)
      real*4        Ftot       (nr__) ! (Jy) Total flux density at 870um
      real*4        amaj       (nr__) ! (arcsec) Major axis (FWHM)
      real*4        amin       (nr__) ! (arcsec) Minor axis (FWHM)
      real*4        PA         (nr__) ! (deg) [-90/90] Position angle (east from north)
      real*4        Size       (nr__) ! (arcsec) ? Mean source size
*                                    (geometrical mean of amaj and amin)
      character*2   Type       (nr__) ! Type of source based on the possible
*                                   associations found in the SIMBAD database (2)
      character*23  SIMBAD     (nr__) ! Source found in the SIMBAD database within the
*                                   FWHM ellipse
      character*3   Stype      (nr__) ! SIMBAD type of the SIMBAD source
      real*4        Dist       (nr__) ! (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    (nr__1) ! (deg) Right Ascension J2000
      real*8        DEdeg_1    (nr__1) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     Seq_1      (nr__1) ! [1/84] Source sequential number (1)
      integer*4     RAh_1      (nr__1) ! (h) Right ascension (J2000)
      integer*4     RAm_1      (nr__1) ! (min) Right ascension (J2000)
      real*4        RAs_1      (nr__1) ! (s) Right ascension (J2000)
      character*1   DE__1      (nr__1) ! Declination sign (J2000)
      integer*4     DEd_1      (nr__1) ! (deg) Declination (J2000)
      integer*4     DEm_1      (nr__1) ! (arcmin) Declination (J2000)
      real*4        DEs_1      (nr__1) ! (arcsec) Declination (J2000)
      real*4        Fp_1       (nr__1) ! (Jy) Peak flux density at 870um (in Jy/21.2"-beam)
      real*4        Ftot_1     (nr__1) ! (Jy) Total flux density at 870um
      real*4        amaj_1     (nr__1) ! (arcsec) Major axis (FWHM)
      real*4        amin_1     (nr__1) ! (arcsec) Minor axis (FWHM)
      real*4        PA_1       (nr__1) ! (deg) [-90/90] Position angle (east from north)
      real*4        Size_1     (nr__1) ! (arcsec) ? Mean source size
*                                    (geometrical mean of amaj and amin)
      character*2   Type_1     (nr__1) ! Type of source based on the possible
*                                   associations found in the SIMBAD database (2)
      character*23  SIMBAD_1   (nr__1) ! Source found in the SIMBAD database within the
*                                   FWHM ellipse
      character*3   Stype_1    (nr__1) ! SIMBAD type of the SIMBAD source
      real*4        Dist_1     (nr__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       (nr__2) ! Name, Cha1-CNN
      integer*4     Seq_2      (nr__2) ! [1/84] Numbering of Gaussclumps sources like
*                                     in Table 2
      real*4        FWHM1      (nr__2) ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHM2      (nr__2) ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHMa1     (nr__2) ! (kAU) Size of the fitted Gaussian (1)
      real*4        FWHMa2     (nr__2) ! (kAU) Size of the fitted Gaussian (1)
      real*4        Ra         (nr__2) ! Aspect ratio (1)
      real*4        Npeak      (nr__2) ! (10+21cm-2) Peak H2 column density (2)
      real*4        AV         (nr__2) ! (mag) Visual extinction derived from 2MASS
      real*4        Mpeak      (nr__2) ! (Msun) Mass in the central beam (HPBW = 21.2") (3)
      real*4        Mtot       (nr__2) ! (Msun) Total mass derived from the Gaussian fit (3)
      real*4        M50_       (nr__2) ! (Msun) Mass computed from the flux measured in an
*                                     aperture of 50" in diameter (3)
      integer*4     CM         (nr__2) ! (%) Mass concentration Mpeak/M50"
      integer*4     e_CM       (nr__2) ! (%) Statistical rms uncertainty on CM
      real*4        npeak_1    (nr__2) ! (10+5cm-3) Beam-averaged free-particle density within
*                                     the central beam (4)
      real*4        n50_       (nr__2) ! (10+5cm-3) Mean free-particle density computed for the
*                                     mass M50" in the aperture of diameter
*                                     50" (4)
      real*4        cn         (nr__2) ! Density contrast npeak/n50"
      real*4        e_cn       (nr__2) ! 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     (nr__3) ! Name, Cha1-SNN
      integer*4     Seq_3      (nr__3) ! [1/87] Numbering of Gaussclumps sources like
*                                       in table2 (table8) or table5 (table9)
      real*4        FWHM1_1    (nr__3) ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHM2_1    (nr__3) ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHMa1_1   (nr__3) ! (kAU) Deconvolved source size (1)
      real*4        FWHMa2_1   (nr__3) ! (kAU) Deconvolved source size (1)
      real*4        SNR        (nr__3) ! Peak signal-to-noise ratio of the fitted
*                                     Gaussian
      real*4        Npeak_2    (nr__3) ! (10+21cm-2) Peak H_2_ column density (2)
      real*4        AV_1       (nr__3) ! (mag) Visual extinction derived from 2MASS
      real*4        Mpeak_1    (nr__3) ! (Msun) Mass in the central beam (HPBW=21.2") (3)
      real*4        Mtot_1     (nr__3) ! (Msun) Total mass derived from the Gaussian fit (3)
      real*4        M50__1     (nr__3) ! (Msun) Mass computed from the flux measured in an
*                                     aperture of 50" in diameter (3)
      integer*4     CM_1       (nr__3) ! (%) Mass concentration Mpeak/M50"
      integer*4     e_CM_1     (nr__3) ! (%) Statistical rms uncertainty on CM
      character*4   Class      (nr__3) ! Infrared class of the nearest YSO listed in
*                                     column SIMBAD of table 2 (table8) or
*                                     table5 (table9)
      integer*4     Ref        (nr__3) ! 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     (nr__4) ! Name, Cha1-SNN
      integer*4     Seq_4      (nr__4) ! [1/87] Numbering of Gaussclumps sources like
*                                       in table2 (table8) or table5 (table9)
      real*4        FWHM1_2    (nr__4) ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHM2_2    (nr__4) ! (arcsec) Size of the fitted Gaussian (1)
      real*4        FWHMa1_2   (nr__4) ! (kAU) Deconvolved source size (1)
      real*4        FWHMa2_2   (nr__4) ! (kAU) Deconvolved source size (1)
      real*4        SNR_1      (nr__4) ! Peak signal-to-noise ratio of the fitted
*                                     Gaussian
      real*4        Npeak_3    (nr__4) ! (10+21cm-2) Peak H_2_ column density (2)
      real*4        AV_2       (nr__4) ! (mag) Visual extinction derived from 2MASS
      real*4        Mpeak_2    (nr__4) ! (Msun) Mass in the central beam (HPBW=21.2") (3)
      real*4        Mtot_2     (nr__4) ! (Msun) Total mass derived from the Gaussian fit (3)
      real*4        M50__2     (nr__4) ! (Msun) Mass computed from the flux measured in an
*                                     aperture of 50" in diameter (3)
      integer*4     CM_2       (nr__4) ! (%) Mass concentration Mpeak/M50"
      integer*4     e_CM_2     (nr__4) ! (%) Statistical rms uncertainty on CM
      character*4   Class_1    (nr__4) ! Infrared class of the nearest YSO listed in
*                                     column SIMBAD of table 2 (table8) or
*                                     table5 (table9)
      integer*4     Ref_1      (nr__4) ! 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(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),Fp(i__),Ftot(i__),amaj(i__),amin(i__),
     +  PA(i__),Size(i__),Type(i__),SIMBAD(i__),Stype(i__),Dist(i__)
        if(ar__(57:60) .EQ. '') Size(i__) = rNULL__
        if(ar__(93:96) .EQ. '') Dist(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  Seq(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),Fp(i__),Ftot(i__),amaj(i__),amin(i__),
     +  PA(i__),Size(i__),Type(i__),SIMBAD(i__),Stype(i__),Dist(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
     +  DEd_1(i__),DEm_1(i__),DEs_1(i__),Fp_1(i__),Ftot_1(i__),
     +  amaj_1(i__),amin_1(i__),PA_1(i__),Size_1(i__),Type_1(i__),
     +  SIMBAD_1(i__),Stype_1(i__),Dist_1(i__)
        if(ar__1(57:60) .EQ. '') Size_1(i__) = rNULL__
        if(ar__1(93:96) .EQ. '') Dist_1(i__) = rNULL__
        RAdeg_1(i__) = rNULL__
        DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
        if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Seq_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
     +  DEd_1(i__),DEm_1(i__),DEs_1(i__),Fp_1(i__),Ftot_1(i__),
     +  amaj_1(i__),amin_1(i__),PA_1(i__),Size_1(i__),Type_1(i__),
     +  SIMBAD_1(i__),Stype_1(i__),Dist_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
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(i__),Seq_2(i__),FWHM1(i__),FWHM2(i__),FWHMa1(i__),
     +  FWHMa2(i__),Ra(i__),Npeak(i__),AV(i__),Mpeak(i__),Mtot(i__),
     +  M50_(i__),CM(i__),e_CM(i__),npeak_1(i__),n50_(i__),cn(i__),
     +  e_cn(i__)
c    ..............Just test output...........
        write(6,3)
     +  Name(i__),Seq_2(i__),FWHM1(i__),FWHM2(i__),FWHMa1(i__),
     +  FWHMa2(i__),Ra(i__),Npeak(i__),AV(i__),Mpeak(i__),Mtot(i__),
     +  M50_(i__),CM(i__),e_CM(i__),npeak_1(i__),n50_(i__),cn(i__),
     +  e_cn(i__)
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(i__),Seq_3(i__),FWHM1_1(i__),FWHM2_1(i__),
     +  FWHMa1_1(i__),FWHMa2_1(i__),SNR(i__),Npeak_2(i__),AV_1(i__),
     +  Mpeak_1(i__),Mtot_1(i__),M50__1(i__),CM_1(i__),e_CM_1(i__),
     +  Class(i__),Ref(i__)
c    ..............Just test output...........
        write(6,4)
     +  Name_1(i__),Seq_3(i__),FWHM1_1(i__),FWHM2_1(i__),
     +  FWHMa1_1(i__),FWHMa2_1(i__),SNR(i__),Npeak_2(i__),AV_1(i__),
     +  Mpeak_1(i__),Mtot_1(i__),M50__1(i__),CM_1(i__),e_CM_1(i__),
     +  Class(i__),Ref(i__)
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(i__),Seq_4(i__),FWHM1_2(i__),FWHM2_2(i__),
     +  FWHMa1_2(i__),FWHMa2_2(i__),SNR_1(i__),Npeak_3(i__),AV_2(i__),
     +  Mpeak_2(i__),Mtot_2(i__),M50__2(i__),CM_2(i__),e_CM_2(i__),
     +  Class_1(i__),Ref_1(i__)
c    ..............Just test output...........
        write(6,5)
     +  Name_2(i__),Seq_4(i__),FWHM1_2(i__),FWHM2_2(i__),
     +  FWHMa1_2(i__),FWHMa2_2(i__),SNR_1(i__),Npeak_3(i__),AV_2(i__),
     +  Mpeak_2(i__),Mtot_2(i__),M50__2(i__),CM_2(i__),e_CM_2(i__),
     +  Class_1(i__),Ref_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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