FORTRAN Generation
(/./ftp/cats/J/AJ/121/3048)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/121/3048 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-19
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/AJ/121/3048       Multicolor photometry in NGC 1512 and NGC 5248 (Maoz+, 2001)
*================================================================================
*An ultraviolet through infrared look at star formation and super star clusters
*in two circumnuclear starburst rings.
*    Maoz D., Barth A.J., Ho L.C., Sternberg A., Filippenko A.V.
*   <Astron. J. 121, 3048 (2001)>
*   =2001AJ....121.3048M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Emission-Line Complexes

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

      character*8   Cluster    (nr__) ! Cluster name
      integer*4     ELC        (nr__) ! Emission-Line Complexe number in the 
*                                       cluster 
*                                       ([MBH2001 NGC NNNN ELC NN in Simbad)
      integer*4     Xpos       (nr__) ! (pix) X position
      integer*4     Ypos       (nr__) ! (pix) Y position
      real*4        DRA        (nr__) ! (arcsec) Offset in right ascension 
*                                       from the galaxy nucleus (1)
      real*4        DDE        (nr__) ! (arcsec) Offset in declination 
*                                       from the galaxy nucleus (1)
      real*4        Rad        (nr__) ! (arcsec) Angular radius of circular aperture 
*                                       used for flux measurements (2)
      real*4        FHa        (nr__) ! (10-19W/m2) Line flux in H{alpha}+[NII]
*                                       (in F658N filter) (3)
      real*4        FPaa       (nr__) ! (10-19W/m2) Line flux in Pa{alpha} (3)
      real*4        F1_6um     (nr__) ! (10-19W/m2) Line flux in 1.6{mu}m (3)
      real*4        FHa_FPaa   (nr__) ! Ratio of H{alpha}+[NII] to Pa{alpha} lines
      integer*4     FPaa_F1_6um(nr__) ! Ratio of Pa{alpha} to 1.6{mu}m lines
*Note (1): Cluster center coordinates:
*            NGC 1512: 04:03:54.29 (J2000) -43:20:56.7 (419.6,  390.7)
*            NGC 5248: 13:37:32.0  (J2000) +08:53:08   (422.16, 396.02)
*Note (2): At a distance of 10 Mpc, 0.1" corresponds to 5pc
*Note (3): At a distance of 10 Mpc this unit corresponds to a luminosity
*           of 1.2x10^36^erg/s.

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

C  Declarations for 'table3.dat'	! Compact sources brightest in V - NGC 1512

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

      integer*4     CS         (nr__1) ! Compact source number
*                                        ([MBH2001 NGC NNNN CS NN in Simbad)
      real*4        Xpix       (nr__1) ! (pix) X pixel PC coordinate in the F547M frame
      real*4        Ypix       (nr__1) ! (pix) Y pixel PC coordinate in the F547M frame
      real*4        oRA        (nr__1) ! (arcsec) Right Ascension offset from
*                                        galaxy nucleus (1)
      real*4        oDE        (nr__1) ! (arcsec) Declination offset from
*                                        galaxy nucleus (1)
      real*4        Vmag       (nr__1) ! (mag) ? The V band magnitude
      real*4        Del1_3     (nr__1) ! (mag) ? Magnitude difference between F547M
*                                          photometry (about Stroemgren y) in
*                                          3-pixel and 1-pixel radius aperture
      character*1   l_F2300    (nr__1) ! Limit flag on F2300 (2)
      real*8        F2300      (nr__1) ! (10-21W/m2/nm) ? The f_{lambda}_ flux at 2300{AA} (3)
      real*4        e_F2300    (nr__1) ! (10-21W/m2/nm) ? The one {sigma} error in F2300 (4)
      character*1   l_F3350    (nr__1) ! Limit flag on F3350 (2)
      real*4        F3350      (nr__1) ! (10-21W/m2/nm) ? The f_{lambda}_ flux at 3350{AA} (3)
      real*4        e_F3350    (nr__1) ! (10-21W/m2/nm) ? The one {sigma} error in F3350 (U) (4)
      character*1   l_F5490    (nr__1) ! Limit flag on F5490 (2)
      real*4        F5490      (nr__1) ! (10-21W/m2/nm) The f_{lambda}_ flux at 5490{AA} (V) (3)
      real*4        e_F5490    (nr__1) ! (10-21W/m2/nm) ? The one {sigma} error in F5490 (4)
      character*1   l_F8040    (nr__1) ! Limit flag on F8040 (2)
      real*4        F8040      (nr__1) ! (10-21W/m2/nm) The f_{lambda}_ flux at 8040{AA} (I) (3)
      real*4        e_F8040    (nr__1) ! (10-21W/m2/nm) ? The one {sigma} error in F8040 (4)
      character*1   l_F1_6     (nr__1) ! Limit flag on F1.6 (2)
      real*4        F1_6       (nr__1) ! (10-21W/m2/nm) ? The f_{lambda}_ flux at 1.6{mu}m (H)
*                                          (3)(4)
      real*4        e_F1_6     (nr__1) ! (10-21W/m2/nm) ? The one {sigma} error in F1.6 (4)
*Note (1): Galaxy nucleus of NGC 1512 (table3.dat) is at pixel coordinates
*           419.6, 390.7.
*          Galaxy nucleus of NGC 5248 (table4.dat) is at pixel coordinates
*           422.16, 396.02
*Note (2): For null detections, the flux limits actually depend on the local
*          background and crowding, and the limits quoted are just those
*          typical for the band.
*Note (3): Fluxes are in units of 10^-19^ergs/s/cm^2^/{AA}
*Note (4): Missing entries are for sources outside the field of view of the
*          FOC or NICMOS.

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

C  Declarations for 'table4.dat'	! Compact sources brightest in V - NGC 5248

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

      integer*4     CS_1       (nr__2) ! Compact source number
*                                        ([MBH2001 NGC NNNN CS NN in Simbad)
      real*4        Xpix_1     (nr__2) ! (pix) X pixel PC coordinate in the F547M frame
      real*4        Ypix_1     (nr__2) ! (pix) Y pixel PC coordinate in the F547M frame
      real*4        oRA_1      (nr__2) ! (arcsec) Right Ascension offset from
*                                        galaxy nucleus (1)
      real*4        oDE_1      (nr__2) ! (arcsec) Declination offset from
*                                        galaxy nucleus (1)
      real*4        Vmag_1     (nr__2) ! (mag) ? The V band magnitude
      real*4        Del1_3_1   (nr__2) ! (mag) ? Magnitude difference between F547M
*                                          photometry (about Stroemgren y) in
*                                          3-pixel and 1-pixel radius aperture
      character*1   l_F2300_1  (nr__2) ! Limit flag on F2300 (2)
      real*8        F2300_1    (nr__2) ! (10-21W/m2/nm) ? The f_{lambda}_ flux at 2300{AA} (3)
      real*4        e_F2300_1  (nr__2) ! (10-21W/m2/nm) ? The one {sigma} error in F2300 (4)
      character*1   l_F3350_1  (nr__2) ! Limit flag on F3350 (2)
      real*4        F3350_1    (nr__2) ! (10-21W/m2/nm) ? The f_{lambda}_ flux at 3350{AA} (3)
      real*4        e_F3350_1  (nr__2) ! (10-21W/m2/nm) ? The one {sigma} error in F3350 (U) (4)
      character*1   l_F5490_1  (nr__2) ! Limit flag on F5490 (2)
      real*4        F5490_1    (nr__2) ! (10-21W/m2/nm) The f_{lambda}_ flux at 5490{AA} (V) (3)
      real*4        e_F5490_1  (nr__2) ! (10-21W/m2/nm) ? The one {sigma} error in F5490 (4)
      character*1   l_F8040_1  (nr__2) ! Limit flag on F8040 (2)
      real*4        F8040_1    (nr__2) ! (10-21W/m2/nm) The f_{lambda}_ flux at 8040{AA} (I) (3)
      real*4        e_F8040_1  (nr__2) ! (10-21W/m2/nm) ? The one {sigma} error in F8040 (4)
      character*1   l_F1_6_1   (nr__2) ! Limit flag on F1.6 (2)
      real*4        F1_6_1     (nr__2) ! (10-21W/m2/nm) ? The f_{lambda}_ flux at 1.6{mu}m (H)
*                                          (3)(4)
      real*4        e_F1_6_1   (nr__2) ! (10-21W/m2/nm) ? The one {sigma} error in F1.6 (4)
*Note (1): Galaxy nucleus of NGC 1512 (table3.dat) is at pixel coordinates
*           419.6, 390.7.
*          Galaxy nucleus of NGC 5248 (table4.dat) is at pixel coordinates
*           422.16, 396.02
*Note (2): For null detections, the flux limits actually depend on the local
*          background and crowding, and the limits quoted are just those
*          typical for the band.
*Note (3): Fluxes are in units of 10^-19^ergs/s/cm^2^/{AA}
*Note (4): Missing entries are for sources outside the field of view of the
*          FOC or NICMOS.

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

C  Declarations for 'table5.dat'	! Cluster ages, extinctions and masses - NGC 1512

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

      integer*4     CS_2       (nr__3) ! Compact source number
*                                      ([MBH2001] NGC NNNN CS NNN in Simbad)
      integer*4     o_CS       (nr__3) ! Number of bands source observed in (1)
      integer*4     B_Age      (nr__3) ! (Myr) Best fitting model age
      real*4        B_AV       (nr__3) ! (mag) Best fitting model V band extinction
      real*4        B_Mass     (nr__3) ! (10+4solMass) Best fitting model mass
      real*4        B_chi2     (nr__3) ! Reduced {chi}^2^ of best fitting model
      integer*4     Y_Age      (nr__3) ! (Myr) ? Youngest acceptable model age
      real*4        Y_AV       (nr__3) ! (mag) ? Youngest acceptable model V band
*                                        extinction
      real*4        Y_Mass     (nr__3) ! (10+4solMass) ? Youngest acceptable model mass
      real*4        Y_chi2     (nr__3) ! ? Reduced {chi}^2^ of youngest acceptable
*                                        model
      integer*4     O_Age      (nr__3) ! (Myr) ? Oldest acceptable  model age
      real*4        O_AV       (nr__3) ! (mag) ? Oldest acceptable model V band extinction
      real*4        O_Mass     (nr__3) ! (10+4solMass) ? Oldest acceptable model mass
      real*4        O_chi2     (nr__3) ! ? Reduced {chi}^2^ of oldest acceptable
*                                        model
*Note (1): Objects with detections in fewer than 4 bands not fitted.  Formally
*          ``acceptable'' models have {chi}^2^ per degree of freedom <2.

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

C  Declarations for 'table6.dat'	! Cluster ages, extinctions and masses - NGC 5248

      integer*4 nr__4
      parameter (nr__4=308)	! Number of records
      character*64 ar__4  	! Full-size record

      integer*4     CS_3       (nr__4) ! Compact source number
*                                      ([MBH2001] NGC NNNN CS NNN in Simbad)
      integer*4     o_CS_1     (nr__4) ! Number of bands source observed in (1)
      integer*4     B_Age_1    (nr__4) ! (Myr) Best fitting model age
      real*4        B_AV_1     (nr__4) ! (mag) Best fitting model V band extinction
      real*4        B_Mass_1   (nr__4) ! (10+4solMass) Best fitting model mass
      real*4        B_chi2_1   (nr__4) ! Reduced {chi}^2^ of best fitting model
      integer*4     Y_Age_1    (nr__4) ! (Myr) ? Youngest acceptable model age
      real*4        Y_AV_1     (nr__4) ! (mag) ? Youngest acceptable model V band
*                                        extinction
      real*4        Y_Mass_1   (nr__4) ! (10+4solMass) ? Youngest acceptable model mass
      real*4        Y_chi2_1   (nr__4) ! ? Reduced {chi}^2^ of youngest acceptable
*                                        model
      integer*4     O_Age_1    (nr__4) ! (Myr) ? Oldest acceptable  model age
      real*4        O_AV_1     (nr__4) ! (mag) ? Oldest acceptable model V band extinction
      real*4        O_Mass_1   (nr__4) ! (10+4solMass) ? Oldest acceptable model mass
      real*4        O_chi2_1   (nr__4) ! ? Reduced {chi}^2^ of oldest acceptable
*                                        model
*Note (1): Objects with detections in fewer than 4 bands not fitted.  Formally
*          ``acceptable'' models have {chi}^2^ per degree of freedom <2.

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

C  Loading file 'table2.dat'	! Emission-Line Complexes

C  Format for file interpretation

    1 format(
     +  A8,1X,I2,1X,I3,1X,I3,1X,F5.2,2X,F5.2,1X,F4.2,1X,F5.1,1X,F4.1,
     +  1X,F5.1,1X,F3.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,32
        read(1,'(A62)')ar__
        read(ar__,1)
     +  Cluster(i__),ELC(i__),Xpos(i__),Ypos(i__),DRA(i__),DDE(i__),
     +  Rad(i__),FHa(i__),FPaa(i__),F1_6um(i__),FHa_FPaa(i__),
     +  FPaa_F1_6um(i__)
c    ..............Just test output...........
        write(6,1)
     +  Cluster(i__),ELC(i__),Xpos(i__),Ypos(i__),DRA(i__),DDE(i__),
     +  Rad(i__),FHa(i__),FPaa(i__),F1_6um(i__),FHa_FPaa(i__),
     +  FPaa_F1_6um(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Compact sources brightest in V - NGC 1512

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,475
        read(1,'(A110)')ar__1
        read(ar__1,2)
     +  CS(i__),Xpix(i__),Ypix(i__),oRA(i__),oDE(i__),Vmag(i__),
     +  Del1_3(i__),l_F2300(i__),F2300(i__),e_F2300(i__),l_F3350(i__),
     +  F3350(i__),e_F3350(i__),l_F5490(i__),F5490(i__),e_F5490(i__),
     +  l_F8040(i__),F8040(i__),e_F8040(i__),l_F1_6(i__),F1_6(i__),
     +  e_F1_6(i__)
        if(ar__1(33:37) .EQ. '') Vmag(i__) = rNULL__
        if(ar__1(39:42) .EQ. '') Del1_3(i__) = rNULL__
        if(ar__1(45:51) .EQ. '') F2300(i__) = rNULL__
        if(ar__1(53:58) .EQ. '') e_F2300(i__) = rNULL__
        if(ar__1(61:66) .EQ. '') F3350(i__) = rNULL__
        if(ar__1(68:72) .EQ. '') e_F3350(i__) = rNULL__
        if(ar__1(82:85) .EQ. '') e_F5490(i__) = rNULL__
        if(ar__1(94:97) .EQ. '') e_F8040(i__) = rNULL__
        if(ar__1(100:104) .EQ. '') F1_6(i__) = rNULL__
        if(ar__1(106:109) .EQ. '') e_F1_6(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  CS(i__),Xpix(i__),Ypix(i__),oRA(i__),oDE(i__),Vmag(i__),
     +  Del1_3(i__),l_F2300(i__),F2300(i__),e_F2300(i__),l_F3350(i__),
     +  F3350(i__),e_F3350(i__),l_F5490(i__),F5490(i__),e_F5490(i__),
     +  l_F8040(i__),F8040(i__),e_F8040(i__),l_F1_6(i__),F1_6(i__),
     +  e_F1_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Compact sources brightest in V - NGC 5248

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,507
        read(1,'(A110)')ar__2
        read(ar__2,3)
     +  CS_1(i__),Xpix_1(i__),Ypix_1(i__),oRA_1(i__),oDE_1(i__),
     +  Vmag_1(i__),Del1_3_1(i__),l_F2300_1(i__),F2300_1(i__),
     +  e_F2300_1(i__),l_F3350_1(i__),F3350_1(i__),e_F3350_1(i__),
     +  l_F5490_1(i__),F5490_1(i__),e_F5490_1(i__),l_F8040_1(i__),
     +  F8040_1(i__),e_F8040_1(i__),l_F1_6_1(i__),F1_6_1(i__),
     +  e_F1_6_1(i__)
        if(ar__2(33:37) .EQ. '') Vmag_1(i__) = rNULL__
        if(ar__2(39:42) .EQ. '') Del1_3_1(i__) = rNULL__
        if(ar__2(45:51) .EQ. '') F2300_1(i__) = rNULL__
        if(ar__2(53:58) .EQ. '') e_F2300_1(i__) = rNULL__
        if(ar__2(61:66) .EQ. '') F3350_1(i__) = rNULL__
        if(ar__2(68:72) .EQ. '') e_F3350_1(i__) = rNULL__
        if(ar__2(82:85) .EQ. '') e_F5490_1(i__) = rNULL__
        if(ar__2(94:97) .EQ. '') e_F8040_1(i__) = rNULL__
        if(ar__2(100:104) .EQ. '') F1_6_1(i__) = rNULL__
        if(ar__2(106:109) .EQ. '') e_F1_6_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  CS_1(i__),Xpix_1(i__),Ypix_1(i__),oRA_1(i__),oDE_1(i__),
     +  Vmag_1(i__),Del1_3_1(i__),l_F2300_1(i__),F2300_1(i__),
     +  e_F2300_1(i__),l_F3350_1(i__),F3350_1(i__),e_F3350_1(i__),
     +  l_F5490_1(i__),F5490_1(i__),e_F5490_1(i__),l_F8040_1(i__),
     +  F8040_1(i__),e_F8040_1(i__),l_F1_6_1(i__),F1_6_1(i__),
     +  e_F1_6_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Cluster ages, extinctions and masses - NGC 1512

C  Format for file interpretation

    4 format(
     +  I3,1X,I1,1X,I3,1X,F3.1,1X,F5.2,1X,F5.2,1X,I3,1X,F3.1,1X,F5.2,
     +  1X,F4.2,1X,I3,1X,F3.1,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,248
        read(1,'(A64)')ar__3
        read(ar__3,4)
     +  CS_2(i__),o_CS(i__),B_Age(i__),B_AV(i__),B_Mass(i__),
     +  B_chi2(i__),Y_Age(i__),Y_AV(i__),Y_Mass(i__),Y_chi2(i__),
     +  O_Age(i__),O_AV(i__),O_Mass(i__),O_chi2(i__)
        if(ar__3(27:29) .EQ. '') Y_Age(i__) = iNULL__
        if(ar__3(31:33) .EQ. '') Y_AV(i__) = rNULL__
        if(ar__3(35:39) .EQ. '') Y_Mass(i__) = rNULL__
        if(ar__3(41:44) .EQ. '') Y_chi2(i__) = rNULL__
        if(ar__3(46:48) .EQ. '') O_Age(i__) = iNULL__
        if(ar__3(50:52) .EQ. '') O_AV(i__) = rNULL__
        if(ar__3(54:58) .EQ. '') O_Mass(i__) = rNULL__
        if(ar__3(60:63) .EQ. '') O_chi2(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  CS_2(i__),o_CS(i__),B_Age(i__),B_AV(i__),B_Mass(i__),
     +  B_chi2(i__),Y_Age(i__),Y_AV(i__),Y_Mass(i__),Y_chi2(i__),
     +  O_Age(i__),O_AV(i__),O_Mass(i__),O_chi2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Cluster ages, extinctions and masses - NGC 5248

C  Format for file interpretation

    5 format(
     +  I3,1X,I1,1X,I3,1X,F3.1,1X,F5.2,1X,F5.2,1X,I3,1X,F3.1,1X,F5.2,
     +  1X,F4.2,1X,I3,1X,F3.1,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,308
        read(1,'(A64)')ar__4
        read(ar__4,5)
     +  CS_3(i__),o_CS_1(i__),B_Age_1(i__),B_AV_1(i__),B_Mass_1(i__),
     +  B_chi2_1(i__),Y_Age_1(i__),Y_AV_1(i__),Y_Mass_1(i__),
     +  Y_chi2_1(i__),O_Age_1(i__),O_AV_1(i__),O_Mass_1(i__),
     +  O_chi2_1(i__)
        if(ar__4(27:29) .EQ. '') Y_Age_1(i__) = iNULL__
        if(ar__4(31:33) .EQ. '') Y_AV_1(i__) = rNULL__
        if(ar__4(35:39) .EQ. '') Y_Mass_1(i__) = rNULL__
        if(ar__4(41:44) .EQ. '') Y_chi2_1(i__) = rNULL__
        if(ar__4(46:48) .EQ. '') O_Age_1(i__) = iNULL__
        if(ar__4(50:52) .EQ. '') O_AV_1(i__) = rNULL__
        if(ar__4(54:58) .EQ. '') O_Mass_1(i__) = rNULL__
        if(ar__4(60:63) .EQ. '') O_chi2_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  CS_3(i__),o_CS_1(i__),B_Age_1(i__),B_AV_1(i__),B_Mass_1(i__),
     +  B_chi2_1(i__),Y_Age_1(i__),Y_AV_1(i__),Y_Mass_1(i__),
     +  Y_chi2_1(i__),O_Age_1(i__),O_AV_1(i__),O_Mass_1(i__),
     +  O_chi2_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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