Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/607/426 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-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 (nr__) ! Object name (G1)
integer*4 Plate (nr__) ! Plate number
integer*4 ObsDate (nr__) ! (d) MJD (JD-2400000.5) date of observation
integer*4 Fiber (nr__) ! Fiber number
real*8 RAdeg (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
character*1 f_umag (nr__) ! [*] Indicates a questionable pipeline
* reduction in umag
real*4 umag (nr__) ! (mag) SDSS PSF u band magnitude
real*4 e_umag (nr__) ! (mag) Uncertainty in umag
character*1 f_gmag (nr__) ! [*] Indicates a questionable pipeline
* reduction in gmag
real*4 gmag (nr__) ! (mag) SDSS PSF g band magnitude
real*4 e_gmag (nr__) ! (mag) Uncertainty in gmag
character*1 f_rmag (nr__) ! [*] Indicates a questionable pipeline
* reduction in rmag
real*4 rmag (nr__) ! (mag) SDSS PSF r band magnitude
real*4 e_rmag (nr__) ! (mag) Uncertainty in rmag
character*1 f_imag (nr__) ! [*] Indicates a questionable pipeline
* reduction in imag
real*4 imag (nr__) ! (mag) SDSS PSF i band magnitude
real*4 e_imag (nr__) ! (mag) Uncertainty in imag
character*1 f_zmag (nr__) ! [*] Indicates a questionable pipeline
* reduction in zmag
real*4 zmag (nr__) ! (mag) SDSS PSF z band magnitude
real*4 e_zmag (nr__) ! (mag) Uncertainty in zmag
real*4 S_N (nr__) ! Final median S/N squared per pixel (2)
integer*4 pmRA (nr__) ! (mas/yr) ? Proper motion in Right Ascension (3)
integer*4 pmDE (nr__) ! (mas/yr) ? Proper motion in Declination (3)
real*4 Ag (nr__) ! (mag) The g band extinction
integer*4 Teff (nr__) ! (K) Model determined effective temperature
integer*4 e_Teff (nr__) ! (K) Uncertainty in Teff
real*4 log_g (nr__) ! ([cm/s2]) Log of the model determined surface gravity
real*4 e_log_g (nr__) ! ([cm/s2]) Uncertainty in log(g)
real*4 Chi2 (nr__) ! The {chi}^2^ fit
character*14 AutoID (nr__) ! The ASCII identification (4)
character*3 Notes (nr__) ! 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 (nr__1) ! Object name (G1)
integer*4 Plate_1 (nr__1) ! Plate number
integer*4 ObsDate_1 (nr__1) ! (d) MJD (JD-2400000.5) date of observation
integer*4 Fiber_1 (nr__1) ! Fiber number
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination in decimal degrees (J2000)
character*1 f_umag_1 (nr__1) ! [*] Indicates a questionable pipeline
* reduction in umag
real*4 umag_1 (nr__1) ! (mag) SDSS PSF u band magnitude
real*4 e_umag_1 (nr__1) ! (mag) Uncertainty in umag
character*1 f_gmag_1 (nr__1) ! [*] Indicates a questionable pipeline
* reduction in gmag
real*4 gmag_1 (nr__1) ! (mag) SDSS PSF g band magnitude
real*4 e_gmag_1 (nr__1) ! (mag) Uncertainty in gmag
character*1 f_rmag_1 (nr__1) ! [*] Indicates a questionable pipeline
* reduction in rmag
real*4 rmag_1 (nr__1) ! (mag) SDSS PSF r band magnitude
real*4 e_rmag_1 (nr__1) ! (mag) Uncertainty in rmag
character*1 f_imag_1 (nr__1) ! [*] Indicates a questionable pipeline
* reduction in imag
real*4 imag_1 (nr__1) ! (mag) SDSS PSF i band magnitude
real*4 e_imag_1 (nr__1) ! (mag) Uncertainty in imag
character*1 f_zmag_1 (nr__1) ! [*] Indicates a questionable pipeline
* reduction in zmag
real*4 zmag_1 (nr__1) ! (mag) SDSS PSF z band magnitude
real*4 e_zmag_1 (nr__1) ! (mag) Uncertainty in zmag
real*4 S_N_1 (nr__1) ! Final median S/N squared per pixel (2)
integer*4 pmRA_1 (nr__1) ! (mas/yr) ? Proper motion in Right Ascension (3)
integer*4 pmDE_1 (nr__1) ! (mas/yr) ? Proper motion in Declination (3)
real*4 Ag_1 (nr__1) ! (mag) The g band extinction
integer*4 Teff_1 (nr__1) ! (K) Model determined effective temperature
integer*4 e_Teff_1 (nr__1) ! (K) Uncertainty in Teff
real*4 log_g_1 (nr__1) ! ([cm/s2]) Log of the model determined surface gravity
real*4 e_log_g_1 (nr__1) ! ([cm/s2]) Uncertainty in log(g)
real*4 Chi2_1 (nr__1) ! The {chi}^2^ fit
character*14 AutoID_1 (nr__1) ! The ASCII identification (4)
character*3 Notes_1 (nr__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 (nr__2) ! Object name (G1)
integer*4 Plate_2 (nr__2) ! Plate number
integer*4 ObsDate_2 (nr__2) ! (d) MJD (JD-2400000.5) date of observation
integer*4 Fiber_2 (nr__2) ! Fiber number
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination in decimal degrees (J2000)
real*8 Epoch (nr__2) ! (yr) Epoch of the observation
character*1 f_umag_2 (nr__2) ! [*] Indicates a questionable pipeline
* reduction in umag
real*4 umag_2 (nr__2) ! (mag) SDSS PSF u band magnitude
real*4 e_umag_2 (nr__2) ! (mag) Uncertainty in umag
character*1 f_gmag_2 (nr__2) ! [*] Indicates a questionable pipeline
* reduction in gmag
real*4 gmag_2 (nr__2) ! (mag) SDSS PSF g band magnitude
real*4 e_gmag_2 (nr__2) ! (mag) Uncertainty in gmag
character*1 f_rmag_2 (nr__2) ! [*] Indicates a questionable pipeline
* reduction in rmag
real*4 rmag_2 (nr__2) ! (mag) SDSS PSF r band magnitude
real*4 e_rmag_2 (nr__2) ! (mag) Uncertainty in rmag
character*1 f_imag_2 (nr__2) ! [*] Indicates a questionable pipeline
* reduction in imag
real*4 imag_2 (nr__2) ! (mag) SDSS PSF i band magnitude
real*4 e_imag_2 (nr__2) ! (mag) Uncertainty in imag
character*1 f_zmag_2 (nr__2) ! [*] Indicates a questionable pipeline
* reduction in zmag
real*4 zmag_2 (nr__2) ! (mag) SDSS PSF z band magnitude
real*4 e_zmag_2 (nr__2) ! (mag) Uncertainty in zmag
real*4 S_N_2 (nr__2) ! Final median S/N squared per pixel (2)
integer*4 pmRA_2 (nr__2) ! (mas/yr) ? Proper motion in Right Ascension (3)
integer*4 pmDE_2 (nr__2) ! (mas/yr) ? Proper motion in Declination (3)
real*4 Ag_2 (nr__2) ! (mag) The g band extinction
character*16 ID (nr__2) ! The human identification of spectral type
character*3 Notes_2 (nr__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 (nr__3) ! Object name (G1)
integer*4 Plate_3 (nr__3) ! Plate number
integer*4 ObsDate_3 (nr__3) ! (d) MJD (JD-2400000.5) date of observation
integer*4 Fiber_3 (nr__3) ! Fiber number
real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_3 (nr__3) ! (deg) Declination in decimal degrees (J2000)
real*8 Epoch_1 (nr__3) ! (yr) Epoch of the observation
character*1 f_umag_3 (nr__3) ! [*] Indicates a questionable pipeline
* reduction in umag
real*4 umag_3 (nr__3) ! (mag) SDSS PSF u band magnitude
real*4 e_umag_3 (nr__3) ! (mag) Uncertainty in umag
character*1 f_gmag_3 (nr__3) ! [*] Indicates a questionable pipeline
* reduction in gmag
real*4 gmag_3 (nr__3) ! (mag) SDSS PSF g band magnitude
real*4 e_gmag_3 (nr__3) ! (mag) Uncertainty in gmag
character*1 f_rmag_3 (nr__3) ! [*] Indicates a questionable pipeline
* reduction in rmag
real*4 rmag_3 (nr__3) ! (mag) SDSS PSF r band magnitude
real*4 e_rmag_3 (nr__3) ! (mag) Uncertainty in rmag
character*1 f_imag_3 (nr__3) ! [*] Indicates a questionable pipeline
* reduction in imag
real*4 imag_3 (nr__3) ! (mag) SDSS PSF i band magnitude
real*4 e_imag_3 (nr__3) ! (mag) Uncertainty in imag
character*1 f_zmag_3 (nr__3) ! [*] Indicates a questionable pipeline
* reduction in zmag
real*4 zmag_3 (nr__3) ! (mag) SDSS PSF z band magnitude
real*4 e_zmag_3 (nr__3) ! (mag) Uncertainty in zmag
real*4 S_N_3 (nr__3) ! Final median S/N squared per pixel (2)
integer*4 pmRA_3 (nr__3) ! (mas/yr) ? Proper motion in Right Ascension (3)
integer*4 pmDE_3 (nr__3) ! (mas/yr) ? Proper motion in Declination (3)
real*4 Ag_3 (nr__3) ! (mag) The g band extinction
character*16 ID_1 (nr__3) ! The human identification of spectral type
character*3 Notes_3 (nr__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 (nr__4) ! Object name (G1)
real*4 gmag_4 (nr__4) ! (mag) SDSS PSF g band magnitude
character*6 SType (nr__4) ! SDSS type
character*12 OName (nr__4) ! Other name
character*7 OType (nr__4) ! Other type
character*1 Note (nr__4) ! 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 (nr__5) ! (deg) Right Ascension J2000
real*8 DEdeg_4 (nr__5) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*10 WD (nr__5) ! Object name from WD catalog (III/235)
character*6 Type (nr__5) ! Object type
real*4 mag (nr__5) ! (mag) ? Observed magnitude
character*1 Filt (nr__5) ! Filter used to obtain mag
integer*4 RAh (nr__5) ! (h) ? Hour of Right Ascension (J2000)
integer*4 RAm (nr__5) ! (min) ? Minute of Right Ascension (J2000)
real*4 RAs (nr__5) ! (s) ? Second of Right Ascension (J2000)
character*1 DE_ (nr__5) ! Sign of the Declination (J2000)
integer*4 DEd (nr__5) ! (deg) ? Degree of Declination (J2000)
integer*4 DEm (nr__5) ! (arcmin) ? Arcminute of Declination (J2000)
real*4 DEs (nr__5) ! (arcsec) ? Arcsecond of Declination (J2000)
character*1 f_umag_4 (nr__5) ! [*] Indicates saturated pixels in umag
real*4 umag_4 (nr__5) ! (mag) ? SDSS PSF u band magnitude
real*4 e_umag_4 (nr__5) ! (mag) ? Uncertainty in umag
character*1 f_gmag_4 (nr__5) ! [*] Indicates saturated pixels in gmag
real*4 gmag_5 (nr__5) ! (mag) ? SDSS PSF g band magnitude
real*4 e_gmag_4 (nr__5) ! (mag) ? Uncertainty in gmag
character*1 f_rmag_4 (nr__5) ! [*] Indicates saturated pixels in rmag
real*4 rmag_4 (nr__5) ! (mag) ? SDSS PSF r band magnitude
real*4 e_rmag_4 (nr__5) ! (mag) ? Uncertainty in rmag
character*1 f_imag_4 (nr__5) ! [*] Indicates saturated pixels in imag
real*4 imag_4 (nr__5) ! (mag) ? SDSS PSF i band magnitude
real*4 e_imag_4 (nr__5) ! (mag) ? Uncertainty in imag
character*1 f_zmag_4 (nr__5) ! [*] Indicates saturated pixels in zmag
real*4 zmag_4 (nr__5) ! (mag) ? SDSS PSF z band magnitude
real*4 e_zmag_4 (nr__5) ! (mag) ? Uncertainty in zmag
character*2 Note_1 (nr__5) ! 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(i__),Plate(i__),ObsDate(i__),Fiber(i__),RAdeg(i__),
+ DEdeg(i__),f_umag(i__),umag(i__),e_umag(i__),f_gmag(i__),
+ gmag(i__),e_gmag(i__),f_rmag(i__),rmag(i__),e_rmag(i__),
+ f_imag(i__),imag(i__),e_imag(i__),f_zmag(i__),zmag(i__),
+ e_zmag(i__),S_N(i__),pmRA(i__),pmDE(i__),Ag(i__),Teff(i__),
+ e_Teff(i__),log_g(i__),e_log_g(i__),Chi2(i__),AutoID(i__),
+ Notes(i__)
if(ar__(127:130) .EQ. '') pmRA(i__) = iNULL__
if(ar__(132:135) .EQ. '') pmDE(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ SDSS(i__),Plate(i__),ObsDate(i__),Fiber(i__),RAdeg(i__),
+ DEdeg(i__),f_umag(i__),umag(i__),e_umag(i__),f_gmag(i__),
+ gmag(i__),e_gmag(i__),f_rmag(i__),rmag(i__),e_rmag(i__),
+ f_imag(i__),imag(i__),e_imag(i__),f_zmag(i__),zmag(i__),
+ e_zmag(i__),S_N(i__),pmRA(i__),pmDE(i__),Ag(i__),Teff(i__),
+ e_Teff(i__),log_g(i__),e_log_g(i__),Chi2(i__),AutoID(i__),
+ Notes(i__)
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(i__),Plate_1(i__),ObsDate_1(i__),Fiber_1(i__),
+ RAdeg_1(i__),DEdeg_1(i__),f_umag_1(i__),umag_1(i__),
+ e_umag_1(i__),f_gmag_1(i__),gmag_1(i__),e_gmag_1(i__),
+ f_rmag_1(i__),rmag_1(i__),e_rmag_1(i__),f_imag_1(i__),
+ imag_1(i__),e_imag_1(i__),f_zmag_1(i__),zmag_1(i__),
+ e_zmag_1(i__),S_N_1(i__),pmRA_1(i__),pmDE_1(i__),Ag_1(i__),
+ Teff_1(i__),e_Teff_1(i__),log_g_1(i__),e_log_g_1(i__),
+ Chi2_1(i__),AutoID_1(i__),Notes_1(i__)
if(ar__1(127:130) .EQ. '') pmRA_1(i__) = iNULL__
if(ar__1(132:135) .EQ. '') pmDE_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ SDSS_1(i__),Plate_1(i__),ObsDate_1(i__),Fiber_1(i__),
+ RAdeg_1(i__),DEdeg_1(i__),f_umag_1(i__),umag_1(i__),
+ e_umag_1(i__),f_gmag_1(i__),gmag_1(i__),e_gmag_1(i__),
+ f_rmag_1(i__),rmag_1(i__),e_rmag_1(i__),f_imag_1(i__),
+ imag_1(i__),e_imag_1(i__),f_zmag_1(i__),zmag_1(i__),
+ e_zmag_1(i__),S_N_1(i__),pmRA_1(i__),pmDE_1(i__),Ag_1(i__),
+ Teff_1(i__),e_Teff_1(i__),log_g_1(i__),e_log_g_1(i__),
+ Chi2_1(i__),AutoID_1(i__),Notes_1(i__)
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(i__),Plate_2(i__),ObsDate_2(i__),Fiber_2(i__),
+ RAdeg_2(i__),DEdeg_2(i__),Epoch(i__),f_umag_2(i__),
+ umag_2(i__),e_umag_2(i__),f_gmag_2(i__),gmag_2(i__),
+ e_gmag_2(i__),f_rmag_2(i__),rmag_2(i__),e_rmag_2(i__),
+ f_imag_2(i__),imag_2(i__),e_imag_2(i__),f_zmag_2(i__),
+ zmag_2(i__),e_zmag_2(i__),S_N_2(i__),pmRA_2(i__),pmDE_2(i__),
+ Ag_2(i__),ID(i__),Notes_2(i__)
if(ar__2(135:138) .EQ. '') pmRA_2(i__) = iNULL__
if(ar__2(140:143) .EQ. '') pmDE_2(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ SDSS_2(i__),Plate_2(i__),ObsDate_2(i__),Fiber_2(i__),
+ RAdeg_2(i__),DEdeg_2(i__),Epoch(i__),f_umag_2(i__),
+ umag_2(i__),e_umag_2(i__),f_gmag_2(i__),gmag_2(i__),
+ e_gmag_2(i__),f_rmag_2(i__),rmag_2(i__),e_rmag_2(i__),
+ f_imag_2(i__),imag_2(i__),e_imag_2(i__),f_zmag_2(i__),
+ zmag_2(i__),e_zmag_2(i__),S_N_2(i__),pmRA_2(i__),pmDE_2(i__),
+ Ag_2(i__),ID(i__),Notes_2(i__)
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(i__),Plate_3(i__),ObsDate_3(i__),Fiber_3(i__),
+ RAdeg_3(i__),DEdeg_3(i__),Epoch_1(i__),f_umag_3(i__),
+ umag_3(i__),e_umag_3(i__),f_gmag_3(i__),gmag_3(i__),
+ e_gmag_3(i__),f_rmag_3(i__),rmag_3(i__),e_rmag_3(i__),
+ f_imag_3(i__),imag_3(i__),e_imag_3(i__),f_zmag_3(i__),
+ zmag_3(i__),e_zmag_3(i__),S_N_3(i__),pmRA_3(i__),pmDE_3(i__),
+ Ag_3(i__),ID_1(i__),Notes_3(i__)
if(ar__3(135:138) .EQ. '') pmRA_3(i__) = iNULL__
if(ar__3(140:143) .EQ. '') pmDE_3(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ SDSS_3(i__),Plate_3(i__),ObsDate_3(i__),Fiber_3(i__),
+ RAdeg_3(i__),DEdeg_3(i__),Epoch_1(i__),f_umag_3(i__),
+ umag_3(i__),e_umag_3(i__),f_gmag_3(i__),gmag_3(i__),
+ e_gmag_3(i__),f_rmag_3(i__),rmag_3(i__),e_rmag_3(i__),
+ f_imag_3(i__),imag_3(i__),e_imag_3(i__),f_zmag_3(i__),
+ zmag_3(i__),e_zmag_3(i__),S_N_3(i__),pmRA_3(i__),pmDE_3(i__),
+ Ag_3(i__),ID_1(i__),Notes_3(i__)
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(i__),gmag_4(i__),SType(i__),OName(i__),OType(i__),
+ Note(i__)
c ..............Just test output...........
write(6,5)
+ SDSS_4(i__),gmag_4(i__),SType(i__),OName(i__),OType(i__),
+ Note(i__)
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(i__),Type(i__),mag(i__),Filt(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),f_umag_4(i__),
+ umag_4(i__),e_umag_4(i__),f_gmag_4(i__),gmag_5(i__),
+ e_gmag_4(i__),f_rmag_4(i__),rmag_4(i__),e_rmag_4(i__),
+ f_imag_4(i__),imag_4(i__),e_imag_4(i__),f_zmag_4(i__),
+ zmag_4(i__),e_zmag_4(i__),Note_1(i__)
if(ar__5(21:25) .EQ. '') mag(i__) = rNULL__
if(ar__5(29:30) .EQ. '') RAh(i__) = iNULL__
if(ar__5(32:33) .EQ. '') RAm(i__) = iNULL__
if(ar__5(35:39) .EQ. '') RAs(i__) = rNULL__
if(ar__5(42:43) .EQ. '') DEd(i__) = iNULL__
if(ar__5(45:46) .EQ. '') DEm(i__) = iNULL__
if(ar__5(48:51) .EQ. '') DEs(i__) = rNULL__
if(ar__5(54:58) .EQ. '') umag_4(i__) = rNULL__
if(ar__5(60:63) .EQ. '') e_umag_4(i__) = rNULL__
if(ar__5(66:70) .EQ. '') gmag_5(i__) = rNULL__
if(ar__5(72:75) .EQ. '') e_gmag_4(i__) = rNULL__
if(ar__5(78:82) .EQ. '') rmag_4(i__) = rNULL__
if(ar__5(84:87) .EQ. '') e_rmag_4(i__) = rNULL__
if(ar__5(90:94) .EQ. '') imag_4(i__) = rNULL__
if(ar__5(96:99) .EQ. '') e_imag_4(i__) = rNULL__
if(ar__5(102:106) .EQ. '') zmag_4(i__) = rNULL__
if(ar__5(108:111) .EQ. '') e_zmag_4(i__) = rNULL__
RAdeg_4(i__) = rNULL__
DEdeg_4(i__) = 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(i__) .GT. -180) RAdeg_4(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg_4(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_4(i__)
c ..............Just test output...........
write(6,6)
+ WD(i__),Type(i__),mag(i__),Filt(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),f_umag_4(i__),
+ umag_4(i__),e_umag_4(i__),f_gmag_4(i__),gmag_5(i__),
+ e_gmag_4(i__),f_rmag_4(i__),rmag_4(i__),e_rmag_4(i__),
+ f_imag_4(i__),imag_4(i__),e_imag_4(i__),f_zmag_4(i__),
+ zmag_4(i__),e_zmag_4(i__),Note_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_4(i__),DEdeg_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end