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