Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/914/58 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/ApJ/914/58 Chandra obs. of Planck ESZ & X-ray clusters (Andrade-Santos+, 2021)
*================================================================================
*Chandra observations of the Planck Early Sunyaev-Zeldovich sample: a
*reexamination of masses and mass proxies.
* Andrade-Santos F., Pratt G.W., Melin J.-B., Arnaud M., Jones C.,
* Forman W.R., Pointecouteau E., Bartalucci I., Vikhlinin A., Murray S.S.,
* Mazzotta P., Borgani S., Lovisari L., van Weeren R.J., Kraft R.P.,
* David L.P., Giacintucci S.
* <Astrophys. J., 914, 58 (2021)>
* =2021ApJ...914...58A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Cluster masses using SZ and X-ray proxies
for the clusters in the Planck ESZ (see VIII/88)
sample
integer*4 nr__
parameter (nr__=147) ! Number of records
character*162 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*13 PLCKESZ ! Cluster identifier (GLLL.ll+BB.bb)
integer*4 RAh ! (h) Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*4 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! Sign of the Declination (J2000)
integer*4 DEd ! (deg) Degree of Declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
real*4 z ! [0.017/0.35] Spectroscopic redshift
real*4 M500YX ! (10+14Msun) [2/17.6] Mass derived using Yx
real*4 e_M500YX ! (10+14Msun) [0.2/2.1] Uncertainty in M500YX
real*4 M500PSX ! (10+14Msun) [2/16.3] Mass re-derived using Ysz with
* X-ray center & size
real*4 e_M500PSX ! (10+14Msun) [0.1/0.6] Uncertainty in M500PSX
real*4 M500PSZ2 ! (10+14Msun) [1.9/16.2] Mass from PSZ2 catalog
real*4 e_M500PSZ2 ! (10+14Msun) [0.05/0.6] Uncertainty in M500PSZ2
real*4 YX_Mpc2 ! (Mpc2) [1/44.7] Yx value in Mpc squared
real*4 e_YX_Mpc2 ! (Mpc2) [0.16/8.5] Lower uncertainty in YX-Mpc2
real*4 E_YX_Mpc2_1 ! (Mpc2) [0.17/9.2] Upper uncertainty in YX-Mpc2
real*8 YX_Arc2 ! (arcmin2) [63/13793] Yx value in arcmin squared
real*8 e_YX_Arc2 ! (arcmin2) [10/2109] Lower uncertainty in YX-Arc2
real*8 E_YX_Arc2_1 ! (arcmin2) [11/2257] Upper uncertainty in YX-Arc2
real*4 YSZ_Mpc2 ! (Mpc2) [0.9/41.4] Ysz value in Mpc squared
real*4 e_YSZ_Mpc2 ! (Mpc2) [0.1/1.6] Lower uncertainty in YSZ-Mpc2
real*4 E_YSZ_Mpc2_1 ! (Mpc2) [0.1/1.6] Upper uncertainty in YSZ-Mpc2
real*8 YSZ_Arc2 ! (arcmin2) [65/8451] Ysz value in arcmin squared
real*4 e_YSZ_Arc2 ! (arcmin2) [8.9/300] Lower uncertainty in YSZ-Arc2
real*4 E_YSZ_Arc2_1 ! (arcmin2) [8.9/300] Upper uncertainty in YSZ-Arc2
real*4 rmax ! [0.45/2.11] Maximum radius, R_500_ (1)
*Note (1): Where the emission integral is computed, in units of R_500_.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Cluster masses using SZ and X-ray proxies
for the clusters in the X-ray (see
Voevodkin & Vikhlinin 2004ApJ...601..610V) sample
integer*4 nr__1
parameter (nr__1=80) ! Number of records
character*163 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*14 Name ! Cluster identifier (GLLL.ll+BB.bb)
integer*4 RAh_1 ! (h) Hour of Right Ascension (J2000)
integer*4 RAm_1 ! (min) Minute of Right Ascension (J2000)
real*4 RAs_1 ! (s) Second of Right Ascension (J2000)
character*1 DE__1 ! Sign of the Declination (J2000)
integer*4 DEd_1 ! (deg) Degree of Declination (J2000)
integer*4 DEm_1 ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs_1 ! (arcsec) Arcsecond of Declination (J2000)
real*4 z_1 ! [0.027/0.22] Spectroscopic redshift
real*4 M500YX_1 ! (10+14Msun) [1.2/17.6] Mass derived using Yx
real*4 e_M500YX_1 ! (10+14Msun) [0.1/2.1] Uncertainty in M500YX
real*4 M500PSX_1 ! (10+14Msun) [1/16.3] Mass re-derived using Ysz with
* X-ray center & size
real*4 e_M500PSX_1 ! (10+14Msun) [0.1/0.5] Uncertainty in M500PSX
real*4 M500PSZ2_1 ! (10+14Msun) [1/16.2] Mass from PSZ2 catalog
real*4 e_M500PSZ2_1 ! (10+14Msun) [0.05/0.5] Uncertainty in M500PSZ2
real*4 YX_Mpc2_1 ! (Mpc2) [0.3/43.8] Yx value in Mpc squared
real*4 e_YX_Mpc2_2 ! (Mpc2) [0.06/8.5] Lower uncertainty in YX-Mpc2
real*4 E_YX_Mpc2_3 ! (Mpc2) [0.07/9.2] Upper uncertainty in YX-Mpc2
real*8 YX_Arc2_1 ! (arcmin2) [96/3657] Yx value in arcmin squared
real*8 e_YX_Arc2_2 ! (arcmin2) [15/1038] Lower uncertainty in YX-Arc2
real*8 E_YX_Arc2_3 ! (arcmin2) [17/1183] Upper uncertainty in YX-Arc2
real*4 YSZ_Mpc2_1 ! (Mpc2) [0.34/41.4] Ysz value in Mpc squared
real*4 e_YSZ_Mpc2_2 ! (Mpc2) [0.07/1.2] Lower uncertainty in YSZ-Mpc2
real*4 E_YSZ_Mpc2_3 ! (Mpc2) [0.07/1.2] Upper uncertainty in YSZ-Mpc2
real*8 YSZ_Arc2_1 ! (arcmin2) [95/2174] Ysz value in arcmin squared
real*4 e_YSZ_Arc2_2 ! (arcmin2) [19/127] Lower uncertainty in YSZ-Arc2
real*4 E_YSZ_Arc2_3 ! (arcmin2) [19/127] Upper uncertainty in YSZ-Arc2
real*4 rmax_1 ! [0.76/2.66] Maximum radius, R_500_ (1)
*Note (1): Where the emission integral is computed, in units of R_500_.
C=============================================================================
C Loading file 'table1.dat' ! Cluster masses using SZ and X-ray proxies
* for the clusters in the Planck ESZ (see VIII/88)
* sample
C Format for file interpretation
1 format(
+ A13,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.3,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F4.2,1X,F8.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,
+ F7.2,1X,F6.2,1X,F6.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,147
read(1,'(A162)')ar__
read(ar__,1)
+ PLCKESZ,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,M500YX,e_M500YX,M500PSX,
+ e_M500PSX,M500PSZ2,e_M500PSZ2,YX_Mpc2,e_YX_Mpc2,E_YX_Mpc2_1,
+ YX_Arc2,e_YX_Arc2,E_YX_Arc2_1,YSZ_Mpc2,e_YSZ_Mpc2,
+ E_YSZ_Mpc2_1,YSZ_Arc2,e_YSZ_Arc2,E_YSZ_Arc2_1,rmax
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ PLCKESZ,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,M500YX,e_M500YX,M500PSX,
+ e_M500PSX,M500PSZ2,e_M500PSZ2,YX_Mpc2,e_YX_Mpc2,E_YX_Mpc2_1,
+ YX_Arc2,e_YX_Arc2,E_YX_Arc2_1,YSZ_Mpc2,e_YSZ_Mpc2,
+ E_YSZ_Mpc2_1,YSZ_Arc2,e_YSZ_Arc2,E_YSZ_Arc2_1,rmax
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Cluster masses using SZ and X-ray proxies
* for the clusters in the X-ray (see
* Voevodkin & Vikhlinin 2004ApJ...601..610V) sample
C Format for file interpretation
2 format(
+ A14,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.3,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,
+ 1X,F4.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,
+ F7.2,1X,F6.2,1X,F6.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,80
read(1,'(A163)')ar__1
read(ar__1,2)
+ Name,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,z_1,M500YX_1,
+ e_M500YX_1,M500PSX_1,e_M500PSX_1,M500PSZ2_1,e_M500PSZ2_1,
+ YX_Mpc2_1,e_YX_Mpc2_2,E_YX_Mpc2_3,YX_Arc2_1,e_YX_Arc2_2,
+ E_YX_Arc2_3,YSZ_Mpc2_1,e_YSZ_Mpc2_2,E_YSZ_Mpc2_3,YSZ_Arc2_1,
+ e_YSZ_Arc2_2,E_YSZ_Arc2_3,rmax_1
RAdeg_1 = rNULL__
DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,2)
+ Name,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,z_1,M500YX_1,
+ e_M500YX_1,M500PSX_1,e_M500PSX_1,M500PSZ2_1,e_M500PSZ2_1,
+ YX_Mpc2_1,e_YX_Mpc2_2,E_YX_Mpc2_3,YX_Arc2_1,e_YX_Arc2_2,
+ E_YX_Arc2_3,YSZ_Mpc2_1,e_YSZ_Mpc2_2,E_YSZ_Mpc2_3,YSZ_Arc2_1,
+ e_YSZ_Arc2_2,E_YSZ_Arc2_3,rmax_1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end