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

Conversion of standardized ReadMe file for file /./ftp/cats/VII/209A 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.8, on 2014-Nov-22
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_        (nr__) ! (deg) Right Ascension 1950
      real*8        Dec        (nr__) ! (deg)     Declination 1950
      integer*4     Vol        (nr__) ! [1,2] Volume number
      integer*4     CGMW       (nr__) ! Running number in volume
      integer*4     Field      (nr__) ! Field number of the ESO/SERC Southern
*                                     Sky Survey
      real*4        Xpos       (nr__) ! (mm) X coordinate on the film (1mm ~ 67")
      real*4        Ypos       (nr__) ! (mm) Y coordinate on the film (1mm ~ 67")
      integer*4     RAh        (nr__) ! (h) Right ascension (1950)
      integer*4     RAm        (nr__) ! (min) Right ascension (1950)
      integer*4     RAs        (nr__) ! (s) [0/60] Right ascension (1950)
      character*1   DE_        (nr__) ! Declination sign
      integer*4     DEd        (nr__) ! (deg) Declination (1950)
      real*4        DEm        (nr__) ! (arcmin) Declination (1950)
      real*8        GLON       (nr__) ! (deg) Galactic longitude
      real*8        GLAT       (nr__) ! (deg) Galactic latitude
      real*4        MajDiam    (nr__) ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam    (nr__) ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*3   MType      (nr__) ! *Morphological type class
      character*10  Feature    (nr__) ! *Feature of the object
      character*59  CrossId    (nr__) ! 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      (nr__1) ! (deg) Right Ascension 1950
      real*8        Dec_1      (nr__1) ! (deg)     Declination 1950
      integer*4     Vol_1      (nr__1) ! [1,2] Volume number
      integer*4     CGMW_1     (nr__1) ! Running number in volume
      integer*4     Field_1    (nr__1) ! Field number of the ESO/SERC Southern
*                                     Sky Survey
      real*4        Xpos_1     (nr__1) ! (mm) X coordinate on the film (1mm ~ 67")
      real*4        Ypos_1     (nr__1) ! (mm) Y coordinate on the film (1mm ~ 67")
      integer*4     RAh_1      (nr__1) ! (h) Right ascension (1950)
      integer*4     RAm_1      (nr__1) ! (min) Right ascension (1950)
      integer*4     RAs_1      (nr__1) ! (s) [0/60] Right ascension (1950)
      character*1   DE__1      (nr__1) ! Declination sign
      integer*4     DEd_1      (nr__1) ! (deg) Declination (1950)
      real*4        DEm_1      (nr__1) ! (arcmin) Declination (1950)
      real*8        GLON_1     (nr__1) ! (deg) Galactic longitude
      real*8        GLAT_1     (nr__1) ! (deg) Galactic latitude
      real*4        MajDiam_1  (nr__1) ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam_1  (nr__1) ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*3   MType_1    (nr__1) ! *Morphological type class
      character*10  Feature_1  (nr__1) ! *Feature of the object
      character*59  CrossId_1  (nr__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      (nr__2) ! (deg) Right Ascension 1950
      real*8        Dec_2      (nr__2) ! (deg)     Declination 1950
      integer*4     Vol_2      (nr__2) ! [3] Volume number
      integer*4     CGMW_2     (nr__2) ! Running number in volume
      character*7   Plate      (nr__2) ! Name of Schmidt Atlas and field number on
*                                   which the search was performed
      character*7   RefSurv    (nr__2) ! *Name of the reference survey
      integer*4     Over1      (nr__2) ! ? Field number of overlapping plate
      integer*4     RAh_2      (nr__2) ! (h) Right ascension (1950)
      integer*4     RAm_2      (nr__2) ! (min) Right ascension (1950)
      real*4        RAs_2      (nr__2) ! (s) [0/60] Right ascension (1950)
      character*1   DE__2      (nr__2) ! Declination sign
      integer*4     DEd_2      (nr__2) ! (deg) Declination (1950)
      real*4        DEm_2      (nr__2) ! (arcmin) Declination (1950)
      real*8        GLON_2     (nr__2) ! (deg) Galactic longitude
      real*8        GLAT_2     (nr__2) ! (deg) Galactic latitude
      real*4        MajDiam_2  (nr__2) ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam_2  (nr__2) ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*13  Feature_2  (nr__2) ! *Feature of the object
      character*12  CrossId_2  (nr__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      (nr__3) ! (deg) Right Ascension J2000
      real*8        Dec_3      (nr__3) ! (deg)     Declination J2000
      integer*4     Vol_3      (nr__3) ! [4] Volume number
      integer*4     CGMW_3     (nr__3) ! Running number in volume
      character*7   Plate_1    (nr__3) ! Name of Schmidt Atlas and field number on
*                                   which the search was performed
      character*7   RefSurv_1  (nr__3) ! *Name of the reference survey
      integer*4     Over1_1    (nr__3) ! ? Field number of overlapping plate
      integer*4     Over2      (nr__3) ! ? Field number of overlapping plate
      integer*4     RAh_3      (nr__3) ! (h) Right ascension (J2000)
      integer*4     RAm_3      (nr__3) ! (min) Right ascension (J2000)
      real*4        RAs_3      (nr__3) ! (s) [0/60] Right ascension (J2000)
      character*1   DE__3      (nr__3) ! Declination sign
      integer*4     DEd_3      (nr__3) ! (deg) Declination (J2000)
      integer*4     DEm_3      (nr__3) ! (arcmin) Declination (J2000)
      integer*4     DEs        (nr__3) ! (arcsec) Declination (J2000)
      real*8        GLON_3     (nr__3) ! (deg) Galactic longitude
      real*8        GLAT_3     (nr__3) ! (deg) Galactic latitude
      real*4        MajDiam_3  (nr__3) ! (mm) Size of apparent major diameter (1mm ~ 67")
      real*4        MinDiam_3  (nr__3) ! (mm) Size of apparent minor diameter (1mm ~ 67")
      character*13  Feature_3  (nr__3) ! Feature of the object
*                                   (see "Note on Feature" section above)
      character*10  CrossId_3  (nr__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      (nr__4) ! (deg) Right Ascension J2000
      real*8        Dec_4      (nr__4) ! (deg)     Declination J2000
      integer*4     Vol_4      (nr__4) ! [5] Volume number
      integer*4     CGMW_4     (nr__4) ! Running number in volume
      integer*4     II_B       (nr__4) ! Field number of the POSS II(B) on which we
*                                  performed the search and detected the
*                                  corresponding object
      character*1   n_II_B     (nr__4) ! [#] #: object also seen on the plate of
*                                         the next field
      character*6   Plate_2    (nr__4) ! Name and field number of the reference plate
*                                  on which the object is also seen (1)
      integer*4     RA1h       (nr__4) ! (h) Right ascension (1950)
      integer*4     RA1m       (nr__4) ! (min) Right ascension (1950)
      real*4        RA1s       (nr__4) ! (s) Right ascension (1950)
      character*1   DE1_       (nr__4) ! Declination sign (1950)
      integer*4     DE1d       (nr__4) ! (deg) Declination (1950)
      integer*4     DE1m       (nr__4) ! (arcmin) Declination (1950)
      integer*4     DE1s       (nr__4) ! (arcsec) Declination (1950)
      integer*4     RAh_4      (nr__4) ! (h) Right ascension (J2000)
      integer*4     RAm_4      (nr__4) ! (min) Right ascension (J2000)
      real*4        RAs_4      (nr__4) ! (s) Right ascension (J2000)
      character*1   DE__4      (nr__4) ! Declination sign (J2000)
      integer*4     DEd_4      (nr__4) ! (deg) Declination (J2000)
      integer*4     DEm_4      (nr__4) ! (arcmin) Declination (J2000)
      integer*4     DEs_1      (nr__4) ! (arcsec) Declination (J2000)
      real*8        GLON_4     (nr__4) ! (deg) Galactic longitude
      real*8        GLAT_4     (nr__4) ! (deg) Galactic latitude
      real*4        MajDiam_4  (nr__4) ! (mm) Diameter of the apparent major axis (1mm~67")
      real*4        MinDiam_4  (nr__4) ! (mm) ? Diameter of the apparent minor axis (1mm~67")
      character*12  Feature_4  (nr__4) ! Feature of the object (2)
      character*22  CrossId_4  (nr__4) ! Optical cross-identifications
      integer*4     cz         (nr__4) ! (km/s) ? Rotational velocity
      integer*4     r_cz       (nr__4) ! [1-10]? Reference for cz (G1)
      character*10  IRAS       (nr__4) ! 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      (nr__5) ! (deg) Right Ascension 1950
      real*8        Dec_5      (nr__5) ! (deg)     Declination 1950
      integer*4     Vol_5      (nr__5) ! [5] Volume number
      integer*4     CGMW_5     (nr__5) ! Running number in volume
      character*10  IRAS_1     (nr__5) ! IRAS (Cat. <II/125>) name
      integer*4     RAh_5      (nr__5) ! (h) Right ascension (1950)
      integer*4     RAm_5      (nr__5) ! (min) Right ascension (1950)
      real*4        RAs_5      (nr__5) ! (s) Right ascension (1950)
      character*1   DE__5      (nr__5) ! Declination sign (1950)
      integer*4     DEd_5      (nr__5) ! (deg) Declination (1950)
      integer*4     DEm_5      (nr__5) ! (arcmin) Declination (1950)
      integer*4     DEs_2      (nr__5) ! (arcsec) Declination (1950)
      real*4        Size       (nr__5) ! (mm) Size (1mm~67")
      integer*4     Sep        (nr__5) ! (arcsec) ? Separation
      real*4        F60um      (nr__5) ! (Jy) Flux density at 60{mu}m
      character*1   n_F60um    (nr__5) ! [:L] Note (1)
      real*4        F100um     (nr__5) ! (Jy) Flux density at 100{mu}m
      character*1   n_F100um   (nr__5) ! [:L] Note (1)
      character*5   Cat        (nr__5) ! Catalog names for known galaxies listed
*                                   in Table 3
      integer*4     cz_1       (nr__5) ! (km/s) ? Measured heliocentric radial velocity
      integer*4     r_cz_1     (nr__5) ! [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(i__),CGMW(i__),Field(i__),Xpos(i__),Ypos(i__),RAh(i__),
     +  RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),GLON(i__),
     +  GLAT(i__),MajDiam(i__),MinDiam(i__),MType(i__),Feature(i__),
     +  CrossId(i__)
c       Derive coordinates RA_ and Dec from input data
c       (RA_ and Dec are set to rNULL_ when unknown)
        RA_(i__) = RAh(i__)
        if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
        if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
        if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
        Dec(i__) = DEd(i__)
        if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
        if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c    ..............Just test output...........
        write(6,1)
     +  Vol(i__),CGMW(i__),Field(i__),Xpos(i__),Ypos(i__),RAh(i__),
     +  RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),GLON(i__),
     +  GLAT(i__),MajDiam(i__),MinDiam(i__),MType(i__),Feature(i__),
     +  CrossId(i__)
        write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
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(i__),CGMW_1(i__),Field_1(i__),Xpos_1(i__),Ypos_1(i__),
     +  RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),
     +  DEm_1(i__),GLON_1(i__),GLAT_1(i__),MajDiam_1(i__),
     +  MinDiam_1(i__),MType_1(i__),Feature_1(i__),CrossId_1(i__)
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(i__) = RAh_1(i__)
        if(RAh_1(i__) .GE. 0) RA__1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAs_1(i__)/240.
        Dec_1(i__) = DEd_1(i__)
        if(DEm_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEm_1(i__)/60.
        if(DE__1(i__).EQ.'-'.AND.Dec_1(i__).GE.0) Dec_1(i__)=-Dec_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Vol_1(i__),CGMW_1(i__),Field_1(i__),Xpos_1(i__),Ypos_1(i__),
     +  RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),
     +  DEm_1(i__),GLON_1(i__),GLAT_1(i__),MajDiam_1(i__),
     +  MinDiam_1(i__),MType_1(i__),Feature_1(i__),CrossId_1(i__)
        write(6,'(6H Pos: 2F8.4)') RA__1(i__),Dec_1(i__)
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(i__),CGMW_2(i__),Plate(i__),RefSurv(i__),Over1(i__),
     +  RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),
     +  DEm_2(i__),GLON_2(i__),GLAT_2(i__),MajDiam_2(i__),
     +  MinDiam_2(i__),Feature_2(i__),CrossId_2(i__)
        if(ar__2(24:27) .EQ. '') Over1(i__) = 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(i__) = RAh_2(i__)
        if(RAh_2(i__) .GE. 0) RA__2(i__)=RAh_2(i__)*15.
        if(RAm_2(i__) .GE. 0) RA__2(i__)=RA__2(i__)+RAm_2(i__)/4.
        if(RAs_2(i__) .GE. 0) RA__2(i__)=RA__2(i__)+RAs_2(i__)/240.
        Dec_2(i__) = DEd_2(i__)
        if(DEm_2(i__) .GE. 0) Dec_2(i__)=Dec_2(i__)+DEm_2(i__)/60.
        if(DE__2(i__).EQ.'-'.AND.Dec_2(i__).GE.0) Dec_2(i__)=-Dec_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  Vol_2(i__),CGMW_2(i__),Plate(i__),RefSurv(i__),Over1(i__),
     +  RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),
     +  DEm_2(i__),GLON_2(i__),GLAT_2(i__),MajDiam_2(i__),
     +  MinDiam_2(i__),Feature_2(i__),CrossId_2(i__)
        write(6,'(6H Pos: 2F8.4)') RA__2(i__),Dec_2(i__)
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(i__),CGMW_3(i__),Plate_1(i__),RefSurv_1(i__),
     +  Over1_1(i__),Over2(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),
     +  DE__3(i__),DEd_3(i__),DEm_3(i__),DEs(i__),GLON_3(i__),
     +  GLAT_3(i__),MajDiam_3(i__),MinDiam_3(i__),Feature_3(i__),
     +  CrossId_3(i__)
        if(ar__3(24:27) .EQ. '') Over1_1(i__) = iNULL__
        if(ar__3(29:32) .EQ. '') Over2(i__) = 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(i__) = RAh_3(i__)
        if(RAh_3(i__) .GE. 0) RA__3(i__)=RAh_3(i__)*15.
        if(RAm_3(i__) .GE. 0) RA__3(i__)=RA__3(i__)+RAm_3(i__)/4.
        if(RAs_3(i__) .GE. 0) RA__3(i__)=RA__3(i__)+RAs_3(i__)/240.
        Dec_3(i__) = DEd_3(i__)
        if(DEm_3(i__) .GE. 0) Dec_3(i__)=Dec_3(i__)+DEm_3(i__)/60.
        if(DEs(i__) .GE. 0) Dec_3(i__)=Dec_3(i__)+DEs(i__)/3600.
        if(DE__3(i__).EQ.'-'.AND.Dec_3(i__).GE.0) Dec_3(i__)=-Dec_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  Vol_3(i__),CGMW_3(i__),Plate_1(i__),RefSurv_1(i__),
     +  Over1_1(i__),Over2(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),
     +  DE__3(i__),DEd_3(i__),DEm_3(i__),DEs(i__),GLON_3(i__),
     +  GLAT_3(i__),MajDiam_3(i__),MinDiam_3(i__),Feature_3(i__),
     +  CrossId_3(i__)
        write(6,'(6H Pos: 2F8.4)') RA__3(i__),Dec_3(i__)
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(i__),CGMW_4(i__),II_B(i__),n_II_B(i__),Plate_2(i__),
     +  RA1h(i__),RA1m(i__),RA1s(i__),DE1_(i__),DE1d(i__),DE1m(i__),
     +  DE1s(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),DE__4(i__),
     +  DEd_4(i__),DEm_4(i__),DEs_1(i__),GLON_4(i__),GLAT_4(i__),
     +  MajDiam_4(i__),MinDiam_4(i__),Feature_4(i__),CrossId_4(i__),
     +  cz(i__),r_cz(i__),IRAS(i__)
        if(ar__4(84:86) .EQ. '') MinDiam_4(i__) = rNULL__
        if(ar__4(125:129) .EQ. '') cz(i__) = iNULL__
        if(ar__4(133:133) .EQ. '') r_cz(i__) = 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(i__) = RAh_4(i__)
        if(RAh_4(i__) .GE. 0) RA__4(i__)=RAh_4(i__)*15.
        if(RAm_4(i__) .GE. 0) RA__4(i__)=RA__4(i__)+RAm_4(i__)/4.
        if(RAs_4(i__) .GE. 0) RA__4(i__)=RA__4(i__)+RAs_4(i__)/240.
        Dec_4(i__) = DEd_4(i__)
        if(DEm_4(i__) .GE. 0) Dec_4(i__)=Dec_4(i__)+DEm_4(i__)/60.
        if(DEs_1(i__) .GE. 0) Dec_4(i__)=Dec_4(i__)+DEs_1(i__)/3600.
        if(DE__4(i__).EQ.'-'.AND.Dec_4(i__).GE.0) Dec_4(i__)=-Dec_4(i__)
c    ..............Just test output...........
        write(6,5)
     +  Vol_4(i__),CGMW_4(i__),II_B(i__),n_II_B(i__),Plate_2(i__),
     +  RA1h(i__),RA1m(i__),RA1s(i__),DE1_(i__),DE1d(i__),DE1m(i__),
     +  DE1s(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),DE__4(i__),
     +  DEd_4(i__),DEm_4(i__),DEs_1(i__),GLON_4(i__),GLAT_4(i__),
     +  MajDiam_4(i__),MinDiam_4(i__),Feature_4(i__),CrossId_4(i__),
     +  cz(i__),r_cz(i__),IRAS(i__)
        write(6,'(6H Pos: 2F8.4)') RA__4(i__),Dec_4(i__)
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(i__),CGMW_5(i__),IRAS_1(i__),RAh_5(i__),RAm_5(i__),
     +  RAs_5(i__),DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_2(i__),
     +  Size(i__),Sep(i__),F60um(i__),n_F60um(i__),F100um(i__),
     +  n_F100um(i__),Cat(i__),cz_1(i__),r_cz_1(i__)
        if(ar__5(45:46) .EQ. '') Sep(i__) = iNULL__
        if(ar__5(69:73) .EQ. '') cz_1(i__) = iNULL__
        if(ar__5(75:76) .EQ. '') r_cz_1(i__) = 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(i__) = RAh_5(i__)
        if(RAh_5(i__) .GE. 0) RA__5(i__)=RAh_5(i__)*15.
        if(RAm_5(i__) .GE. 0) RA__5(i__)=RA__5(i__)+RAm_5(i__)/4.
        if(RAs_5(i__) .GE. 0) RA__5(i__)=RA__5(i__)+RAs_5(i__)/240.
        Dec_5(i__) = DEd_5(i__)
        if(DEm_5(i__) .GE. 0) Dec_5(i__)=Dec_5(i__)+DEm_5(i__)/60.
        if(DEs_2(i__) .GE. 0) Dec_5(i__)=Dec_5(i__)+DEs_2(i__)/3600.
        if(DE__5(i__).EQ.'-'.AND.Dec_5(i__).GE.0) Dec_5(i__)=-Dec_5(i__)
c    ..............Just test output...........
        write(6,6)
     +  Vol_5(i__),CGMW_5(i__),IRAS_1(i__),RAh_5(i__),RAm_5(i__),
     +  RAs_5(i__),DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_2(i__),
     +  Size(i__),Sep(i__),F60um(i__),n_F60um(i__),F100um(i__),
     +  n_F100um(i__),Cat(i__),cz_1(i__),r_cz_1(i__)
        write(6,'(6H Pos: 2F8.4)') RA__5(i__),Dec_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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