FORTRAN Generation
(/./ftp/cats/J/ApJ/659/162)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/659/162 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 2024-Apr-20
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/659/162     CAS parameters of 199 nearby galaxies    (Taylor-Mager+, 2007)
*================================================================================
*Dependence of galaxy structure on rest-frame wavelength and galaxy type.
*    Taylor-Mager V.A., Conselice C.J., Windhorst R.A., Jansen R.A.
*   <Astrophys. J., 659, 162-187 (2007)>
*   =2007ApJ...659..162T
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! Observed galaxy list

      integer*4 nr__
      parameter (nr__=199)	! Number of records
      character*71 ar__   	! Full-size record

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID         (nr__) ! Running identification number
      character*11  Name       (nr__) ! Galaxy name
      integer*4     RAh        (nr__) ! (h) Hour of Right Ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Minute of Right Ascension (J2000)
      real*4        RAs        (nr__) ! (s) Second of Right Ascension (J2000)
      character*1   DE_        (nr__) ! Sign of the Declination (J2000)
      integer*4     DEd        (nr__) ! (deg) Degree of Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Arcminute of Declination (J2000)
      real*4        DEs        (nr__) ! (arcsec) Arcsecond of Declination (J2000)
      real*4        Type       (nr__) ! Classification (1)
      real*4        B_V        (nr__) ! (mag) ? Total (B-V) color index
      real*4        e_B_V      (nr__) ! (mag) ? Uncertainty in B-V
      integer*4     VGSR       (nr__) ! (km/s) ? Galactic standard of rest velocity
      character*7   Com        (nr__) ! Special comments for the galaxy
*Note (1): de Vaucouleurs numerical types were used for normal galaxies,
*     with the following types assigned to peculiar/merging galaxies:
*     13.0=pre-merger, 13.5=minor merger, 14.0=major merger, 
*     14.5=merger remnant, 15.0=peculiar

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! Concentration Index (C)

      integer*4 nr__1
      parameter (nr__1=199)	! Number of records
      character*139 ar__1  	! Full-size record

      integer*4     ID_1       (nr__1) ! Running identification number
      real*4        FUV        (nr__1) ! ? The GALEX FUV (120-177nm) filter index (1)
      character*1   f_FUV      (nr__1) ! [12] S/N or resolution too low in FUV (1)
      real*4        e_FUV      (nr__1) ! ? Uncertainty in FUV
      real*4        NUV        (nr__1) ! ? The GALEX NUV (177-300nm) filter index (1)
      character*1   f_NUV      (nr__1) ! [12] S/N or resolution too low in NUV (1)
      real*4        e_NUV      (nr__1) ! ? Uncertainty in NUV
      real*4        F255W      (nr__1) ! ? The HST F225W filter index (1)
      character*1   f_F255W    (nr__1) ! [12] S/N or resolution too low in F255W (1)
      real*4        e_F255W    (nr__1) ! ? Uncertainty in F255W
      real*4        F300W      (nr__1) ! ? The HST F300W filter index (1)
      character*1   f_F300W    (nr__1) ! [12] S/N or resolution too low in F300W (1)
      real*4        e_F300W    (nr__1) ! ? Uncertainty in F300W
      real*4        U          (nr__1) ! ? The U filter index (1)
      character*1   f_U        (nr__1) ! [12] S/N or resolution too low in U (1)
      real*4        e_U        (nr__1) ! ? Uncertainty in U
      real*4        B          (nr__1) ! ? The B filter index (1)
      character*1   f_B        (nr__1) ! [12] S/N or resolution too low in B (1)
      real*4        e_B        (nr__1) ! ? Uncertainty in B
      real*4        V          (nr__1) ! ? The V filter index (1)
      character*1   f_V        (nr__1) ! [12] S/N or resolution too low in V (1)
      real*4        e_V        (nr__1) ! ? Uncertainty in V
      real*4        R          (nr__1) ! ? The R filter index (1)
      character*1   f_R        (nr__1) ! [12] S/N or resolution too low in R (1)
      real*4        e_R        (nr__1) ! ? Uncertainty in R
      real*4        F814W      (nr__1) ! ? The HST F814W filter index (1)
      character*1   f_F814W    (nr__1) ! [12] S/N or resolution too low in F814W (2)
      real*4        e_F814W    (nr__1) ! ? Uncertainty in F814W
*Note (1): Concentration (C) index in table3, Asymmetry (A) index in table4
*     and Clumpiness (S) index in table5.
*Note (2): Flag as follows:
*      1 = Galaxy does not meet the S/N requirements for accurate CAS
*          parameter measurements in this image (S/N<75)
*      2 = Galaxy does not meet the resolution requirements for accurate CAS
*          parameter measurements in this image (R>1.25kpc)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table4.dat'	! Asymmetry Index (A)

      integer*4 nr__2
      parameter (nr__2=199)	! Number of records
      character*139 ar__2  	! Full-size record

      integer*4     ID_2       (nr__2) ! Running identification number
      real*4        FUV_1      (nr__2) ! ? The GALEX FUV (120-177nm) filter index (1)
      character*1   f_FUV_1    (nr__2) ! [12] S/N or resolution too low in FUV (1)
      real*4        e_FUV_1    (nr__2) ! ? Uncertainty in FUV
      real*4        NUV_1      (nr__2) ! ? The GALEX NUV (177-300nm) filter index (1)
      character*1   f_NUV_1    (nr__2) ! [12] S/N or resolution too low in NUV (1)
      real*4        e_NUV_1    (nr__2) ! ? Uncertainty in NUV
      real*4        F255W_1    (nr__2) ! ? The HST F225W filter index (1)
      character*1   f_F255W_1  (nr__2) ! [12] S/N or resolution too low in F255W (1)
      real*4        e_F255W_1  (nr__2) ! ? Uncertainty in F255W
      real*4        F300W_1    (nr__2) ! ? The HST F300W filter index (1)
      character*1   f_F300W_1  (nr__2) ! [12] S/N or resolution too low in F300W (1)
      real*4        e_F300W_1  (nr__2) ! ? Uncertainty in F300W
      real*4        U_1        (nr__2) ! ? The U filter index (1)
      character*1   f_U_1      (nr__2) ! [12] S/N or resolution too low in U (1)
      real*4        e_U_1      (nr__2) ! ? Uncertainty in U
      real*4        B_1        (nr__2) ! ? The B filter index (1)
      character*1   f_B_1      (nr__2) ! [12] S/N or resolution too low in B (1)
      real*4        e_B_1      (nr__2) ! ? Uncertainty in B
      real*4        V_1        (nr__2) ! ? The V filter index (1)
      character*1   f_V_1      (nr__2) ! [12] S/N or resolution too low in V (1)
      real*4        e_V_1      (nr__2) ! ? Uncertainty in V
      real*4        R_1        (nr__2) ! ? The R filter index (1)
      character*1   f_R_1      (nr__2) ! [12] S/N or resolution too low in R (1)
      real*4        e_R_1      (nr__2) ! ? Uncertainty in R
      real*4        F814W_1    (nr__2) ! ? The HST F814W filter index (1)
      character*1   f_F814W_1  (nr__2) ! [12] S/N or resolution too low in F814W (2)
      real*4        e_F814W_1  (nr__2) ! ? Uncertainty in F814W
*Note (1): Concentration (C) index in table3, Asymmetry (A) index in table4
*     and Clumpiness (S) index in table5.
*Note (2): Flag as follows:
*      1 = Galaxy does not meet the S/N requirements for accurate CAS
*          parameter measurements in this image (S/N<75)
*      2 = Galaxy does not meet the resolution requirements for accurate CAS
*          parameter measurements in this image (R>1.25kpc)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table5.dat'	! Clumpiness Index (S)

      integer*4 nr__3
      parameter (nr__3=199)	! Number of records
      character*139 ar__3  	! Full-size record

      integer*4     ID_3       (nr__3) ! Running identification number
      real*4        FUV_2      (nr__3) ! ? The GALEX FUV (120-177nm) filter index (1)
      character*1   f_FUV_2    (nr__3) ! [12] S/N or resolution too low in FUV (1)
      real*4        e_FUV_2    (nr__3) ! ? Uncertainty in FUV
      real*4        NUV_2      (nr__3) ! ? The GALEX NUV (177-300nm) filter index (1)
      character*1   f_NUV_2    (nr__3) ! [12] S/N or resolution too low in NUV (1)
      real*4        e_NUV_2    (nr__3) ! ? Uncertainty in NUV
      real*4        F255W_2    (nr__3) ! ? The HST F225W filter index (1)
      character*1   f_F255W_2  (nr__3) ! [12] S/N or resolution too low in F255W (1)
      real*4        e_F255W_2  (nr__3) ! ? Uncertainty in F255W
      real*4        F300W_2    (nr__3) ! ? The HST F300W filter index (1)
      character*1   f_F300W_2  (nr__3) ! [12] S/N or resolution too low in F300W (1)
      real*4        e_F300W_2  (nr__3) ! ? Uncertainty in F300W
      real*4        U_2        (nr__3) ! ? The U filter index (1)
      character*1   f_U_2      (nr__3) ! [12] S/N or resolution too low in U (1)
      real*4        e_U_2      (nr__3) ! ? Uncertainty in U
      real*4        B_2        (nr__3) ! ? The B filter index (1)
      character*1   f_B_2      (nr__3) ! [12] S/N or resolution too low in B (1)
      real*4        e_B_2      (nr__3) ! ? Uncertainty in B
      real*4        V_2        (nr__3) ! ? The V filter index (1)
      character*1   f_V_2      (nr__3) ! [12] S/N or resolution too low in V (1)
      real*4        e_V_2      (nr__3) ! ? Uncertainty in V
      real*4        R_2        (nr__3) ! ? The R filter index (1)
      character*1   f_R_2      (nr__3) ! [12] S/N or resolution too low in R (1)
      real*4        e_R_2      (nr__3) ! ? Uncertainty in R
      real*4        F814W_2    (nr__3) ! ? The HST F814W filter index (1)
      character*1   f_F814W_2  (nr__3) ! [12] S/N or resolution too low in F814W (2)
      real*4        e_F814W_2  (nr__3) ! ? Uncertainty in F814W
*Note (1): Concentration (C) index in table3, Asymmetry (A) index in table4
*     and Clumpiness (S) index in table5.
*Note (2): Flag as follows:
*      1 = Galaxy does not meet the S/N requirements for accurate CAS
*          parameter measurements in this image (S/N<75)
*      2 = Galaxy does not meet the resolution requirements for accurate CAS
*          parameter measurements in this image (R>1.25kpc)

C=============================================================================

C  Loading file 'table1.dat'	! Observed galaxy list

C  Format for file interpretation

    1 format(
     +  I3,1X,A11,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.1,
     +  1X,F6.3,1X,F5.3,1X,I5,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,199
        read(1,'(A71)')ar__
        read(ar__,1)
     +  ID(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Type(i__),B_V(i__),e_B_V(i__),
     +  VGSR(i__),Com(i__)
        if(ar__(46:51) .EQ. '') B_V(i__) = rNULL__
        if(ar__(53:57) .EQ. '') e_B_V(i__) = rNULL__
        if(ar__(59:63) .EQ. '') VGSR(i__) = iNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Type(i__),B_V(i__),e_B_V(i__),
     +  VGSR(i__),Com(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3.dat'	! Concentration Index (C)

C  Format for file interpretation

    2 format(
     +  I3,1X,F6.3,1X,A1,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F6.3,1X,A1,
     +  1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,
     +  1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,
     +  F6.3,1X,A1,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,199
        read(1,'(A139)')ar__1
        read(ar__1,2)
     +  ID_1(i__),FUV(i__),f_FUV(i__),e_FUV(i__),NUV(i__),f_NUV(i__),
     +  e_NUV(i__),F255W(i__),f_F255W(i__),e_F255W(i__),F300W(i__),
     +  f_F300W(i__),e_F300W(i__),U(i__),f_U(i__),e_U(i__),B(i__),
     +  f_B(i__),e_B(i__),V(i__),f_V(i__),e_V(i__),R(i__),f_R(i__),
     +  e_R(i__),F814W(i__),f_F814W(i__),e_F814W(i__)
        if(ar__1(5:10) .EQ. '') FUV(i__) = rNULL__
        if(ar__1(14:19) .EQ. '') e_FUV(i__) = rNULL__
        if(ar__1(21:25) .EQ. '') NUV(i__) = rNULL__
        if(ar__1(29:34) .EQ. '') e_NUV(i__) = rNULL__
        if(ar__1(36:41) .EQ. '') F255W(i__) = rNULL__
        if(ar__1(45:49) .EQ. '') e_F255W(i__) = rNULL__
        if(ar__1(51:56) .EQ. '') F300W(i__) = rNULL__
        if(ar__1(60:64) .EQ. '') e_F300W(i__) = rNULL__
        if(ar__1(66:71) .EQ. '') U(i__) = rNULL__
        if(ar__1(75:79) .EQ. '') e_U(i__) = rNULL__
        if(ar__1(81:86) .EQ. '') B(i__) = rNULL__
        if(ar__1(90:94) .EQ. '') e_B(i__) = rNULL__
        if(ar__1(96:101) .EQ. '') V(i__) = rNULL__
        if(ar__1(105:109) .EQ. '') e_V(i__) = rNULL__
        if(ar__1(111:116) .EQ. '') R(i__) = rNULL__
        if(ar__1(120:124) .EQ. '') e_R(i__) = rNULL__
        if(ar__1(126:131) .EQ. '') F814W(i__) = rNULL__
        if(ar__1(135:139) .EQ. '') e_F814W(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  ID_1(i__),FUV(i__),f_FUV(i__),e_FUV(i__),NUV(i__),f_NUV(i__),
     +  e_NUV(i__),F255W(i__),f_F255W(i__),e_F255W(i__),F300W(i__),
     +  f_F300W(i__),e_F300W(i__),U(i__),f_U(i__),e_U(i__),B(i__),
     +  f_B(i__),e_B(i__),V(i__),f_V(i__),e_V(i__),R(i__),f_R(i__),
     +  e_R(i__),F814W(i__),f_F814W(i__),e_F814W(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table4.dat'	! Asymmetry Index (A)

C  Format for file interpretation

    3 format(
     +  I3,1X,F6.3,1X,A1,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F6.3,1X,A1,
     +  1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,
     +  1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,
     +  F6.3,1X,A1,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,199
        read(1,'(A139)')ar__2
        read(ar__2,3)
     +  ID_2(i__),FUV_1(i__),f_FUV_1(i__),e_FUV_1(i__),NUV_1(i__),
     +  f_NUV_1(i__),e_NUV_1(i__),F255W_1(i__),f_F255W_1(i__),
     +  e_F255W_1(i__),F300W_1(i__),f_F300W_1(i__),e_F300W_1(i__),
     +  U_1(i__),f_U_1(i__),e_U_1(i__),B_1(i__),f_B_1(i__),e_B_1(i__),
     +  V_1(i__),f_V_1(i__),e_V_1(i__),R_1(i__),f_R_1(i__),e_R_1(i__),
     +  F814W_1(i__),f_F814W_1(i__),e_F814W_1(i__)
        if(ar__2(5:10) .EQ. '') FUV_1(i__) = rNULL__
        if(ar__2(14:19) .EQ. '') e_FUV_1(i__) = rNULL__
        if(ar__2(21:25) .EQ. '') NUV_1(i__) = rNULL__
        if(ar__2(29:34) .EQ. '') e_NUV_1(i__) = rNULL__
        if(ar__2(36:41) .EQ. '') F255W_1(i__) = rNULL__
        if(ar__2(45:49) .EQ. '') e_F255W_1(i__) = rNULL__
        if(ar__2(51:56) .EQ. '') F300W_1(i__) = rNULL__
        if(ar__2(60:64) .EQ. '') e_F300W_1(i__) = rNULL__
        if(ar__2(66:71) .EQ. '') U_1(i__) = rNULL__
        if(ar__2(75:79) .EQ. '') e_U_1(i__) = rNULL__
        if(ar__2(81:86) .EQ. '') B_1(i__) = rNULL__
        if(ar__2(90:94) .EQ. '') e_B_1(i__) = rNULL__
        if(ar__2(96:101) .EQ. '') V_1(i__) = rNULL__
        if(ar__2(105:109) .EQ. '') e_V_1(i__) = rNULL__
        if(ar__2(111:116) .EQ. '') R_1(i__) = rNULL__
        if(ar__2(120:124) .EQ. '') e_R_1(i__) = rNULL__
        if(ar__2(126:131) .EQ. '') F814W_1(i__) = rNULL__
        if(ar__2(135:139) .EQ. '') e_F814W_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  ID_2(i__),FUV_1(i__),f_FUV_1(i__),e_FUV_1(i__),NUV_1(i__),
     +  f_NUV_1(i__),e_NUV_1(i__),F255W_1(i__),f_F255W_1(i__),
     +  e_F255W_1(i__),F300W_1(i__),f_F300W_1(i__),e_F300W_1(i__),
     +  U_1(i__),f_U_1(i__),e_U_1(i__),B_1(i__),f_B_1(i__),e_B_1(i__),
     +  V_1(i__),f_V_1(i__),e_V_1(i__),R_1(i__),f_R_1(i__),e_R_1(i__),
     +  F814W_1(i__),f_F814W_1(i__),e_F814W_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table5.dat'	! Clumpiness Index (S)

C  Format for file interpretation

    4 format(
     +  I3,1X,F6.3,1X,A1,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F6.3,1X,A1,
     +  1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,
     +  1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,F6.3,1X,A1,1X,F5.3,1X,
     +  F6.3,1X,A1,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,199
        read(1,'(A139)')ar__3
        read(ar__3,4)
     +  ID_3(i__),FUV_2(i__),f_FUV_2(i__),e_FUV_2(i__),NUV_2(i__),
     +  f_NUV_2(i__),e_NUV_2(i__),F255W_2(i__),f_F255W_2(i__),
     +  e_F255W_2(i__),F300W_2(i__),f_F300W_2(i__),e_F300W_2(i__),
     +  U_2(i__),f_U_2(i__),e_U_2(i__),B_2(i__),f_B_2(i__),e_B_2(i__),
     +  V_2(i__),f_V_2(i__),e_V_2(i__),R_2(i__),f_R_2(i__),e_R_2(i__),
     +  F814W_2(i__),f_F814W_2(i__),e_F814W_2(i__)
        if(ar__3(5:10) .EQ. '') FUV_2(i__) = rNULL__
        if(ar__3(14:19) .EQ. '') e_FUV_2(i__) = rNULL__
        if(ar__3(21:25) .EQ. '') NUV_2(i__) = rNULL__
        if(ar__3(29:34) .EQ. '') e_NUV_2(i__) = rNULL__
        if(ar__3(36:41) .EQ. '') F255W_2(i__) = rNULL__
        if(ar__3(45:49) .EQ. '') e_F255W_2(i__) = rNULL__
        if(ar__3(51:56) .EQ. '') F300W_2(i__) = rNULL__
        if(ar__3(60:64) .EQ. '') e_F300W_2(i__) = rNULL__
        if(ar__3(66:71) .EQ. '') U_2(i__) = rNULL__
        if(ar__3(75:79) .EQ. '') e_U_2(i__) = rNULL__
        if(ar__3(81:86) .EQ. '') B_2(i__) = rNULL__
        if(ar__3(90:94) .EQ. '') e_B_2(i__) = rNULL__
        if(ar__3(96:101) .EQ. '') V_2(i__) = rNULL__
        if(ar__3(105:109) .EQ. '') e_V_2(i__) = rNULL__
        if(ar__3(111:116) .EQ. '') R_2(i__) = rNULL__
        if(ar__3(120:124) .EQ. '') e_R_2(i__) = rNULL__
        if(ar__3(126:131) .EQ. '') F814W_2(i__) = rNULL__
        if(ar__3(135:139) .EQ. '') e_F814W_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  ID_3(i__),FUV_2(i__),f_FUV_2(i__),e_FUV_2(i__),NUV_2(i__),
     +  f_NUV_2(i__),e_NUV_2(i__),F255W_2(i__),f_F255W_2(i__),
     +  e_F255W_2(i__),F300W_2(i__),f_F300W_2(i__),e_F300W_2(i__),
     +  U_2(i__),f_U_2(i__),e_U_2(i__),B_2(i__),f_B_2(i__),e_B_2(i__),
     +  V_2(i__),f_V_2(i__),e_V_2(i__),R_2(i__),f_R_2(i__),e_R_2(i__),
     +  F814W_2(i__),f_F814W_2(i__),e_F814W_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end