Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/392/135 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.81 (2015-09-23), on 2026-Oct-12
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. J/MNRAS/392/135 Galaxy clusters from SDSS (Lopes+, 2009)
*================================================================================
*NoSOCS in SDSS. I. Sample definition and comparison of mass estimates.
* Lopes P.A.A., De Carvalho R.R., Kohl-Moreira J.L., Jones C.
* <Mon. Not. R. Astron. Soc., 392, 135-152 (2009)>
* =2009MNRAS.392..135L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Main properties of the 127 NoSOCS clusters at
z<=0.01 and with at least 10 member galaxies
integer*4 nr__
parameter (nr__=127) ! Number of records
character*108 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 Name ! Name (1)
real*8 RAdeg ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000)
real*4 zhel ! Heliocentric redshift (from CIRS for table3)
real*4 zhel2 ! ? New heliocentric redshift, from SDSS,
* (only for table3)
real*4 Ngal ! Number of galaxies (richness)
real*4 e_Ngal ! rms uncertainty on Ngal
real*4 Lopt ! (10+12Lsun) Optical luminosity
real*4 e_Lopt ! (10+12Lsun) rms uncertainty on Lopt
real*4 beta ! ?=-10.0 Asymetry (beta parameter) (2)
real*8 e_beta ! ?=-99.000 rms uncertainty on beta
real*8 delta ! ?=-9.999 Cumulative deviation
* ({Delta} parameter) (2)
real*4 e_delta ! ?=-9.999 rms uncertainty on delta
*Note (1): NSCS JHHMMSS+DDMMSS in table1.
*Note (2): Parameters obtained from the substructure tests and their
* significance (see Section 5 of paper)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Velocity dispersion, characteristic radii and
masses of 127 NoSOCS clusters
integer*4 nr__1
parameter (nr__1=127) ! Number of records
character*108 ar__1 ! Full-size record
character*19 Name_1 ! Name
real*4 sig ! (km/s) Line-of-sight velocity dispersion
real*4 e_sig ! (km/s) Error on sig (lower limit)
real*4 E_sig_1 ! (km/s) Error on sig (upper limit)
integer*4 Nm ! Number of cluster members used for that
* estimate
integer*4 N200 ! Number of cluster members within R200
real*4 R500 ! (Mpc) Radius R500 (1)
real*4 e_R500 ! (Mpc) Error on R500 (lower limit)
real*4 E_R500_1 ! (Mpc) Error on R500 (upper limit)
real*4 M500 ! (10+14Msun) Mass for a radius of R500
real*4 e_M500 ! (10+14Msun) Error on M500 (lower limit)
real*4 E_M500_1 ! (10+14Msun) Error on M500 (upper limit)
real*4 R200 ! (Mpc) Radius R200 (1)
real*4 e_R200 ! (Mpc) Error on R200 (lower limit)
real*4 E_R200_1 ! (Mpc) Error on R200 (upper limit)
real*4 M200 ! (10+14Msun) Mass for a radius of R200
real*4 e_M200 ! (10+14Msun) Error on M200 (lower limit)
real*4 E_M200_1 ! (10+14Msun) Error on M200 (upper limit)
*Note (1): Radius defined as:
* R200 = R_A_({rho}[200{rho}_c_(z)])^1/2^ and
* R500 = R_A_({rho}_V_[200{rho}_c_(z)])^1/2^ where
* * R_A_ is the Radius
* * {rho}_V_=3M_V_(3{pi}R_A_^3^)
* * M_V_ is the virial mass in a volume of radius R_A_
* * {rho}_c_(z) is the critical density at redshift z
* * M500 is the virial mass for R_A_=R500,
* * M200 is the virial mass for R_A_=R200.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Main properties of the 56 CIRS clusters
(J/AJ/132/1275) at z<=0.01 and with at least
10 member galaxies
integer*4 nr__2
parameter (nr__2=56) ! Number of records
character*108 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 Name_2 ! Name (1)
real*8 RAdeg_1 ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg_1 ! (deg) Declination in decimal degrees (J2000)
real*4 zhel_1 ! Heliocentric redshift (from CIRS for table3)
real*4 zhel2_1 ! ? New heliocentric redshift, from SDSS,
* (only for table3)
real*4 Ngal_1 ! Number of galaxies (richness)
real*4 e_Ngal_1 ! rms uncertainty on Ngal
real*4 Lopt_1 ! (10+12Lsun) Optical luminosity
real*4 e_Lopt_1 ! (10+12Lsun) rms uncertainty on Lopt
real*4 beta_1 ! ?=-10.0 Asymetry (beta parameter) (2)
real*8 e_beta_1 ! ?=-99.000 rms uncertainty on beta
real*8 delta_1 ! ?=-9.999 Cumulative deviation
* ({Delta} parameter) (2)
real*4 e_delta_1 ! ?=-9.999 rms uncertainty on delta
*Note (1): NSCS JHHMMSS+DDMMSS in table1.
*Note (2): Parameters obtained from the substructure tests and their
* significance (see Section 5 of paper)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Velocity dispersion, characteristic radii and
masses of the 56 CIRS clusters
integer*4 nr__3
parameter (nr__3=56) ! Number of records
character*108 ar__3 ! Full-size record
character*19 Name_3 ! Name
real*4 sig_1 ! (km/s) Line-of-sight velocity dispersion
real*4 e_sig_2 ! (km/s) Error on sig (lower limit)
real*4 E_sig_3 ! (km/s) Error on sig (upper limit)
integer*4 Nm_1 ! Number of cluster members used for that
* estimate
integer*4 N200_1 ! Number of cluster members within R200
real*4 R500_1 ! (Mpc) Radius R500 (1)
real*4 e_R500_2 ! (Mpc) Error on R500 (lower limit)
real*4 E_R500_3 ! (Mpc) Error on R500 (upper limit)
real*4 M500_1 ! (10+14Msun) Mass for a radius of R500
real*4 e_M500_2 ! (10+14Msun) Error on M500 (lower limit)
real*4 E_M500_3 ! (10+14Msun) Error on M500 (upper limit)
real*4 R200_1 ! (Mpc) Radius R200 (1)
real*4 e_R200_2 ! (Mpc) Error on R200 (lower limit)
real*4 E_R200_3 ! (Mpc) Error on R200 (upper limit)
real*4 M200_1 ! (10+14Msun) Mass for a radius of R200
real*4 e_M200_2 ! (10+14Msun) Error on M200 (lower limit)
real*4 E_M200_3 ! (10+14Msun) Error on M200 (upper limit)
*Note (1): Radius defined as:
* R200 = R_A_({rho}[200{rho}_c_(z)])^1/2^ and
* R500 = R_A_({rho}_V_[200{rho}_c_(z)])^1/2^ where
* * R_A_ is the Radius
* * {rho}_V_=3M_V_(3{pi}R_A_^3^)
* * M_V_ is the virial mass in a volume of radius R_A_
* * {rho}_c_(z) is the critical density at redshift z
* * M500 is the virial mass for R_A_=R500,
* * M200 is the virial mass for R_A_=R200.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! *X-ray luminosity and temperature of the
NoSOCS plus CIRS clusters
integer*4 nr__4
parameter (nr__4=183) ! Number of records
character*87 ar__4 ! Full-size record
character*19 Name_4 ! Cluster name
real*4 LXa ! (10+37W) ?=99.999 Annulus background X-ray luminosity
real*4 e_LXa ! (10+37W) ?=99.999 rms uncertainty on LXa
real*4 LXf ! (10+37W) ?=99.999 Frame background X-ray luminosity
real*4 e_LXf ! (10+37W) ?=99.999 rms uncertainty on LXf
real*4 LXb ! (10+37W) ?=99.999 Box background X-ray luminosity
real*4 e_LXb ! (10+37W) ?=99.999 rms uncertainty on LXb
real*4 TXB ! (keV) ?=- BAX X-ray temperature (1)
real*4 e_TXB ! (keV) ?=- Error on TXB (lower limit)
real*4 E_TXB_1 ! (keV) ?=- Error on TXB (upper limit)
real*4 TX ! (keV) ?=- Interpolated RASS X-ray temperature
* obtained from the L_X_-T_X_ relation
character*2 Q ! Detection flag (2)
*Note (1): BAX: X-ray galaxy clusters database, http://bax.ast.obs-mip.fr/
*Note (2): Detection flag as follows:
* SD = significant detection
* UL = upper limit
* XX = measurements should be taken with concern
C=============================================================================
C Loading file 'table1.dat' ! Main properties of the 127 NoSOCS clusters at
* z<=0.01 and with at least 10 member galaxies
C Format for file interpretation
1 format(
+ A19,1X,F9.5,1X,F9.5,1X,F6.4,1X,F6.4,1X,F5.1,1X,F4.1,1X,F5.3,
+ 1X,F5.3,1X,F6.1,3X,F7.3,1X,F7.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,127
read(1,'(A108)')ar__
read(ar__,1)
+ Name,RAdeg,DEdeg,zhel,zhel2,Ngal,e_Ngal,Lopt,e_Lopt,beta,
+ e_beta,delta,e_delta
if(ar__(48:53) .EQ. '') zhel2 = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,RAdeg,DEdeg,zhel,zhel2,Ngal,e_Ngal,Lopt,e_Lopt,beta,
+ e_beta,delta,e_delta
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Velocity dispersion, characteristic radii and
* masses of 127 NoSOCS clusters
C Format for file interpretation
2 format(
+ A19,1X,F6.1,1X,F5.1,1X,F5.1,1X,I3,1X,I3,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,127
read(1,'(A108)')ar__1
read(ar__1,2)
+ Name_1,sig,e_sig,E_sig_1,Nm,N200,R500,e_R500,E_R500_1,M500,
+ e_M500,E_M500_1,R200,e_R200,E_R200_1,M200,e_M200,E_M200_1
c ..............Just test output...........
write(6,2)
+ Name_1,sig,e_sig,E_sig_1,Nm,N200,R500,e_R500,E_R500_1,M500,
+ e_M500,E_M500_1,R200,e_R200,E_R200_1,M200,e_M200,E_M200_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Main properties of the 56 CIRS clusters
* (J/AJ/132/1275) at z<=0.01 and with at least
* 10 member galaxies
C Format for file interpretation
3 format(
+ A19,1X,F9.5,1X,F9.5,1X,F6.4,1X,F6.4,1X,F5.1,1X,F4.1,1X,F5.3,
+ 1X,F5.3,1X,F6.1,3X,F7.3,1X,F7.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,56
read(1,'(A108)')ar__2
read(ar__2,3)
+ Name_2,RAdeg_1,DEdeg_1,zhel_1,zhel2_1,Ngal_1,e_Ngal_1,Lopt_1,
+ e_Lopt_1,beta_1,e_beta_1,delta_1,e_delta_1
if(ar__2(48:53) .EQ. '') zhel2_1 = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,RAdeg_1,DEdeg_1,zhel_1,zhel2_1,Ngal_1,e_Ngal_1,Lopt_1,
+ e_Lopt_1,beta_1,e_beta_1,delta_1,e_delta_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Velocity dispersion, characteristic radii and
* masses of the 56 CIRS clusters
C Format for file interpretation
4 format(
+ A19,1X,F6.1,1X,F5.1,1X,F5.1,1X,I3,1X,I3,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,56
read(1,'(A108)')ar__3
read(ar__3,4)
+ Name_3,sig_1,e_sig_2,E_sig_3,Nm_1,N200_1,R500_1,e_R500_2,
+ E_R500_3,M500_1,e_M500_2,E_M500_3,R200_1,e_R200_2,E_R200_3,
+ M200_1,e_M200_2,E_M200_3
c ..............Just test output...........
write(6,4)
+ Name_3,sig_1,e_sig_2,E_sig_3,Nm_1,N200_1,R500_1,e_R500_2,
+ E_R500_3,M500_1,e_M500_2,E_M500_3,R200_1,e_R200_2,E_R200_3,
+ M200_1,e_M200_2,E_M200_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! *X-ray luminosity and temperature of the
* NoSOCS plus CIRS clusters
C Format for file interpretation
5 format(
+ A19,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.2,
+ 2X,F4.2,2X,F4.2,1X,F5.2,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,183
read(1,'(A87)')ar__4
read(ar__4,5)
+ Name_4,LXa,e_LXa,LXf,e_LXf,LXb,e_LXb,TXB,e_TXB,E_TXB_1,TX,Q
if (idig(ar__4(63:67)).EQ.0) TXB = rNULL__
if (idig(ar__4(70:73)).EQ.0) e_TXB = rNULL__
if (idig(ar__4(76:79)).EQ.0) E_TXB_1 = rNULL__
if (idig(ar__4(81:85)).EQ.0) TX = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4,LXa,e_LXa,LXf,e_LXf,LXb,e_LXb,TXB,e_TXB,E_TXB_1,TX,Q
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end