FORTRAN Generation
(/./ftp/cats/VII/209A)

Conversion of standardized ReadMe file for file /./ftp/cats/VII/209A 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.5, on 2014-Jul-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__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. VII/209A  Galaxies Behind the Milky Way (Saito+ 1990-91; Roman+ 1996-98)
*================================================================================
*A search for galaxies behind the Milky Way between l=210deg and 230deg. (Vol. 1)
*    Saito M., Ohtani H., Asonuma A., Kashikawa N., Maki T., Nishida S.,
*    Watanabe T.
*   <Publ. Astron. Soc. Japan, Vol. 42, 603 (1990)>
*   =1990PASJ...42..603S
*Catalog of Galaxies Behind the Milky Way, l-230 to 250 degrees. (Vol. 2)
*    Saito M., Ohtani H., Baba A., Hotta N., Kameno S., Kurosu S.,
*    Nakada K., Takata T.
*   <Publ. Astron. Soc. Japan, Vol. 43, 449 (1991)>
*   =1991PASJ...43..449S
*A search for galaxies behind the Milky Way at Aquila and Sagittarius. (Vol. 3)
*    Roman A.T., Nakanishi K., Tomita A., Saito M.
*   <Publ. Astron. Soc. Jap., 48, 679 (1996)
*   =1996PASJ...48..679R
*A systematic search for galaxies behind the Milky Way at Sagittarius. (Vol. 4)
*    Roman A.T., Nakanishi K., Saito M.
*   <Publ. Astron. Soc. Jap., 50, 37-46 (1998)>
*   =1998PASJ...50...37R
*A visual search for galaxies in a Milky way region around the
*North Supergalactic Pole. (Vol. 5)
*    Roman A.T., Iwata I., Saito M.
*   <Astrophys. J., Suppl. Ser., 127, 27 (2000)>
*   =2000ApJS..127...27R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'cgmw1.dat'	! Volume 1: l=210 to 230 deg

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

      real*8        RA_         ! (deg) Right Ascension 1950
      real*8        Dec         ! (deg)     Declination 1950
      integer*4     Vol         ! [1,2] Volume number
      integer*4     CGMW        ! Running number in volume
      integer*4     Field       ! Field number of the ESO/SERC Southern
*                                     Sky Survey
      real*4        Xpos        ! (mm) X coordinate on the film (1mm ~ 67")
      real*4        Ypos        ! (mm) Y coordinate on the film (1mm ~ 67")
      integer*4     RAh         ! (h) Right ascension (1950)
      integer*4     RAm         ! (min) Right ascension (1950)
      integer*4     RAs         ! (s) [0/60] Right ascension (1950)
      character*1   DE_         ! Declination sign
      integer*4     DEd         ! (deg) Declination (1950)
      real*4        DEm         ! (arcmin) Declination (1950)
      real*8        GLON        ! (deg) Galactic longitude
      real*8        GLAT        ! (deg) Galactic latitude
      real*4        MajDiam     ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam     ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*3   MType       ! *Morphological type class
      character*10  Feature     ! *Feature of the object
      character*59  CrossId     ! Cross identifications in other catalogs
*Note on MType: the classes are:
*      E : Elliptical
*      I : Irregular
*      S : Spiral
*     SB : Barred Spiral
*Note on Feature: the following features are included:
*     Bright, Faint : for surface brightness
*     Asymm, Ring   : for morphology
*     Pair, Triple  : for close association (If accompanying galaxies are
*           almost in the same position, they are not listed in the catalog).
*     PN+galaxy     : a foreground planetary nebula overlaps a galaxy

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

C  Declarations for 'cgmw2.dat'	! Volume 2: l=230 to 250 deg

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

      real*8        RA__1       ! (deg) Right Ascension 1950
      real*8        Dec_1       ! (deg)     Declination 1950
      integer*4     Vol_1       ! [1,2] Volume number
      integer*4     CGMW_1      ! Running number in volume
      integer*4     Field_1     ! Field number of the ESO/SERC Southern
*                                     Sky Survey
      real*4        Xpos_1      ! (mm) X coordinate on the film (1mm ~ 67")
      real*4        Ypos_1      ! (mm) Y coordinate on the film (1mm ~ 67")
      integer*4     RAh_1       ! (h) Right ascension (1950)
      integer*4     RAm_1       ! (min) Right ascension (1950)
      integer*4     RAs_1       ! (s) [0/60] Right ascension (1950)
      character*1   DE__1       ! Declination sign
      integer*4     DEd_1       ! (deg) Declination (1950)
      real*4        DEm_1       ! (arcmin) Declination (1950)
      real*8        GLON_1      ! (deg) Galactic longitude
      real*8        GLAT_1      ! (deg) Galactic latitude
      real*4        MajDiam_1   ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam_1   ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*3   MType_1     ! *Morphological type class
      character*10  Feature_1   ! *Feature of the object
      character*59  CrossId_1   ! Cross identifications in other catalogs
*Note on MType: the classes are:
*      E : Elliptical
*      I : Irregular
*      S : Spiral
*     SB : Barred Spiral
*Note on Feature: the following features are included:
*     Bright, Faint : for surface brightness
*     Asymm, Ring   : for morphology
*     Pair, Triple  : for close association (If accompanying galaxies are
*           almost in the same position, they are not listed in the catalog).
*     PN+galaxy     : a foreground planetary nebula overlaps a galaxy

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

C  Declarations for 'cgmw3.dat'	! Volume 3: l=  8 to  43 deg

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

      real*8        RA__2       ! (deg) Right Ascension 1950
      real*8        Dec_2       ! (deg)     Declination 1950
      integer*4     Vol_2       ! [3] Volume number
      integer*4     CGMW_2      ! Running number in volume
      character*7   Plate       ! Name of Schmidt Atlas and field number on
*                                   which the search was performed
      character*7   RefSurv     ! *Name of the reference survey
      integer*4     Over1       ! ? Field number of overlapping plate
      integer*4     RAh_2       ! (h) Right ascension (1950)
      integer*4     RAm_2       ! (min) Right ascension (1950)
      real*4        RAs_2       ! (s) [0/60] Right ascension (1950)
      character*1   DE__2       ! Declination sign
      integer*4     DEd_2       ! (deg) Declination (1950)
      real*4        DEm_2       ! (arcmin) Declination (1950)
      real*8        GLON_2      ! (deg) Galactic longitude
      real*8        GLAT_2      ! (deg) Galactic latitude
      real*4        MajDiam_2   ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam_2   ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*13  Feature_2   ! *Feature of the object
      character*12  CrossId_2   ! Cross identifications in other catalogs
*Note on RefSurv: the survey are designated by:
*     SR = UK-SR
*      I = UK-I
*      O = POSS I-O
*     EJ = SERC-EJ
*Note on Feature: the following features are included:
*     for the surface brightness: LSB (low) or bright
*     for morphological type: S (spiral), E (elliptical), pec (peculiar)
*             and ND (nuclear dominant)
*     for multiplicity: pair or triplet, which indicates the presence of
*             one galaxy or two galaxies within a few arcminutes.
*         If an accompanying galaxy is nearly at the same position, the
*             companion is not listed in the catalog and is merely
*             indicated as 'pair'.
*         If the object seems to be close to a star or on a star, we
*             denote this as 'near *' or 'with *', respectively.
*     for aspect:
*         If the object is a possible star (the image has circular shape
*             or the faint extended component around a bright star-like
*             object), we denote it as 'star:', where the colon means
*             possible feature.
*         'PN:' means a possible planetary nebulae, i.e., the image is
*             rather round and the surface brightness distribution seems
*             to be homogeneous or the outer edge of the image is
*             relatively sharp compared with that of a galaxy.

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

C  Declarations for 'cgmw4.dat'	! Volume 4: l= -7 to  16 deg

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

      real*8        RA__3       ! (deg) Right Ascension J2000
      real*8        Dec_3       ! (deg)     Declination J2000
      integer*4     Vol_3       ! [4] Volume number
      integer*4     CGMW_3      ! Running number in volume
      character*7   Plate_1     ! Name of Schmidt Atlas and field number on
*                                   which the search was performed
      character*7   RefSurv_1   ! *Name of the reference survey
      integer*4     Over1_1     ! ? Field number of overlapping plate
      integer*4     Over2       ! ? Field number of overlapping plate
      integer*4     RAh_3       ! (h) Right ascension (J2000)
      integer*4     RAm_3       ! (min) Right ascension (J2000)
      real*4        RAs_3       ! (s) [0/60] Right ascension (J2000)
      character*1   DE__3       ! Declination sign
      integer*4     DEd_3       ! (deg) Declination (J2000)
      integer*4     DEm_3       ! (arcmin) Declination (J2000)
      integer*4     DEs         ! (arcsec) Declination (J2000)
      real*8        GLON_3      ! (deg) Galactic longitude
      real*8        GLAT_3      ! (deg) Galactic latitude
      real*4        MajDiam_3   ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam_3   ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*13  Feature_3   ! Feature of the object
*                                   (see "Note on Feature" section above)
      character*10  CrossId_3   ! Cross identifications in other catalogs
*Note on RefSurv: the survey are designated by:
*     SR = UK-SR
*      R = ESO/SRC-R

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

C  Declarations for 'cgmw5.dat'	! Volume 5: l= 32 to  68 deg

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

      real*8        RA__4       ! (deg) Right Ascension J2000
      real*8        Dec_4       ! (deg)     Declination J2000
      integer*4     Vol_4       ! [5] Volume number
      integer*4     CGMW_4      ! Running number in volume
      integer*4     II_B        ! Field number of the POSS II(B) on which we
*                                  performed the search and detected the
*                                  corresponding object
      character*1   n_II_B      ! [#] #: object also seen on the plate of
*                                         the next field
      character*6   Plate_2     ! Name and field number of the reference plate
*                                  on which the object is also seen (1)
      integer*4     RA1h        ! (h) Right ascension (1950)
      integer*4     RA1m        ! (min) Right ascension (1950)
      real*4        RA1s        ! (s) Right ascension (1950)
      character*1   DE1_        ! Declination sign (1950)
      integer*4     DE1d        ! (deg) Declination (1950)
      integer*4     DE1m        ! (arcmin) Declination (1950)
      integer*4     DE1s        ! (arcsec) Declination (1950)
      integer*4     RAh_4       ! (h) Right ascension (J2000)
      integer*4     RAm_4       ! (min) Right ascension (J2000)
      real*4        RAs_4       ! (s) Right ascension (J2000)
      character*1   DE__4       ! Declination sign (J2000)
      integer*4     DEd_4       ! (deg) Declination (J2000)
      integer*4     DEm_4       ! (arcmin) Declination (J2000)
      integer*4     DEs_1       ! (arcsec) Declination (J2000)
      real*8        GLON_4      ! (deg) Galactic longitude
      real*8        GLAT_4      ! (deg) Galactic latitude
      real*4        MajDiam_4   ! (mm) Diameter of the apparent major axis (1mm~67")
      real*4        MinDiam_4   ! (mm) ? Diameter of the apparent minor axis (1mm~67")
      character*12  Feature_4   ! Feature of the object (2)
      character*22  CrossId_4   ! Optical cross-identifications
      integer*4     cz          ! (km/s) ? Rotational velocity
      integer*4     r_cz        ! [1-10]? Reference for cz (G1)
      character*10  IRAS        ! IRAS (Cat. <II/125>) cross-identification
*Note (1): The reference plates, POSS II(I) and POSS I(E), are indicated
*     as I and E, respectively. About 70% of the detected objects appear
*     in the reference plates.
*Note (2): Feature of the object in surface brightness, morphology, and
*    multiplicity: low (LSB) and high (HSB) for the surface brightness;
*    spiral (S), elliptical (E), peculiar (pec), and nuclear dominant (ND)
*    for the morphological type; and pair, triplet, or group for the
*    multiplicity, which indicates the presence of one galaxy, two
*    galaxies, or a group of galaxies within a few arcminutes of the
*    object. If an accompanying galaxy is nearly at the same position, the
*    companion is not listed in the catalog and is merely indicated as
*    "pair." If the object seems to be close to a star or on a star, we
*    denote this as "near*" or "with*," respectively. If the object is a
*    possible star, we denote it as "star:"; "PN:" means a possible
*    planetary nebula and "neb:" a possible Galactic nebula
*--------------------------------------------------------------------------------
*
*List of Surveyed POSS II(B) Fields (Vol. 5)
*-------------------------------------------
* Field number    RA     DE
*-------------------------------------------
*   458         18 47   +30
*   459         19 10   +30
*   460         19 33   +30
*   523         18 42   +25
*   524         19 04   +25
*   525         19 26   +25
*   526         19 48   +25
*   591         18 33   +20
*   592         18 54   +20
*   593         19 15   +20
*   594         19 36   +20
*   662         18 20   +15
*   663         18 40   +15
*   664         19 00   +15
*   665         19 20   +15
*   733         18 00   +10
*   734         18 20   +10
*   735         18 40   +10
*   736         19 00   +10

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

C  Declarations for 'cgmw5ir.dat'	! Cross identification of the CGMW5 detected objects
                               with IRAS point sources

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

      real*8        RA__5       ! (deg) Right Ascension 1950
      real*8        Dec_5       ! (deg)     Declination 1950
      integer*4     Vol_5       ! [5] Volume number
      integer*4     CGMW_5      ! Running number in volume
      character*10  IRAS_1      ! IRAS (Cat. <II/125>) name
      integer*4     RAh_5       ! (h) Right ascension (1950)
      integer*4     RAm_5       ! (min) Right ascension (1950)
      real*4        RAs_5       ! (s) Right ascension (1950)
      character*1   DE__5       ! Declination sign (1950)
      integer*4     DEd_5       ! (deg) Declination (1950)
      integer*4     DEm_5       ! (arcmin) Declination (1950)
      integer*4     DEs_2       ! (arcsec) Declination (1950)
      real*4        Size        ! (mm) Size (1mm~67")
      integer*4     Sep         ! (arcsec) ? Separation
      real*4        F60um       ! (Jy) Flux density at 60{mu}m
      character*1   n_F60um     ! [:L] Note (1)
      real*4        F100um      ! (Jy) Flux density at 100{mu}m
      character*1   n_F100um    ! [:L] Note (1)
      character*5   Cat         ! Catalog names for known galaxies listed
*                                   in Table 3
      integer*4     cz_1        ! (km/s) ? Measured heliocentric radial velocity
      integer*4     r_cz_1      ! [1-11]? Reference for cz (G1)
*Note (1): ":": Uncertainty flag, "L": Point source correlation coefficient

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

C  Loading file 'cgmw1.dat'	! Volume 1: l=210 to 230 deg

C  Format for file interpretation

    1 format(
     +  I1,1X,I4,1X,I3,1X,F6.1,1X,F6.1,1X,I2,1X,I2,1X,I2,3X,A1,I2,1X,
     +  F4.1,2X,F7.3,1X,F7.3,1X,F3.1,1X,F3.1,1X,A3,1X,A10,1X,A59)

C  Effective file loading

      open(unit=1,file='cgmw1.dat', status='old')
      write(6,*) '....Loading file: cgmw1.dat'
      do i__=1,2411
        read(1,'(A144)')ar__
        read(ar__,1)
     +  Vol,CGMW,Field,Xpos,Ypos,RAh,RAm,RAs,DE_,DEd,DEm,GLON,GLAT,
     +  MajDiam,MinDiam,MType,Feature,CrossId
c       Derive coordinates RA_ and Dec from input data
c       (RA_ and Dec are set to rNULL_ when unknown)
        RA_ = RAh
        if(RAh .GE. 0) RA_=RAh*15.
        if(RAm .GE. 0) RA_=RA_+RAm/4.
        if(RAs .GE. 0) RA_=RA_+RAs/240.
        Dec = DEd
        if(DEm .GE. 0) Dec=Dec+DEm/60.
        if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c    ..............Just test output...........
        write(6,1)
     +  Vol,CGMW,Field,Xpos,Ypos,RAh,RAm,RAs,DE_,DEd,DEm,GLON,GLAT,
     +  MajDiam,MinDiam,MType,Feature,CrossId
        write(6,'(6H Pos: 2F8.4)') RA_,Dec
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cgmw2.dat'	! Volume 2: l=230 to 250 deg

C  Format for file interpretation

    2 format(
     +  I1,1X,I4,1X,I3,1X,F6.1,1X,F6.1,1X,I2,1X,I2,1X,I2,3X,A1,I2,1X,
     +  F4.1,2X,F7.3,1X,F7.3,1X,F3.1,1X,F3.1,1X,A3,1X,A10,1X,A59)

C  Effective file loading

      open(unit=1,file='cgmw2.dat', status='old')
      write(6,*) '....Loading file: cgmw2.dat'
      do i__=1,4633
        read(1,'(A144)')ar__1
        read(ar__1,2)
     +  Vol_1,CGMW_1,Field_1,Xpos_1,Ypos_1,RAh_1,RAm_1,RAs_1,DE__1,
     +  DEd_1,DEm_1,GLON_1,GLAT_1,MajDiam_1,MinDiam_1,MType_1,
     +  Feature_1,CrossId_1
c       Derive coordinates RA__1 and Dec_1 from input data
c       (RA__1 and Dec_1 are set to rNULL_ when unknown)
        RA__1 = RAh_1
        if(RAh_1 .GE. 0) RA__1=RAh_1*15.
        if(RAm_1 .GE. 0) RA__1=RA__1+RAm_1/4.
        if(RAs_1 .GE. 0) RA__1=RA__1+RAs_1/240.
        Dec_1 = DEd_1
        if(DEm_1 .GE. 0) Dec_1=Dec_1+DEm_1/60.
        if(DE__1.EQ.'-'.AND.Dec_1.GE.0) Dec_1=-Dec_1
c    ..............Just test output...........
        write(6,2)
     +  Vol_1,CGMW_1,Field_1,Xpos_1,Ypos_1,RAh_1,RAm_1,RAs_1,DE__1,
     +  DEd_1,DEm_1,GLON_1,GLAT_1,MajDiam_1,MinDiam_1,MType_1,
     +  Feature_1,CrossId_1
        write(6,'(6H Pos: 2F8.4)') RA__1,Dec_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cgmw3.dat'	! Volume 3: l=  8 to  43 deg

C  Format for file interpretation

    3 format(
     +  I1,1X,I4,1X,A7,1X,A7,1X,I4,6X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,
     +  F4.1,2X,F7.3,1X,F7.3,1X,F3.1,1X,F3.1,1X,A13,1X,A12)

C  Effective file loading

      open(unit=1,file='cgmw3.dat', status='old')
      write(6,*) '....Loading file: cgmw3.dat'
      do i__=1,5332
        read(1,'(A104)')ar__2
        read(ar__2,3)
     +  Vol_2,CGMW_2,Plate,RefSurv,Over1,RAh_2,RAm_2,RAs_2,DE__2,
     +  DEd_2,DEm_2,GLON_2,GLAT_2,MajDiam_2,MinDiam_2,Feature_2,
     +  CrossId_2
        if(ar__2(24:27) .EQ. '') Over1 = iNULL__
c       Derive coordinates RA__2 and Dec_2 from input data
c       (RA__2 and Dec_2 are set to rNULL_ when unknown)
        RA__2 = RAh_2
        if(RAh_2 .GE. 0) RA__2=RAh_2*15.
        if(RAm_2 .GE. 0) RA__2=RA__2+RAm_2/4.
        if(RAs_2 .GE. 0) RA__2=RA__2+RAs_2/240.
        Dec_2 = DEd_2
        if(DEm_2 .GE. 0) Dec_2=Dec_2+DEm_2/60.
        if(DE__2.EQ.'-'.AND.Dec_2.GE.0) Dec_2=-Dec_2
c    ..............Just test output...........
        write(6,3)
     +  Vol_2,CGMW_2,Plate,RefSurv,Over1,RAh_2,RAm_2,RAs_2,DE__2,
     +  DEd_2,DEm_2,GLON_2,GLAT_2,MajDiam_2,MinDiam_2,Feature_2,
     +  CrossId_2
        write(6,'(6H Pos: 2F8.4)') RA__2,Dec_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cgmw4.dat'	! Volume 4: l= -7 to  16 deg

C  Format for file interpretation

    4 format(
     +  I1,1X,I4,1X,A7,1X,A7,1X,I4,1X,I4,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,
     +  1X,I2,1X,I2,1X,F7.3,1X,F7.3,1X,F3.1,1X,F3.1,1X,A13,1X,A10)

C  Effective file loading

      open(unit=1,file='cgmw4.dat', status='old')
      write(6,*) '....Loading file: cgmw4.dat'
      do i__=1,7154
        read(1,'(A102)')ar__3
        read(ar__3,4)
     +  Vol_3,CGMW_3,Plate_1,RefSurv_1,Over1_1,Over2,RAh_3,RAm_3,
     +  RAs_3,DE__3,DEd_3,DEm_3,DEs,GLON_3,GLAT_3,MajDiam_3,MinDiam_3,
     +  Feature_3,CrossId_3
        if(ar__3(24:27) .EQ. '') Over1_1 = iNULL__
        if(ar__3(29:32) .EQ. '') Over2 = iNULL__
c       Derive coordinates RA__3 and Dec_3 from input data
c       (RA__3 and Dec_3 are set to rNULL_ when unknown)
        RA__3 = RAh_3
        if(RAh_3 .GE. 0) RA__3=RAh_3*15.
        if(RAm_3 .GE. 0) RA__3=RA__3+RAm_3/4.
        if(RAs_3 .GE. 0) RA__3=RA__3+RAs_3/240.
        Dec_3 = DEd_3
        if(DEm_3 .GE. 0) Dec_3=Dec_3+DEm_3/60.
        if(DEs .GE. 0) Dec_3=Dec_3+DEs/3600.
        if(DE__3.EQ.'-'.AND.Dec_3.GE.0) Dec_3=-Dec_3
c    ..............Just test output...........
        write(6,4)
     +  Vol_3,CGMW_3,Plate_1,RefSurv_1,Over1_1,Over2,RAh_3,RAm_3,
     +  RAs_3,DE__3,DEd_3,DEm_3,DEs,GLON_3,GLAT_3,MajDiam_3,MinDiam_3,
     +  Feature_3,CrossId_3
        write(6,'(6H Pos: 2F8.4)') RA__3,Dec_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cgmw5.dat'	! Volume 5: l= 32 to  68 deg

C  Format for file interpretation

    5 format(
     +  I1,1X,I5,1X,I3,A1,1X,A6,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,
     +  I2,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F7.3,1X,F7.3,
     +  1X,F3.1,1X,F3.1,1X,A12,1X,A22,2X,I5,3X,I1,1X,A10)

C  Effective file loading

      open(unit=1,file='cgmw5.dat', status='old')
      write(6,*) '....Loading file: cgmw5.dat'
      do i__=1,11310
        read(1,'(A145)')ar__4
        read(ar__4,5)
     +  Vol_4,CGMW_4,II_B,n_II_B,Plate_2,RA1h,RA1m,RA1s,DE1_,DE1d,
     +  DE1m,DE1s,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_1,GLON_4,
     +  GLAT_4,MajDiam_4,MinDiam_4,Feature_4,CrossId_4,cz,r_cz,IRAS
        if(ar__4(84:86) .EQ. '') MinDiam_4 = rNULL__
        if(ar__4(125:129) .EQ. '') cz = iNULL__
        if(ar__4(133:133) .EQ. '') r_cz = iNULL__
c       Derive coordinates RA__4 and Dec_4 from input data
c       (RA__4 and Dec_4 are set to rNULL_ when unknown)
        RA__4 = RAh_4
        if(RAh_4 .GE. 0) RA__4=RAh_4*15.
        if(RAm_4 .GE. 0) RA__4=RA__4+RAm_4/4.
        if(RAs_4 .GE. 0) RA__4=RA__4+RAs_4/240.
        Dec_4 = DEd_4
        if(DEm_4 .GE. 0) Dec_4=Dec_4+DEm_4/60.
        if(DEs_1 .GE. 0) Dec_4=Dec_4+DEs_1/3600.
        if(DE__4.EQ.'-'.AND.Dec_4.GE.0) Dec_4=-Dec_4
c    ..............Just test output...........
        write(6,5)
     +  Vol_4,CGMW_4,II_B,n_II_B,Plate_2,RA1h,RA1m,RA1s,DE1_,DE1d,
     +  DE1m,DE1s,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_1,GLON_4,
     +  GLAT_4,MajDiam_4,MinDiam_4,Feature_4,CrossId_4,cz,r_cz,IRAS
        write(6,'(6H Pos: 2F8.4)') RA__4,Dec_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'cgmw5ir.dat'	! Cross identification of the CGMW5 detected objects
*                               with IRAS point sources

C  Format for file interpretation

    6 format(
     +  I1,1X,I5,1X,A10,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,
     +  F3.1,1X,I2,1X,F5.2,A1,1X,F5.2,A1,3X,A5,I5,1X,I2)

C  Effective file loading

      open(unit=1,file='cgmw5ir.dat', status='old')
      write(6,*) '....Loading file: cgmw5ir.dat'
      do i__=1,159
        read(1,'(A76)')ar__5
        read(ar__5,6)
     +  Vol_5,CGMW_5,IRAS_1,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,DEm_5,DEs_2,
     +  Size,Sep,F60um,n_F60um,F100um,n_F100um,Cat,cz_1,r_cz_1
        if(ar__5(45:46) .EQ. '') Sep = iNULL__
        if(ar__5(69:73) .EQ. '') cz_1 = iNULL__
        if(ar__5(75:76) .EQ. '') r_cz_1 = iNULL__
c       Derive coordinates RA__5 and Dec_5 from input data
c       (RA__5 and Dec_5 are set to rNULL_ when unknown)
        RA__5 = RAh_5
        if(RAh_5 .GE. 0) RA__5=RAh_5*15.
        if(RAm_5 .GE. 0) RA__5=RA__5+RAm_5/4.
        if(RAs_5 .GE. 0) RA__5=RA__5+RAs_5/240.
        Dec_5 = DEd_5
        if(DEm_5 .GE. 0) Dec_5=Dec_5+DEm_5/60.
        if(DEs_2 .GE. 0) Dec_5=Dec_5+DEs_2/3600.
        if(DE__5.EQ.'-'.AND.Dec_5.GE.0) Dec_5=-Dec_5
c    ..............Just test output...........
        write(6,6)
     +  Vol_5,CGMW_5,IRAS_1,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,DEm_5,DEs_2,
     +  Size,Sep,F60um,n_F60um,F100um,n_F100um,Cat,cz_1,r_cz_1
        write(6,'(6H Pos: 2F8.4)') RA__5,Dec_5
c    .......End.of.Just test output...........
      end do
      close(1)

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