FORTRAN Generation
(/./ftp/cats/J/ApJ/607/426)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/607/426 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-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/ApJ/607/426       Catalog of white dwarfs in SDSS-DR1       (Kleinman+, 2004)
*================================================================================
*A catalog of spectroscopically identified white dwarf stars in the
*First Data Release of the Sloan Digital Sky Survey.
*    Kleinman S.J., Harris H.C., Eisenstein D.J., Liebert J., Nitta A.,
*    Krzesinski J., Munn J.A., Dahn C.C., Hawley S.L., Pier J.R., Schmidt G.,
*    Silvestri N.M., Smith J.A., Szkody P., Strauss M.A., Knapp G.R.,
*    Collinge M.J., Mukadam A.S., Koester D., Uomoto A., Schlegel D.J.,
*    Anderson S.F., Brinkmann J., Lamb D.Q., Schneider D.P., York D.G.
*   <Astrophys. J., 607, 426-444 (2004)>
*   =2004ApJ...607..426K
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! SDSS DR1 spectroscopically identified DA stars

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

C  J2000 position composed of: RAdeg DEdeg
      character*19  SDSS        ! Object name (G1)
      integer*4     Plate       ! Plate number
      integer*4     ObsDate     ! (d) MJD (JD-2400000.5) date of observation
      integer*4     Fiber       ! Fiber number
      real*8        RAdeg       ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg       ! (deg) Declination in decimal degrees (J2000)
      character*1   f_umag      ! [*] Indicates a questionable pipeline
*                                         reduction in umag
      real*4        umag        ! (mag) SDSS PSF u band magnitude
      real*4        e_umag      ! (mag) Uncertainty in umag
      character*1   f_gmag      ! [*] Indicates a questionable pipeline
*                                         reduction in gmag
      real*4        gmag        ! (mag) SDSS PSF g band magnitude
      real*4        e_gmag      ! (mag) Uncertainty in gmag
      character*1   f_rmag      ! [*] Indicates a questionable pipeline
*                                         reduction in rmag
      real*4        rmag        ! (mag) SDSS PSF r band magnitude
      real*4        e_rmag      ! (mag) Uncertainty in rmag
      character*1   f_imag      ! [*] Indicates a questionable pipeline
*                                         reduction in imag
      real*4        imag        ! (mag) SDSS PSF i band magnitude
      real*4        e_imag      ! (mag) Uncertainty in imag
      character*1   f_zmag      ! [*] Indicates a questionable pipeline
*                                         reduction in zmag
      real*4        zmag        ! (mag) SDSS PSF z band magnitude
      real*4        e_zmag      ! (mag) Uncertainty in zmag
      real*4        S_N         ! Final median S/N squared per pixel (2)
      integer*4     pmRA        ! (mas/yr) ? Proper motion in Right Ascension (3)
      integer*4     pmDE        ! (mas/yr) ? Proper motion in Declination (3)
      real*4        Ag          ! (mag) The g band extinction
      integer*4     Teff        ! (K) Model determined effective temperature
      integer*4     e_Teff      ! (K) Uncertainty in Teff
      real*4        log_g       ! ([cm/s2]) Log of the model determined surface gravity
      real*4        e_log_g     ! ([cm/s2]) Uncertainty in log(g)
      real*4        Chi2        ! The {chi}^2^ fit
      character*14  AutoID      ! The ASCII identification (4)
      character*3   Notes       ! Additional notes (G2)
*Note (2): Of the spectrum from the DR1 best database.
*Note (3): NULL value indicates N/A.
*Note (4): The format is simply DA or DB followed by the Teff subtype
*          (50400K/Teff) and the log(g)s, separated by an underscore.
*    The CH: prefix and colon suffix are quality assessment checks which
*    are discussed in Section 5. Where our check on the auto fitting deemed
*    the fit invalid, we replace the fit parameters with zeroes and the fit
*    AutoID with N/A.

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

C  Declarations for 'table4.dat'	! SDSS DR1 spectroscopically identified DB stars

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

C  J2000 position composed of: RAdeg DEdeg
      character*19  SDSS_1      ! Object name (G1)
      integer*4     Plate_1     ! Plate number
      integer*4     ObsDate_1   ! (d) MJD (JD-2400000.5) date of observation
      integer*4     Fiber_1     ! Fiber number
      real*8        RAdeg_1     ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_1     ! (deg) Declination in decimal degrees (J2000)
      character*1   f_umag_1    ! [*] Indicates a questionable pipeline
*                                         reduction in umag
      real*4        umag_1      ! (mag) SDSS PSF u band magnitude
      real*4        e_umag_1    ! (mag) Uncertainty in umag
      character*1   f_gmag_1    ! [*] Indicates a questionable pipeline
*                                         reduction in gmag
      real*4        gmag_1      ! (mag) SDSS PSF g band magnitude
      real*4        e_gmag_1    ! (mag) Uncertainty in gmag
      character*1   f_rmag_1    ! [*] Indicates a questionable pipeline
*                                         reduction in rmag
      real*4        rmag_1      ! (mag) SDSS PSF r band magnitude
      real*4        e_rmag_1    ! (mag) Uncertainty in rmag
      character*1   f_imag_1    ! [*] Indicates a questionable pipeline
*                                         reduction in imag
      real*4        imag_1      ! (mag) SDSS PSF i band magnitude
      real*4        e_imag_1    ! (mag) Uncertainty in imag
      character*1   f_zmag_1    ! [*] Indicates a questionable pipeline
*                                         reduction in zmag
      real*4        zmag_1      ! (mag) SDSS PSF z band magnitude
      real*4        e_zmag_1    ! (mag) Uncertainty in zmag
      real*4        S_N_1       ! Final median S/N squared per pixel (2)
      integer*4     pmRA_1      ! (mas/yr) ? Proper motion in Right Ascension (3)
      integer*4     pmDE_1      ! (mas/yr) ? Proper motion in Declination (3)
      real*4        Ag_1        ! (mag) The g band extinction
      integer*4     Teff_1      ! (K) Model determined effective temperature
      integer*4     e_Teff_1    ! (K) Uncertainty in Teff
      real*4        log_g_1     ! ([cm/s2]) Log of the model determined surface gravity
      real*4        e_log_g_1   ! ([cm/s2]) Uncertainty in log(g)
      real*4        Chi2_1      ! The {chi}^2^ fit
      character*14  AutoID_1    ! The ASCII identification (4)
      character*3   Notes_1     ! Additional notes (G2)
*Note (2): Of the spectrum from the DR1 best database.
*Note (3): NULL value indicates N/A.
*Note (4): The format is simply DA or DB followed by the Teff subtype
*          (50400K/Teff) and the log(g)s, separated by an underscore.
*    The CH: prefix and colon suffix are quality assessment checks which
*    are discussed in Section 5. Where our check on the auto fitting deemed
*    the fit invalid, we replace the fit parameters with zeroes and the fit
*    AutoID with N/A.

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

C  Declarations for 'table5.dat'	! All SDSS DR1 spectroscopically identified white dwarf
                           and subdwarf stars

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

C  J2000 position composed of: RAdeg DEdeg
      character*19  SDSS_2      ! Object name (G1)
      integer*4     Plate_2     ! Plate number
      integer*4     ObsDate_2   ! (d) MJD (JD-2400000.5) date of observation
      integer*4     Fiber_2     ! Fiber number
      real*8        RAdeg_2     ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_2     ! (deg) Declination in decimal degrees (J2000)
      real*8        Epoch       ! (yr) Epoch of the observation
      character*1   f_umag_2    ! [*] Indicates a questionable pipeline
*                                       reduction in umag
      real*4        umag_2      ! (mag) SDSS PSF u band magnitude
      real*4        e_umag_2    ! (mag) Uncertainty in umag
      character*1   f_gmag_2    ! [*] Indicates a questionable pipeline
*                                       reduction in gmag
      real*4        gmag_2      ! (mag) SDSS PSF g band magnitude
      real*4        e_gmag_2    ! (mag) Uncertainty in gmag
      character*1   f_rmag_2    ! [*] Indicates a questionable pipeline
*                                       reduction in rmag
      real*4        rmag_2      ! (mag) SDSS PSF r band magnitude
      real*4        e_rmag_2    ! (mag) Uncertainty in rmag
      character*1   f_imag_2    ! [*] Indicates a questionable pipeline
*                                       reduction in imag
      real*4        imag_2      ! (mag) SDSS PSF i band magnitude
      real*4        e_imag_2    ! (mag) Uncertainty in imag
      character*1   f_zmag_2    ! [*] Indicates a questionable pipeline
*                                       reduction in zmag
      real*4        zmag_2      ! (mag) SDSS PSF z band magnitude
      real*4        e_zmag_2    ! (mag) Uncertainty in zmag
      real*4        S_N_2       ! Final median S/N squared per pixel (2)
      integer*4     pmRA_2      ! (mas/yr) ? Proper motion in Right Ascension (3)
      integer*4     pmDE_2      ! (mas/yr) ? Proper motion in Declination (3)
      real*4        Ag_2        ! (mag) The g band extinction
      character*16  ID          ! The human identification of spectral type
      character*3   Notes_2     ! Additional notes (G2)
*Note (2): Of the spectrum from the DR1 best database.
*Note (3): NULL value indicates N/A.

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

C  Declarations for 'table6.dat'	! All SDSS DR1 spectroscopically identified white dwarf
                           and subdwarf stars with uncertain identifications

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

C  J2000 position composed of: RAdeg DEdeg
      character*19  SDSS_3      ! Object name (G1)
      integer*4     Plate_3     ! Plate number
      integer*4     ObsDate_3   ! (d) MJD (JD-2400000.5) date of observation
      integer*4     Fiber_3     ! Fiber number
      real*8        RAdeg_3     ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_3     ! (deg) Declination in decimal degrees (J2000)
      real*8        Epoch_1     ! (yr) Epoch of the observation
      character*1   f_umag_3    ! [*] Indicates a questionable pipeline
*                                       reduction in umag
      real*4        umag_3      ! (mag) SDSS PSF u band magnitude
      real*4        e_umag_3    ! (mag) Uncertainty in umag
      character*1   f_gmag_3    ! [*] Indicates a questionable pipeline
*                                       reduction in gmag
      real*4        gmag_3      ! (mag) SDSS PSF g band magnitude
      real*4        e_gmag_3    ! (mag) Uncertainty in gmag
      character*1   f_rmag_3    ! [*] Indicates a questionable pipeline
*                                       reduction in rmag
      real*4        rmag_3      ! (mag) SDSS PSF r band magnitude
      real*4        e_rmag_3    ! (mag) Uncertainty in rmag
      character*1   f_imag_3    ! [*] Indicates a questionable pipeline
*                                       reduction in imag
      real*4        imag_3      ! (mag) SDSS PSF i band magnitude
      real*4        e_imag_3    ! (mag) Uncertainty in imag
      character*1   f_zmag_3    ! [*] Indicates a questionable pipeline
*                                       reduction in zmag
      real*4        zmag_3      ! (mag) SDSS PSF z band magnitude
      real*4        e_zmag_3    ! (mag) Uncertainty in zmag
      real*4        S_N_3       ! Final median S/N squared per pixel (2)
      integer*4     pmRA_3      ! (mas/yr) ? Proper motion in Right Ascension (3)
      integer*4     pmDE_3      ! (mas/yr) ? Proper motion in Declination (3)
      real*4        Ag_3        ! (mag) The g band extinction
      character*16  ID_1        ! The human identification of spectral type
      character*3   Notes_3     ! Additional notes (G2)
*Note (2): Of the spectrum from the DR1 best database.
*Note (3): NULL value indicates N/A.

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

C  Declarations for 'table11.dat'	! Previously known white dwarf stars recovered in
                           this work

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

      character*19  SDSS_4      ! Object name (G1)
      real*4        gmag_4      ! (mag) SDSS PSF g band magnitude
      character*6   SType       ! SDSS type
      character*12  OName       ! Other name
      character*7   OType       ! Other type
      character*1   Note        ! Additional notes (2)
*Note (2): Additional notes as follows:
*      1 = SDSS position different from previous by more than 10 arcsec.
*      2 = Unpublished white dwarf, in Hamburg Quasar Survey
*          (Homeier 2002, priv. comm.)
*      3 = White dwarf discovered in the SDSS EDR
*          (Gansicke et al., 2002A&A...394..957G)
*      4 = White dwarf discovered in the SDSS
*          (Harris et al., 2001ApJ...549L.109H)

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

C  Declarations for 'table12.dat'	! Previously known white dwarf stars not
                           spectroscopically recovered in this work

      integer*4 nr__5
      parameter (nr__5=218)	! Number of records
      character*115 ar__5  	! Full-size record

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

      character*10  WD          ! Object name from WD catalog (III/235)
      character*6   Type        ! Object type
      real*4        mag         ! (mag) ? Observed magnitude
      character*1   Filt        ! Filter used to obtain mag
      integer*4     RAh         ! (h) ? Hour of Right Ascension (J2000)
      integer*4     RAm         ! (min) ? Minute of Right Ascension (J2000)
      real*4        RAs         ! (s) ? Second of Right Ascension (J2000)
      character*1   DE_         ! Sign of the Declination (J2000)
      integer*4     DEd         ! (deg) ? Degree of Declination (J2000)
      integer*4     DEm         ! (arcmin) ? Arcminute of Declination (J2000)
      real*4        DEs         ! (arcsec) ? Arcsecond of Declination (J2000)
      character*1   f_umag_4    ! [*] Indicates saturated pixels in umag
      real*4        umag_4      ! (mag) ? SDSS PSF u band magnitude
      real*4        e_umag_4    ! (mag) ? Uncertainty in umag
      character*1   f_gmag_4    ! [*] Indicates saturated pixels in gmag
      real*4        gmag_5      ! (mag) ? SDSS PSF g band magnitude
      real*4        e_gmag_4    ! (mag) ? Uncertainty in gmag
      character*1   f_rmag_4    ! [*] Indicates saturated pixels in rmag
      real*4        rmag_4      ! (mag) ? SDSS PSF r band magnitude
      real*4        e_rmag_4    ! (mag) ? Uncertainty in rmag
      character*1   f_imag_4    ! [*] Indicates saturated pixels in imag
      real*4        imag_4      ! (mag) ? SDSS PSF i band magnitude
      real*4        e_imag_4    ! (mag) ? Uncertainty in imag
      character*1   f_zmag_4    ! [*] Indicates saturated pixels in zmag
      real*4        zmag_4      ! (mag) ? SDSS PSF z band magnitude
      real*4        e_zmag_4    ! (mag) ? Uncertainty in zmag
      character*2   Note_1      ! Additional note (1)
*Note (1): Additional notes as follows:
*      1 = No white dwarf found near this position.
*      2 = SDSS spectrum 0410-51816-565 shows WD0255+009.1 is a QSO.
*      3 = No SDSS imaging data, too near bright star.
*      4 = Probably same star as WD0330-009.
*      5 = WD0820+021 probably not detected, too faint.
*      6 = Spectrum not in DR1, but is given in
*          Initial Survey paper (Harris et al., 2003, Cat. <J/AJ/126/1023>).
*      7 = Falls in small gap in SDSS imaging data.
*      8 = Unresolved with WD1330+015.1.
*      9 = The colors of WD1401+005 are red and indicate it is not a
*          white dwarf.
*     10 = Probably same star as WD1422+028.
*     11 = Colors are very red and indicate WD1449+003 is not a white dwarf.
*     12 = Colors indicate WD1451--004, WD1455+019, and WD1500+006 are
*          horizontal-branch stars, not white dwarf stars.
*     13 = Probably same star as WD1544+009.

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

C  Loading file 'table3.dat'	! SDSS DR1 spectroscopically identified DA stars

C  Format for file interpretation

    1 format(
     +  5X,A19,1X,I4,1X,I5,1X,I3,1X,F9.5,1X,F9.5,1X,A1,F5.2,1X,F4.2,
     +  1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,
     +  A1,F5.2,1X,F4.2,1X,F5.2,1X,I4,1X,I4,1X,F5.3,1X,I5,1X,I5,1X,
     +  F4.2,1X,F4.2,1X,F4.2,1X,A14,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,1888
        read(1,'(A188)')ar__
        read(ar__,1)
     +  SDSS,Plate,ObsDate,Fiber,RAdeg,DEdeg,f_umag,umag,e_umag,
     +  f_gmag,gmag,e_gmag,f_rmag,rmag,e_rmag,f_imag,imag,e_imag,
     +  f_zmag,zmag,e_zmag,S_N,pmRA,pmDE,Ag,Teff,e_Teff,log_g,e_log_g,
     +  Chi2,AutoID,Notes
        if(ar__(127:130) .EQ. '') pmRA = iNULL__
        if(ar__(132:135) .EQ. '') pmDE = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  SDSS,Plate,ObsDate,Fiber,RAdeg,DEdeg,f_umag,umag,e_umag,
     +  f_gmag,gmag,e_gmag,f_rmag,rmag,e_rmag,f_imag,imag,e_imag,
     +  f_zmag,zmag,e_zmag,S_N,pmRA,pmDE,Ag,Teff,e_Teff,log_g,e_log_g,
     +  Chi2,AutoID,Notes
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! SDSS DR1 spectroscopically identified DB stars

C  Format for file interpretation

    2 format(
     +  5X,A19,1X,I4,1X,I5,1X,I3,1X,F9.5,1X,F9.5,1X,A1,F5.2,1X,F4.2,
     +  1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,
     +  A1,F5.2,1X,F4.2,1X,F5.2,1X,I4,1X,I4,1X,F5.3,1X,I5,1X,I5,1X,
     +  F4.2,1X,F4.2,1X,F4.2,1X,A14,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,171
        read(1,'(A189)')ar__1
        read(ar__1,2)
     +  SDSS_1,Plate_1,ObsDate_1,Fiber_1,RAdeg_1,DEdeg_1,f_umag_1,
     +  umag_1,e_umag_1,f_gmag_1,gmag_1,e_gmag_1,f_rmag_1,rmag_1,
     +  e_rmag_1,f_imag_1,imag_1,e_imag_1,f_zmag_1,zmag_1,e_zmag_1,
     +  S_N_1,pmRA_1,pmDE_1,Ag_1,Teff_1,e_Teff_1,log_g_1,e_log_g_1,
     +  Chi2_1,AutoID_1,Notes_1
        if(ar__1(127:130) .EQ. '') pmRA_1 = iNULL__
        if(ar__1(132:135) .EQ. '') pmDE_1 = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  SDSS_1,Plate_1,ObsDate_1,Fiber_1,RAdeg_1,DEdeg_1,f_umag_1,
     +  umag_1,e_umag_1,f_gmag_1,gmag_1,e_gmag_1,f_rmag_1,rmag_1,
     +  e_rmag_1,f_imag_1,imag_1,e_imag_1,f_zmag_1,zmag_1,e_zmag_1,
     +  S_N_1,pmRA_1,pmDE_1,Ag_1,Teff_1,e_Teff_1,log_g_1,e_log_g_1,
     +  Chi2_1,AutoID_1,Notes_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! All SDSS DR1 spectroscopically identified white dwarf
*                           and subdwarf stars

C  Format for file interpretation

    3 format(
     +  5X,A19,1X,I4,1X,I5,1X,I3,1X,F9.5,1X,F9.5,1X,F7.2,1X,A1,F5.2,
     +  1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
     +  F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,I4,1X,I4,1X,F5.3,1X,A16,1X,
     +  A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,2791
        read(1,'(A175)')ar__2
        read(ar__2,3)
     +  SDSS_2,Plate_2,ObsDate_2,Fiber_2,RAdeg_2,DEdeg_2,Epoch,
     +  f_umag_2,umag_2,e_umag_2,f_gmag_2,gmag_2,e_gmag_2,f_rmag_2,
     +  rmag_2,e_rmag_2,f_imag_2,imag_2,e_imag_2,f_zmag_2,zmag_2,
     +  e_zmag_2,S_N_2,pmRA_2,pmDE_2,Ag_2,ID,Notes_2
        if(ar__2(135:138) .EQ. '') pmRA_2 = iNULL__
        if(ar__2(140:143) .EQ. '') pmDE_2 = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  SDSS_2,Plate_2,ObsDate_2,Fiber_2,RAdeg_2,DEdeg_2,Epoch,
     +  f_umag_2,umag_2,e_umag_2,f_gmag_2,gmag_2,e_gmag_2,f_rmag_2,
     +  rmag_2,e_rmag_2,f_imag_2,imag_2,e_imag_2,f_zmag_2,zmag_2,
     +  e_zmag_2,S_N_2,pmRA_2,pmDE_2,Ag_2,ID,Notes_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! All SDSS DR1 spectroscopically identified white dwarf
*                           and subdwarf stars with uncertain identifications

C  Format for file interpretation

    4 format(
     +  5X,A19,1X,I4,1X,I5,1X,I3,1X,F9.5,1X,F9.5,1X,F7.2,1X,A1,F5.2,
     +  1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
     +  F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,I4,1X,I4,1X,F5.3,1X,A16,1X,
     +  A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,144
        read(1,'(A170)')ar__3
        read(ar__3,4)
     +  SDSS_3,Plate_3,ObsDate_3,Fiber_3,RAdeg_3,DEdeg_3,Epoch_1,
     +  f_umag_3,umag_3,e_umag_3,f_gmag_3,gmag_3,e_gmag_3,f_rmag_3,
     +  rmag_3,e_rmag_3,f_imag_3,imag_3,e_imag_3,f_zmag_3,zmag_3,
     +  e_zmag_3,S_N_3,pmRA_3,pmDE_3,Ag_3,ID_1,Notes_3
        if(ar__3(135:138) .EQ. '') pmRA_3 = iNULL__
        if(ar__3(140:143) .EQ. '') pmDE_3 = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  SDSS_3,Plate_3,ObsDate_3,Fiber_3,RAdeg_3,DEdeg_3,Epoch_1,
     +  f_umag_3,umag_3,e_umag_3,f_gmag_3,gmag_3,e_gmag_3,f_rmag_3,
     +  rmag_3,e_rmag_3,f_imag_3,imag_3,e_imag_3,f_zmag_3,zmag_3,
     +  e_zmag_3,S_N_3,pmRA_3,pmDE_3,Ag_3,ID_1,Notes_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table11.dat'	! Previously known white dwarf stars recovered in
*                           this work

C  Format for file interpretation

    5 format(5X,A19,1X,F5.2,1X,A6,1X,A12,1X,A7,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table11.dat')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,109
        read(1,'(A61)')ar__4
        read(ar__4,5)SDSS_4,gmag_4,SType,OName,OType,Note
c    ..............Just test output...........
        write(6,5)SDSS_4,gmag_4,SType,OName,OType,Note
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table12.dat'	! Previously known white dwarf stars not
*                           spectroscopically recovered in this work

C  Format for file interpretation

    6 format(
     +  2X,A10,1X,A6,1X,F5.2,1X,A1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,
     +  1X,F4.1,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
     +  F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table12.dat')
      write(6,*) '....Loading file: table12.dat'
      do i__=1,218
        read(1,'(A115)')ar__5
        read(ar__5,6)
     +  WD,Type,mag,Filt,RAh,RAm,RAs,DE_,DEd,DEm,DEs,f_umag_4,umag_4,
     +  e_umag_4,f_gmag_4,gmag_5,e_gmag_4,f_rmag_4,rmag_4,e_rmag_4,
     +  f_imag_4,imag_4,e_imag_4,f_zmag_4,zmag_4,e_zmag_4,Note_1
        if(ar__5(21:25) .EQ. '') mag = rNULL__
        if(ar__5(29:30) .EQ. '') RAh = iNULL__
        if(ar__5(32:33) .EQ. '') RAm = iNULL__
        if(ar__5(35:39) .EQ. '') RAs = rNULL__
        if(ar__5(42:43) .EQ. '') DEd = iNULL__
        if(ar__5(45:46) .EQ. '') DEm = iNULL__
        if(ar__5(48:51) .EQ. '') DEs = rNULL__
        if(ar__5(54:58) .EQ. '') umag_4 = rNULL__
        if(ar__5(60:63) .EQ. '') e_umag_4 = rNULL__
        if(ar__5(66:70) .EQ. '') gmag_5 = rNULL__
        if(ar__5(72:75) .EQ. '') e_gmag_4 = rNULL__
        if(ar__5(78:82) .EQ. '') rmag_4 = rNULL__
        if(ar__5(84:87) .EQ. '') e_rmag_4 = rNULL__
        if(ar__5(90:94) .EQ. '') imag_4 = rNULL__
        if(ar__5(96:99) .EQ. '') e_imag_4 = rNULL__
        if(ar__5(102:106) .EQ. '') zmag_4 = rNULL__
        if(ar__5(108:111) .EQ. '') e_zmag_4 = rNULL__
        RAdeg_4 = rNULL__
        DEdeg_4 = rNULL__
c  Derive coordinates RAdeg_4 and DEdeg_4 from input data
c  (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg_4=RAh*15.
        if(RAm .GT. -180) RAdeg_4=RAdeg_4+RAm/4.
        if(RAs .GT. -180) RAdeg_4=RAdeg_4+RAs/240.
        if(DEd .GE. 0) DEdeg_4=DEd
        if(DEm .GE. 0) DEdeg_4=DEdeg_4+DEm/60.
        if(DEs .GE. 0) DEdeg_4=DEdeg_4+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg_4.GE.0) DEdeg_4=-DEdeg_4
c    ..............Just test output...........
        write(6,6)
     +  WD,Type,mag,Filt,RAh,RAm,RAs,DE_,DEd,DEm,DEs,f_umag_4,umag_4,
     +  e_umag_4,f_gmag_4,gmag_5,e_gmag_4,f_rmag_4,rmag_4,e_rmag_4,
     +  f_imag_4,imag_4,e_imag_4,f_zmag_4,zmag_4,e_zmag_4,Note_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_4,DEdeg_4
c    .......End.of.Just test output...........
      end do
      close(1)

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