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.81 (2015-09-23), on 2024-Apr-20
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. 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
C 1950 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg (nr__) ! (deg) Right Ascension 1950
real*8 DEdeg (nr__) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
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
C 1950 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension 1950
real*8 DEdeg_1 (nr__1) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
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
C 1950 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension 1950
real*8 DEdeg_2 (nr__2) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
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
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension J2000
real*8 DEdeg_3 (nr__3) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
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
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_4 (nr__4) ! (deg) Right Ascension J2000
real*8 DEdeg_4 (nr__4) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
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
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_5 (nr__5) ! (deg) Right Ascension 1950
real*8 DEdeg_5 (nr__5) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
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,status='old',file=
+'cgmw1.dat')
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__)
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(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(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)') RAdeg(i__),DEdeg(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,status='old',file=
+'cgmw2.dat')
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__)
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(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_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)') RAdeg_1(i__),DEdeg_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,status='old',file=
+'cgmw3.dat')
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__
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(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_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)') RAdeg_2(i__),DEdeg_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,status='old',file=
+'cgmw4.dat')
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__
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(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs(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,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)') RAdeg_3(i__),DEdeg_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,status='old',file=
+'cgmw5.dat')
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__
RAdeg_4(i__) = rNULL__
DEdeg_4(i__) = rNULL__
c Derive coordinates RAdeg_4 and DEdeg_4 from input data
c (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown)
if(RAh_4(i__) .GT. -180) RAdeg_4(i__)=RAh_4(i__)*15.
if(RAm_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAm_4(i__)/4.
if(RAs_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAs_4(i__)/240.
if(DEd_4(i__) .GE. 0) DEdeg_4(i__)=DEd_4(i__)
if(DEm_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm_4(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs_1(i__)/3600.
if(DE__4(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_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)') RAdeg_4(i__),DEdeg_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,status='old',file=
+'cgmw5ir.dat')
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__
RAdeg_5(i__) = rNULL__
DEdeg_5(i__) = rNULL__
c Derive coordinates RAdeg_5 and DEdeg_5 from input data
c (RAdeg_5 and DEdeg_5 are set to rNULL__ when unknown)
if(RAh_5(i__) .GT. -180) RAdeg_5(i__)=RAh_5(i__)*15.
if(RAm_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAm_5(i__)/4.
if(RAs_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAs_5(i__)/240.
if(DEd_5(i__) .GE. 0) DEdeg_5(i__)=DEd_5(i__)
if(DEm_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEm_5(i__)/60.
if(DEs_2(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEs_2(i__)/3600.
if(DE__5(i__).EQ.'-'.AND.DEdeg_5(i__).GE.0) DEdeg_5(i__)=-DEdeg_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)') RAdeg_5(i__),DEdeg_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end