FORTRAN Generation
(/./ftp/cats/J/ApJ/702/178)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/702/178 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-28
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/702/178      YSOs in the central 400pc of the Galaxy  (Yusef-Zadeh+, 2009)
*================================================================================
*Star formation in the central 400 pc of the Milky Way:
*evidence for a population of massive young stellar objects.
*    Yusef-Zadeh F., Hewitt J.W., Arendt R.G., Whitney B., Rieke G., Wardle M.,
*    Hinz J.L., Stolovy S., Lang C.C., Burton M.G., Ramirez S.
*   <Astrophys. J., 702, 178-225 (2009)>
*   =2009ApJ...702..178Y
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Galactic coordinates and flux density of
                             sources in G359.43+0.02

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

C  Position composed of: GLON GLAT
      integer*4     v_YHA2009_ (nr__) ! Source number
      character*1   m__YHA2009_(nr__) ! [e] e: Extended component
      character*1   n__YHA2009_(nr__) ! [abc] Individual notes (1)
      real*8        GLON       (nr__) ! (deg) Galactic longitude
      real*8        GLAT       (nr__) ! (deg) Galactic latitude
      real*4        S6cm       (nr__) ! ? VLA 6cm flux density
      real*4        e_S6cm     (nr__) ! ? S6cm uncertainty
      real*4        S20cm      (nr__) ! ? VLA 20cm flux density
      real*4        e_S20cm    (nr__) ! ? S20cm uncertainty
      real*4        amaj       (nr__) ! (arcsec) ? Major axis
      real*4        amin       (nr__) ! (arcsec) ? Minor axis
      integer*4     Surf       (nr__) ! (arcsec2) ? Surface when circular
*Note (1): Flag as follows:
*   a = Infrared counterpart 300758.
*   b = Infrared counterpart 284291.
*   c = Source 11 has a central compact source and an extended component.

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

C  Declarations for 'table3.dat'	! parameters of the SED fits to candidate YSOs
                             in l<0{deg}

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

C  Position composed of: GLON GLAT
      integer*4     SSTGC      (nr__1) ! SST (Spitzer Space Telescope) Galactic
*                                    center (GC) number (1)
      character*3   n_SSTGC    (nr__1) ! [abc, ] Notes for individual sources (2)
      real*8        GLON_1     (nr__1) ! (deg) Galactic longitude
      real*8        GLAT_1     (nr__1) ! (deg) Galactic latitude
      real*4        Chi2       (nr__1) ! {chi}^2^ value
      integer*4     Nfits      (nr__1) ! Number of acceptably fitting models
      real*4        Av         (nr__1) ! (mag) Visual extinction A_v_
      real*4        Mass       (nr__1) ! (Msun) Average mass of the protostar
      real*4        e_Mass     (nr__1) ! (Msun) Mass uncertainty
      real*4        Lum        (nr__1) ! (Lsun) Average luminosity of the protostar
      real*4        e_Lum      (nr__1) ! (Lsun) Lum uncertainty
      real*4        Menv       (nr__1) ! (Msun/yr) Average mass accretion rate
      real*4        e_Menv     (nr__1) ! (Msun/yr) Menv 1{sigma} standard deviation
      character*3   Stage      (nr__1) ! Evolutionary stage of the fitted source
*                                    (I, II or III)
*Note (1): Based on the IRAC catalogs (Ramirez et al. 2008, Cat. II/295;
*          <SSTGC NNNNNNN> in Simbad)
*Note (2): Flag as follows:
*    a = radio counterpart 6
*    b = radio counterpart 5
*    c = Candidate YSO in the G359.43+0.02 cluster (table 2 of the paper; note
*        that "GLAT" for SSTGC 307175 in table 2 is not 0.02637 but, indeed,
*        -0.00026 as in table 3).

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

C  Declarations for 'table4.dat'	! parameters of the SED fits to candidate YSOs
                             in l>0{deg}

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

C  Position composed of: GLON GLAT
      integer*4     SSTGC_1    (nr__2) ! SST (Spitzer Space Telescope) Galactic
*                                    center (GC) number (1)
      character*3   n_SSTGC_1  (nr__2) ! [abc, ] Notes for individual sources (2)
      real*8        GLON_2     (nr__2) ! (deg) Galactic longitude
      real*8        GLAT_2     (nr__2) ! (deg) Galactic latitude
      real*4        Chi2_1     (nr__2) ! {chi}^2^ value
      integer*4     Nfits_1    (nr__2) ! Number of acceptably fitting models
      real*4        Av_1       (nr__2) ! (mag) Visual extinction A_v_
      real*4        Mass_1     (nr__2) ! (Msun) Average mass of the protostar
      real*4        e_Mass_1   (nr__2) ! (Msun) Mass uncertainty
      real*4        Lum_1      (nr__2) ! (Lsun) Average luminosity of the protostar
      real*4        e_Lum_1    (nr__2) ! (Lsun) Lum uncertainty
      real*4        Menv_1     (nr__2) ! (Msun/yr) Average mass accretion rate
      real*4        e_Menv_1   (nr__2) ! (Msun/yr) Menv 1{sigma} standard deviation
      character*3   Stage_1    (nr__2) ! Evolutionary stage of the fitted source
*                                    (I, II or III)
*Note (1): Based on the IRAC catalogs (Ramirez et al. 2008, Cat. II/295;
*          <SSTGC NNNNNNN> in Simbad)
*Note (2): Flag as follows:
*    a = radio counterpart 6
*    b = radio counterpart 5
*    c = Candidate YSO in the G359.43+0.02 cluster (table 2 of the paper; note
*        that "GLAT" for SSTGC 307175 in table 2 is not 0.02637 but, indeed,
*        -0.00026 as in table 3).

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

C  Declarations for 'table5.dat'	! Table of sources with 4.5um Excess Emission

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

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

      integer*4     ID         (nr__3) ! Green source identification
*                                  (<[YHA2009] gNN> in Simbad) (1)
      real*8        GLON_3     (nr__3) ! (deg) Galactic longitude
      real*4        GLAT_3     (nr__3) ! (deg) Galactic latitude
      integer*4     RAh        (nr__3) ! (h) Hour of right ascension (J2000)
      integer*4     RAm        (nr__3) ! (min) Minute of right ascension (J2000)
      real*4        RAs        (nr__3) ! (s) Second of right ascension (J2000)
      character*1   DE_        (nr__3) ! Sign of declination (J2000)
      integer*4     DEd        (nr__3) ! (deg) Degree of declination (J2000)
      integer*4     DEm        (nr__3) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs        (nr__3) ! (arcsec) Arcsecond of declination (J2000)
      character*27  Comm       (nr__3) ! Description
*Note (1): Visually identified from IRAC data. Many point sources may be
*     multiple or extended sources at higher resolution.

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

C  Declarations for 'table6.dat'	! Table of 70um Sources

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

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

      integer*4     ID_1       (nr__4) ! Source identification from table 5, or from
*                                  table8 when n_ID is 'a' (1)
      character*1   n_ID       (nr__4) ! [a] Methanol class II maser source (2)
      real*8        GLON_4     (nr__4) ! (deg) Galactic longitude
      real*8        GLAT_4     (nr__4) ! (deg) Galactic latitude
      integer*4     RAh_1      (nr__4) ! (h) Hour of right ascension (J2000)
      integer*4     RAm_1      (nr__4) ! (min) Minute of right ascension (J2000)
      real*4        RAs_1      (nr__4) ! (s) Second of right ascension (J2000)
      character*1   DE__1      (nr__4) ! Sign of declination (J2000)
      integer*4     DEd_1      (nr__4) ! (deg) Degree of declination (J2000)
      integer*4     DEm_1      (nr__4) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs_1      (nr__4) ! (arcsec) Arcsecond of declination (J2000)
      integer*4     Flux       (nr__4) ! (mJy) ? Flux
      real*4        S_N        (nr__4) ! ? Signal to noise ratio
*Note (1): Subset of 4.5um excess sources that were observed with MIPS at 70um.
*Note (2): Source No. given is for table of methanol class II maser sources
*          (table 8).

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

C  Declarations for 'table7.dat'	! Table of Methanol class I maser sources

      integer*4 nr__5
      parameter (nr__5=21)	! Number of records
      character*55 ar__5  	! Full-size record

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

      integer*4     ID_2       (nr__5) ! Source identification
      real*8        GLON_5     (nr__5) ! (deg) Galactic longitude
      real*4        GLAT_5     (nr__5) ! (deg) Galactic latitude
      integer*4     RAh_2      (nr__5) ! (h) Hour of right ascension (J2000)
      integer*4     RAm_2      (nr__5) ! (min) Minute of right ascension (J2000)
      real*4        RAs_2      (nr__5) ! (s) Second of right ascension (J2000)
      character*1   DE__2      (nr__5) ! Sign of declination (J2000)
      integer*4     DEd_2      (nr__5) ! (deg) Degree of declination (J2000)
      integer*4     DEm_2      (nr__5) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs_2      (nr__5) ! (arcsec) Arcsecond of declination (J2000)
      real*4        Vlsr       (nr__5) ! (km/s) LSR velocity
      character*7   Ref        (nr__5) ! Reference(s) (1)
*Note (1): References as follows:
*   H90 = Haschick et al. 1990ApJ...354..556H (44GHz);
*   S94 = Slysh et al. 1994MNRAS.268..464S (44GHz);
*   M97 = Mehringer & Menten 1997ApJ...474..346M (44GHz, multiple sources in
*         Sgr B2);
*   V00 = Val'tts et al. 2000MNRAS.317..315V (95GHz).
*   C96 = Caswell 1996MNRAS.283..606C;
*   W98 = Walsh et al. 1998, Cat. J/MNRAS/301/640

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

C  Declarations for 'table8.dat'	! Table of Methanol class II maser sources

      integer*4 nr__6
      parameter (nr__6=23)	! Number of records
      character*55 ar__6  	! Full-size record

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

      integer*4     ID_3       (nr__6) ! Source identification
      real*8        GLON_6     (nr__6) ! (deg) Galactic longitude
      real*4        GLAT_6     (nr__6) ! (deg) Galactic latitude
      integer*4     RAh_3      (nr__6) ! (h) Hour of right ascension (J2000)
      integer*4     RAm_3      (nr__6) ! (min) Minute of right ascension (J2000)
      real*4        RAs_3      (nr__6) ! (s) Second of right ascension (J2000)
      character*1   DE__3      (nr__6) ! Sign of declination (J2000)
      integer*4     DEd_3      (nr__6) ! (deg) Degree of declination (J2000)
      integer*4     DEm_3      (nr__6) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs_3      (nr__6) ! (arcsec) Arcsecond of declination (J2000)
      real*4        Vlsr_1     (nr__6) ! (km/s) LSR velocity
      character*7   Ref_1      (nr__6) ! Reference(s) (1)
*Note (1): References as follows:
*   H90 = Haschick et al. 1990ApJ...354..556H (44GHz);
*   S94 = Slysh et al. 1994MNRAS.268..464S (44GHz);
*   M97 = Mehringer & Menten 1997ApJ...474..346M (44GHz, multiple sources in
*         Sgr B2);
*   V00 = Val'tts et al. 2000MNRAS.317..315V (95GHz).
*   C96 = Caswell 1996MNRAS.283..606C;
*   W98 = Walsh et al. 1998, Cat. J/MNRAS/301/640

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

C  Declarations for 'table9.dat'	! Parameters of SED fits to 4.5um excess sources
                             near the Galactic Center

      integer*4 nr__7
      parameter (nr__7=14)	! Number of records
      character*90 ar__7  	! Full-size record

C  Position composed of: GLON GLAT
      character*4   v_YHA2009__1(nr__7) ! Green number (gNNA)
      integer*4     SSTGC_2    (nr__7) ! SSTGC number
      real*8        GLON_7     (nr__7) ! (deg) Galactic longitude
      real*8        GLAT_7     (nr__7) ! (deg) Galactic latitude
      real*4        Chi2_2     (nr__7) ! {chi}^2^ value
      integer*4     Nfits_2    (nr__7) ! Number of acceptably fitting models
      real*4        Av_2       (nr__7) ! (mag) Visual extinction A_v_
      real*4        Mass_2     (nr__7) ! (Msun) Average mass of the protostar
      real*4        e_Mass_2   (nr__7) ! (Msun) Mass uncertainty
      real*4        Lum_2      (nr__7) ! (Lsun) Average luminosity of the protostar
      real*4        e_Lum_2    (nr__7) ! (Lsun) Lum uncertainty
      real*4        Menv_2     (nr__7) ! (Msun/yr) Average mass accretion rate
      real*4        e_Menv_2   (nr__7) ! (Msun/yr) Menv 1{sigma} standard deviation
      character*3   Stage_2    (nr__7) ! Evolutionary stage of the fitted source

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

C  Declarations for 'table10.dat'	! parameters of SED fits to 4.5um excess sources
                             foreground to the Galactic Center

      integer*4 nr__8
      parameter (nr__8=19)	! Number of records
      character*90 ar__8  	! Full-size record

C  Position composed of: GLON GLAT
      character*4   v_YHA2009__2(nr__8) ! Green number (gNNA)
      integer*4     SSTGC_3    (nr__8) ! SSTGC number
      real*8        GLON_8     (nr__8) ! (deg) Galactic longitude
      real*8        GLAT_8     (nr__8) ! (deg) Galactic latitude
      real*4        Chi2_3     (nr__8) ! {chi}^2^ value
      integer*4     Nfits_3    (nr__8) ! Number of acceptably fitting models
      real*4        Av_3       (nr__8) ! (mag) Visual extinction A_v_
      real*4        Mass_3     (nr__8) ! (Msun) Average mass of the protostar
      real*4        e_Mass_3   (nr__8) ! (Msun) Mass uncertainty
      real*4        Lum_3      (nr__8) ! (Lsun) Average luminosity of the protostar
      real*4        e_Lum_3    (nr__8) ! (Lsun) Lum uncertainty
      real*4        Menv_3     (nr__8) ! (Msun/yr) Average mass accretion rate
      real*4        e_Menv_3   (nr__8) ! (Msun/yr) Menv 1{sigma} standard deviation
      character*3   Stage_3    (nr__8) ! Evolutionary stage of the fitted source

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

C  Loading file 'table1.dat'	! Galactic coordinates and flux density of
*                             sources in G359.43+0.02

C  Format for file interpretation

    1 format(
     +  I2,A1,1X,A1,1X,F10.6,1X,F9.6,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,
     +  1X,F4.1,1X,F4.1,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,18
        read(1,'(A61)')ar__
        read(ar__,1)
     +  v_YHA2009_(i__),m__YHA2009_(i__),n__YHA2009_(i__),GLON(i__),
     +  GLAT(i__),S6cm(i__),e_S6cm(i__),S20cm(i__),e_S20cm(i__),
     +  amaj(i__),amin(i__),Surf(i__)
        if(ar__(28:31) .EQ. '') S6cm(i__) = rNULL__
        if(ar__(33:36) .EQ. '') e_S6cm(i__) = rNULL__
        if(ar__(38:41) .EQ. '') S20cm(i__) = rNULL__
        if(ar__(43:46) .EQ. '') e_S20cm(i__) = rNULL__
        if(ar__(48:51) .EQ. '') amaj(i__) = rNULL__
        if(ar__(53:56) .EQ. '') amin(i__) = rNULL__
        if(ar__(58:61) .EQ. '') Surf(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  v_YHA2009_(i__),m__YHA2009_(i__),n__YHA2009_(i__),GLON(i__),
     +  GLAT(i__),S6cm(i__),e_S6cm(i__),S20cm(i__),e_S20cm(i__),
     +  amaj(i__),amin(i__),Surf(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! parameters of the SED fits to candidate YSOs
*                             in l<0{deg}

C  Format for file interpretation

    2 format(
     +  I6,A3,1X,F8.4,1X,F8.5,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X,F3.1,
     +  1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,223
        read(1,'(A88)')ar__1
        read(ar__1,2)
     +  SSTGC(i__),n_SSTGC(i__),GLON_1(i__),GLAT_1(i__),Chi2(i__),
     +  Nfits(i__),Av(i__),Mass(i__),e_Mass(i__),Lum(i__),e_Lum(i__),
     +  Menv(i__),e_Menv(i__),Stage(i__)
c    ..............Just test output...........
        write(6,2)
     +  SSTGC(i__),n_SSTGC(i__),GLON_1(i__),GLAT_1(i__),Chi2(i__),
     +  Nfits(i__),Av(i__),Mass(i__),e_Mass(i__),Lum(i__),e_Lum(i__),
     +  Menv(i__),e_Menv(i__),Stage(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! parameters of the SED fits to candidate YSOs
*                             in l>0{deg}

C  Format for file interpretation

    3 format(
     +  I6,A3,1X,F8.4,1X,F8.5,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X,F3.1,
     +  1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,137
        read(1,'(A88)')ar__2
        read(ar__2,3)
     +  SSTGC_1(i__),n_SSTGC_1(i__),GLON_2(i__),GLAT_2(i__),
     +  Chi2_1(i__),Nfits_1(i__),Av_1(i__),Mass_1(i__),e_Mass_1(i__),
     +  Lum_1(i__),e_Lum_1(i__),Menv_1(i__),e_Menv_1(i__),Stage_1(i__)
c    ..............Just test output...........
        write(6,3)
     +  SSTGC_1(i__),n_SSTGC_1(i__),GLON_2(i__),GLAT_2(i__),
     +  Chi2_1(i__),Nfits_1(i__),Av_1(i__),Mass_1(i__),e_Mass_1(i__),
     +  Lum_1(i__),e_Lum_1(i__),Menv_1(i__),e_Menv_1(i__),Stage_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Table of sources with 4.5um Excess Emission

C  Format for file interpretation

    4 format(
     +  I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,
     +  1X,A27)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,33
        read(1,'(A69)')ar__3
        read(ar__3,4)
     +  ID(i__),GLON_3(i__),GLAT_3(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),Comm(i__)
        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,4)
     +  ID(i__),GLON_3(i__),GLAT_3(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),Comm(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 'table6.dat'	! Table of 70um Sources

C  Format for file interpretation

    5 format(
     +  I2,A1,1X,F9.5,1X,F8.5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,
     +  F4.1,1X,I5,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,9
        read(1,'(A58)')ar__4
        read(ar__4,5)
     +  ID_1(i__),n_ID(i__),GLON_4(i__),GLAT_4(i__),RAh_1(i__),
     +  RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),Flux(i__),S_N(i__)
        if(ar__4(48:52) .EQ. '') Flux(i__) = iNULL__
        if(ar__4(54:58) .EQ. '') S_N(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,5)
     +  ID_1(i__),n_ID(i__),GLON_4(i__),GLAT_4(i__),RAh_1(i__),
     +  RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),Flux(i__),S_N(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 'table7.dat'	! Table of Methanol class I maser sources

C  Format for file interpretation

    6 format(
     +  I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,
     +  1X,F5.1,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,21
        read(1,'(A55)')ar__5
        read(ar__5,6)
     +  ID_2(i__),GLON_5(i__),GLAT_5(i__),RAh_2(i__),RAm_2(i__),
     +  RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
     +  Vlsr(i__),Ref(i__)
        RAdeg_2(i__) = rNULL__
        DEdeg_2(i__) = rNULL__
c  Derive coordinates RAdeg_2 and DEdeg_2 from input data
c  (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
        if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
        if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
        if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
        if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
        if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
        if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
        if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c    ..............Just test output...........
        write(6,6)
     +  ID_2(i__),GLON_5(i__),GLAT_5(i__),RAh_2(i__),RAm_2(i__),
     +  RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
     +  Vlsr(i__),Ref(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Table of Methanol class II maser sources

C  Format for file interpretation

    7 format(
     +  I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,
     +  1X,F5.1,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,23
        read(1,'(A55)')ar__6
        read(ar__6,7)
     +  ID_3(i__),GLON_6(i__),GLAT_6(i__),RAh_3(i__),RAm_3(i__),
     +  RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),
     +  Vlsr_1(i__),Ref_1(i__)
        RAdeg_3(i__) = rNULL__
        DEdeg_3(i__) = rNULL__
c  Derive coordinates RAdeg_3 and DEdeg_3 from input data
c  (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
        if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
        if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
        if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
        if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
        if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
        if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
        if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c    ..............Just test output...........
        write(6,7)
     +  ID_3(i__),GLON_6(i__),GLAT_6(i__),RAh_3(i__),RAm_3(i__),
     +  RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),
     +  Vlsr_1(i__),Ref_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Parameters of SED fits to 4.5um excess sources
*                             near the Galactic Center

C  Format for file interpretation

    8 format(
     +  A4,1X,I7,1X,F8.4,1X,F7.4,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X,
     +  F3.1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,14
        read(1,'(A90)')ar__7
        read(ar__7,8)
     +  v_YHA2009__1(i__),SSTGC_2(i__),GLON_7(i__),GLAT_7(i__),
     +  Chi2_2(i__),Nfits_2(i__),Av_2(i__),Mass_2(i__),e_Mass_2(i__),
     +  Lum_2(i__),e_Lum_2(i__),Menv_2(i__),e_Menv_2(i__),Stage_2(i__)
c    ..............Just test output...........
        write(6,8)
     +  v_YHA2009__1(i__),SSTGC_2(i__),GLON_7(i__),GLAT_7(i__),
     +  Chi2_2(i__),Nfits_2(i__),Av_2(i__),Mass_2(i__),e_Mass_2(i__),
     +  Lum_2(i__),e_Lum_2(i__),Menv_2(i__),e_Menv_2(i__),Stage_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! parameters of SED fits to 4.5um excess sources
*                             foreground to the Galactic Center

C  Format for file interpretation

    9 format(
     +  A4,1X,I7,1X,F8.4,1X,F7.4,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X,
     +  F3.1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,19
        read(1,'(A90)')ar__8
        read(ar__8,9)
     +  v_YHA2009__2(i__),SSTGC_3(i__),GLON_8(i__),GLAT_8(i__),
     +  Chi2_3(i__),Nfits_3(i__),Av_3(i__),Mass_3(i__),e_Mass_3(i__),
     +  Lum_3(i__),e_Lum_3(i__),Menv_3(i__),e_Menv_3(i__),Stage_3(i__)
c    ..............Just test output...........
        write(6,9)
     +  v_YHA2009__2(i__),SSTGC_3(i__),GLON_8(i__),GLAT_8(i__),
     +  Chi2_3(i__),Nfits_3(i__),Av_3(i__),Mass_3(i__),e_Mass_3(i__),
     +  Lum_3(i__),e_Lum_3(i__),Menv_3(i__),e_Menv_3(i__),Stage_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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