Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/831/102 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 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/ApJ/831/102 Swift and NuSTAR obs. of the BL Lac Mrk 421 (Kapanadze+, 2016)
*================================================================================
*X-ray flaring activity of Mrk 421 in the first half of 2013.
* Kapanadze B., Dorner D., Vercellone S., Romano P., Aller H., Aller M.,
* Hughes P., Reynolds M., Kapanadze S., Tabagari L.
* <Astrophys. J., 831, 102-102 (2016)>
* =2016ApJ...831..102K (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Summary of the XRT and UVOT observations
integer*4 nr__
parameter (nr__=68) ! Number of records
character*132 ar__ ! Full-size record
integer*4 obsID ! Observation identifier (1)
character*10 St_date ! ("Y/M/D") Date of the observation start; UTC
* (year 2013 added by CDS)
character*8 St_time ! ("h:m:s") Time of the observation start; UTC
character*10 End_date ! ("Y/M/D") Date of the observation end; UTC
* (year 2013 added by CDS)
character*8 End_time ! ("h:m:s") Time of the observation end; UTC
integer*4 Exp ! (s) [84/10263] Exposure time
real*4 CR ! (ct/s) [7/163] Swift/XRT count rate
real*4 e_CR ! (ct/s) [0.04/0.8] Uncertainty in CR
real*4 uvw1mag ! (mag) [10.5/11.8]? De-reddened Swift/UVOT uvw1 band
* magnitude
real*4 e_uvw1mag ! (mag) [0.03/0.08]? Uncertainty in uvw1mag
real*4 Fuvw1 ! (mJy) [16.9/52]? De-reddened Swift/UVOT uvw1 band
* flux density
real*4 e_Fuvw1 ! (mJy) [0.5/3]? Uncertainty in Fuvw1
real*4 uvm2mag ! (mag) [10.5/12]? De-reddened Swift/UVOT uvm2 band
* magnitude
real*4 e_uvm2mag ! (mag) [0.03/0.08]? Uncertainty in uvm2mag
real*4 Fuvm2 ! (mJy) [15.3/47]? De-reddened Swift/UVOT uvm2 band
* flux density
real*4 e_Fuvm2 ! (mJy) [0.2/2]? Uncertainty in Fuvm2
real*4 uvw2mag ! (mag) [10.5/12]? De-reddened Swift/UVOT uvw2 band
* magnitude
real*4 e_uvw2mag ! (mag) [0.03/0.08]? Uncertainty in uvw2mag
real*4 Fuvw2 ! (mJy) [15/46.5]? De-reddened Swift/UVOT uvw2 band
* flux density
real*4 e_Fuvw2 ! (mJy) [0.4/3]? Uncertainty in Fuvw2
*Note (1): The three leading zeroes are omitted.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Summary of the NuSTAR observations
integer*4 nr__1
parameter (nr__1=23) ! Number of records
character*68 ar__1 ! Full-size record
integer*8 obsID_1 ! Observation identifier
character*10 St_date_1 ! ("Y/M/D") Date of the observation start; UTC
* (year 2013 added by CDS)
character*8 St_time_1 ! ("h:m:s") Time of the observation start; UTC
character*10 End_date_1 ! ("Y/M/D") Date of the observation end; UTC
* (year 2013 added by CDS)
character*8 End_time_1 ! ("h:m:s") Time of the observation end; UTC
integer*4 Exp_1 ! (s) Exposure time
real*4 CR_1 ! (ct/s) NuSTAR count rate
real*4 e_CR_1 ! (ct/s) Uncertainty in CR
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Summary of intra-day X-ray variabilities (IDVs) from
the XRT and NuSTAR observations
integer*4 nr__2
parameter (nr__2=47) ! Number of records
character*131 ar__2 ! Full-size record
character*6 Fac ! Facility identifier (NuSTAR or XRT)
character*22 obsID_2 ! Observation identifier
real*4 Dur ! (h) [0.4/23.5] Observed duration
real*4 chi2 ! [1/460] Reduced {Chi}^2^ of fit
integer*4 DOF ! [4/392] Degrees of Freedom in fit
integer*4 Bin ! (s) [30/120] Bin width
real*4 Fvar ! (%) [1/47.5] Fraction rms variability amplitude
real*4 e_Fvar ! (%) [0.06/5.2] Uncertainty in Fvar
real*4 Dela1 ! [1.6/3.2]? Range of photon index at 1keV;
* first orbit
real*4 e_Dela1 ! [0.01/0.2]? Uncertainty in Dela1
real*4 Dela2 ! [0.2/3.1]? Range of photon index at 1keV;
* second orbit
real*4 e_Dela2 ! [0.01/0.05]? Uncertainty in Dela2
real*4 Delb1 ! [-0.3/0.5]? Curvature parameter; first orbit
real*4 e_Delb1 ! [0.01/0.2]? Uncertainty in Delb1
real*4 Delb2 ! [0.1/0.6]? Curvature parameter; second orbit
real*4 e_Delb2 ! [0.02/0.2]? Uncertainty in Delb2
real*4 DelEp1 ! (keV) [0.001/10]? Position of synchrotron SED;
* first orbit
real*4 e_DelEp1 ! (keV) [0.002/2.2]? Uncertainty in DelEp1
real*4 DelEp2 ! (keV) [0.04/10]? Position of synchrotron SED;
* second orbit
real*4 e_DelEp2 ! (keV) [0.01/2.4]? Uncertainty in DelEp2
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Summary of sub-hour X-ray variability
integer*4 nr__3
parameter (nr__3=93) ! Number of records
character*131 ar__3 ! Full-size record
character*6 Fac_1 ! Facility identifier (NuSTAR or XRT)
character*22 obsID_3 ! Observation identifier
real*4 Dur_1 ! (h) [0.4/23.5] Observed duration
real*4 chi2_1 ! [1/460] Reduced {Chi}^2^ of fit
integer*4 DOF_1 ! [4/392] Degrees of Freedom in fit
integer*4 Bin_1 ! (s) [30/120] Bin width
real*4 Fvar_1 ! (%) [1/47.5] Fraction rms variability amplitude
real*4 e_Fvar_1 ! (%) [0.06/5.2] Uncertainty in Fvar
real*4 Dela1_1 ! [1.6/3.2]? Range of photon index at 1keV;
* first orbit
real*4 e_Dela1_1 ! [0.01/0.2]? Uncertainty in Dela1
real*4 Dela2_1 ! [0.2/3.1]? Range of photon index at 1keV;
* second orbit
real*4 e_Dela2_1 ! [0.01/0.05]? Uncertainty in Dela2
real*4 Delb1_1 ! [-0.3/0.5]? Curvature parameter; first orbit
real*4 e_Delb1_1 ! [0.01/0.2]? Uncertainty in Delb1
real*4 Delb2_1 ! [0.1/0.6]? Curvature parameter; second orbit
real*4 e_Delb2_1 ! [0.02/0.2]? Uncertainty in Delb2
real*4 DelEp1_1 ! (keV) [0.001/10]? Position of synchrotron SED;
* first orbit
real*4 e_DelEp1_1 ! (keV) [0.002/2.2]? Uncertainty in DelEp1
real*4 DelEp2_1 ! (keV) [0.04/10]? Position of synchrotron SED;
* second orbit
real*4 e_DelEp2_1 ! (keV) [0.01/2.4]? Uncertainty in DelEp2
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Summary of the XRT spectral analysis with the
log-parabolic (LP) model
integer*4 nr__4
parameter (nr__4=328) ! Number of records
character*136 ar__4 ! Full-size record
character*24 obsID_4 ! Observation identifier
real*4 a ! [1/2.9] Photon index at 1kev
real*4 e_a ! [0.01/0.05] Uncertainty in a
real*4 b ! [0.09/0.5] Curvature parameter
real*4 e_b ! [0.02/0.2] Uncertainty in b
real*4 Ep ! (keV) [0.003/10] Position of synchrotron SED
real*4 e_Ep ! (keV) [0.001/2.2] Uncertainty in Ep
real*4 K ! [0.01/1.1] Normalization
real*4 e_K ! [0.0009/0.02] Uncertainty in K
real*8 chi2_2 ! [0.8/1.4] Reduced {Chi}^2^ of fit
integer*4 DOF_2 ! [69/475] Degrees of Freedom of fit
character*1 u_DOF ! [?] Uncertainty flag in DOF
real*4 logFS ! ([mW/m2]) [-10/-8.5] Log unabsorbed soft (0.3-2keV)
* flux (in units of erg/s/cm^2^)
real*4 e_logFS ! ([mW/m2]) [0.002/0.05] Uncertainty in logFS
real*8 logFH ! ([mW/m2]) [-10.6/9] Log unabsorbed hard (2-10keV)
* flux (in units of erg/s/cm^2^)
real*4 e_logFH ! ([mW/m2]) [0.005/0.2] Uncertainty in logFH
real*4 logFT ! ([mW/m2]) [-9.9/-8.2] Log unabsorbed total (0.3-10keV)
* flux (in units of erg/s/cm^2^)
real*4 e_logFT ! ([mW/m2]) [0.002/0.09] Uncertainty in logFT
real*4 HR ! [0.1/1.3] Hardness ratio (1)
real*4 e_HR ! [0.005/0.04] Uncertainty in HR
*Note (1): Where HR is the ratio of the unabsorbed 2-10keV to 0.3-2keV fluxes.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Summary of the XRT spectral analysis with a simple
power-law (PL) model
integer*4 nr__5
parameter (nr__5=59) ! Number of records
character*115 ar__5 ! Full-size record
character*28 obsID_5 ! Observation identifier
real*4 Gamma ! Photon index
real*4 e_Gamma ! Uncertainty in Gamma
real*4 K_1 ! Normalization
real*4 e_K_1 ! Uncertainty in K
real*4 chi2_3 ! Reduced {Chi}^2^ of fit
integer*4 DOF_3 ! Degrees of Freedom of fit
real*8 logFS_1 ! ([mW/m2]) Log unabsorbed soft (0.3-2 keV) flux (1)
real*4 e_logFS_1 ! ([mW/m2]) Uncertainty in logFS
real*8 logFH_1 ! ([mW/m2]) Log unabsorbed hard (2-10 keV) flux (1)
real*4 e_logFH_1 ! ([mW/m2]) Uncertainty in logFH
real*8 logFT_1 ! ([mW/m2]) Log unabsorbed total (0.3-10 keV) flux (1)
real*4 e_logFT_1 ! ([mW/m2]) Uncertainty in logFT
real*4 HR_1 ! Hardness ratio (2)
real*4 e_HR_1 ! Uncertainty in HR
*Note (1): In units of erg/s/cm^2^.
*Note (2): Where HR is the ratio of the unabsorbed 2-10keV to 0.3-2keV fluxes.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table10.dat' ! Summary of the NuSTAR spectral analysis with
the LP model
integer*4 nr__6
parameter (nr__6=251) ! Number of records
character*129 ar__6 ! Full-size record
character*19 obsID_6 ! Observation identifier
real*4 a_1 ! [2.2/3.2] Photon index at 1kev
real*4 e_a_1 ! [0.006/0.06] Uncertainty in a
real*4 b_1 ! [0.1/0.6] Curvature parameter
real*4 e_b_1 ! [0/0.2] Uncertainty in b
real*4 Ep_1 ! (keV) [0.001/4.4] Position of synchrotron SED
real*4 e_Ep_1 ! (keV) [0.002/0.9]? Uncertainty in Ep
real*4 K_2 ! [0.03/11] Normalization
real*4 e_K_2 ! [0.01/0.4] Uncertainty in K
real*4 chi2_4 ! [0.8/11] Reduced {Chi}^2^ of fit
integer*4 DOF_4 ! [158/952] Degrees of Freedom of fit
real*8 logFS_2 ! ([mW/m2]) [-10.2/10.1] Log unabsorbed soft (3-10keV)
* flux (in units of erg/s/cm^2^)
real*4 e_logFS_2 ! ([mW/m2]) [0.001/0.04] Uncertainty in logFS
real*8 logFH_2 ! ([mW/m2]) [-10.8/-8.6] Log unabsorbed hard (10-79keV)
* flux (in units of erg/s/cm^2^)
real*4 e_logFH_2 ! ([mW/m2]) [0.004/0.04] Uncertainty in logFH
real*8 logFT_2 ! ([mW/m2]) [-10.1/-8.3] Log unabsorbed total (3-79keV)
* flux (in units of erg/s/cm^2^)
real*4 e_logFT_2 ! ([mW/m2]) [0.002/0.02] Uncertainty in logFT
real*4 HR_2 ! [0.2/1.1] Hardness ratio (1)
real*4 e_HR_2 ! [0.006/0.06] Uncertainty in HR
*Note (1): Where HR is the ratio of the unabsorbed 10-79keV to 3-10keV fluxes.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Summary of the NuSTAR spectral analysis with a simple
PL model
integer*4 nr__7
parameter (nr__7=22) ! Number of records
character*115 ar__7 ! Full-size record
character*28 obsID_7 ! Observation identifier
real*4 Gamma_1 ! Photon index
real*4 e_Gamma_1 ! Uncertainty in Gamma
real*4 K_3 ! Normalization
real*4 e_K_3 ! Uncertainty in K
real*4 chi2_5 ! Reduced {Chi}^2^ of fit
integer*4 DOF_5 ! Degrees of Freedom of fit
real*8 logFS_3 ! ([mW/m2]) Log unabsorbed soft (3-10 keV) flux (1)
real*4 e_logFS_3 ! ([mW/m2]) Uncertainty in logFS
real*8 logFH_3 ! ([mW/m2]) Log unabsorbed hard (10-79) keV flux (1)
real*4 e_logFH_3 ! ([mW/m2]) Uncertainty in logFH
real*8 logFT_3 ! ([mW/m2]) Log unabsorbed total (3-79 keV) flux (1)
real*4 e_logFT_3 ! ([mW/m2]) Uncertainty in logFT
real*4 HR_3 ! Hardness ratio (2)
real*4 e_HR_3 ! Uncertainty in HR
*Note (1): In units of erg/s/cm^2^.
*Note (2): Where HR is the ratio of the unabsorbed 10-79keV to 3-10keV fluxes.
C=============================================================================
C Loading file 'table1.dat' ! Summary of the XRT and UVOT observations
C Format for file interpretation
1 format(
+ I8,1X,A10,1X,A8,1X,A10,1X,A8,1X,I5,1X,F6.2,1X,F4.2,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,F5.2,
+ 1X,F4.2,1X,F5.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,68
read(1,'(A132)')ar__
read(ar__,1)
+ obsID,St_date,St_time,End_date,End_time,Exp,CR,e_CR,uvw1mag,
+ e_uvw1mag,Fuvw1,e_Fuvw1,uvm2mag,e_uvm2mag,Fuvm2,e_Fuvm2,
+ uvw2mag,e_uvw2mag,Fuvw2,e_Fuvw2
if(ar__(68:72) .EQ. '') uvw1mag = rNULL__
if(ar__(74:77) .EQ. '') e_uvw1mag = rNULL__
if(ar__(79:83) .EQ. '') Fuvw1 = rNULL__
if(ar__(85:88) .EQ. '') e_Fuvw1 = rNULL__
if(ar__(90:94) .EQ. '') uvm2mag = rNULL__
if(ar__(96:99) .EQ. '') e_uvm2mag = rNULL__
if(ar__(101:105) .EQ. '') Fuvm2 = rNULL__
if(ar__(107:110) .EQ. '') e_Fuvm2 = rNULL__
if(ar__(112:116) .EQ. '') uvw2mag = rNULL__
if(ar__(118:121) .EQ. '') e_uvw2mag = rNULL__
if(ar__(123:127) .EQ. '') Fuvw2 = rNULL__
if(ar__(129:132) .EQ. '') e_Fuvw2 = rNULL__
c ..............Just test output...........
write(6,1)
+ obsID,St_date,St_time,End_date,End_time,Exp,CR,e_CR,uvw1mag,
+ e_uvw1mag,Fuvw1,e_Fuvw1,uvm2mag,e_uvm2mag,Fuvm2,e_Fuvm2,
+ uvw2mag,e_uvw2mag,Fuvw2,e_Fuvw2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Summary of the NuSTAR observations
C Format for file interpretation
2 format(I11,1X,A10,1X,A8,1X,A10,1X,A8,1X,I5,1X,F5.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,23
read(1,'(A68)')ar__1
read(ar__1,2)
+ obsID_1,St_date_1,St_time_1,End_date_1,End_time_1,Exp_1,CR_1,
+ e_CR_1
c ..............Just test output...........
write(6,2)
+ obsID_1,St_date_1,St_time_1,End_date_1,End_time_1,Exp_1,CR_1,
+ e_CR_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Summary of intra-day X-ray variabilities (IDVs) from
* the XRT and NuSTAR observations
C Format for file interpretation
3 format(
+ A6,1X,A22,1X,F5.2,1X,F6.2,1X,I3,1X,I3,1X,F6.3,1X,F5.3,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.3,1X,F4.2,1X,F4.2,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,47
read(1,'(A131)')ar__2
read(ar__2,3)
+ Fac,obsID_2,Dur,chi2,DOF,Bin,Fvar,e_Fvar,Dela1,e_Dela1,Dela2,
+ e_Dela2,Delb1,e_Delb1,Delb2,e_Delb2,DelEp1,e_DelEp1,DelEp2,
+ e_DelEp2
if(ar__2(65:68) .EQ. '') Dela1 = rNULL__
if(ar__2(70:73) .EQ. '') e_Dela1 = rNULL__
if(ar__2(75:78) .EQ. '') Dela2 = rNULL__
if(ar__2(80:83) .EQ. '') e_Dela2 = rNULL__
if(ar__2(85:89) .EQ. '') Delb1 = rNULL__
if(ar__2(91:95) .EQ. '') e_Delb1 = rNULL__
if(ar__2(97:100) .EQ. '') Delb2 = rNULL__
if(ar__2(102:105) .EQ. '') e_Delb2 = rNULL__
if(ar__2(107:112) .EQ. '') DelEp1 = rNULL__
if(ar__2(114:118) .EQ. '') e_DelEp1 = rNULL__
if(ar__2(120:125) .EQ. '') DelEp2 = rNULL__
if(ar__2(127:131) .EQ. '') e_DelEp2 = rNULL__
c ..............Just test output...........
write(6,3)
+ Fac,obsID_2,Dur,chi2,DOF,Bin,Fvar,e_Fvar,Dela1,e_Dela1,Dela2,
+ e_Dela2,Delb1,e_Delb1,Delb2,e_Delb2,DelEp1,e_DelEp1,DelEp2,
+ e_DelEp2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Summary of sub-hour X-ray variability
C Format for file interpretation
4 format(
+ A6,1X,A22,1X,F5.2,1X,F6.2,1X,I3,1X,I3,1X,F6.3,1X,F5.3,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.3,1X,F4.2,1X,F4.2,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,93
read(1,'(A131)')ar__3
read(ar__3,4)
+ Fac_1,obsID_3,Dur_1,chi2_1,DOF_1,Bin_1,Fvar_1,e_Fvar_1,
+ Dela1_1,e_Dela1_1,Dela2_1,e_Dela2_1,Delb1_1,e_Delb1_1,Delb2_1,
+ e_Delb2_1,DelEp1_1,e_DelEp1_1,DelEp2_1,e_DelEp2_1
if(ar__3(65:68) .EQ. '') Dela1_1 = rNULL__
if(ar__3(70:73) .EQ. '') e_Dela1_1 = rNULL__
if(ar__3(75:78) .EQ. '') Dela2_1 = rNULL__
if(ar__3(80:83) .EQ. '') e_Dela2_1 = rNULL__
if(ar__3(85:89) .EQ. '') Delb1_1 = rNULL__
if(ar__3(91:95) .EQ. '') e_Delb1_1 = rNULL__
if(ar__3(97:100) .EQ. '') Delb2_1 = rNULL__
if(ar__3(102:105) .EQ. '') e_Delb2_1 = rNULL__
if(ar__3(107:112) .EQ. '') DelEp1_1 = rNULL__
if(ar__3(114:118) .EQ. '') e_DelEp1_1 = rNULL__
if(ar__3(120:125) .EQ. '') DelEp2_1 = rNULL__
if(ar__3(127:131) .EQ. '') e_DelEp2_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ Fac_1,obsID_3,Dur_1,chi2_1,DOF_1,Bin_1,Fvar_1,e_Fvar_1,
+ Dela1_1,e_Dela1_1,Dela2_1,e_Dela2_1,Delb1_1,e_Delb1_1,Delb2_1,
+ e_Delb2_1,DelEp1_1,e_DelEp1_1,DelEp2_1,e_DelEp2_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Summary of the XRT spectral analysis with the
* log-parabolic (LP) model
C Format for file interpretation
5 format(
+ A24,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.3,1X,F5.3,1X,F6.4,
+ 1X,F6.4,1X,F7.4,1X,I3,A1,1X,F6.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,328
read(1,'(A136)')ar__4
read(ar__4,5)
+ obsID_4,a,e_a,b,e_b,Ep,e_Ep,K,e_K,chi2_2,DOF_2,u_DOF,logFS,
+ e_logFS,logFH,e_logFH,logFT,e_logFT,HR,e_HR
c ..............Just test output...........
write(6,5)
+ obsID_4,a,e_a,b,e_b,Ep,e_Ep,K,e_K,chi2_2,DOF_2,u_DOF,logFS,
+ e_logFS,logFH,e_logFH,logFT,e_logFT,HR,e_HR
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Summary of the XRT spectral analysis with a simple
* power-law (PL) model
C Format for file interpretation
6 format(
+ A28,1X,F4.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I3,1X,F7.3,1X,
+ F5.3,1X,F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,59
read(1,'(A115)')ar__5
read(ar__5,6)
+ obsID_5,Gamma,e_Gamma,K_1,e_K_1,chi2_3,DOF_3,logFS_1,
+ e_logFS_1,logFH_1,e_logFH_1,logFT_1,e_logFT_1,HR_1,e_HR_1
c ..............Just test output...........
write(6,6)
+ obsID_5,Gamma,e_Gamma,K_1,e_K_1,chi2_3,DOF_3,logFS_1,
+ e_logFS_1,logFH_1,e_logFH_1,logFT_1,e_logFT_1,HR_1,e_HR_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table10.dat' ! Summary of the NuSTAR spectral analysis with
* the LP model
C Format for file interpretation
7 format(
+ A19,1X,F4.2,1X,F5.3,1X,F4.2,1X,F4.2,1X,F5.3,1X,F5.3,1X,F6.3,
+ 1X,F4.2,1X,F6.3,1X,I3,1X,F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F7.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table10.dat')
write(6,*) '....Loading file: table10.dat'
do i__=1,251
read(1,'(A129)')ar__6
read(ar__6,7)
+ obsID_6,a_1,e_a_1,b_1,e_b_1,Ep_1,e_Ep_1,K_2,e_K_2,chi2_4,
+ DOF_4,logFS_2,e_logFS_2,logFH_2,e_logFH_2,logFT_2,e_logFT_2,
+ HR_2,e_HR_2
if(ar__6(48:52) .EQ. '') e_Ep_1 = rNULL__
c ..............Just test output...........
write(6,7)
+ obsID_6,a_1,e_a_1,b_1,e_b_1,Ep_1,e_Ep_1,K_2,e_K_2,chi2_4,
+ DOF_4,logFS_2,e_logFS_2,logFH_2,e_logFH_2,logFT_2,e_logFT_2,
+ HR_2,e_HR_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Summary of the NuSTAR spectral analysis with a simple
* PL model
C Format for file interpretation
8 format(
+ A28,1X,F4.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I3,1X,F7.3,1X,
+ F5.3,1X,F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,22
read(1,'(A115)')ar__7
read(ar__7,8)
+ obsID_7,Gamma_1,e_Gamma_1,K_3,e_K_3,chi2_5,DOF_5,logFS_3,
+ e_logFS_3,logFH_3,e_logFH_3,logFT_3,e_logFT_3,HR_3,e_HR_3
c ..............Just test output...........
write(6,8)
+ obsID_7,Gamma_1,e_Gamma_1,K_3,e_K_3,chi2_5,DOF_5,logFS_3,
+ e_logFS_3,logFH_3,e_logFH_3,logFT_3,e_logFT_3,HR_3,e_HR_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end