Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-25
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. J/MNRAS/269/151 APM cluster redshift survey (Dalton+, 1994)
*================================================================================
*The APM galaxy survey. IV. Redshifts of rich clusters of galaxies.
* Dalton G.B., Efstathiou G., Maddox S.J., Sutherland W.J.
* <Mon. Not. R. Astron. Soc. 269, 151 (1994)>
* =1994MNRAS.269..151D (SIMBAD/NED Reference)
* +1995MNRAS.273..258D (Erratum)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Clusters for which redshifts were found in
Andernach's compilation
integer*4 nr__
parameter (nr__=36) ! Number of records
character*49 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension 1950
real*8 Dec (nr__) ! (deg) Declination 1950
integer*4 RAh (nr__) ! (h) Right ascension 1950
integer*4 RAm (nr__) ! (min) Right ascension 1950
real*4 RAs (nr__) ! (s) Right ascension 1950
character*1 DE_ (nr__) ! Declination sign
integer*4 DEd (nr__) ! (deg) Declination 1950
integer*4 DEm (nr__) ! (arcmin) Declination 1950
integer*4 DEs (nr__) ! (arcsec) Declination 1950
real*4 mX (nr__) ! (mag) []? Characteristic magnitude (1)
real*4 R (nr__) ! []? Richness (1)
real*4 z (nr__) ! Redshift
character*5 Abell (nr__) ! Abell name
*Note (1): Those clusters with no entries in mX and R are nearby systems which
* could be missed in the percolation stage of our cluster selection
* procedure
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Clusters observed using AUTOFIB
integer*4 nr__1
parameter (nr__1=23) ! Number of records
character*55 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 RAh_1 (nr__1) ! (h) Right ascension 1950
integer*4 RAm_1 (nr__1) ! (min) Right ascension 1950
real*4 RAs_1 (nr__1) ! (s) Right ascension 1950
character*1 DE__1 (nr__1) ! Declination sign
integer*4 DEd_1 (nr__1) ! (deg) Declination 1950
integer*4 DEm_1 (nr__1) ! (arcmin) Declination 1950
integer*4 DEs_1 (nr__1) ! (arcsec) Declination 1950
real*4 mX_1 (nr__1) ! (mag) []? Characteristic magnitude (1)
real*4 R_1 (nr__1) ! []? Richness (1)
real*4 z_1 (nr__1) ! Redshift
integer*4 o_z (nr__1) ! Number of galaxies for which redshifts
* were obtained
character*5 Abell_1 (nr__1) ! Abell name
character*1 Sep (nr__1) ! [/] when Abell2 exists
character*5 Abell2 (nr__1) ! Second Abell name
*Note (1): Those clusters with no entries in mX and R were included in the
* original selection for the AUTOFIB observations and are included in
* the final list as nearby systems
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Redshifts for galaxies in clusters observed
using a single slit
integer*4 nr__2
parameter (nr__2=365) ! Number of records
character*61 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 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) Right ascension 1950
character*1 DE__2 (nr__2) ! Declination sign
integer*4 DEd_2 (nr__2) ! (deg) Declination 1950
integer*4 DEm_2 (nr__2) ! (arcmin) Declination 1950
real*4 DEs_2 (nr__2) ! (arcsec) Declination 1950
real*4 zObs (nr__2) ! Observed redshift
real*4 zL (nr__2) ! Maximum likelihood redshift estimate
real*4 Proba (nr__2) ! (%) Probability to be a cluster member
character*1 n_Proba (nr__2) ! [*] Note on Proba (1)
real*4 z_2 (nr__2) ! []? Adopted maximum likelihood redshift
real*4 Proba2 (nr__2) ! (%) []? Adopted probability of cluster membership
character*1 n_Proba2 (nr__2) ! [t*] Note on Pint (1)
*Note (1): * : This galaxy is rejected on the basis of the likelihood ratio test
* t : The likelihood function for this cluster is unstable. If the
* galaxy redshifts agree then this is adopted as the cluster
* redshift.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Cluster redshifts obtained using the
likelihood ratio test
integer*4 nr__3
parameter (nr__3=205) ! Number of records
character*55 ar__3 ! Full-size record
real*8 RA__3 (nr__3) ! (deg) Right Ascension 1950
real*8 Dec_3 (nr__3) ! (deg) Declination 1950
integer*4 RAh_3 (nr__3) ! (h) Right ascension 1950
integer*4 RAm_3 (nr__3) ! (min) Right ascension 1950
real*4 RAs_3 (nr__3) ! (s) Right ascension 1950
character*1 DE__3 (nr__3) ! Declination sign
integer*4 DEd_3 (nr__3) ! (deg) Declination 1950
integer*4 DEm_3 (nr__3) ! (arcmin) Declination 1950
integer*4 DEs_3 (nr__3) ! (arcsec) Declination 1950
real*4 mX_2 (nr__3) ! (mag) []? Characteristic magnitude (1)
real*4 R_2 (nr__3) ! []? Richness (1)
real*4 z_3 (nr__3) ! Redshift
integer*4 r_z (nr__3) ! Source of the redshift
character*5 Abell_2 (nr__3) ! Abell name
character*1 Sep_1 (nr__3) ! [/] when Abell2 exists
character*5 Abell2_1 (nr__3) ! Second Abell name
*Note (1): 1: This work (AAT long-slit)
* 2: This work (A.N.U. 2.3m)
* 3: Data taken form Andernach (1991)
* 4: Obtained from Huchra (private communication)
C=============================================================================
C Loading file 'table1.dat' ! Clusters for which redshifts were found in
* Andernach's compilation
C Format for file interpretation
1 format(
+ I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,F5.2,1X,F4.1,2X,F5.3,
+ 5X,A5)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,36
read(1,'(A49)')ar__
read(ar__,1)
+ RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),mX(i__),R(i__),z(i__),Abell(i__)
if(ar__(23:27) .EQ. '') mX(i__) = rNULL__
if(ar__(29:32) .EQ. '') R(i__) = rNULL__
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(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,1)
+ RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),mX(i__),R(i__),z(i__),Abell(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 'table2.dat' ! Clusters observed using AUTOFIB
C Format for file interpretation
2 format(
+ I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,F5.2,1X,F4.1,2X,F5.3,
+ 2X,I2,1X,A5,A1,A5)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,23
read(1,'(A55)')ar__1
read(ar__1,2)
+ RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),
+ DEm_1(i__),DEs_1(i__),mX_1(i__),R_1(i__),z_1(i__),o_z(i__),
+ Abell_1(i__),Sep(i__),Abell2(i__)
if(ar__1(23:27) .EQ. '') mX_1(i__) = rNULL__
if(ar__1(29:32) .EQ. '') R_1(i__) = rNULL__
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(DEs_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.Dec_1(i__).GE.0) Dec_1(i__)=-Dec_1(i__)
c ..............Just test output...........
write(6,2)
+ RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),
+ DEm_1(i__),DEs_1(i__),mX_1(i__),R_1(i__),z_1(i__),o_z(i__),
+ Abell_1(i__),Sep(i__),Abell2(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 'table3.dat' ! Redshifts for galaxies in clusters observed
* using a single slit
C Format for file interpretation
3 format(
+ I2,I2,F5.2,1X,A1,I2,I2,F4.1,3X,F6.4,3X,F6.4,5X,F4.2,A1,1X,
+ F6.4,2X,F4.2,A1)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,365
read(1,'(A61)')ar__2
read(ar__2,3)
+ RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),
+ DEm_2(i__),DEs_2(i__),zObs(i__),zL(i__),Proba(i__),
+ n_Proba(i__),z_2(i__),Proba2(i__),n_Proba2(i__)
if(ar__2(49:54) .EQ. '') z_2(i__) = rNULL__
if(ar__2(57:60) .EQ. '') Proba2(i__) = rNULL__
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(DEs_2(i__) .GE. 0) Dec_2(i__)=Dec_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.Dec_2(i__).GE.0) Dec_2(i__)=-Dec_2(i__)
c ..............Just test output...........
write(6,3)
+ RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),
+ DEm_2(i__),DEs_2(i__),zObs(i__),zL(i__),Proba(i__),
+ n_Proba(i__),z_2(i__),Proba2(i__),n_Proba2(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 'table4.dat' ! Cluster redshifts obtained using the
* likelihood ratio test
C Format for file interpretation
4 format(
+ I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,F5.2,1X,F4.1,2X,F5.3,
+ 2X,I2,1X,A5,A1,A5)
C Effective file loading
open(unit=1,file='table4.dat', status='old')
write(6,*) '....Loading file: table4.dat'
do i__=1,205
read(1,'(A55)')ar__3
read(ar__3,4)
+ RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),DEd_3(i__),
+ DEm_3(i__),DEs_3(i__),mX_2(i__),R_2(i__),z_3(i__),r_z(i__),
+ Abell_2(i__),Sep_1(i__),Abell2_1(i__)
if(ar__3(23:27) .EQ. '') mX_2(i__) = rNULL__
if(ar__3(29:32) .EQ. '') R_2(i__) = rNULL__
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_3(i__) .GE. 0) Dec_3(i__)=Dec_3(i__)+DEs_3(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)
+ RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),DEd_3(i__),
+ DEm_3(i__),DEs_3(i__),mX_2(i__),R_2(i__),z_3(i__),r_z(i__),
+ Abell_2(i__),Sep_1(i__),Abell2_1(i__)
write(6,'(6H Pos: 2F8.4)') RA__3(i__),Dec_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end