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 2013-May-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__
ter (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