Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/179/1 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 2026-Oct-12
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/ApJS/179/1 The OPTX project: CLANS, CLASXS and CDF-N (Trouille+, 2008)
*================================================================================
*The OPTX project.
*I. The flux and redshift catalogs for the CLANS, CLASXS, and CDF-N fields.
* Trouille L., Barger A.J., Cowie L.L., Yang Y., Mushotzky R.F.
* <Astrophys. J. Suppl. Ser., 179, 1-18 (2008)>
* =2008ApJS..179....1T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! CLANS X-ray catalog: basic source properties
integer*4 nr__
parameter (nr__=761) ! Number of records
character*146 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 CLANS (nr__) ! Identification number (in RA order)
real*8 RAdeg (nr__) ! (deg) X-ray right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) X-ray declination (J2000)
real*4 e_RAdeg (nr__) ! (arcsec) Wavdetect 95% confidence error in RAdeg
real*4 e_DEdeg (nr__) ! (arcsec) Wavdetect 95% confidence error in DEdeg
character*1 l_F0_5_2 (nr__) ! Limit flag on F0.5-2 (2)
real*4 F0_5_2 (nr__) ! (aW/m2) The 0.5-2keV flux (3)
real*4 e_F0_5_2 (nr__) ! (aW/m2) ?=-99 1{sigma} Poisson lower error (F0.5-2) (4)
real*4 E_F0_5_2_1 (nr__) ! (aW/m2) ?=-99 1{sigma} Poisson upper error (F0.5-2) (4)
character*1 l_F2_8 (nr__) ! Limit flag on F2-8 (2)
real*8 F2_8 (nr__) ! (aW/m2) The 2-8 keV flux (3)
real*4 e_F2_8 (nr__) ! (aW/m2) ?=-99 1{sigma} Poisson lower error (F2-8) (4)
real*4 E_F2_8_1 (nr__) ! (aW/m2) ?=-99 1{sigma} Poisson upper error (F2-8) (4)
real*4 F0_5_8 (nr__) ! (aW/m2) The 0.5-8 keV flux (3)
real*4 e_F0_5_8 (nr__) ! (aW/m2) ?=-99 1{sigma} Poisson lower error (F0.5-8) (4)
real*4 E_F0_5_8_1 (nr__) ! (aW/m2) ?=-99 1{sigma} Poisson upper error (F0.5-8) (4)
character*1 l_HR (nr__) ! Limit Flag on HR (5)
real*8 HR (nr__) ! Hardness ratio (6)
real*8 e_HR (nr__) ! ?=-99 1{sigma} lower error in HR
real*8 E_HR_1 (nr__) ! ?=-99 1{sigma} upper error in HR
real*4 Gamma (nr__) ! ?=-99 Power law photon index
real*4 e_Gamma (nr__) ! ?=-99 1{sigma} lower error in Gamma
real*4 E_Gamma_1 (nr__) ! ?=-99 1{sigma} upper error in Gamma
*Note (2): For sources detected in one band but with a very weak signal in the
* other, we flag the flux value with a < symbol.
*Note (3): For any source detected in one band but with a very weak signal in
* the other, the background-subtracted flux could be negative. In this
* case, only the upper Poisson error is quoted.
*Note (4): Using the approximations from Gehrels (1986). They do not include the
* uncertainty in the flux conversion factor; however, the errors are
* generally dominated by the Poisson errors.
*Note (5): For sources with a very weak signal in the 0.5-2keV (2-8keV) band,
* we flag the HR value with a > (<) symbol.
*Note (6): HR=C(2-8keV)/C(0.5-2keV) where C(2-8keV) and C(0.5-2keV) are
* the count rates in the given bands. For sources with a very weak
* signal in the 0.5-2keV (2-8keV) band, we list the lower
* (upper) limit.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! CLANS X-ray catalog: additional source properties
integer*4 nr__1
parameter (nr__1=761) ! Number of records
character*81 ar__1 ! Full-size record
integer*4 CLANS_1 (nr__1) ! Identification number (in RA order)
real*4 n0_5_2 (nr__1) ! (ct) Net 0.5-2 keV counts (2)
real*4 e_n0_5_2 (nr__1) ! (ct) 1{sigma} Poisson lower error in n0.5-2 (3)
real*4 E_n0_5_2_1 (nr__1) ! (ct) 1{sigma} Poisson upper error in n0.5-2 (3)
real*4 n2_8 (nr__1) ! (ct) Net 2-8 keV counts (2)
real*4 e_n2_8 (nr__1) ! (ct) 1{sigma} Poisson lower error in n2-8 (3)
real*4 E_n2_8_1 (nr__1) ! (ct) 1{sigma} Poisson upper error in n2-8 (3)
real*8 n0_5_8 (nr__1) ! (ct) Net 5-8 counts (2)
real*4 e_n0_5_8 (nr__1) ! (ct) 1{sigma} Poisson lower error in n0.5-8 (3)
real*4 E_n0_5_8_1 (nr__1) ! (ct) 1{sigma} Poisson upper error in n0.5-8 (3)
real*4 t0_5_2 (nr__1) ! (10+4s) Effective 0.5-2 exposure time
real*4 t2_8 (nr__1) ! (10+4s) Effective 2-8 exposure time
real*4 t0_5_8 (nr__1) ! (10+4s) Effective 0.5-8 exposure time
real*4 R (nr__1) ! (arcmin) Source distance from pointing center
*Note (2): The eight-piece segmented annulus region we use to estimate the
* background for each source is four times as large as the source cell
* area. Therefore, the net counts for the majority of the sources are
* multiples of 0.25. For a few sources we exclude background segments
* containing more counts than the 3{sigma} Poisson upper limit in order
* to avoid nearby sources. Therefore the net counts for these sources
* are not multiples of 0.25.
*Note (3): Using the approximations from Gehrels (1986).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! CLANS Optical/IR counterparts catalog
integer*4 nr__2
parameter (nr__2=761) ! Number of records
character*99 ar__2 ! Full-size record
C J2000 position composed of: RAodeg DEodeg
integer*4 CLANS_2 (nr__2) ! Identification number (in RA order)
real*8 RAodeg (nr__2) ! (deg) ?=-99 Opt/NIR Right Ascension in decimal degrees
* (J2000)
real*8 DEodeg (nr__2) ! (deg) ?=-99 Opt/NIR Declination in decimal degrees
* (J2000)
real*4 g_mag (nr__2) ! (mag) ?=-99 g' band AB magnitude (2)
real*4 r_mag (nr__2) ! (mag) ?=-99 r' band AB magnitude (2)
real*4 i_mag (nr__2) ! (mag) ?=-99 i' band AB magnitude (2)
real*4 z_mag (nr__2) ! (mag) ?=-99 z' band AB magnitude (2)
real*4 Jmag (nr__2) ! (mag) ?=-99 J band AB magnitude (2)
real*4 Hmag (nr__2) ! (mag) ?=-99 H band AB magnitude (2)
real*4 Kmag (nr__2) ! (mag) ?=-99 K band AB magnitude (2)
real*4 v3_6mag (nr__2) ! (mag) ?=-99 Spitzer 3.6{mu}m band AB magnitude (2)
real*4 v24mag (nr__2) ! (mag) ?=-99 Spitzer 24{mu}m band AB magnitude (2)
real*4 zsp (nr__2) ! ?=- Spectroscopic redshift (G3)
real*4 zph (nr__2) ! ?=- Photometric redshift (G2)
integer*4 Cl (nr__2) ! ?=- Optical spectral classification (G1)
*Note (2): Typical photometric uncertainties are given in Section 3.6.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! CLASXS Optical/IR counterparts catalog
integer*4 nr__3
parameter (nr__3=525) ! Number of records
character*134 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__3) ! Identification number (1)
real*8 RAdeg_1 (nr__3) ! (deg) ?=-99 Opt/NIR Right Ascension in decimal
* degrees (J2000)
real*8 DEdeg_1 (nr__3) ! (deg) ?=-99 Opt/NIR Declination in decimal degrees
* (J2000)
real*4 Umag (nr__3) ! (mag) ? U band AB magnitude, for table 13 (2)
real*4 umag_1 (nr__3) ! (mag) ? u band AB magnitude, for table 12 (2)
real*4 Bmag (nr__3) ! (mag) ? B band AB magnitude (2)
real*4 g_mag_1 (nr__3) ! (mag) ? g' band AB magnitude, for table 12 (2)
real*4 Vmag (nr__3) ! (mag) ? V band AB magnitude (2)
real*4 Rmag (nr__3) ! (mag) ? R band AB magnitude (2)
real*4 i_mag_1 (nr__3) ! (mag) ? i' band AB magnitude, for table 12 (2)
real*4 Imag (nr__3) ! (mag) ? I band AB magnitude, for table 13 (2)
real*4 z_mag_1 (nr__3) ! (mag) ? z' band AB magnitude (2)
real*4 Jmag_1 (nr__3) ! (mag) ? J band AB magnitude (2)
real*4 Hmag_1 (nr__3) ! (mag) ? H band AB magnitude (2)
real*4 Kmag_1 (nr__3) ! (mag) ? K band AB magnitude (2), (3)
real*4 v3_6mag_1 (nr__3) ! (mag) ? Spitzer 3.6{mu}m band AB magnitude (2)
real*4 v24mag_1 (nr__3) ! (mag) ? Spitzer 24{mu}m band AB magnitude (2)
real*4 zsp_1 (nr__3) ! Spectroscopic redshift (G3)
character*1 r_zsp (nr__3) ! [a-s] Reference for zsp (5)
real*4 zph_1 (nr__3) ! Photometric redshift (G2)
integer*4 Cl_1 (nr__3) ! Optical spectral classification (G1)
*Note (1): Nomenclature note:
* * in Table12: CLASXS NNN, corresponds to objects in
* Tables 2 and 3 of Yang et al. (2004, Cat. J/AJ/128/1501)
* * in Table13: CDFN NNN, corresponds to objects in Table A1 of
* Barger et al. (2003, Cat. J/AJ/126/632), [ABB2003] CDFN NNN in Simbad
*Note (2): Typical photometric uncertainties are given in Section 3.6.
* -99.000 = Not detected (2{sigma} significance).
*Note (3): K band magnitude for table 12 but K_s_ band magnitude for table 13.
*Note (5): References as follows:
* a = Barger et al., 2003AJ....126..632B, Cat. J/AJ/126/632
* b = Swinbank et al., 2004ApJ...617...64S
* c = Chapman et al., 2005ApJ...622..772C, Cat. J/ApJ/622/772
* d = Reddy et al., 2006ApJ...653.1004R, Cat. J/ApJ/653/1004
* s = Steffen et al., 2004AJ....128.1483S, Cat. J/AJ/128/1483
* All other spectroscopic redshifts presented here for the first time.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! CDF-N Optical/IR counterparts catalog
integer*4 nr__4
parameter (nr__4=503) ! Number of records
character*134 ar__4 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq_1 (nr__4) ! Identification number (1)
real*8 RAdeg_2 (nr__4) ! (deg) ?=-99 Opt/NIR Right Ascension in decimal
* degrees (J2000)
real*8 DEdeg_2 (nr__4) ! (deg) ?=-99 Opt/NIR Declination in decimal degrees
* (J2000)
real*4 Umag_2 (nr__4) ! (mag) ? U band AB magnitude, for table 13 (2)
real*4 umag_3 (nr__4) ! (mag) ? u band AB magnitude, for table 12 (2)
real*4 Bmag_1 (nr__4) ! (mag) ? B band AB magnitude (2)
real*4 g_mag_2 (nr__4) ! (mag) ? g' band AB magnitude, for table 12 (2)
real*4 Vmag_1 (nr__4) ! (mag) ? V band AB magnitude (2)
real*4 Rmag_1 (nr__4) ! (mag) ? R band AB magnitude (2)
real*4 i_mag_2 (nr__4) ! (mag) ? i' band AB magnitude, for table 12 (2)
real*4 Imag_1 (nr__4) ! (mag) ? I band AB magnitude, for table 13 (2)
real*4 z_mag_2 (nr__4) ! (mag) ? z' band AB magnitude (2)
real*4 Jmag_2 (nr__4) ! (mag) ? J band AB magnitude (2)
real*4 Hmag_2 (nr__4) ! (mag) ? H band AB magnitude (2)
real*4 Kmag_2 (nr__4) ! (mag) ? K band AB magnitude (2), (3)
real*4 v3_6mag_2 (nr__4) ! (mag) ? Spitzer 3.6{mu}m band AB magnitude (2)
real*4 v24mag_2 (nr__4) ! (mag) ? Spitzer 24{mu}m band AB magnitude (2)
real*4 zsp_2 (nr__4) ! Spectroscopic redshift (G3)
character*1 r_zsp_1 (nr__4) ! [a-s] Reference for zsp (5)
real*4 zph_2 (nr__4) ! Photometric redshift (G2)
integer*4 Cl_2 (nr__4) ! Optical spectral classification (G1)
*Note (1): Nomenclature note:
* * in Table12: CLASXS NNN, corresponds to objects in
* Tables 2 and 3 of Yang et al. (2004, Cat. J/AJ/128/1501)
* * in Table13: CDFN NNN, corresponds to objects in Table A1 of
* Barger et al. (2003, Cat. J/AJ/126/632), [ABB2003] CDFN NNN in Simbad
*Note (2): Typical photometric uncertainties are given in Section 3.6.
* -99.000 = Not detected (2{sigma} significance).
*Note (3): K band magnitude for table 12 but K_s_ band magnitude for table 13.
*Note (5): References as follows:
* a = Barger et al., 2003AJ....126..632B, Cat. J/AJ/126/632
* b = Swinbank et al., 2004ApJ...617...64S
* c = Chapman et al., 2005ApJ...622..772C, Cat. J/ApJ/622/772
* d = Reddy et al., 2006ApJ...653.1004R, Cat. J/ApJ/653/1004
* s = Steffen et al., 2004AJ....128.1483S, Cat. J/AJ/128/1483
* All other spectroscopic redshifts presented here for the first time.
C=============================================================================
C Loading file 'table4.dat' ! CLANS X-ray catalog: basic source properties
C Format for file interpretation
1 format(
+ I3,1X,F9.5,1X,F8.5,1X,F5.3,1X,F5.3,1X,A1,F6.2,1X,F6.2,1X,F6.2,
+ 1X,A1,F7.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,A1,F7.2,
+ 1X,F7.2,1X,F7.2,1X,F6.2,1X,F6.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,761
read(1,'(A146)')ar__
read(ar__,1)
+ CLANS(i__),RAdeg(i__),DEdeg(i__),e_RAdeg(i__),e_DEdeg(i__),
+ l_F0_5_2(i__),F0_5_2(i__),e_F0_5_2(i__),E_F0_5_2_1(i__),
+ l_F2_8(i__),F2_8(i__),e_F2_8(i__),E_F2_8_1(i__),F0_5_8(i__),
+ e_F0_5_8(i__),E_F0_5_8_1(i__),l_HR(i__),HR(i__),e_HR(i__),
+ E_HR_1(i__),Gamma(i__),e_Gamma(i__),E_Gamma_1(i__)
c ..............Just test output...........
write(6,1)
+ CLANS(i__),RAdeg(i__),DEdeg(i__),e_RAdeg(i__),e_DEdeg(i__),
+ l_F0_5_2(i__),F0_5_2(i__),e_F0_5_2(i__),E_F0_5_2_1(i__),
+ l_F2_8(i__),F2_8(i__),e_F2_8(i__),E_F2_8_1(i__),F0_5_8(i__),
+ e_F0_5_8(i__),E_F0_5_8_1(i__),l_HR(i__),HR(i__),e_HR(i__),
+ E_HR_1(i__),Gamma(i__),e_Gamma(i__),E_Gamma_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! CLANS X-ray catalog: additional source properties
C Format for file interpretation
2 format(
+ I3,1X,F6.2,1X,F5.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,F5.2,1X,F7.2,1X,
+ F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,761
read(1,'(A81)')ar__1
read(ar__1,2)
+ CLANS_1(i__),n0_5_2(i__),e_n0_5_2(i__),E_n0_5_2_1(i__),
+ n2_8(i__),e_n2_8(i__),E_n2_8_1(i__),n0_5_8(i__),e_n0_5_8(i__),
+ E_n0_5_8_1(i__),t0_5_2(i__),t2_8(i__),t0_5_8(i__),R(i__)
c ..............Just test output...........
write(6,2)
+ CLANS_1(i__),n0_5_2(i__),e_n0_5_2(i__),E_n0_5_2_1(i__),
+ n2_8(i__),e_n2_8(i__),E_n2_8_1(i__),n0_5_8(i__),e_n0_5_8(i__),
+ E_n0_5_8_1(i__),t0_5_2(i__),t2_8(i__),t0_5_8(i__),R(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! CLANS Optical/IR counterparts catalog
C Format for file interpretation
3 format(
+ I3,1X,F7.3,1X,F7.3,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,761
read(1,'(A99)')ar__2
read(ar__2,3)
+ CLANS_2(i__),RAodeg(i__),DEodeg(i__),g_mag(i__),r_mag(i__),
+ i_mag(i__),z_mag(i__),Jmag(i__),Hmag(i__),Kmag(i__),
+ v3_6mag(i__),v24mag(i__),zsp(i__),zph(i__),Cl(i__)
if (idig(ar__2(84:88)).EQ.0) zsp(i__) = rNULL__
if (idig(ar__2(90:95)).EQ.0) zph(i__) = rNULL__
if (Cl(i__) .EQ. 45) Cl(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ CLANS_2(i__),RAodeg(i__),DEodeg(i__),g_mag(i__),r_mag(i__),
+ i_mag(i__),z_mag(i__),Jmag(i__),Hmag(i__),Kmag(i__),
+ v3_6mag(i__),v24mag(i__),zsp(i__),zph(i__),Cl(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table12.dat' ! CLASXS Optical/IR counterparts catalog
C Format for file interpretation
4 format(
+ I3,1X,F7.3,1X,F7.3,1X,F5.1,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F5.1,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F5.2,1X,A1,1X,F6.2,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,525
read(1,'(A134)')ar__3
read(ar__3,4)
+ Seq(i__),RAdeg_1(i__),DEdeg_1(i__),Umag(i__),umag_1(i__),
+ Bmag(i__),g_mag_1(i__),Vmag(i__),Rmag(i__),i_mag_1(i__),
+ Imag(i__),z_mag_1(i__),Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),
+ v3_6mag_1(i__),v24mag_1(i__),zsp_1(i__),r_zsp(i__),zph_1(i__),
+ Cl_1(i__)
if(ar__3(21:25) .EQ. '') Umag(i__) = rNULL__
if(ar__3(27:32) .EQ. '') umag_1(i__) = rNULL__
if(ar__3(34:39) .EQ. '') Bmag(i__) = rNULL__
if(ar__3(41:46) .EQ. '') g_mag_1(i__) = rNULL__
if(ar__3(48:53) .EQ. '') Vmag(i__) = rNULL__
if(ar__3(55:60) .EQ. '') Rmag(i__) = rNULL__
if(ar__3(62:67) .EQ. '') i_mag_1(i__) = rNULL__
if(ar__3(69:73) .EQ. '') Imag(i__) = rNULL__
if(ar__3(75:80) .EQ. '') z_mag_1(i__) = rNULL__
if(ar__3(82:87) .EQ. '') Jmag_1(i__) = rNULL__
if(ar__3(89:94) .EQ. '') Hmag_1(i__) = rNULL__
if(ar__3(96:101) .EQ. '') Kmag_1(i__) = rNULL__
if(ar__3(103:108) .EQ. '') v3_6mag_1(i__) = rNULL__
if(ar__3(110:115) .EQ. '') v24mag_1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Seq(i__),RAdeg_1(i__),DEdeg_1(i__),Umag(i__),umag_1(i__),
+ Bmag(i__),g_mag_1(i__),Vmag(i__),Rmag(i__),i_mag_1(i__),
+ Imag(i__),z_mag_1(i__),Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),
+ v3_6mag_1(i__),v24mag_1(i__),zsp_1(i__),r_zsp(i__),zph_1(i__),
+ Cl_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! CDF-N Optical/IR counterparts catalog
C Format for file interpretation
5 format(
+ I3,1X,F7.3,1X,F7.3,1X,F5.1,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F5.1,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F5.2,1X,A1,1X,F6.2,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,503
read(1,'(A134)')ar__4
read(ar__4,5)
+ Seq_1(i__),RAdeg_2(i__),DEdeg_2(i__),Umag_2(i__),umag_3(i__),
+ Bmag_1(i__),g_mag_2(i__),Vmag_1(i__),Rmag_1(i__),i_mag_2(i__),
+ Imag_1(i__),z_mag_2(i__),Jmag_2(i__),Hmag_2(i__),Kmag_2(i__),
+ v3_6mag_2(i__),v24mag_2(i__),zsp_2(i__),r_zsp_1(i__),
+ zph_2(i__),Cl_2(i__)
if(ar__4(21:25) .EQ. '') Umag_2(i__) = rNULL__
if(ar__4(27:32) .EQ. '') umag_3(i__) = rNULL__
if(ar__4(34:39) .EQ. '') Bmag_1(i__) = rNULL__
if(ar__4(41:46) .EQ. '') g_mag_2(i__) = rNULL__
if(ar__4(48:53) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__4(55:60) .EQ. '') Rmag_1(i__) = rNULL__
if(ar__4(62:67) .EQ. '') i_mag_2(i__) = rNULL__
if(ar__4(69:73) .EQ. '') Imag_1(i__) = rNULL__
if(ar__4(75:80) .EQ. '') z_mag_2(i__) = rNULL__
if(ar__4(82:87) .EQ. '') Jmag_2(i__) = rNULL__
if(ar__4(89:94) .EQ. '') Hmag_2(i__) = rNULL__
if(ar__4(96:101) .EQ. '') Kmag_2(i__) = rNULL__
if(ar__4(103:108) .EQ. '') v3_6mag_2(i__) = rNULL__
if(ar__4(110:115) .EQ. '') v24mag_2(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Seq_1(i__),RAdeg_2(i__),DEdeg_2(i__),Umag_2(i__),umag_3(i__),
+ Bmag_1(i__),g_mag_2(i__),Vmag_1(i__),Rmag_1(i__),i_mag_2(i__),
+ Imag_1(i__),z_mag_2(i__),Jmag_2(i__),Hmag_2(i__),Kmag_2(i__),
+ v3_6mag_2(i__),v24mag_2(i__),zsp_2(i__),r_zsp_1(i__),
+ zph_2(i__),Cl_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end