Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/292/419 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 2024-Apr-18
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/292/419 Cooling flows in 207 clusters of galaxies (White+, 1997)
*================================================================================
*An investigation of cooling flows and general cluster properties from an
*X-ray image deprojection analysis of 207 clusters of galaxies.
* White D.A., Jones C., Forman W.
* <Mon. Not. R. Astron. Soc., 292, 419-467 (1997)>
* =1997MNRAS.292..419W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Cluster reference data
integer*4 nr__
parameter (nr__=260) ! Number of records
character*99 ar__ ! Full-size record
character*12 Name (nr__) ! Cluster name (K38 = Klem 38)
real*4 NH (nr__) ! (10+20cm-2) ? Line-of-sight Galactic hydrogen column
* density
real*4 z (nr__) ! ? Redshift
character*5 r_z (nr__) ! Redshift reference
integer*4 sigma (nr__) ! (km/s) Optical velocity dispersion
character*1 u_sigma (nr__) ! Uncertainty flag on sigma
integer*4 E_sigma (nr__) ! (km/s) ? Error in sigma (upper limit)
integer*4 e_sigma_1 (nr__) ! (km/s) ? Error in sigma (lower limit)
character*5 r_sigma (nr__) ! Sigma reference
real*4 TX (nr__) ! (keV) Average cluster X-ray temperature
character*1 u_TX (nr__) ! Uncertainty flag on TX
real*4 E_TX (nr__) ! (keV) ? Error in TX (upper limit)
real*4 e_TX_1 (nr__) ! (keV) ? Error in TX (lower limit)
character*4 r_TX (nr__) ! Reference for TX (Vopt = sigma+opt)
real*8 LX (nr__) ! (10+37W) X-ray luminosity
character*1 u_LX (nr__) ! Uncertainty flag on LX
real*4 E_LX (nr__) ! (10+37W) ? Error in LX (upper limit)
real*4 e_LX_1 (nr__) ! (10+37W) ? Error in LX (lower limit)
character*4 r_LX (nr__) ! Reference for LX
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__1
parameter (nr__1=53) ! Number of records
character*136 ar__1 ! Full-size record
character*5 Ref (nr__1) ! Reference code
character*19 BibCode (nr__1) ! BibCode
character*19 Aut (nr__1) ! Author's name
character*89 Com (nr__1) ! Comments
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Deprojection input data
integer*4 nr__2
parameter (nr__2=228) ! Number of records
character*128 ar__2 ! Full-size record
integer*4 Seq (nr__2) ! Sequential number
character*12 Name_1 (nr__2) ! Cluster name
character*22 OName (nr__2) ! Alternative names
character*6 Sequence (nr__2) ! Einstein Observatory Sequence (G1)
real*4 IPC1 (nr__2) ! ? IPC energy range of the extracted data
real*4 IPC2 (nr__2) ! ? IPC energy range of the extracted data
real*4 Gain (nr__2) ! ? IPC gain
integer*4 ExpTime (nr__2) ! (s) Exposure time
real*4 z_1 (nr__2) ! Optical redshift
integer*4 DL (nr__2) ! (Mpc) Luminosity distance (H=50km/s/Mpc, q0=0.5)
* of the cluster
real*4 NH_1 (nr__2) ! (10+21cm-2) Galactic hydrogen column density
real*4 kT (nr__2) ! (keV) Spatial average X-ray temperature
integer*4 sigma_1 (nr__2) ! (km/s) Cluster optical velocity dispersion
integer*4 e_sigma_2 (nr__2) ! (km/s) ? rms uncertainty on sigma
real*4 Rcore (nr__2) ! (Mpc) Core radius
integer*4 Nbins (nr__2) ! Total number of bins
integer*4 dR (nr__2) ! (arcsec) Bin size, in arcsec, of the projection
real*4 dRMpc (nr__2) ! (Mpc) Bin size, in Mpc, of the projection
real*4 Rout (nr__2) ! (Mpc) Outer radius of the projection
real*4 Pout (nr__2) ! (10-4K/cm+3) Gas pressure, at the outer radius
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Deprojection results
integer*4 nr__3
parameter (nr__3=228) ! Number of records
character*162 ar__3 ! Full-size record
integer*4 Seq_1 (nr__3) ! Sequential number
character*12 Name_2 (nr__3) ! Cluster name
character*6 Sequence_1 (nr__3) ! Einstein Observatory Sequence (G1)
real*4 Rcore_1 (nr__3) ! (Mpc) Core radius
real*4 Rout_1 (nr__3) ! (Mpc) Outer radius of the deprojection
real*4 kT_1 (nr__3) ! (K) Spatially averaged, emission-weighted
* reprojected temperature from the
* deprojected profile
real*4 E_kT (nr__3) ! (K) Error in kT (upper limit)
real*4 e_kT_1 (nr__3) ! (K) Error in kT (lower limit)
real*4 kTRef (nr__3) ! (K) ? Reference temperature
real*4 LX_1 (nr__3) ! (10+37W) Spatially averaged bolometric X-ray
* luminosity from the deprojection
real*4 LXRef (nr__3) ! (10+37W) ? Reference X-ray luminosity
real*4 Mgas (nr__3) ! (10+12solMass) Integrated gas mass to the outer radius
real*4 e_Mgas (nr__3) ! (10+12solMass) rms uncertainty on Mgas
real*4 Mgrav (nr__3) ! (10+12solMass) Intergeted gravitational mass at the
* outer radius
real*4 Mgas_Mgrav (nr__3) ! ? gas to gravitational masses ratio
real*4 e_Mgas_Mgrav(nr__3) ! ? rms uncertainty on Mgas/Mgrav
real*4 tcool (nr__3) ! Cooling time (in H0^-1^)
real*4 E_tcool (nr__3) ! Error in tcool (upper limit)
real*4 e_tcool_1 (nr__3) ! Error in tcool (lower limit)
character*1 l_Rcool (nr__3) ! Limit flag on Rcool
integer*4 Rcool (nr__3) ! (kpc) ? Cooling radius
character*1 l_E_Rcool (nr__3) ! Limit flag on E_Rcool
integer*4 E_Rcool (nr__3) ! (kpc) ? Error in Rcool (upper limit)
character*1 l_e_Rcool_1(nr__3) ! Limit flag on e_Rcool
integer*4 e_Rcool_1 (nr__3) ! (kpc) ? Error in Rcool (lower limit)
character*1 l_dM_dt (nr__3) ! Limit flag on dM/dt
real*4 dM_dt (nr__3) ! (solMass/yr) Mass loss of the cooling flow
character*1 l_E_dM_dt (nr__3) ! Limit flag on E_dM/dt
real*4 E_dM_dt (nr__3) ! (solMass/yr) Error on dM/dt (upper limit)
character*1 l_e_dM_dt_1(nr__3) ! Limit flag on e_dM/dt
real*4 e_dM_dt_1 (nr__3) ! (solMass/yr) Error on dM/dt (lower limit)
character*4 CF (nr__3) ! Cooling flow class (1)
real*4 S_Z_dT (nr__3) ! (mK) Projected Sunyaev-Zel'dovich microwave
* decrement from the over 6arcmin radius
real*4 T_Depth (nr__3) ! (%) Thompson depth from Rout to the centre
* of the cluster
*Note (1): Cooling flow classes are defined as follows:
* LCF = large cooling flows, dM/dt>50M_{sun}_/yr
* MCF = moderate cooling flows, 10<dM/dt<=50M_{sun}_/yr
* SCF = small cooling flows, 0<dM/dt<=10M_{sun}_/yr
* PCF = possible cooling flows, {Delta}R>Rcool
* (Rcool is the average size of the cooling flow)
* PLCF = possible large cooling flows
* PMCF = possible moderate cooling flows
* PSCF = possible small cooling flows
* XCF = excluded cooling flows, {Delta}R>Rcool
* (Rcool is the average size of the cooling flow)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Deprojection results summed to different radii
integer*4 nr__4
parameter (nr__4=228) ! Number of records
character*104 ar__4 ! Full-size record
integer*4 Seq_2 (nr__4) ! ? Sequential number
character*12 Name_3 (nr__4) ! Cluster name
character*6 Sequence_2 (nr__4) ! Einstein Observatory Sequence (G1)
real*4 R0_2 (nr__4) ! (Mpc) ? Half-light radius, at the
* reference radius R<0.2Mpc (1)
real*4 R0_5 (nr__4) ! (Mpc) ? Half-light radius, at the
* reference radius R<0.5Mpc
real*4 R1_0 (nr__4) ! (Mpc) ? Half-light radius, at the
* reference radius R<1.0Mpc
real*4 LX0_5 (nr__4) ! (10+37W) ? Integrated value of X-ray
* luminosity, within R<0.5Mpc
real*4 LX1_0 (nr__4) ! (10+37W) ? Integrated value of X-ray
* luminosity, within R<1.0Mpc
real*4 Mgas0_5 (nr__4) ! (10+12solMass) ? Integrated value of the gas mass,
* within R<0.5Mpc
real*4 e_Mgas0_5 (nr__4) ! (10+12solMass) ? rms uncertainty on Mgas0.5
real*4 Mgas1_0 (nr__4) ! (10+12solMass) ? Integrated value of the gas mass,
* within R<1.0Mpc
real*4 e_Mgas1_0 (nr__4) ! (10+12solMass) ? rms uncertainty on Mgas1.0
real*4 Mgrav0_5 (nr__4) ! (10+12solMass) ? Gravitational mass integrated
* value, within R<0.5Mpc
real*4 Mgrav1_0 (nr__4) ! (10+12solMass) ? Gravitational mass integrated
* value, within R<1.0Mpc
real*4 Mgas_Mgrav0_5(nr__4) ! (%) ? Gas-to-gravitational masses ratio,
* within R<0.5Mpc
real*4 e_Mgas_Mgrav0_5(nr__4) ! (%) ? rms uncertainty on Mgas/Mgrav0.5
real*4 Mgas_Mgrav1_0(nr__4) ! (%) ? Gas-to-gravitational masses ratio,
* within R<1.0Mpc
real*4 e_Mgas_Mgrav1_0(nr__4) ! (%) ? rms uncertainty on Mgas/Mgrav1.0
*Note (1): The half-light radius is the radius which contains half of
* the X-ray luminosity at the reference radius 0.2, 0.5 or 1.0Mpc.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Catalogue of cooling flows detected in the
Einstein Observatory deprojections
integer*4 nr__5
parameter (nr__5=97) ! Number of records
character*70 ar__5 ! Full-size record
integer*4 Rank (nr__5) ! Rank of the flow, ordered by
* decreasing dM/dt (1)
character*1 n_Rank (nr__5) ! [b] Note on Rank (1)
integer*4 Seq_3 (nr__5) ! Sequential number
character*12 Name_4 (nr__5) ! Cluster name
character*6 Sequence_3 (nr__5) ! Einstein Observatory Sequence (G1)
character*1 l_Rcool_1 (nr__5) ! Limit flag on Rcool
real*4 Rcool_1 (nr__5) ! (kpc) Cooling radius
character*1 l_E_Rcool_2(nr__5) ! Limit flag on E_Rcool (2)
real*4 E_Rcool_2 (nr__5) ! (kpc) Error on Rcool (upper limit)
character*1 l_e_Rcool_3(nr__5) ! Limit flag on E_Rcool (2)
real*4 e_Rcool_3 (nr__5) ! (kpc) Error on Rcool (lower limit)
character*1 l_dM_dt_1 (nr__5) ! Limit flag on dM/dt (2)
real*4 dM_dt_1 (nr__5) ! (solMass/yr) Mass loss of the cooling flow
character*1 l_E_dM_dt_2(nr__5) ! Limit flag on E_dM/dt (2)
real*4 E_dM_dt_2 (nr__5) ! (solMass/yr) Error on dM/dt (upper limit)
character*1 l_e_dM_dt_3(nr__5) ! Limit flag on e_dM/dt (2)
real*4 e_dM_dt_3 (nr__5) ! (solMass/yr) Error on dM/dt (lower limit)
*Note (1): The entries are ordered by decreasing dM/dt (Rank number),
* and may contain multiple entries for a single cluster if both IPC and
* HRI data are available, and indicate a cooling flow (then n_Rank = b).
*Note (2): The occurence of '>' symbol indicates (usually for HRI results)
* that the outer radius of the deprojection is not large enough that
* tcool increase above t0, and therefore the numbers are quoted at the
* radius Rout.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Supplementary reference data from the literature
integer*4 nr__6
parameter (nr__6=58) ! Number of records
character*53 ar__6 ! Full-size record
character*12 Name_5 (nr__6) ! Cluster name
real*8 LHa (nr__6) ! (10+33W) ? H{alpha} luminosity
real*4 tcool_1 (nr__6) ! ? Cooling time, in Ho^-1^ unit
real*4 E_tcool_2 (nr__6) ! ? Error in tcool (upper limit)
real*4 e_tcool_3 (nr__6) ! ? Error in tcool (lower limit)
integer*4 dM_dt_2 (nr__6) ! (solMass/yr) ? Mass loss
integer*4 E_dM_dt_4 (nr__6) ! (solMass/yr) ? Error in dM/dt (upper limit)
integer*4 e_dM_dt_5 (nr__6) ! (solMass/yr) ? Error in dM/dt (lower limit)
C=============================================================================
C Loading file 'table1.dat' ! Cluster reference data
C Format for file interpretation
1 format(
+ A12,2X,F5.2,1X,F5.4,2X,A5,1X,I4,A1,2X,I3,1X,I3,1X,A5,1X,F4.1,
+ A1,1X,F3.1,1X,F3.1,1X,A4,2X,F7.3,A1,1X,F5.3,1X,F5.3,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,260
read(1,'(A99)')ar__
read(ar__,1)
+ Name(i__),NH(i__),z(i__),r_z(i__),sigma(i__),u_sigma(i__),
+ E_sigma(i__),e_sigma_1(i__),r_sigma(i__),TX(i__),u_TX(i__),
+ E_TX(i__),e_TX_1(i__),r_TX(i__),LX(i__),u_LX(i__),E_LX(i__),
+ e_LX_1(i__),r_LX(i__)
if(ar__(15:19) .EQ. '') NH(i__) = rNULL__
if(ar__(21:25) .EQ. '') z(i__) = rNULL__
if(ar__(41:43) .EQ. '') E_sigma(i__) = iNULL__
if(ar__(45:47) .EQ. '') e_sigma_1(i__) = iNULL__
if(ar__(61:63) .EQ. '') E_TX(i__) = rNULL__
if(ar__(65:67) .EQ. '') e_TX_1(i__) = rNULL__
if(ar__(84:88) .EQ. '') E_LX(i__) = rNULL__
if(ar__(90:94) .EQ. '') e_LX_1(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),NH(i__),z(i__),r_z(i__),sigma(i__),u_sigma(i__),
+ E_sigma(i__),e_sigma_1(i__),r_sigma(i__),TX(i__),u_TX(i__),
+ E_TX(i__),e_TX_1(i__),r_TX(i__),LX(i__),u_LX(i__),E_LX(i__),
+ e_LX_1(i__),r_LX(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
2 format(A5,2X,A19,1X,A19,1X,A89)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,53
read(1,'(A136)')ar__1
read(ar__1,2)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,2)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Deprojection input data
C Format for file interpretation
3 format(
+ I3,1X,A12,2X,A22,1X,A6,1X,F3.1,1X,F3.1,1X,F4.1,1X,I5,1X,F5.3,
+ 1X,I4,1X,F4.2,1X,F4.1,1X,I4,1X,I4,1X,F4.2,1X,I2,1X,I3,1X,F5.3,
+ 1X,F5.3,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,228
read(1,'(A128)')ar__2
read(ar__2,3)
+ Seq(i__),Name_1(i__),OName(i__),Sequence(i__),IPC1(i__),
+ IPC2(i__),Gain(i__),ExpTime(i__),z_1(i__),DL(i__),NH_1(i__),
+ kT(i__),sigma_1(i__),e_sigma_2(i__),Rcore(i__),Nbins(i__),
+ dR(i__),dRMpc(i__),Rout(i__),Pout(i__)
if(ar__2(49:51) .EQ. '') IPC1(i__) = rNULL__
if(ar__2(53:55) .EQ. '') IPC2(i__) = rNULL__
if(ar__2(57:60) .EQ. '') Gain(i__) = rNULL__
if(ar__2(94:97) .EQ. '') e_sigma_2(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Seq(i__),Name_1(i__),OName(i__),Sequence(i__),IPC1(i__),
+ IPC2(i__),Gain(i__),ExpTime(i__),z_1(i__),DL(i__),NH_1(i__),
+ kT(i__),sigma_1(i__),e_sigma_2(i__),Rcore(i__),Nbins(i__),
+ dR(i__),dRMpc(i__),Rout(i__),Pout(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Deprojection results
C Format for file interpretation
4 format(
+ I3,1X,A12,2X,A6,1X,F4.2,1X,F5.3,1X,F4.1,1X,F3.1,1X,F3.1,1X,
+ F4.1,1X,F5.2,1X,F6.2,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F3.1,
+ 1X,F4.2,1X,F5.2,1X,F4.2,2X,A1,I3,1X,A1,I3,1X,A1,I3,1X,A1,F5.1,
+ 1X,A1,F5.1,1X,A1,F5.1,1X,A4,F5.2,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,228
read(1,'(A162)')ar__3
read(ar__3,4)
+ Seq_1(i__),Name_2(i__),Sequence_1(i__),Rcore_1(i__),
+ Rout_1(i__),kT_1(i__),E_kT(i__),e_kT_1(i__),kTRef(i__),
+ LX_1(i__),LXRef(i__),Mgas(i__),e_Mgas(i__),Mgrav(i__),
+ Mgas_Mgrav(i__),e_Mgas_Mgrav(i__),tcool(i__),E_tcool(i__),
+ e_tcool_1(i__),l_Rcool(i__),Rcool(i__),l_E_Rcool(i__),
+ E_Rcool(i__),l_e_Rcool_1(i__),e_Rcool_1(i__),l_dM_dt(i__),
+ dM_dt(i__),l_E_dM_dt(i__),E_dM_dt(i__),l_e_dM_dt_1(i__),
+ e_dM_dt_1(i__),CF(i__),S_Z_dT(i__),T_Depth(i__)
if(ar__3(50:53) .EQ. '') kTRef(i__) = rNULL__
if(ar__3(61:66) .EQ. '') LXRef(i__) = rNULL__
if(ar__3(86:89) .EQ. '') Mgas_Mgrav(i__) = rNULL__
if(ar__3(91:93) .EQ. '') e_Mgas_Mgrav(i__) = rNULL__
if(ar__3(113:115) .EQ. '') Rcool(i__) = iNULL__
if(ar__3(118:120) .EQ. '') E_Rcool(i__) = iNULL__
if(ar__3(123:125) .EQ. '') e_Rcool_1(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ Seq_1(i__),Name_2(i__),Sequence_1(i__),Rcore_1(i__),
+ Rout_1(i__),kT_1(i__),E_kT(i__),e_kT_1(i__),kTRef(i__),
+ LX_1(i__),LXRef(i__),Mgas(i__),e_Mgas(i__),Mgrav(i__),
+ Mgas_Mgrav(i__),e_Mgas_Mgrav(i__),tcool(i__),E_tcool(i__),
+ e_tcool_1(i__),l_Rcool(i__),Rcool(i__),l_E_Rcool(i__),
+ E_Rcool(i__),l_e_Rcool_1(i__),e_Rcool_1(i__),l_dM_dt(i__),
+ dM_dt(i__),l_E_dM_dt(i__),E_dM_dt(i__),l_e_dM_dt_1(i__),
+ e_dM_dt_1(i__),CF(i__),S_Z_dT(i__),T_Depth(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Deprojection results summed to different radii
C Format for file interpretation
5 format(
+ I3,1X,A12,1X,A6,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.2,1X,F5.2,1X,
+ F4.1,1X,F3.1,1X,F5.1,1X,F4.1,1X,F5.1,1X,F6.1,1X,F4.1,1X,F3.1,
+ 1X,F4.1,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,228
read(1,'(A104)')ar__4
read(ar__4,5)
+ Seq_2(i__),Name_3(i__),Sequence_2(i__),R0_2(i__),R0_5(i__),
+ R1_0(i__),LX0_5(i__),LX1_0(i__),Mgas0_5(i__),e_Mgas0_5(i__),
+ Mgas1_0(i__),e_Mgas1_0(i__),Mgrav0_5(i__),Mgrav1_0(i__),
+ Mgas_Mgrav0_5(i__),e_Mgas_Mgrav0_5(i__),Mgas_Mgrav1_0(i__),
+ e_Mgas_Mgrav1_0(i__)
if(ar__4(1:3) .EQ. '') Seq_2(i__) = iNULL__
if(ar__4(25:29) .EQ. '') R0_2(i__) = rNULL__
if(ar__4(31:35) .EQ. '') R0_5(i__) = rNULL__
if(ar__4(37:41) .EQ. '') R1_0(i__) = rNULL__
if(ar__4(43:47) .EQ. '') LX0_5(i__) = rNULL__
if(ar__4(49:53) .EQ. '') LX1_0(i__) = rNULL__
if(ar__4(55:58) .EQ. '') Mgas0_5(i__) = rNULL__
if(ar__4(60:62) .EQ. '') e_Mgas0_5(i__) = rNULL__
if(ar__4(64:68) .EQ. '') Mgas1_0(i__) = rNULL__
if(ar__4(70:73) .EQ. '') e_Mgas1_0(i__) = rNULL__
if(ar__4(75:79) .EQ. '') Mgrav0_5(i__) = rNULL__
if(ar__4(81:86) .EQ. '') Mgrav1_0(i__) = rNULL__
if(ar__4(88:91) .EQ. '') Mgas_Mgrav0_5(i__) = rNULL__
if(ar__4(93:95) .EQ. '') e_Mgas_Mgrav0_5(i__) = rNULL__
if(ar__4(97:100) .EQ. '') Mgas_Mgrav1_0(i__) = rNULL__
if(ar__4(102:104) .EQ. '') e_Mgas_Mgrav1_0(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Seq_2(i__),Name_3(i__),Sequence_2(i__),R0_2(i__),R0_5(i__),
+ R1_0(i__),LX0_5(i__),LX1_0(i__),Mgas0_5(i__),e_Mgas0_5(i__),
+ Mgas1_0(i__),e_Mgas1_0(i__),Mgrav0_5(i__),Mgrav1_0(i__),
+ Mgas_Mgrav0_5(i__),e_Mgas_Mgrav0_5(i__),Mgas_Mgrav1_0(i__),
+ e_Mgas_Mgrav1_0(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Catalogue of cooling flows detected in the
* Einstein Observatory deprojections
C Format for file interpretation
6 format(
+ I2,A1,1X,I3,1X,A12,2X,A6,1X,A1,F5.1,1X,A1,F5.1,1X,A1,F5.1,1X,
+ A1,F5.1,1X,A1,F5.1,1X,A1,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,97
read(1,'(A70)')ar__5
read(ar__5,6)
+ Rank(i__),n_Rank(i__),Seq_3(i__),Name_4(i__),Sequence_3(i__),
+ l_Rcool_1(i__),Rcool_1(i__),l_E_Rcool_2(i__),E_Rcool_2(i__),
+ l_e_Rcool_3(i__),e_Rcool_3(i__),l_dM_dt_1(i__),dM_dt_1(i__),
+ l_E_dM_dt_2(i__),E_dM_dt_2(i__),l_e_dM_dt_3(i__),
+ e_dM_dt_3(i__)
c ..............Just test output...........
write(6,6)
+ Rank(i__),n_Rank(i__),Seq_3(i__),Name_4(i__),Sequence_3(i__),
+ l_Rcool_1(i__),Rcool_1(i__),l_E_Rcool_2(i__),E_Rcool_2(i__),
+ l_e_Rcool_3(i__),e_Rcool_3(i__),l_dM_dt_1(i__),dM_dt_1(i__),
+ l_E_dM_dt_2(i__),E_dM_dt_2(i__),l_e_dM_dt_3(i__),
+ e_dM_dt_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Supplementary reference data from the literature
C Format for file interpretation
7 format(A12,1X,F7.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,58
read(1,'(A53)')ar__6
read(ar__6,7)
+ Name_5(i__),LHa(i__),tcool_1(i__),E_tcool_2(i__),
+ e_tcool_3(i__),dM_dt_2(i__),E_dM_dt_4(i__),e_dM_dt_5(i__)
if(ar__6(14:20) .EQ. '') LHa(i__) = rNULL__
if(ar__6(22:27) .EQ. '') tcool_1(i__) = rNULL__
if(ar__6(29:34) .EQ. '') E_tcool_2(i__) = rNULL__
if(ar__6(36:41) .EQ. '') e_tcool_3(i__) = rNULL__
if(ar__6(43:45) .EQ. '') dM_dt_2(i__) = iNULL__
if(ar__6(47:49) .EQ. '') E_dM_dt_4(i__) = iNULL__
if(ar__6(51:53) .EQ. '') e_dM_dt_5(i__) = iNULL__
c ..............Just test output...........
write(6,7)
+ Name_5(i__),LHa(i__),tcool_1(i__),E_tcool_2(i__),
+ e_tcool_3(i__),dM_dt_2(i__),E_dM_dt_4(i__),e_dM_dt_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end