Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/116/584 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-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/116/584 Cl J0023+0423 and Cl J1604+4304 morphology (Lubin+ 1998)
*================================================================================
*A study of nine high-redshift clusters of galaxies.
*III. Hubble space telescope morphology of clusters 0023+0423 and 1604+4304.
* Lubin L.M., Postman M., Oke J.B., Ratnatunga K.U., Gunn J.E.,
* Hoessel J.G., Schneider D.P.
* <Astron. J. 116, 584 (1998)>
* =1998AJ....116..584L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1a.dat' ! Cl J0023+0423 field MDS parameters from the
automated galaxy classifications
integer*4 nr__
parameter (nr__=674) ! Number of records
character*150 ar__ ! Full-size record
integer*4 v_OPL98_ ! *MDS identification number
integer*4 Nmod ! *Model number
integer*4 Ntot ! *Total number of pixels in selected region
integer*4 Npix ! *Total number of pixels in selected region
* above 1{sigma}
integer*4 Nfit ! *Number of parameters in MLE model
integer*4 Chip ! *WFPC2 CCD number
real*8 Xpos ! (pix) *X position of centroid of MLE model
real*8 Ypos ! (pix) *Y position of centroid of MLE model
real*4 Sky ! (mag) *MLE model Sky magnitude
real*4 mtot ! (mag) *Total MLE model magnitude
real*4 e_mtot ! (mag) rms uncertainty on mtot
real*4 logR1_2 ! ([arcsec]) *Logarithmic half-light radius of MLE model
real*4 e_logR1_2 ! ([arcsec]) rms uncertainty on logR1/2
real*4 P_A_ ! (rad) *[-1.57/1.57] Orientation of MLE model
real*4 e_P_A_ ! (rad) rms uncertainty on P.A.
real*4 b_aD ! *Disk axis ratio of MLE model
real*4 e_b_aD ! rms uncertainty on e_b/aD
real*4 b_aB ! *Bulge axis ratio of MLE model
real*4 e_b_aB ! rms uncertainty on e_b/aB
real*4 B__B_D ! *Bulge-to-(disk+bulge) luminosity ratio
real*4 e_B__B_D ! rms uncertainty on B/(B+D)
real*4 log_B_D_R1_2 ! log ratio of bulge/disk half-light radius
* of MLE model
real*4 e_log_B_D_R1_2 ! rms uncertainty on (B/D)R1/2
real*4 SNRIL ! log integrated signal-to-noise ratio
character*6 Class ! *Name of MLE classification
real*4 Asymm ! ?=99.99 Asymmetry parameter
*Note on [OPL98]:
* Original identification number from the MDS reduction pipeline. Those numbers
* that are excluded indicate detected objects that were classified as stars or
* which did not reach the required signal-to-noise ratio.
* [OPL98] MDS J0023+0423 NNN in Simbad for table1a.dat
* [OPL98] MDS J1604+4304 NNN in Simbad for table1b.dat
*Note on Nmod:
* MDS model number: 1, disk; 2, bulge; 3, disk+bulge
*Note on Ntot, Npix:
* Ntot is the total number of usable pixels in the selected 64 or 128 pixel
* square region around each object; Npix is the number of Ntot pixels that are
* more than 1{sigma} above the estimated sky background.
*Note on Nfit:
* The number of parameters in the fit that actually vary. A maximum of 12
* parameters are fit, though in most cases one or more parameters are held
* fixed.
*Note on Chip:
* CCD chip number: PC=1; WFC=2,3,4
*Note on Xpos, Ypos:
* The mean error between the model centroid and the actual centroid of the
* object image is 0.2 pixels. Coordinates are relative to each individual chip.
*Note on Sky:
* A maximum likelihood estimate for the local sky background is determined
* simultaneously with the other model parameters. Sufficient pixels are used
* to ensure that the sky level is determined to an accuracy of 0.5%. The sky
* background is assumed to be flat over the small region selected for analysis.
*Note on mtot:
* Small differences between the analytic total magnitude and the true total
* magnitude may arise because the true galaxy is not smooth, and the model
* may not average properly over bright regions of star formation.
*Note on logR1/2:
* The radius within which half of the light of the unconvolved model would
* be contained if it were radially symmetric (an axis ratio of unity). Lower
* and upper limits of 0.1 pixel and one-half of the maximum radius of the
* region selected for analysis have been imposed. The half-light radius of
* the individual components can be derived used the bulge/disk half-radius
* ratio.
*Note on P.A.:
* The adopted position angle (measured east from north) is that of the axis
* of symmetry of the model. It is set equal to zero when the source is
* assumed to be azimuthally symmetric with an axis ratio of unity. For
* disk+bulge models, the orientation is assumed to be the same for both
* components.
*Note on b/aD, b/aB:
* The ratio of the minor axis half-light ratio to that of the major axis.
* For disk+bulge models it is defined independently for each component.
* If the axis is not significantly different from unity, it is fixed at 1.
* A ratio of zero indicates that this component was not fitted to the data.
*Note on B/(B+D):
* The luminosity of the individual component can be derived from this
* parameter and the total magnitude mtot.
*Note on Class:
* Object classification: object, galaxy, disk, bulge or disk+bulge.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1b.dat' ! Cl J1604+4304 field MDS parameters from the
automated galaxy classifications
integer*4 nr__1
parameter (nr__1=559) ! Number of records
character*150 ar__1 ! Full-size record
integer*4 v_OPL98__1 ! *MDS identification number
integer*4 Nmod_1 ! *Model number
integer*4 Ntot_1 ! *Total number of pixels in selected region
integer*4 Npix_1 ! *Total number of pixels in selected region
* above 1{sigma}
integer*4 Nfit_1 ! *Number of parameters in MLE model
integer*4 Chip_1 ! *WFPC2 CCD number
real*8 Xpos_1 ! (pix) *X position of centroid of MLE model
real*8 Ypos_1 ! (pix) *Y position of centroid of MLE model
real*4 Sky_1 ! (mag) *MLE model Sky magnitude
real*4 mtot_1 ! (mag) *Total MLE model magnitude
real*4 e_mtot_1 ! (mag) rms uncertainty on mtot
real*4 logR1_2_1 ! ([arcsec]) *Logarithmic half-light radius of MLE model
real*4 e_logR1_2_1 ! ([arcsec]) rms uncertainty on logR1/2
real*4 P_A__1 ! (rad) *[-1.57/1.57] Orientation of MLE model
real*4 e_P_A__1 ! (rad) rms uncertainty on P.A.
real*4 b_aD_1 ! *Disk axis ratio of MLE model
real*4 e_b_aD_1 ! rms uncertainty on e_b/aD
real*4 b_aB_1 ! *Bulge axis ratio of MLE model
real*4 e_b_aB_1 ! rms uncertainty on e_b/aB
real*4 B__B_D_1 ! *Bulge-to-(disk+bulge) luminosity ratio
real*4 e_B__B_D_1 ! rms uncertainty on B/(B+D)
real*4 log_B_D_R1_2_1 ! log ratio of bulge/disk half-light radius
* of MLE model
real*4 e_log_B_D_R1_2_1 ! rms uncertainty on (B/D)R1/2
real*4 SNRIL_1 ! log integrated signal-to-noise ratio
character*6 Class_1 ! *Name of MLE classification
real*4 Asymm_1 ! ?=99.99 Asymmetry parameter
*Note on [OPL98]:
* Original identification number from the MDS reduction pipeline. Those numbers
* that are excluded indicate detected objects that were classified as stars or
* which did not reach the required signal-to-noise ratio.
* [OPL98] MDS J0023+0423 NNN in Simbad for table1a.dat
* [OPL98] MDS J1604+4304 NNN in Simbad for table1b.dat
*Note on Nmod:
* MDS model number: 1, disk; 2, bulge; 3, disk+bulge
*Note on Ntot, Npix:
* Ntot is the total number of usable pixels in the selected 64 or 128 pixel
* square region around each object; Npix is the number of Ntot pixels that are
* more than 1{sigma} above the estimated sky background.
*Note on Nfit:
* The number of parameters in the fit that actually vary. A maximum of 12
* parameters are fit, though in most cases one or more parameters are held
* fixed.
*Note on Chip:
* CCD chip number: PC=1; WFC=2,3,4
*Note on Xpos, Ypos:
* The mean error between the model centroid and the actual centroid of the
* object image is 0.2 pixels. Coordinates are relative to each individual chip.
*Note on Sky:
* A maximum likelihood estimate for the local sky background is determined
* simultaneously with the other model parameters. Sufficient pixels are used
* to ensure that the sky level is determined to an accuracy of 0.5%. The sky
* background is assumed to be flat over the small region selected for analysis.
*Note on mtot:
* Small differences between the analytic total magnitude and the true total
* magnitude may arise because the true galaxy is not smooth, and the model
* may not average properly over bright regions of star formation.
*Note on logR1/2:
* The radius within which half of the light of the unconvolved model would
* be contained if it were radially symmetric (an axis ratio of unity). Lower
* and upper limits of 0.1 pixel and one-half of the maximum radius of the
* region selected for analysis have been imposed. The half-light radius of
* the individual components can be derived used the bulge/disk half-radius
* ratio.
*Note on P.A.:
* The adopted position angle (measured east from north) is that of the axis
* of symmetry of the model. It is set equal to zero when the source is
* assumed to be azimuthally symmetric with an axis ratio of unity. For
* disk+bulge models, the orientation is assumed to be the same for both
* components.
*Note on b/aD, b/aB:
* The ratio of the minor axis half-light ratio to that of the major axis.
* For disk+bulge models it is defined independently for each component.
* If the axis is not significantly different from unity, it is fixed at 1.
* A ratio of zero indicates that this component was not fitted to the data.
*Note on B/(B+D):
* The luminosity of the individual component can be derived from this
* parameter and the total magnitude mtot.
*Note on Class:
* Object classification: object, galaxy, disk, bulge or disk+bulge.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2a.dat' ! Cl J0023+0423 field visually classified galaxies
integer*4 nr__2
parameter (nr__2=209) ! Number of records
character*122 ar__2 ! Full-size record
integer*4 v_OPL98__2 ! *MDS object identifier
real*4 mtot_2 ! (mag) *Analytic total magnitude
character*7 Class_2 ! *Hubble class
integer*4 D ! *Disturbance index
character*6 n_D ! *Interpretation of the disturbance
character*95 Comm ! *Description of object morphology
*Note on [OPL98]:
* [OPL98] MDS J0023+0423 NNN in Simbad for table2a.dat
* [OPL98] MDS J1604+4304 NNN in Simbad for table2b.dat
*Note on mtot:
* Total magnitude of best-fit analytic model in either R_702_ (Cl 0023+0423)
* or I_814_ (Cl 1604+4304) as given in table 1.
*Note on Class:
* The standard Hubble classification scheme (e.g. E, S0, Sa, Sab, Sb) with the
* addition of the following: E/S0 or S0/E, cannot distinguish between E or S0;
* X, compact object (likely nonstellar, but too compact to see structure);
* , peculiar or unclassifiable.
*Note on D:
* 0, normal; 1, moderate asymmetry; 2, strong asymmetric;
* 3, moderate distortion; 4, strong distortion.
*Note on n_D:
* M, merger; I, tidal interaction with neighbor; T, tidal feature; C, chaotic.
*Note on Comm:
* Here, comp. is short for "companion"; @ indicates position relative to the
* galaxy going clockwise, e.g. @7 = "at 7 o'clock".
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2b.dat' ! Cl J1604+4304 field visually classified galaxies
integer*4 nr__3
parameter (nr__3=205) ! Number of records
character*122 ar__3 ! Full-size record
integer*4 v_OPL98__3 ! *MDS object identifier
real*4 mtot_3 ! (mag) *Analytic total magnitude
character*7 Class_3 ! *Hubble class
integer*4 D_1 ! *Disturbance index
character*6 n_D_1 ! *Interpretation of the disturbance
character*95 Comm_1 ! *Description of object morphology
*Note on [OPL98]:
* [OPL98] MDS J0023+0423 NNN in Simbad for table2a.dat
* [OPL98] MDS J1604+4304 NNN in Simbad for table2b.dat
*Note on mtot:
* Total magnitude of best-fit analytic model in either R_702_ (Cl 0023+0423)
* or I_814_ (Cl 1604+4304) as given in table 1.
*Note on Class:
* The standard Hubble classification scheme (e.g. E, S0, Sa, Sab, Sb) with the
* addition of the following: E/S0 or S0/E, cannot distinguish between E or S0;
* X, compact object (likely nonstellar, but too compact to see structure);
* , peculiar or unclassifiable.
*Note on D:
* 0, normal; 1, moderate asymmetry; 2, strong asymmetric;
* 3, moderate distortion; 4, strong distortion.
*Note on n_D:
* M, merger; I, tidal interaction with neighbor; T, tidal feature; C, chaotic.
*Note on Comm:
* Here, comp. is short for "companion"; @ indicates position relative to the
* galaxy going clockwise, e.g. @7 = "at 7 o'clock".
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4a.dat' ! Cl J0023+0423 field photometry of the visually
classified galaxies
integer*4 nr__4
parameter (nr__4=209) ! Number of records
character*52 ar__4 ! Full-size record
integer*4 v_OPL98__4 ! *MDS identification number
integer*4 Keck ! *Keck Number
real*4 Bmag ! (mag) *B magnitude
real*4 Vmag ! (mag) *V magnitude
real*4 Rmag ! (mag) *R magnitude
real*4 Imag ! (mag) *I magnitude
real*4 z ! ? Redshift
character*1 l_Age ! limit flag on Age
real*4 Age ! (Gyr) *? Galaxy age
real*4 e_Age ! (Gyr) ? rms uncertainty on Age
*Note on [OPL98]:
* [OPL98] MDS J0023+0423 NNN in Simbad for table4a.dat
* [OPL98] MDS J1604+4304 NNN in Simbad for table4b.dat
*Note on Keck:
* Corresponding number from the Keck photometry tables (see Paper II,
* Cat. <J/AJ/116/560>). Because of the better resolution of the HST
* observations, an individual galaxy in the Keck image may be associated with
* more than one galaxy in the corresponding HST image. Therefore, these HST
* galaxies have the same Keck number.
*Note on Bmag, Vmag, Rmag, Imag:
* Keck BVRI aperture photometry. A value of -1.00 indicates no detection in
* that band
*Note on Age:
* Galaxy age determined from a comparison between the {tau}=0.6Gyr Bruzual &
* Charlot (1993ApJ...405..538B) population synthesis model and the broadband
* AB magnitudes derived from the Keck photometry.
* Upper limits are indicated as values with no errors.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4b.dat' ! Cl J1604+4304 field photometry of the visually
classified galaxies
integer*4 nr__5
parameter (nr__5=205) ! Number of records
character*52 ar__5 ! Full-size record
integer*4 v_OPL98__5 ! *MDS identification number
integer*4 Keck_1 ! *Keck Number
real*4 Bmag_1 ! (mag) *B magnitude
real*4 Vmag_1 ! (mag) *V magnitude
real*4 Rmag_1 ! (mag) *R magnitude
real*4 Imag_1 ! (mag) *I magnitude
real*4 z_1 ! ? Redshift
character*1 l_Age_1 ! limit flag on Age
real*4 Age_1 ! (Gyr) *? Galaxy age
real*4 e_Age_1 ! (Gyr) ? rms uncertainty on Age
*Note on [OPL98]:
* [OPL98] MDS J0023+0423 NNN in Simbad for table4a.dat
* [OPL98] MDS J1604+4304 NNN in Simbad for table4b.dat
*Note on Keck:
* Corresponding number from the Keck photometry tables (see Paper II,
* Cat. <J/AJ/116/560>). Because of the better resolution of the HST
* observations, an individual galaxy in the Keck image may be associated with
* more than one galaxy in the corresponding HST image. Therefore, these HST
* galaxies have the same Keck number.
*Note on Bmag, Vmag, Rmag, Imag:
* Keck BVRI aperture photometry. A value of -1.00 indicates no detection in
* that band
*Note on Age:
* Galaxy age determined from a comparison between the {tau}=0.6Gyr Bruzual &
* Charlot (1993ApJ...405..538B) population synthesis model and the broadband
* AB magnitudes derived from the Keck photometry.
* Upper limits are indicated as values with no errors.
C=============================================================================
C Loading file 'table1a.dat' ! Cl J0023+0423 field MDS parameters from the
* automated galaxy classifications
C Format for file interpretation
1 format(
+ I3,1X,I1,1X,I4,1X,I4,1X,I2,1X,I1,1X,F7.3,1X,F7.3,1X,F6.3,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F4.2,1X,A6,
+ 1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'table1a.dat')
write(6,*) '....Loading file: table1a.dat'
do i__=1,674
read(1,'(A150)')ar__
read(ar__,1)
+ v_OPL98_,Nmod,Ntot,Npix,Nfit,Chip,Xpos,Ypos,Sky,mtot,e_mtot,
+ logR1_2,e_logR1_2,P_A_,e_P_A_,b_aD,e_b_aD,b_aB,e_b_aB,B__B_D,
+ e_B__B_D,log_B_D_R1_2,e_log_B_D_R1_2,SNRIL,Class,Asymm
c ..............Just test output...........
write(6,1)
+ v_OPL98_,Nmod,Ntot,Npix,Nfit,Chip,Xpos,Ypos,Sky,mtot,e_mtot,
+ logR1_2,e_logR1_2,P_A_,e_P_A_,b_aD,e_b_aD,b_aB,e_b_aB,B__B_D,
+ e_B__B_D,log_B_D_R1_2,e_log_B_D_R1_2,SNRIL,Class,Asymm
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1b.dat' ! Cl J1604+4304 field MDS parameters from the
* automated galaxy classifications
C Format for file interpretation
2 format(
+ I3,1X,I1,1X,I4,1X,I4,1X,I2,1X,I1,1X,F7.3,1X,F7.3,1X,F6.3,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F4.2,1X,A6,
+ 1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'table1b.dat')
write(6,*) '....Loading file: table1b.dat'
do i__=1,559
read(1,'(A150)')ar__1
read(ar__1,2)
+ v_OPL98__1,Nmod_1,Ntot_1,Npix_1,Nfit_1,Chip_1,Xpos_1,Ypos_1,
+ Sky_1,mtot_1,e_mtot_1,logR1_2_1,e_logR1_2_1,P_A__1,e_P_A__1,
+ b_aD_1,e_b_aD_1,b_aB_1,e_b_aB_1,B__B_D_1,e_B__B_D_1,
+ log_B_D_R1_2_1,e_log_B_D_R1_2_1,SNRIL_1,Class_1,Asymm_1
c ..............Just test output...........
write(6,2)
+ v_OPL98__1,Nmod_1,Ntot_1,Npix_1,Nfit_1,Chip_1,Xpos_1,Ypos_1,
+ Sky_1,mtot_1,e_mtot_1,logR1_2_1,e_logR1_2_1,P_A__1,e_P_A__1,
+ b_aD_1,e_b_aD_1,b_aB_1,e_b_aB_1,B__B_D_1,e_B__B_D_1,
+ log_B_D_R1_2_1,e_log_B_D_R1_2_1,SNRIL_1,Class_1,Asymm_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2a.dat' ! Cl J0023+0423 field visually classified galaxies
C Format for file interpretation
3 format(I3,1X,F5.2,1X,A7,1X,I1,1X,A6,1X,A95)
C Effective file loading
open(unit=1,status='old',file=
+'table2a.dat')
write(6,*) '....Loading file: table2a.dat'
do i__=1,209
read(1,'(A122)')ar__2
read(ar__2,3)v_OPL98__2,mtot_2,Class_2,D,n_D,Comm
c ..............Just test output...........
write(6,3)v_OPL98__2,mtot_2,Class_2,D,n_D,Comm
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2b.dat' ! Cl J1604+4304 field visually classified galaxies
C Format for file interpretation
4 format(I3,1X,F5.2,1X,A7,1X,I1,1X,A6,1X,A95)
C Effective file loading
open(unit=1,status='old',file=
+'table2b.dat')
write(6,*) '....Loading file: table2b.dat'
do i__=1,205
read(1,'(A122)')ar__3
read(ar__3,4)v_OPL98__3,mtot_3,Class_3,D_1,n_D_1,Comm_1
c ..............Just test output...........
write(6,4)v_OPL98__3,mtot_3,Class_3,D_1,n_D_1,Comm_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4a.dat' ! Cl J0023+0423 field photometry of the visually
* classified galaxies
C Format for file interpretation
5 format(
+ I3,1X,I4,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F6.4,1X,A1,F4.2,
+ 1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4a.dat')
write(6,*) '....Loading file: table4a.dat'
do i__=1,209
read(1,'(A52)')ar__4
read(ar__4,5)
+ v_OPL98__4,Keck,Bmag,Vmag,Rmag,Imag,z,l_Age,Age,e_Age
if(ar__4(34:39) .EQ. '') z = rNULL__
if(ar__4(42:45) .EQ. '') Age = rNULL__
if(ar__4(47:50) .EQ. '') e_Age = rNULL__
c ..............Just test output...........
write(6,5)
+ v_OPL98__4,Keck,Bmag,Vmag,Rmag,Imag,z,l_Age,Age,e_Age
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4b.dat' ! Cl J1604+4304 field photometry of the visually
* classified galaxies
C Format for file interpretation
6 format(
+ I3,1X,I4,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F6.4,1X,A1,F4.2,
+ 1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4b.dat')
write(6,*) '....Loading file: table4b.dat'
do i__=1,205
read(1,'(A52)')ar__5
read(ar__5,6)
+ v_OPL98__5,Keck_1,Bmag_1,Vmag_1,Rmag_1,Imag_1,z_1,l_Age_1,
+ Age_1,e_Age_1
if(ar__5(34:39) .EQ. '') z_1 = rNULL__
if(ar__5(42:45) .EQ. '') Age_1 = rNULL__
if(ar__5(47:50) .EQ. '') e_Age_1 = rNULL__
c ..............Just test output...........
write(6,6)
+ v_OPL98__5,Keck_1,Bmag_1,Vmag_1,Rmag_1,Imag_1,z_1,l_Age_1,
+ Age_1,e_Age_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end