Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/392/135 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 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 (nr__) ! Name (1)
real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
real*4 zhel (nr__) ! Heliocentric redshift (from CIRS for table3)
real*4 zhel2 (nr__) ! ? New heliocentric redshift, from SDSS,
* (only for table3)
real*4 Ngal (nr__) ! Number of galaxies (richness)
real*4 e_Ngal (nr__) ! rms uncertainty on Ngal
real*4 Lopt (nr__) ! (10+12Lsun) Optical luminosity
real*4 e_Lopt (nr__) ! (10+12Lsun) rms uncertainty on Lopt
real*4 beta (nr__) ! ?=-10.0 Asymetry (beta parameter) (2)
real*8 e_beta (nr__) ! ?=-99.000 rms uncertainty on beta
real*8 delta (nr__) ! ?=-9.999 Cumulative deviation
* ({Delta} parameter) (2)
real*4 e_delta (nr__) ! ?=-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 (nr__1) ! Name
real*4 sig (nr__1) ! (km/s) Line-of-sight velocity dispersion
real*4 e_sig (nr__1) ! (km/s) Error on sig (lower limit)
real*4 E_sig_1 (nr__1) ! (km/s) Error on sig (upper limit)
integer*4 Nm (nr__1) ! Number of cluster members used for that
* estimate
integer*4 N200 (nr__1) ! Number of cluster members within R200
real*4 R500 (nr__1) ! (Mpc) Radius R500 (1)
real*4 e_R500 (nr__1) ! (Mpc) Error on R500 (lower limit)
real*4 E_R500_1 (nr__1) ! (Mpc) Error on R500 (upper limit)
real*4 M500 (nr__1) ! (10+14Msun) Mass for a radius of R500
real*4 e_M500 (nr__1) ! (10+14Msun) Error on M500 (lower limit)
real*4 E_M500_1 (nr__1) ! (10+14Msun) Error on M500 (upper limit)
real*4 R200 (nr__1) ! (Mpc) Radius R200 (1)
real*4 e_R200 (nr__1) ! (Mpc) Error on R200 (lower limit)
real*4 E_R200_1 (nr__1) ! (Mpc) Error on R200 (upper limit)
real*4 M200 (nr__1) ! (10+14Msun) Mass for a radius of R200
real*4 e_M200 (nr__1) ! (10+14Msun) Error on M200 (lower limit)
real*4 E_M200_1 (nr__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 (nr__2) ! Name (1)
real*8 RAdeg_1 (nr__2) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__2) ! (deg) Declination in decimal degrees (J2000)
real*4 zhel_1 (nr__2) ! Heliocentric redshift (from CIRS for table3)
real*4 zhel2_1 (nr__2) ! ? New heliocentric redshift, from SDSS,
* (only for table3)
real*4 Ngal_1 (nr__2) ! Number of galaxies (richness)
real*4 e_Ngal_1 (nr__2) ! rms uncertainty on Ngal
real*4 Lopt_1 (nr__2) ! (10+12Lsun) Optical luminosity
real*4 e_Lopt_1 (nr__2) ! (10+12Lsun) rms uncertainty on Lopt
real*4 beta_1 (nr__2) ! ?=-10.0 Asymetry (beta parameter) (2)
real*8 e_beta_1 (nr__2) ! ?=-99.000 rms uncertainty on beta
real*8 delta_1 (nr__2) ! ?=-9.999 Cumulative deviation
* ({Delta} parameter) (2)
real*4 e_delta_1 (nr__2) ! ?=-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 (nr__3) ! Name
real*4 sig_1 (nr__3) ! (km/s) Line-of-sight velocity dispersion
real*4 e_sig_2 (nr__3) ! (km/s) Error on sig (lower limit)
real*4 E_sig_3 (nr__3) ! (km/s) Error on sig (upper limit)
integer*4 Nm_1 (nr__3) ! Number of cluster members used for that
* estimate
integer*4 N200_1 (nr__3) ! Number of cluster members within R200
real*4 R500_1 (nr__3) ! (Mpc) Radius R500 (1)
real*4 e_R500_2 (nr__3) ! (Mpc) Error on R500 (lower limit)
real*4 E_R500_3 (nr__3) ! (Mpc) Error on R500 (upper limit)
real*4 M500_1 (nr__3) ! (10+14Msun) Mass for a radius of R500
real*4 e_M500_2 (nr__3) ! (10+14Msun) Error on M500 (lower limit)
real*4 E_M500_3 (nr__3) ! (10+14Msun) Error on M500 (upper limit)
real*4 R200_1 (nr__3) ! (Mpc) Radius R200 (1)
real*4 e_R200_2 (nr__3) ! (Mpc) Error on R200 (lower limit)
real*4 E_R200_3 (nr__3) ! (Mpc) Error on R200 (upper limit)
real*4 M200_1 (nr__3) ! (10+14Msun) Mass for a radius of R200
real*4 e_M200_2 (nr__3) ! (10+14Msun) Error on M200 (lower limit)
real*4 E_M200_3 (nr__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 (nr__4) ! Cluster name
real*4 LXa (nr__4) ! (10+37W) ?=99.999 Annulus background X-ray luminosity
real*4 e_LXa (nr__4) ! (10+37W) ?=99.999 rms uncertainty on LXa
real*4 LXf (nr__4) ! (10+37W) ?=99.999 Frame background X-ray luminosity
real*4 e_LXf (nr__4) ! (10+37W) ?=99.999 rms uncertainty on LXf
real*4 LXb (nr__4) ! (10+37W) ?=99.999 Box background X-ray luminosity
real*4 e_LXb (nr__4) ! (10+37W) ?=99.999 rms uncertainty on LXb
real*4 TXB (nr__4) ! (keV) ?=- BAX X-ray temperature (1)
real*4 e_TXB (nr__4) ! (keV) ?=- Error on TXB (lower limit)
real*4 E_TXB_1 (nr__4) ! (keV) ?=- Error on TXB (upper limit)
real*4 TX (nr__4) ! (keV) ?=- Interpolated RASS X-ray temperature
* obtained from the L_X_-T_X_ relation
character*2 Q (nr__4) ! 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(i__),RAdeg(i__),DEdeg(i__),zhel(i__),zhel2(i__),
+ Ngal(i__),e_Ngal(i__),Lopt(i__),e_Lopt(i__),beta(i__),
+ e_beta(i__),delta(i__),e_delta(i__)
if(ar__(48:53) .EQ. '') zhel2(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),zhel(i__),zhel2(i__),
+ Ngal(i__),e_Ngal(i__),Lopt(i__),e_Lopt(i__),beta(i__),
+ e_beta(i__),delta(i__),e_delta(i__)
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(i__),sig(i__),e_sig(i__),E_sig_1(i__),Nm(i__),
+ N200(i__),R500(i__),e_R500(i__),E_R500_1(i__),M500(i__),
+ e_M500(i__),E_M500_1(i__),R200(i__),e_R200(i__),E_R200_1(i__),
+ M200(i__),e_M200(i__),E_M200_1(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),sig(i__),e_sig(i__),E_sig_1(i__),Nm(i__),
+ N200(i__),R500(i__),e_R500(i__),E_R500_1(i__),M500(i__),
+ e_M500(i__),E_M500_1(i__),R200(i__),e_R200(i__),E_R200_1(i__),
+ M200(i__),e_M200(i__),E_M200_1(i__)
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(i__),RAdeg_1(i__),DEdeg_1(i__),zhel_1(i__),
+ zhel2_1(i__),Ngal_1(i__),e_Ngal_1(i__),Lopt_1(i__),
+ e_Lopt_1(i__),beta_1(i__),e_beta_1(i__),delta_1(i__),
+ e_delta_1(i__)
if(ar__2(48:53) .EQ. '') zhel2_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),RAdeg_1(i__),DEdeg_1(i__),zhel_1(i__),
+ zhel2_1(i__),Ngal_1(i__),e_Ngal_1(i__),Lopt_1(i__),
+ e_Lopt_1(i__),beta_1(i__),e_beta_1(i__),delta_1(i__),
+ e_delta_1(i__)
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(i__),sig_1(i__),e_sig_2(i__),E_sig_3(i__),Nm_1(i__),
+ N200_1(i__),R500_1(i__),e_R500_2(i__),E_R500_3(i__),
+ M500_1(i__),e_M500_2(i__),E_M500_3(i__),R200_1(i__),
+ e_R200_2(i__),E_R200_3(i__),M200_1(i__),e_M200_2(i__),
+ E_M200_3(i__)
c ..............Just test output...........
write(6,4)
+ Name_3(i__),sig_1(i__),e_sig_2(i__),E_sig_3(i__),Nm_1(i__),
+ N200_1(i__),R500_1(i__),e_R500_2(i__),E_R500_3(i__),
+ M500_1(i__),e_M500_2(i__),E_M500_3(i__),R200_1(i__),
+ e_R200_2(i__),E_R200_3(i__),M200_1(i__),e_M200_2(i__),
+ E_M200_3(i__)
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(i__),LXa(i__),e_LXa(i__),LXf(i__),e_LXf(i__),LXb(i__),
+ e_LXb(i__),TXB(i__),e_TXB(i__),E_TXB_1(i__),TX(i__),Q(i__)
if (idig(ar__4(63:67)).EQ.0) TXB(i__) = rNULL__
if (idig(ar__4(70:73)).EQ.0) e_TXB(i__) = rNULL__
if (idig(ar__4(76:79)).EQ.0) E_TXB_1(i__) = rNULL__
if (idig(ar__4(81:85)).EQ.0) TX(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4(i__),LXa(i__),e_LXa(i__),LXf(i__),e_LXf(i__),LXb(i__),
+ e_LXb(i__),TXB(i__),e_TXB(i__),E_TXB_1(i__),TX(i__),Q(i__)
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