Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/356/568 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-Aug-11
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/356/568 Deep Chandra survey of the Groth Strip (Nandra+, 2005)
*================================================================================
*A deep Chandra survey of the Groth Strip - I. The X-ray data.
* Nandra K., Laird E.S., Adelberger K., Gardner J.P., Mushotzky R.F.,
* Rhodes J., Steidel C.C., Teplitz H.I., Arnaud K.A.
* <Mon. Not. R. Astron. Soc., 356, 568-586 (2005)>
* =2005MNRAS.356..568N
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Chandra GWS X-ray catalogue
integer*4 nr__
parameter (nr__=158) ! Number of records
character*188 ar__ ! Full-size record
character*4 Seq ! Sequential number (1)
character*16 CXO ! CXO GWS designation (JHHMMSS.s+DDMMSS) (2)
character*1 l_FB ! Limit flag on FB
real*8 FB ! (ct) Full-band (0.5-7keV) counts
real*4 E_FB ! (ct) ? Error on FB, upper limit
real*4 e_FB_1 ! (ct) ? Error on FB, lower limit
character*1 l_SB ! Limit flag on SB
real*8 SB ! (ct) Soft-band, 0.5-2keV) counts
real*4 E_SB ! (ct) ? Error on SB, upper limit
real*4 e_SB_1 ! (ct) ? Error on SB, lower limit
character*1 l_HB ! Limit flag on HB
real*4 HB ! (ct) Hard-band (2-7keV) counts
real*4 E_HB ! (ct) ? Error on HB, upper limit
real*4 e_HB_1 ! (ct) ? Error on HB, lower limit
character*1 l_UB ! Limit flag on UB
real*4 UB ! (ct) Ultra-hard-band (4-7keV) counts
real*4 E_UB ! (ct) ? Error on UB, upper limit
real*4 e_UB_1 ! (ct) ? Error on UB, lower limit
character*1 l_F0_5_10keV ! Limit flag on F0.5-10keV
real*4 F0_5_10keV ! (10-18W/m2) Flux at 0.5-10keV
real*4 E_F0_5_10keV ! (10-18W/m2) ? Error on F0.5-10keV, upper limit
real*4 e_F0_5_10keV_1 ! (10-18W/m2) ? Error on F0.5-10keV, lower limit
character*1 l_F0_5_2keV ! Limit flag on F0.5-2keV
real*4 F0_5_2keV ! (10-18W/m2) Flux at 0.5-2keV
real*4 E_F0_5_2keV ! (10-18W/m2) ? Error on F0.5-2keV, upper limit
real*4 e_F0_5_2keV_1 ! (10-18W/m2) ? Error on F0.5-2keV, lower limit
character*1 l_F2_10keV ! Limit flag on F2-10keV
real*4 F2_10keV ! (10-18W/m2) Flux at 2-10keV
real*4 E_F2_10keV ! (10-18W/m2) ? Error on F2-10keV, upper limit
real*4 e_F2_10keV_1 ! (10-18W/m2) ? Error on F2-10keV, lower limit
character*1 l_F5_10keV ! Limit flag on F5-10keV
real*4 F5_10keV ! (10-18W/m2) Flux at 5-10keV
real*4 E_F5_10keV ! (10-18W/m2) ? rms uncertainty on F5-10keV, upper limit
real*4 e_F5_10keV_1 ! (10-18W/m2) ? rms uncertainty on F5-10keV, lower limit
real*4 log_Pmin ! ([---]) Decimal log of the lowest false detection
* probability found for the four bands (3)
real*4 OAA ! (arcmin) Off-axis angle
real*4 HR ! ? Hardness ratio (4)
character*4 Flags ! [fshu] Flags (5)
*Note (1): Sequential number:
* cNNN in table3 for Chandra GWS X-ray catalogue
* aNNN in tablea1 for additional sources detected in ELAIS-N1
* fNNN in tablea2 for additional sources detected in HDF-N
*Note (2): CXO nomenclature:
* table3: CXO GWS JHHMMSS.s+DDMMSS for Chandra GWS X-ray catalogue
* tablea1: CXO EN1 JHHMMSS.s+DDMMSS for additional sources detected in ELAIS-N1
* tablea2: CXO HDF-N JHHMMSS.s+DDMMSS for additional sources detected in HDF-N
*Note (3): A probability of 1E-8 is assigned if the Poisson probability
* is less than this value
*Note (4): Hardness ratio, HR=(H-S)/(H+S), where H and S are the 2-7keV
* and 0.5-2keV counts, corrected to on-axis values
*Note (5): fshu = source detected at <4x10^-6^ probability in this band,
* where the bands are full (f), soft (s), hard (h) and ultra-hard (u) -
* the first band quoted is the one with the lowest probability
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Additional sources detected in ELAIS-N1
integer*4 nr__1
parameter (nr__1=21) ! Number of records
character*188 ar__1 ! Full-size record
character*4 Seq_1 ! Sequential number (1)
character*16 CXO_1 ! CXO GWS designation (JHHMMSS.s+DDMMSS) (2)
character*1 l_FB_1 ! Limit flag on FB
real*8 FB_1 ! (ct) Full-band (0.5-7keV) counts
real*4 E_FB_2 ! (ct) ? Error on FB, upper limit
real*4 e_FB_3 ! (ct) ? Error on FB, lower limit
character*1 l_SB_1 ! Limit flag on SB
real*8 SB_1 ! (ct) Soft-band, 0.5-2keV) counts
real*4 E_SB_2 ! (ct) ? Error on SB, upper limit
real*4 e_SB_3 ! (ct) ? Error on SB, lower limit
character*1 l_HB_1 ! Limit flag on HB
real*4 HB_1 ! (ct) Hard-band (2-7keV) counts
real*4 E_HB_2 ! (ct) ? Error on HB, upper limit
real*4 e_HB_3 ! (ct) ? Error on HB, lower limit
character*1 l_UB_1 ! Limit flag on UB
real*4 UB_1 ! (ct) Ultra-hard-band (4-7keV) counts
real*4 E_UB_2 ! (ct) ? Error on UB, upper limit
real*4 e_UB_3 ! (ct) ? Error on UB, lower limit
character*1 l_F0_5_10keV_1 ! Limit flag on F0.5-10keV
real*4 F0_5_10keV_1 ! (10-18W/m2) Flux at 0.5-10keV
real*4 E_F0_5_10keV_2 ! (10-18W/m2) ? Error on F0.5-10keV, upper limit
real*4 e_F0_5_10keV_3 ! (10-18W/m2) ? Error on F0.5-10keV, lower limit
character*1 l_F0_5_2keV_1 ! Limit flag on F0.5-2keV
real*4 F0_5_2keV_1 ! (10-18W/m2) Flux at 0.5-2keV
real*4 E_F0_5_2keV_2 ! (10-18W/m2) ? Error on F0.5-2keV, upper limit
real*4 e_F0_5_2keV_3 ! (10-18W/m2) ? Error on F0.5-2keV, lower limit
character*1 l_F2_10keV_1 ! Limit flag on F2-10keV
real*4 F2_10keV_1 ! (10-18W/m2) Flux at 2-10keV
real*4 E_F2_10keV_2 ! (10-18W/m2) ? Error on F2-10keV, upper limit
real*4 e_F2_10keV_3 ! (10-18W/m2) ? Error on F2-10keV, lower limit
character*1 l_F5_10keV_1 ! Limit flag on F5-10keV
real*4 F5_10keV_1 ! (10-18W/m2) Flux at 5-10keV
real*4 E_F5_10keV_2 ! (10-18W/m2) ? rms uncertainty on F5-10keV, upper limit
real*4 e_F5_10keV_3 ! (10-18W/m2) ? rms uncertainty on F5-10keV, lower limit
real*4 log_Pmin_1 ! ([---]) Decimal log of the lowest false detection
* probability found for the four bands (3)
real*4 OAA_1 ! (arcmin) Off-axis angle
real*4 HR_1 ! ? Hardness ratio (4)
character*4 Flags_1 ! [fshu] Flags (5)
*Note (1): Sequential number:
* cNNN in table3 for Chandra GWS X-ray catalogue
* aNNN in tablea1 for additional sources detected in ELAIS-N1
* fNNN in tablea2 for additional sources detected in HDF-N
*Note (2): CXO nomenclature:
* table3: CXO GWS JHHMMSS.s+DDMMSS for Chandra GWS X-ray catalogue
* tablea1: CXO EN1 JHHMMSS.s+DDMMSS for additional sources detected in ELAIS-N1
* tablea2: CXO HDF-N JHHMMSS.s+DDMMSS for additional sources detected in HDF-N
*Note (3): A probability of 1E-8 is assigned if the Poisson probability
* is less than this value
*Note (4): Hardness ratio, HR=(H-S)/(H+S), where H and S are the 2-7keV
* and 0.5-2keV counts, corrected to on-axis values
*Note (5): fshu = source detected at <4x10^-6^ probability in this band,
* where the bands are full (f), soft (s), hard (h) and ultra-hard (u) -
* the first band quoted is the one with the lowest probability
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Additional sources detected in HDF-N
integer*4 nr__2
parameter (nr__2=75) ! Number of records
character*188 ar__2 ! Full-size record
character*4 Seq_2 ! Sequential number (1)
character*16 CXO_2 ! CXO GWS designation (JHHMMSS.s+DDMMSS) (2)
character*1 l_FB_2 ! Limit flag on FB
real*8 FB_2 ! (ct) Full-band (0.5-7keV) counts
real*4 E_FB_4 ! (ct) ? Error on FB, upper limit
real*4 e_FB_5 ! (ct) ? Error on FB, lower limit
character*1 l_SB_2 ! Limit flag on SB
real*8 SB_2 ! (ct) Soft-band, 0.5-2keV) counts
real*4 E_SB_4 ! (ct) ? Error on SB, upper limit
real*4 e_SB_5 ! (ct) ? Error on SB, lower limit
character*1 l_HB_2 ! Limit flag on HB
real*4 HB_2 ! (ct) Hard-band (2-7keV) counts
real*4 E_HB_4 ! (ct) ? Error on HB, upper limit
real*4 e_HB_5 ! (ct) ? Error on HB, lower limit
character*1 l_UB_2 ! Limit flag on UB
real*4 UB_2 ! (ct) Ultra-hard-band (4-7keV) counts
real*4 E_UB_4 ! (ct) ? Error on UB, upper limit
real*4 e_UB_5 ! (ct) ? Error on UB, lower limit
character*1 l_F0_5_10keV_2 ! Limit flag on F0.5-10keV
real*4 F0_5_10keV_2 ! (10-18W/m2) Flux at 0.5-10keV
real*4 E_F0_5_10keV_4 ! (10-18W/m2) ? Error on F0.5-10keV, upper limit
real*4 e_F0_5_10keV_5 ! (10-18W/m2) ? Error on F0.5-10keV, lower limit
character*1 l_F0_5_2keV_2 ! Limit flag on F0.5-2keV
real*4 F0_5_2keV_2 ! (10-18W/m2) Flux at 0.5-2keV
real*4 E_F0_5_2keV_4 ! (10-18W/m2) ? Error on F0.5-2keV, upper limit
real*4 e_F0_5_2keV_5 ! (10-18W/m2) ? Error on F0.5-2keV, lower limit
character*1 l_F2_10keV_2 ! Limit flag on F2-10keV
real*4 F2_10keV_2 ! (10-18W/m2) Flux at 2-10keV
real*4 E_F2_10keV_4 ! (10-18W/m2) ? Error on F2-10keV, upper limit
real*4 e_F2_10keV_5 ! (10-18W/m2) ? Error on F2-10keV, lower limit
character*1 l_F5_10keV_2 ! Limit flag on F5-10keV
real*4 F5_10keV_2 ! (10-18W/m2) Flux at 5-10keV
real*4 E_F5_10keV_4 ! (10-18W/m2) ? rms uncertainty on F5-10keV, upper limit
real*4 e_F5_10keV_5 ! (10-18W/m2) ? rms uncertainty on F5-10keV, lower limit
real*4 log_Pmin_1 ! ([---]) Decimal log of the lowest false detection
* probability found for the four bands (3)
real*4 OAA_2 ! (arcmin) Off-axis angle
real*4 HR_2 ! ? Hardness ratio (4)
character*4 Flags_2 ! [fshu] Flags (5)
*Note (1): Sequential number:
* cNNN in table3 for Chandra GWS X-ray catalogue
* aNNN in tablea1 for additional sources detected in ELAIS-N1
* fNNN in tablea2 for additional sources detected in HDF-N
*Note (2): CXO nomenclature:
* table3: CXO GWS JHHMMSS.s+DDMMSS for Chandra GWS X-ray catalogue
* tablea1: CXO EN1 JHHMMSS.s+DDMMSS for additional sources detected in ELAIS-N1
* tablea2: CXO HDF-N JHHMMSS.s+DDMMSS for additional sources detected in HDF-N
*Note (3): A probability of 1E-8 is assigned if the Poisson probability
* is less than this value
*Note (4): Hardness ratio, HR=(H-S)/(H+S), where H and S are the 2-7keV
* and 0.5-2keV counts, corrected to on-axis values
*Note (5): fshu = source detected at <4x10^-6^ probability in this band,
* where the bands are full (f), soft (s), hard (h) and ultra-hard (u) -
* the first band quoted is the one with the lowest probability
C=============================================================================
C Loading file 'table3.dat' ! Chandra GWS X-ray catalogue
C Format for file interpretation
1 format(
+ A4,1X,A16,1X,A1,F7.2,1X,F4.1,1X,F4.1,1X,A1,F7.2,1X,F4.1,1X,
+ F4.1,1X,A1,F6.2,1X,F4.1,1X,F4.1,1X,A1,F6.2,1X,F4.1,1X,F5.2,A1,
+ F6.2,1X,F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,A1,F6.2,1X,
+ F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,3X,F4.1,1X,F5.2,1X,
+ F4.1,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,158
read(1,'(A188)')ar__
read(ar__,1)
+ Seq,CXO,l_FB,FB,E_FB,e_FB_1,l_SB,SB,E_SB,e_SB_1,l_HB,HB,E_HB,
+ e_HB_1,l_UB,UB,E_UB,e_UB_1,l_F0_5_10keV,F0_5_10keV,
+ E_F0_5_10keV,e_F0_5_10keV_1,l_F0_5_2keV,F0_5_2keV,E_F0_5_2keV,
+ e_F0_5_2keV_1,l_F2_10keV,F2_10keV,E_F2_10keV,e_F2_10keV_1,
+ l_F5_10keV,F5_10keV,E_F5_10keV,e_F5_10keV_1,log_Pmin,OAA,HR,
+ Flags
if(ar__(32:35) .EQ. '') E_FB = rNULL__
if(ar__(37:40) .EQ. '') e_FB_1 = rNULL__
if(ar__(51:54) .EQ. '') E_SB = rNULL__
if(ar__(56:59) .EQ. '') e_SB_1 = rNULL__
if(ar__(69:72) .EQ. '') E_HB = rNULL__
if(ar__(74:77) .EQ. '') e_HB_1 = rNULL__
if(ar__(87:90) .EQ. '') E_UB = rNULL__
if(ar__(92:96) .EQ. '') e_UB_1 = rNULL__
if(ar__(105:108) .EQ. '') E_F0_5_10keV = rNULL__
if(ar__(110:113) .EQ. '') e_F0_5_10keV_1 = rNULL__
if(ar__(122:125) .EQ. '') E_F0_5_2keV = rNULL__
if(ar__(127:130) .EQ. '') e_F0_5_2keV_1 = rNULL__
if(ar__(140:143) .EQ. '') E_F2_10keV = rNULL__
if(ar__(145:148) .EQ. '') e_F2_10keV_1 = rNULL__
if(ar__(157:160) .EQ. '') E_F5_10keV = rNULL__
if(ar__(162:165) .EQ. '') e_F5_10keV_1 = rNULL__
if(ar__(180:183) .EQ. '') HR = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq,CXO,l_FB,FB,E_FB,e_FB_1,l_SB,SB,E_SB,e_SB_1,l_HB,HB,E_HB,
+ e_HB_1,l_UB,UB,E_UB,e_UB_1,l_F0_5_10keV,F0_5_10keV,
+ E_F0_5_10keV,e_F0_5_10keV_1,l_F0_5_2keV,F0_5_2keV,E_F0_5_2keV,
+ e_F0_5_2keV_1,l_F2_10keV,F2_10keV,E_F2_10keV,e_F2_10keV_1,
+ l_F5_10keV,F5_10keV,E_F5_10keV,e_F5_10keV_1,log_Pmin,OAA,HR,
+ Flags
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea1.dat' ! Additional sources detected in ELAIS-N1
C Format for file interpretation
2 format(
+ A4,1X,A16,1X,A1,F7.2,1X,F4.1,1X,F4.1,1X,A1,F7.2,1X,F4.1,1X,
+ F4.1,1X,A1,F6.2,1X,F4.1,1X,F4.1,1X,A1,F6.2,1X,F4.1,1X,F5.2,A1,
+ F6.2,1X,F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,A1,F6.2,1X,
+ F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,3X,F4.1,1X,F5.2,1X,
+ F4.1,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,21
read(1,'(A188)')ar__1
read(ar__1,2)
+ Seq_1,CXO_1,l_FB_1,FB_1,E_FB_2,e_FB_3,l_SB_1,SB_1,E_SB_2,
+ e_SB_3,l_HB_1,HB_1,E_HB_2,e_HB_3,l_UB_1,UB_1,E_UB_2,e_UB_3,
+ l_F0_5_10keV_1,F0_5_10keV_1,E_F0_5_10keV_2,e_F0_5_10keV_3,
+ l_F0_5_2keV_1,F0_5_2keV_1,E_F0_5_2keV_2,e_F0_5_2keV_3,
+ l_F2_10keV_1,F2_10keV_1,E_F2_10keV_2,e_F2_10keV_3,
+ l_F5_10keV_1,F5_10keV_1,E_F5_10keV_2,e_F5_10keV_3,log_Pmin_1,
+ OAA_1,HR_1,Flags_1
if(ar__1(32:35) .EQ. '') E_FB_2 = rNULL__
if(ar__1(37:40) .EQ. '') e_FB_3 = rNULL__
if(ar__1(51:54) .EQ. '') E_SB_2 = rNULL__
if(ar__1(56:59) .EQ. '') e_SB_3 = rNULL__
if(ar__1(69:72) .EQ. '') E_HB_2 = rNULL__
if(ar__1(74:77) .EQ. '') e_HB_3 = rNULL__
if(ar__1(87:90) .EQ. '') E_UB_2 = rNULL__
if(ar__1(92:96) .EQ. '') e_UB_3 = rNULL__
if(ar__1(105:108) .EQ. '') E_F0_5_10keV_2 = rNULL__
if(ar__1(110:113) .EQ. '') e_F0_5_10keV_3 = rNULL__
if(ar__1(122:125) .EQ. '') E_F0_5_2keV_2 = rNULL__
if(ar__1(127:130) .EQ. '') e_F0_5_2keV_3 = rNULL__
if(ar__1(140:143) .EQ. '') E_F2_10keV_2 = rNULL__
if(ar__1(145:148) .EQ. '') e_F2_10keV_3 = rNULL__
if(ar__1(157:160) .EQ. '') E_F5_10keV_2 = rNULL__
if(ar__1(162:165) .EQ. '') e_F5_10keV_3 = rNULL__
if(ar__1(180:183) .EQ. '') HR_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1,CXO_1,l_FB_1,FB_1,E_FB_2,e_FB_3,l_SB_1,SB_1,E_SB_2,
+ e_SB_3,l_HB_1,HB_1,E_HB_2,e_HB_3,l_UB_1,UB_1,E_UB_2,e_UB_3,
+ l_F0_5_10keV_1,F0_5_10keV_1,E_F0_5_10keV_2,e_F0_5_10keV_3,
+ l_F0_5_2keV_1,F0_5_2keV_1,E_F0_5_2keV_2,e_F0_5_2keV_3,
+ l_F2_10keV_1,F2_10keV_1,E_F2_10keV_2,e_F2_10keV_3,
+ l_F5_10keV_1,F5_10keV_1,E_F5_10keV_2,e_F5_10keV_3,log_Pmin_1,
+ OAA_1,HR_1,Flags_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Additional sources detected in HDF-N
C Format for file interpretation
3 format(
+ A4,1X,A16,1X,A1,F7.2,1X,F4.1,1X,F4.1,1X,A1,F7.2,1X,F4.1,1X,
+ F4.1,1X,A1,F6.2,1X,F4.1,1X,F4.1,1X,A1,F6.2,1X,F4.1,1X,F5.2,A1,
+ F6.2,1X,F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,A1,F6.2,1X,
+ F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,3X,F4.1,1X,F5.2,1X,
+ F4.1,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,75
read(1,'(A188)')ar__2
read(ar__2,3)
+ Seq_2,CXO_2,l_FB_2,FB_2,E_FB_4,e_FB_5,l_SB_2,SB_2,E_SB_4,
+ e_SB_5,l_HB_2,HB_2,E_HB_4,e_HB_5,l_UB_2,UB_2,E_UB_4,e_UB_5,
+ l_F0_5_10keV_2,F0_5_10keV_2,E_F0_5_10keV_4,e_F0_5_10keV_5,
+ l_F0_5_2keV_2,F0_5_2keV_2,E_F0_5_2keV_4,e_F0_5_2keV_5,
+ l_F2_10keV_2,F2_10keV_2,E_F2_10keV_4,e_F2_10keV_5,
+ l_F5_10keV_2,F5_10keV_2,E_F5_10keV_4,e_F5_10keV_5,log_Pmin_1,
+ OAA_2,HR_2,Flags_2
if(ar__2(32:35) .EQ. '') E_FB_4 = rNULL__
if(ar__2(37:40) .EQ. '') e_FB_5 = rNULL__
if(ar__2(51:54) .EQ. '') E_SB_4 = rNULL__
if(ar__2(56:59) .EQ. '') e_SB_5 = rNULL__
if(ar__2(69:72) .EQ. '') E_HB_4 = rNULL__
if(ar__2(74:77) .EQ. '') e_HB_5 = rNULL__
if(ar__2(87:90) .EQ. '') E_UB_4 = rNULL__
if(ar__2(92:96) .EQ. '') e_UB_5 = rNULL__
if(ar__2(105:108) .EQ. '') E_F0_5_10keV_4 = rNULL__
if(ar__2(110:113) .EQ. '') e_F0_5_10keV_5 = rNULL__
if(ar__2(122:125) .EQ. '') E_F0_5_2keV_4 = rNULL__
if(ar__2(127:130) .EQ. '') e_F0_5_2keV_5 = rNULL__
if(ar__2(140:143) .EQ. '') E_F2_10keV_4 = rNULL__
if(ar__2(145:148) .EQ. '') e_F2_10keV_5 = rNULL__
if(ar__2(157:160) .EQ. '') E_F5_10keV_4 = rNULL__
if(ar__2(162:165) .EQ. '') e_F5_10keV_5 = rNULL__
if(ar__2(180:183) .EQ. '') HR_2 = rNULL__
c ..............Just test output...........
write(6,3)
+ Seq_2,CXO_2,l_FB_2,FB_2,E_FB_4,e_FB_5,l_SB_2,SB_2,E_SB_4,
+ e_SB_5,l_HB_2,HB_2,E_HB_4,e_HB_5,l_UB_2,UB_2,E_UB_4,e_UB_5,
+ l_F0_5_10keV_2,F0_5_10keV_2,E_F0_5_10keV_4,e_F0_5_10keV_5,
+ l_F0_5_2keV_2,F0_5_2keV_2,E_F0_5_2keV_4,e_F0_5_2keV_5,
+ l_F2_10keV_2,F2_10keV_2,E_F2_10keV_4,e_F2_10keV_5,
+ l_F5_10keV_2,F5_10keV_2,E_F5_10keV_4,e_F5_10keV_5,log_Pmin_1,
+ OAA_2,HR_2,Flags_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end