Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/808/19 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-Sep-09
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/808/19 Megasecond Chandra XVP obs. of NGC3115. II. (Lin+, 2015)
*================================================================================
*The Megasecond Chandra X-ray Visionary Project observation of NGC 3115.
*II. Properties of point sources.
* Lin D., Irwin J.A., Wong K.-W., Jennings Z.G., Homan J., Romanowsky A.J.,
* Strader J., Sivakoff G.R., Brodie J.P., Remillard R.A.
* <Astrophys. J., 808, 19 (2015)>
* =2015ApJ...808...19L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! The globular cluster LMXBs and candidates
integer*4 nr__
parameter (nr__=48) ! Number of records
character*92 ar__ ! Full-size record
character*3 Type ! Object type (1)
integer*4 MSID ! [4/451] The master source unique index
character*4 GC ! GC in Jennings et al. (2014, J/AJ/148/32;
* <[JSR2014] ANNN> in Simbad)
real*4 dXO ! (arcsec) [0.01/1] The offset between the GC center
* and our X-ray source
real*4 gmag ! (mag) [19.9/24.4]? HST/ACS g-band magnitude
real*4 e_gmag ! (mag) [0.002/0.07]? gmag uncertainty
real*4 zmag ! (mag) [18.7/23.2]? HST/ACS z-band magnitude
real*4 e_zmag ! (mag) [0.002/0.07]? zmag uncertainty
real*4 gSmag ! (mag) [19.5/24.6]? Subaru/Suprime-Cam g-band magnitude
real*4 e_gSmag ! (mag) [0.003/0.3]? gSCmag uncertainty
real*4 rSmag ! (mag) [18.8/23.6]? Subaru/Suprime-Cam r-band magnitude
real*4 e_rSmag ! (mag) [0.003/0.2]? rSCmag uncertainty
real*4 iSmag ! (mag) [18.4/23.6]? Subaru/Suprime-Cam i-band magnitude
real*4 e_iSmag ! (mag) [0.003/0.3]? iSCmag uncertainty
real*4 Rh ! (pc) [0.2/4.8]? Half-light radius
integer*4 Vel ! (km/s) [238/1123]? Heliocentric velocity from
* Pota et al. (2013, J/MNRAS/428/389)
*Note (1): Type as follows:
* GC = LMXBs coincident with HST/ACS GCs (37 sources).
* GCc = LMXBs coincident with GCs detected/covered only in the
* Subaru/Suprime-Cam images, but not in the HST/ACS images (7 sources).
* GCp = GC LMXB candidates whose optical matches were not classified as GCs by
* Jennings et al. (2014, J/AJ/148/32) but were assumed to be GCs by us
* (there is another similar source, S65, which is not included in the
* table because the photometry is not available due to being too close to
* the bright galaxy center; see the text for details).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! The Master Source Catalog
integer*4 nr__1
parameter (nr__1=482) ! Number of records
character*172 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 MSID_1 ! [1/482] The master source unique index
character*16 CXOU ! CXOU Name (JHHMMSS.s+DDMMSS)
real*8 RAdeg ! (deg) Right Ascension in decimal degrees (J2000; ICRS)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000; ICRS)
real*4 PU ! (arcsec) [0.03/8] The 1{sigma} statistical positional
* uncertainty in each coordinate
real*4 Sep ! (arcmin) [0.003/22] Angular separation from galaxy center
real*4 SepR ! Ratio of angular separation to D25 isophote
* elliptical radius
real*4 psf50 ! (arcsec) The 50% PSF radius at the source position in
* merged observation
real*4 psf90 ! (arcsec) The 90% PSF radius at the source position in
* merged observation
integer*4 ObsID ! Observation ID in which the source is
* detected; 0 for the merged observation
real*4 SNR ! [2/95.2] Signal-to-noise ratio (1)
real*4 Vvar ! [1/127] Flux variation factor in 0.5-7keV band
real*4 Svar ! [1/46] Significance of flux variation; equation 1
real*4 Lmax ! (10-7W) Maximum 0.5-7keV luminosity (2)
real*4 b_Lmax ! (10-7W) Lower limit on Lmax at 68% confidence level (2)
real*4 B_Lmax_1 ! (10-7W) Upper limit of Lmax at 68% confidence level (2)
character*5 OLmax ! Corresponding observation ID of Lmax
real*4 Lmin ! (10-7W) Minimum 0.5-7 keV luminosity (2)
real*4 b_Lmin ! (10-7W) Lower limit on Lmin at 68% confidence level (2)
real*4 B_Lmin_1 ! (10-7W) Upper limit of Lmin at 68% confidence level (2)
character*5 OLmin ! Corresponding observation ID of Lmin
integer*4 GLmax ! Maximum Gregory-Loredo short-term
* variability index (6-10 = variable) (3)
character*6 Type_1 ! Source type (4)
*Note (1): The 0.5-7keV net counts divided by the error in the observation the
* source was detected.
*Note (2): Assuming a source distance of 9.7Mpc. In units of erg/s
*Note (3): A higher value indicates variability of a higher degree of
* confidence; sources with indices of 6-10 are definitely variable.
*Note (4): Source type as follows:
* BH = black hole X-ray binary;
* GC = low-mass X-ray binary in a (candidate) globular cluster; those not
* designated as BH are most likely neutron star low-mass
* X-ray binaries;
* F = low-mass X-ray binary in the field; those not designated as BH are
* most likely neutron star low-mass X-ray binaries;
* SSS = supersoft X-ray source;
* AGN = background active galactic nucleus;
* Star = coronally active star in our Galaxy;
* G = hot gas emission in a galaxy
* T = Transient
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! The source counts, flux, and hardness ratio
in the merged observation for the 482 master source
integer*4 nr__2
parameter (nr__2=496) ! Number of records
character*306 ar__2 ! Full-size record
integer*4 MSID_2 ! [1/482] The master source unique index
character*1 Obs ! Observation (G1)
real*4 Exp ! (ks) [35/1131] Exposure time
real*4 Cb ! Counts in the broad band; 0.5-7.0keV
real*4 b_Cb ! Lower limit of Cb at 68% confidence level
real*4 B_Cb_1 ! Upper limit of Cb at 68% confidence level
real*4 Fb ! (mW/m2) Flux in the broad band; 0.5-7.0keV (2)
real*4 b_Fb ! (mW/m2) Lower limit of Fb at 68% confidence level
real*4 B_Fb_1 ! (mW/m2) Upper limit of Fb at 68% confidence level
real*4 Cu ! Counts in the ultra-soft band; 0.2-0.5keV
real*4 b_Cu ! Lower limit of Cu at 68% confidence level
real*4 B_Cu_1 ! Upper limit of Cu at 68% confidence level
real*4 Fu ! (mW/m2) Flux in the ultra-soft band; 0.2-0.5keV (2)
real*4 b_Fu ! (mW/m2) Lower limit of Fu at 68% confidence level
real*4 B_Fu_1 ! (mW/m2) Upper limit of Fu at 68% confidence level
real*4 Cs ! Counts in the soft band; 0.5-1.2keV
real*4 b_Cs ! Lower limit of Cs at 68% confidence level
real*4 B_Cs_1 ! Upper limit of Cs at 68% confidence level
real*4 Fs ! (mW/m2) Flux in the soft band; 0.5-1.2keV (2)
real*4 b_Fs ! (mW/m2) Lower limit of Fs at 68% confidence level
real*4 B_Fs_1 ! (mW/m2) Upper limit of Fs at 68% confidence level
real*4 Cm ! Counts in the medium band; 1.2-2.0keV
real*4 b_Cm ! Lower limit of Cm at 68% confidence level
real*4 B_Cm_1 ! Upper limit of Cm at 68% confidence level
real*4 Fm ! (mW/m2) Flux in the medium band; 1.2-2.0keV (2)
real*4 b_Fm ! (mW/m2) Lower limit of Fm at 68% confidence level
real*4 B_Fm_1 ! (mW/m2) Upper limit of Fm at 68% confidence level
real*4 Ch ! Counts in the hard band; 2.0-7.0keV
real*4 b_Ch ! Lower limit of Ch at 68% confidence level
real*4 B_Ch_1 ! Upper limit of Ch at 68% confidence level
real*4 Fh ! (mW/m2) Flux in the hard band; 2.0-7.0keV (2)
real*4 b_Fh ! (mW/m2) Lower limit of Fh at 68% confidence level
real*4 B_Fh_1 ! (mW/m2) Upper limit of Fh at 68% confidence level
real*4 HR1 ! ? Hardness Ratio using the 0.5-1.2keV and
* the 1.2-2.0 keV fluxes (3)
real*4 bHR1 ! ? Lower limit of HR1 at 68% confidence level
real*4 BHR1_1 ! ? Upper limit of HR1 at 68% confidence level
real*4 HR2 ! Hardness Ratio using the 1.2-2.0keV and
* the 2.0-7.0 keV fluxes (3)
real*4 bHR2 ! Lower limit of HR2 at 68% confidence level
real*4 BHR2_1 ! Upper limit of HR2 at 68% confidence level
real*4 HR3 ! Hardness Ratio using the 0.5-2.0keV and
* the 2.0-7.0 keV fluxes (3)
real*4 bHR3 ! Lower limit of HR3 at 68% confidence level
real*4 BHR3_1 ! Upper limit of HR3 at 68% confidence level
*Note (2): In units of erg/s/cm^2^.
*Note (3): HR is defined as (H-S)/(H+S), where S and H are the energy fluxes
* in the soft and hard energy bands, respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! The source counts and flux in individual observations
for the 482 master sources
integer*4 nr__3
parameter (nr__3=4188) ! Number of records
character*237 ar__3 ! Full-size record
integer*4 MSID_3 ! [1/482] The master source unique index
integer*4 ObsID_1 ! The observation ID
real*4 Cb_1 ! Counts in the broad band; 0.5-7.0keV
real*4 b_Cb_2 ! Lower limit of Cb at 68% confidence level
real*4 B_Cb_3 ! Upper limit of Cb at 68% confidence level
real*4 Fb_1 ! (mW/m2) Flux in the broad band; 0.5-7.0keV (1)
real*4 b_Fb_2 ! (mW/m2) Lower limit of Fb at 68% confidence level
real*4 B_Fb_3 ! (mW/m2) Upper limit of Fb at 68% confidence level
real*4 Cu_1 ! Counts in the ultra-soft band; 0.2-0.5keV
real*4 b_Cu_2 ! Lower limit of Cu at 68% confidence level
real*4 B_Cu_3 ! Upper limit of Cu at 68% confidence level
real*4 Fu_1 ! (mW/m2) Flux in the ultra-soft band; 0.2-0.5keV (1)
real*4 b_Fu_2 ! (mW/m2) Lower limit of Fu at 68% confidence level
real*4 B_Fu_3 ! (mW/m2) Upper limit of Fu at 68% confidence level
real*4 Cs_1 ! Counts in the soft band; 0.5-1.2keV
real*4 b_Cs_2 ! Lower limit of Cs at 68% confidence level
real*4 B_Cs_3 ! Upper limit of Cs at 68% confidence level
real*4 Fs_1 ! (mW/m2) Flux in the soft band; 0.5-1.2keV (1)
real*4 b_Fs_2 ! (mW/m2) Lower limit of Fs at 68% confidence level
real*4 B_Fs_3 ! (mW/m2) Upper limit of Fs at 68% confidence level
real*4 Cm_1 ! Counts in the medium band; 1.2-2.0keV
real*4 b_Cm_2 ! Lower limit of Cm at 68% confidence level
real*4 B_Cm_3 ! Upper limit of Cm at 68% confidence level
real*4 Fm_1 ! (mW/m2) Flux in the medium band; 1.2-2.0keV (1)
real*4 b_Fm_2 ! (mW/m2) Lower limit of Fm at 68% confidence level
real*4 B_Fm_3 ! (mW/m2) Upper limit of Fm at 68% confidence level
real*4 Ch_1 ! Counts in the hard band; 2.0-7.0keV
real*4 b_Ch_2 ! Lower limit of Ch at 68% confidence level
real*4 B_Ch_3 ! Upper limit of Ch at 68% confidence level
real*4 Fh_1 ! (mW/m2) Flux in the hard band; 2.0-7.0keV (1)
real*4 b_Fh_2 ! (mW/m2) Lower limit of Fh at 68% confidence level
real*4 B_Fh_3 ! (mW/m2) Upper limit of Fh at 68% confidence level
*Note (1): In units of erg/s/cm^2^.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! *Spectral fit results
integer*4 nr__4
parameter (nr__4=396) ! Number of records
character*169 ar__4 ! Full-size record
integer*4 MSID_4 ! [1/481] Master source unique index
character*1 Obs_1 ! Observation (G1)
real*4 NH_PL ! (10+20/cm2) Intrinsic column density from PL fit
real*4 e_NH_PL ! (10+20/cm2) ? Lower uncertainty on NH at 90%
* confidence level
real*4 E_NH_PL_1 ! (10+20/cm2) ? Upper uncertainty on NH at 90%
* confidence level
real*4 Gamma_PL ! Photon index
real*4 e_Gamma_PL ! Lower uncertainty on Gamma_PL at 90%
* confidence level
real*4 E_Gamma_PL_1 ! Upper uncertainty on Gamma_PL at 90%
* confidence level
real*4 N_PL ! (10-6) Normalization of the PL fit
real*4 e_N_PL ! (10-6) Lower uncertainty on N-PL at 90%
* confidence level
real*4 E_N_PL_1 ! (10-6) Upper uncertainty on N-PL at 90%
* confidence level
real*4 L_PL ! (10+30W) Luminosity in 0.3-8.0keV from PL fit (2)
real*4 e_L_PL ! (10+30W) Lower uncertainty on L-PL at 90%
* confidence level (2)
real*4 E_L_PL_1 ! (10+30W) Upper uncertainty on L-PL at 90%
* confidence level (2)
real*4 NH_MCD ! (10+20/cm2) Intrinsic column density from MCD fit
real*4 e_NH_MCD ! (10+20/cm2) ? Lower uncertainty on NH-MCD at 90%
* confidence level
real*4 E_NH_MCD_1 ! (10+20/cm2) ? Upper uncertainty on NH-MCD at 90%
* confidence level
real*4 kT_MCD ! (keV) MCD fit apparent inner disk temperature
real*4 e_kT_MCD ! (keV) Lower uncertainty on kT-MCD at 90%
* confidence level
real*4 E_kT_MCD_1 ! (keV) Upper uncertainty on kT-MCD at 90%
* confidence level
real*8 R_MCD ! (km) Apparent inner disk radius from MCD fit (3)
real*8 e_R_MCD ! (km) Lower uncertainty on R-MCD at 90%
* confidence level
real*8 E_R_MCD_1 ! (km) Upper uncertainty on R-MCD at 90%
* confidence level
real*4 L_MCD ! (10+30W) Luminosity in 0.3-8.0keV from MCD fit (2)
real*4 e_L_MCD ! (10+30W) Lower uncertainty on L-MCD at 90%
* confidence level (2)
real*4 E_L_MCD_1 ! (10+30W) Upper uncertainty on L-MCD at 90%
* confidence level (2)
*Note (2): Assume a source distance of 9.7Mpc and isotropic emission for
* the PL fit or an inclination of 60 degree for the MCD fit.
* The units are 1e37erg/s.
*Note (3): From the MCD normalization (R-MCD/km)/(D/10kpc))^2^ cos(theta), where
* D is the source distance assumed to be 9.7Mpc and theta is the disk
* inclination assumed to be 60 degrees.
C=============================================================================
C Loading file 'table2.dat' ! The globular cluster LMXBs and candidates
C Format for file interpretation
1 format(
+ A3,1X,I3,1X,A4,1X,F4.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F4.2,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,48
read(1,'(A92)')ar__
read(ar__,1)
+ Type,MSID,GC,dXO,gmag,e_gmag,zmag,e_zmag,gSmag,e_gSmag,rSmag,
+ e_rSmag,iSmag,e_iSmag,Rh,Vel
if(ar__(19:24) .EQ. '') gmag = rNULL__
if(ar__(26:30) .EQ. '') e_gmag = rNULL__
if(ar__(32:37) .EQ. '') zmag = rNULL__
if(ar__(39:43) .EQ. '') e_zmag = rNULL__
if(ar__(45:50) .EQ. '') gSmag = rNULL__
if(ar__(52:56) .EQ. '') e_gSmag = rNULL__
if(ar__(58:63) .EQ. '') rSmag = rNULL__
if(ar__(65:69) .EQ. '') e_rSmag = rNULL__
if(ar__(71:76) .EQ. '') iSmag = rNULL__
if(ar__(78:82) .EQ. '') e_iSmag = rNULL__
if(ar__(84:87) .EQ. '') Rh = rNULL__
if(ar__(89:92) .EQ. '') Vel = iNULL__
c ..............Just test output...........
write(6,1)
+ Type,MSID,GC,dXO,gmag,e_gmag,zmag,e_zmag,gSmag,e_gSmag,rSmag,
+ e_rSmag,iSmag,e_iSmag,Rh,Vel
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! The Master Source Catalog
C Format for file interpretation
2 format(
+ I3,1X,A16,1X,F10.6,1X,F10.6,1X,F4.2,1X,F6.3,1X,F6.3,1X,F5.2,
+ 1X,F5.2,1X,I5,1X,F6.1,1X,F5.1,1X,F4.1,1X,E8.2,1X,E8.2,1X,E8.2,
+ 1X,A5,1X,E8.2,1X,E8.2,1X,E8.2,1X,A5,1X,I1,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,482
read(1,'(A172)')ar__1
read(ar__1,2)
+ MSID_1,CXOU,RAdeg,DEdeg,PU,Sep,SepR,psf50,psf90,ObsID,SNR,
+ Vvar,Svar,Lmax,b_Lmax,B_Lmax_1,OLmax,Lmin,b_Lmin,B_Lmin_1,
+ OLmin,GLmax,Type_1
c ..............Just test output...........
write(6,2)
+ MSID_1,CXOU,RAdeg,DEdeg,PU,Sep,SepR,psf50,psf90,ObsID,SNR,
+ Vvar,Svar,Lmax,b_Lmax,B_Lmax_1,OLmax,Lmin,b_Lmin,B_Lmin_1,
+ OLmin,GLmax,Type_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! The source counts, flux, and hardness ratio
* in the merged observation for the 482 master source
C Format for file interpretation
3 format(
+ I3,1X,A1,1X,F6.1,1X,F6.1,1X,F6.1,1X,F6.1,1X,E8.2,1X,E8.2,1X,
+ E8.2,1X,F6.1,1X,F6.1,1X,F6.1,1X,E8.2,1X,E8.2,1X,E8.2,1X,F6.1,
+ 1X,F6.1,1X,F6.1,1X,E8.2,1X,E8.2,1X,E8.2,1X,F6.1,1X,F6.1,1X,
+ F6.1,1X,E8.2,1X,E8.2,1X,E8.2,1X,F6.1,1X,F6.1,1X,F6.1,1X,E8.2,
+ 1X,E8.2,1X,E8.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,496
read(1,'(A306)')ar__2
read(ar__2,3)
+ MSID_2,Obs,Exp,Cb,b_Cb,B_Cb_1,Fb,b_Fb,B_Fb_1,Cu,b_Cu,B_Cu_1,
+ Fu,b_Fu,B_Fu_1,Cs,b_Cs,B_Cs_1,Fs,b_Fs,B_Fs_1,Cm,b_Cm,B_Cm_1,
+ Fm,b_Fm,B_Fm_1,Ch,b_Ch,B_Ch_1,Fh,b_Fh,B_Fh_1,HR1,bHR1,BHR1_1,
+ HR2,bHR2,BHR2_1,HR3,bHR3,BHR3_1
if(ar__2(254:258) .EQ. '') HR1 = rNULL__
if(ar__2(260:264) .EQ. '') bHR1 = rNULL__
if(ar__2(266:270) .EQ. '') BHR1_1 = rNULL__
c ..............Just test output...........
write(6,3)
+ MSID_2,Obs,Exp,Cb,b_Cb,B_Cb_1,Fb,b_Fb,B_Fb_1,Cu,b_Cu,B_Cu_1,
+ Fu,b_Fu,B_Fu_1,Cs,b_Cs,B_Cs_1,Fs,b_Fs,B_Fs_1,Cm,b_Cm,B_Cm_1,
+ Fm,b_Fm,B_Fm_1,Ch,b_Ch,B_Ch_1,Fh,b_Fh,B_Fh_1,HR1,bHR1,BHR1_1,
+ HR2,bHR2,BHR2_1,HR3,bHR3,BHR3_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! The source counts and flux in individual observations
* for the 482 master sources
C Format for file interpretation
4 format(
+ I3,1X,I5,1X,F6.1,1X,F6.1,1X,F6.1,1X,E8.2,1X,E8.2,1X,E8.2,1X,
+ F5.1,1X,F5.1,1X,F5.1,1X,E8.2,1X,E8.2,1X,E8.2,1X,F5.1,1X,F5.1,
+ 1X,F5.1,1X,E8.2,1X,E8.2,1X,E8.2,1X,F5.1,1X,F5.1,1X,F5.1,1X,
+ E8.2,1X,E8.2,1X,E8.2,1X,F5.1,1X,F5.1,1X,F5.1,1X,E8.2,1X,E8.2,
+ 1X,E8.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,4188
read(1,'(A237)')ar__3
read(ar__3,4)
+ MSID_3,ObsID_1,Cb_1,b_Cb_2,B_Cb_3,Fb_1,b_Fb_2,B_Fb_3,Cu_1,
+ b_Cu_2,B_Cu_3,Fu_1,b_Fu_2,B_Fu_3,Cs_1,b_Cs_2,B_Cs_3,Fs_1,
+ b_Fs_2,B_Fs_3,Cm_1,b_Cm_2,B_Cm_3,Fm_1,b_Fm_2,B_Fm_3,Ch_1,
+ b_Ch_2,B_Ch_3,Fh_1,b_Fh_2,B_Fh_3
c ..............Just test output...........
write(6,4)
+ MSID_3,ObsID_1,Cb_1,b_Cb_2,B_Cb_3,Fb_1,b_Fb_2,B_Fb_3,Cu_1,
+ b_Cu_2,B_Cu_3,Fu_1,b_Fu_2,B_Fu_3,Cs_1,b_Cs_2,B_Cs_3,Fs_1,
+ b_Fs_2,B_Fs_3,Cm_1,b_Cm_2,B_Cm_3,Fm_1,b_Fm_2,B_Fm_3,Ch_1,
+ b_Ch_2,B_Ch_3,Fh_1,b_Fh_2,B_Fh_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! *Spectral fit results
C Format for file interpretation
5 format(
+ I3,1X,A1,1X,F6.1,1X,F6.1,1X,F6.1,1X,F4.2,1X,F5.2,1X,F4.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F5.1,1X,F5.1,1X,F5.1,1X,F6.1,1X,F6.1,
+ 1X,F6.1,1X,F4.2,1X,F5.2,1X,F4.2,1X,F10.2,1X,F10.2,1X,F10.2,1X,
+ F5.1,1X,F5.1,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,396
read(1,'(A169)')ar__4
read(ar__4,5)
+ MSID_4,Obs_1,NH_PL,e_NH_PL,E_NH_PL_1,Gamma_PL,e_Gamma_PL,
+ E_Gamma_PL_1,N_PL,e_N_PL,E_N_PL_1,L_PL,e_L_PL,E_L_PL_1,NH_MCD,
+ e_NH_MCD,E_NH_MCD_1,kT_MCD,e_kT_MCD,E_kT_MCD_1,R_MCD,e_R_MCD,
+ E_R_MCD_1,L_MCD,e_L_MCD,E_L_MCD_1
if(ar__4(14:19) .EQ. '') e_NH_PL = rNULL__
if(ar__4(21:26) .EQ. '') E_NH_PL_1 = rNULL__
if(ar__4(90:95) .EQ. '') e_NH_MCD = rNULL__
if(ar__4(97:102) .EQ. '') E_NH_MCD_1 = rNULL__
c ..............Just test output...........
write(6,5)
+ MSID_4,Obs_1,NH_PL,e_NH_PL,E_NH_PL_1,Gamma_PL,e_Gamma_PL,
+ E_Gamma_PL_1,N_PL,e_N_PL,E_N_PL_1,L_PL,e_L_PL,E_L_PL_1,NH_MCD,
+ e_NH_MCD,E_NH_MCD_1,kT_MCD,e_kT_MCD,E_kT_MCD_1,R_MCD,e_R_MCD,
+ E_R_MCD_1,L_MCD,e_L_MCD,E_L_MCD_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end