Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/914/58 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/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 (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*13 PLCKESZ (nr__) ! Cluster identifier (GLLL.ll+BB.bb)
integer*4 RAh (nr__) ! (h) Hour of Right Ascension (J2000)
integer*4 RAm (nr__) ! (min) Minute of Right Ascension (J2000)
real*4 RAs (nr__) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__) ! Sign of the Declination (J2000)
integer*4 DEd (nr__) ! (deg) Degree of Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Arcsecond of Declination (J2000)
real*4 z (nr__) ! [0.017/0.35] Spectroscopic redshift
real*4 M500YX (nr__) ! (10+14Msun) [2/17.6] Mass derived using Yx
real*4 e_M500YX (nr__) ! (10+14Msun) [0.2/2.1] Uncertainty in M500YX
real*4 M500PSX (nr__) ! (10+14Msun) [2/16.3] Mass re-derived using Ysz with
* X-ray center & size
real*4 e_M500PSX (nr__) ! (10+14Msun) [0.1/0.6] Uncertainty in M500PSX
real*4 M500PSZ2 (nr__) ! (10+14Msun) [1.9/16.2] Mass from PSZ2 catalog
real*4 e_M500PSZ2 (nr__) ! (10+14Msun) [0.05/0.6] Uncertainty in M500PSZ2
real*4 YX_Mpc2 (nr__) ! (Mpc2) [1/44.7] Yx value in Mpc squared
real*4 e_YX_Mpc2 (nr__) ! (Mpc2) [0.16/8.5] Lower uncertainty in YX-Mpc2
real*4 E_YX_Mpc2_1(nr__) ! (Mpc2) [0.17/9.2] Upper uncertainty in YX-Mpc2
real*8 YX_Arc2 (nr__) ! (arcmin2) [63/13793] Yx value in arcmin squared
real*8 e_YX_Arc2 (nr__) ! (arcmin2) [10/2109] Lower uncertainty in YX-Arc2
real*8 E_YX_Arc2_1(nr__) ! (arcmin2) [11/2257] Upper uncertainty in YX-Arc2
real*4 YSZ_Mpc2 (nr__) ! (Mpc2) [0.9/41.4] Ysz value in Mpc squared
real*4 e_YSZ_Mpc2 (nr__) ! (Mpc2) [0.1/1.6] Lower uncertainty in YSZ-Mpc2
real*4 E_YSZ_Mpc2_1(nr__) ! (Mpc2) [0.1/1.6] Upper uncertainty in YSZ-Mpc2
real*8 YSZ_Arc2 (nr__) ! (arcmin2) [65/8451] Ysz value in arcmin squared
real*4 e_YSZ_Arc2 (nr__) ! (arcmin2) [8.9/300] Lower uncertainty in YSZ-Arc2
real*4 E_YSZ_Arc2_1(nr__) ! (arcmin2) [8.9/300] Upper uncertainty in YSZ-Arc2
real*4 rmax (nr__) ! [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 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*14 Name (nr__1) ! Cluster identifier (GLLL.ll+BB.bb)
integer*4 RAh_1 (nr__1) ! (h) Hour of Right Ascension (J2000)
integer*4 RAm_1 (nr__1) ! (min) Minute of Right Ascension (J2000)
real*4 RAs_1 (nr__1) ! (s) Second of Right Ascension (J2000)
character*1 DE__1 (nr__1) ! Sign of the Declination (J2000)
integer*4 DEd_1 (nr__1) ! (deg) Degree of Declination (J2000)
integer*4 DEm_1 (nr__1) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs_1 (nr__1) ! (arcsec) Arcsecond of Declination (J2000)
real*4 z_1 (nr__1) ! [0.027/0.22] Spectroscopic redshift
real*4 M500YX_1 (nr__1) ! (10+14Msun) [1.2/17.6] Mass derived using Yx
real*4 e_M500YX_1 (nr__1) ! (10+14Msun) [0.1/2.1] Uncertainty in M500YX
real*4 M500PSX_1 (nr__1) ! (10+14Msun) [1/16.3] Mass re-derived using Ysz with
* X-ray center & size
real*4 e_M500PSX_1(nr__1) ! (10+14Msun) [0.1/0.5] Uncertainty in M500PSX
real*4 M500PSZ2_1 (nr__1) ! (10+14Msun) [1/16.2] Mass from PSZ2 catalog
real*4 e_M500PSZ2_1(nr__1) ! (10+14Msun) [0.05/0.5] Uncertainty in M500PSZ2
real*4 YX_Mpc2_1 (nr__1) ! (Mpc2) [0.3/43.8] Yx value in Mpc squared
real*4 e_YX_Mpc2_2(nr__1) ! (Mpc2) [0.06/8.5] Lower uncertainty in YX-Mpc2
real*4 E_YX_Mpc2_3(nr__1) ! (Mpc2) [0.07/9.2] Upper uncertainty in YX-Mpc2
real*8 YX_Arc2_1 (nr__1) ! (arcmin2) [96/3657] Yx value in arcmin squared
real*8 e_YX_Arc2_2(nr__1) ! (arcmin2) [15/1038] Lower uncertainty in YX-Arc2
real*8 E_YX_Arc2_3(nr__1) ! (arcmin2) [17/1183] Upper uncertainty in YX-Arc2
real*4 YSZ_Mpc2_1 (nr__1) ! (Mpc2) [0.34/41.4] Ysz value in Mpc squared
real*4 e_YSZ_Mpc2_2(nr__1) ! (Mpc2) [0.07/1.2] Lower uncertainty in YSZ-Mpc2
real*4 E_YSZ_Mpc2_3(nr__1) ! (Mpc2) [0.07/1.2] Upper uncertainty in YSZ-Mpc2
real*8 YSZ_Arc2_1 (nr__1) ! (arcmin2) [95/2174] Ysz value in arcmin squared
real*4 e_YSZ_Arc2_2(nr__1) ! (arcmin2) [19/127] Lower uncertainty in YSZ-Arc2
real*4 E_YSZ_Arc2_3(nr__1) ! (arcmin2) [19/127] Upper uncertainty in YSZ-Arc2
real*4 rmax_1 (nr__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(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),z(i__),M500YX(i__),e_M500YX(i__),
+ M500PSX(i__),e_M500PSX(i__),M500PSZ2(i__),e_M500PSZ2(i__),
+ YX_Mpc2(i__),e_YX_Mpc2(i__),E_YX_Mpc2_1(i__),YX_Arc2(i__),
+ e_YX_Arc2(i__),E_YX_Arc2_1(i__),YSZ_Mpc2(i__),e_YSZ_Mpc2(i__),
+ E_YSZ_Mpc2_1(i__),YSZ_Arc2(i__),e_YSZ_Arc2(i__),
+ E_YSZ_Arc2_1(i__),rmax(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ PLCKESZ(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),z(i__),M500YX(i__),e_M500YX(i__),
+ M500PSX(i__),e_M500PSX(i__),M500PSZ2(i__),e_M500PSZ2(i__),
+ YX_Mpc2(i__),e_YX_Mpc2(i__),E_YX_Mpc2_1(i__),YX_Arc2(i__),
+ e_YX_Arc2(i__),E_YX_Arc2_1(i__),YSZ_Mpc2(i__),e_YSZ_Mpc2(i__),
+ E_YSZ_Mpc2_1(i__),YSZ_Arc2(i__),e_YSZ_Arc2(i__),
+ E_YSZ_Arc2_1(i__),rmax(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),z_1(i__),M500YX_1(i__),
+ e_M500YX_1(i__),M500PSX_1(i__),e_M500PSX_1(i__),
+ M500PSZ2_1(i__),e_M500PSZ2_1(i__),YX_Mpc2_1(i__),
+ e_YX_Mpc2_2(i__),E_YX_Mpc2_3(i__),YX_Arc2_1(i__),
+ e_YX_Arc2_2(i__),E_YX_Arc2_3(i__),YSZ_Mpc2_1(i__),
+ e_YSZ_Mpc2_2(i__),E_YSZ_Mpc2_3(i__),YSZ_Arc2_1(i__),
+ e_YSZ_Arc2_2(i__),E_YSZ_Arc2_3(i__),rmax_1(i__)
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ Name(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),z_1(i__),M500YX_1(i__),
+ e_M500YX_1(i__),M500PSX_1(i__),e_M500PSX_1(i__),
+ M500PSZ2_1(i__),e_M500PSZ2_1(i__),YX_Mpc2_1(i__),
+ e_YX_Mpc2_2(i__),E_YX_Mpc2_3(i__),YX_Arc2_1(i__),
+ e_YX_Arc2_2(i__),E_YX_Arc2_3(i__),YSZ_Mpc2_1(i__),
+ e_YSZ_Mpc2_2(i__),E_YSZ_Mpc2_3(i__),YSZ_Arc2_1(i__),
+ e_YSZ_Arc2_2(i__),E_YSZ_Arc2_3(i__),rmax_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end