Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/402/620 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-Oct-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/MNRAS/402/620 SDSS WD main-sequence binaries (Rebassa-Mansergas+, 2010)
*================================================================================
*Post-common envelope binaries from SDSS.
*VII. A catalogue of white dwarf-main sequence binaries.
* Rebassa-Mansergas A., Gansicke B.T., Schreiber M.R., Koester D.,
* Rodriguez-Gil P.
* <Mon. Not. R. Astron. Soc., 402, 620-640 (2010)>
* =2010MNRAS.402..620R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Updated classification of the 204 objects from
Silvestri et al. (2007, Cat <J/AJ/134/741>)
which are not considered as WDMS binaries by us
integer*4 nr__
parameter (nr__=204) ! Number of records
character*66 ar__ ! Full-size record
character*19 SDSS (nr__) ! SDSS designation (JHHMMSS.ss+DDMMSS.s)
character*42 Com (nr__) ! Comments
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Updated classification of the 81 objects from
Heller et al. (2009, Cat. <J/A+A/496/191>)
which are not considered as WDMS binaries by us
integer*4 nr__1
parameter (nr__1=81) ! Number of records
character*66 ar__1 ! Full-size record
character*19 SDSS_1 (nr__1) ! SDSS designation (JHHMMSS.ss+DDMMSS.s)
character*42 Com_1 (nr__1) ! Comments
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! The complete catalogue. Coordinates and
GALEX-SDSS-UKIDSS magnitudes for the 1602 WDMS
binaries and candidates are also included
integer*4 nr__2
parameter (nr__2=1602) ! Number of records
character*242 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*19 SDSS_2 (nr__2) ! SDSS designation (JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg (nr__2) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__2) ! (deg) Declination in decimal degrees (J2000)
real*8 NUV (nr__2) ! (mag) ?=0 GALEX NUV magnitude (177-300nm)
real*8 e_NUV (nr__2) ! (mag) ?=0 rms uncertainty on NUV
real*8 FUV (nr__2) ! (mag) ?=0 GALEX FUV magnitude (135-175nm)
real*8 e_FUV (nr__2) ! (mag) ?=0 rms uncertainty on FUV
real*8 umag (nr__2) ! (mag) ?=0 SDSS u magnitude
real*8 e_umag (nr__2) ! (mag) ?=0 rms uncertainty on umag
real*8 gmag (nr__2) ! (mag) ?=0 SDSS g magnitude
real*8 e_gmag (nr__2) ! (mag) ?=0 rms uncertainty on gmag
real*8 rmag (nr__2) ! (mag) ?=0 SDSS r magnitude
real*8 e_rmag (nr__2) ! (mag) ?=0 rms uncertainty on rmag
real*8 imag (nr__2) ! (mag) ?=0 SDSS i magnitude
real*8 e_imag (nr__2) ! (mag) ?=0 rms uncertainty on imag
real*8 zmag (nr__2) ! (mag) ?=0 SDSS z magnitude (1)
real*8 e_zmag (nr__2) ! (mag) ?=0 rms uncertainty on zmag
real*8 Ymag (nr__2) ! (mag) ?=0 UKIDSS Y magnitude (1)
real*8 e_Ymag (nr__2) ! (mag) ?=0 rms Yncertainty on Ymag
real*8 Jmag (nr__2) ! (mag) ?=0 UKIDSS J magnitude (1)
real*8 e_Jmag (nr__2) ! (mag) ?=0 rms uncertainty on Jmag
real*8 Hmag (nr__2) ! (mag) ?=0 UKIDSS H magnitude (1)
real*8 e_Hmag (nr__2) ! (mag) ?=0 rms uncertainty on Hmag
real*8 Kmag (nr__2) ! (mag) ?=0 UKIDSS K magnitude (1)
real*8 e_Kmag (nr__2) ! (mag) ?=0 rms uncertainty on Kmag
character*6 Type (nr__2) ! Object type
*Note (1): UKIDSS filters 50% cut-off wavelengths ({mu}m):
* -------------------
* Z 0.83 0.925
* Y 0.97 1.07
* J 1.17 1.33
* H 1.49 1.78
* K 2.03 2.37
* -------------------
* (details at http://www.ukidss.org/technical/instrument/filters.html)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! *White-dwarf masses, effective temperatures,
surface gravities, spectral types and distances
of the 1602 WDMS binaries in our catalogue, as
determined from spectral modelling
integer*4 nr__3
parameter (nr__3=3592) ! Number of records
character*118 ar__3 ! Full-size record
integer*4 Sol (nr__3) ! [0/1]? Adopted (1) or ruled-out (0) solution
character*19 SDSS_3 (nr__3) ! SDSS designation (JHHMMSS.ss+DDMMSS.s)
character*6 Type_1 (nr__3) ! Object type
integer*4 MJD (nr__3) ! SDSS MJD identifier
integer*4 Plate (nr__3) ! SDSS Plate identifier
integer*4 Fiber (nr__3) ! SDSS Fiber identifier
integer*4 Teff (nr__3) ! (K) ?=0 Effective temperature
integer*4 e_Teff (nr__3) ! (K) ?=0 rms uncertainty on Teff
real*4 logg (nr__3) ! ([cm/s2]) ?=0 Surface gravity
real*4 e_logg (nr__3) ! ([cm/s2]) ?=0 rms uncertainty on logg
real*4 Mwd (nr__3) ! (solMass) ?=0.00 Mass of the white dwarf
real*4 e_Mwd (nr__3) ! (solMass) ?=0.00 rms uncertainty on Mwd
integer*4 dwd (nr__3) ! (pc) ?=0 Distance of the white dwarf
integer*4 e_dwd (nr__3) ! (pc) ?=0 rms uncertainty on dwd
integer*4 MType (nr__3) ! ?=-1 Sub-type of the M dwarf
integer*4 dsec (nr__3) ! (pc) ?=0 Distance of the secondary
integer*4 e_dsec (nr__3) ! (pc) ?=0 rms uncertainty on dsec
character*9 Flag (nr__3) ! Flag(s) for previously studied systems and
* resolved WDMS binary pairs (2)
character*4 Notes (nr__3) ! ? Individual notes (3)
*Note (2): Flags as follows:
* e = Eisenstein et al. (2006, Cat. <J/ApJS/167/40>)
* s = Silvestri et al. (2007, Cat <J/AJ/134/741>)
* a = Augusteijn et al. (2008, Cat <J/A+A/486/843>)
* h = Heller et al. (2009, Cat. <J/A+A/496/191>)
* re = resolved WDMS binary pairs
*Note (3): Notes as follows:
* 1 = Nearby Galaxy
* 2 = GALEX magnitudes changes the solution obtained from the fit
* 3 = No SDSS image available
* 4 = Star/s nearby
* 5 = PCEB
* 6 = One of the spectra comes from a SEGUE plate
* 7 = Irradiated secondary star
* 8 = Broken spectrum
* 9 = Star superimposed in the SDSS image,triple?
* 10 = Wrong flux calibration in at least one spectrum
* 11 = Cataclysmic Variable?
* 12 = Wide binary?
* 13 = Peculiar secondary star
* 14 = Detached magnetic Cataclysmic Variable
* 15 = Bad SDSS image
* 16 = Low accretion rate polar
* 17 = PG1159
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Radial velocities measured from the NaI
{lambda}{lambda}8183.27,8194.81 doublet and the
H{alpha} emission for 1068 systems in our catalogue
integer*4 nr__4
parameter (nr__4=1225) ! Number of records
character*68 ar__4 ! Full-size record
character*19 SDSS_4 (nr__4) ! SDSS designation (JHHMMSS.ss+DDMMSS.s)
integer*4 m_SDSS (nr__4) ! [1/5]? Multiplicity index on SDSS
real*8 HJD (nr__4) ! (d) Heliocentric Julian date
real*8 RVNa (nr__4) ! (km/s) ?=0 Radial velocities measured from the
* NaI {lambda}{lambda}8183.27,8194.81 doublet
real*4 e_RVNa (nr__4) ! (km/s) ?=0 rms uncertainty on RVNa
real*8 RVHa (nr__4) ! (km/s) ?=0 Radial velocities measured from the H{alpha}
* emission
real*4 e_RVHa (nr__4) ! (km/s) ?=0 rms uncertainty on RVHa
character*1 Rel (nr__4) ! [y/n] Radial velocity values obtained from
* spectra that are (y) and that are not (n)
* combined from individual exposures taken on
* different nights
C=============================================================================
C Loading file 'table2.dat' ! Updated classification of the 204 objects from
* Silvestri et al. (2007, Cat <J/AJ/134/741>)
* which are not considered as WDMS binaries by us
C Format for file interpretation
1 format(4X,A19,1X,A42)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,204
read(1,'(A66)')ar__
read(ar__,1)SDSS(i__),Com(i__)
c ..............Just test output...........
write(6,1)SDSS(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Updated classification of the 81 objects from
* Heller et al. (2009, Cat. <J/A+A/496/191>)
* which are not considered as WDMS binaries by us
C Format for file interpretation
2 format(4X,A19,1X,A42)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,81
read(1,'(A66)')ar__1
read(ar__1,2)SDSS_1(i__),Com_1(i__)
c ..............Just test output...........
write(6,2)SDSS_1(i__),Com_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! The complete catalogue. Coordinates and
* GALEX-SDSS-UKIDSS magnitudes for the 1602 WDMS
* binaries and candidates are also included
C Format for file interpretation
3 format(
+ 4X,A19,1X,F10.6,1X,F10.6,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,
+ F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F8.5,
+ 1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.6,1X,F8.5,1X,F7.5,1X,F8.5,1X,
+ F7.5,1X,F8.5,1X,F7.5,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,1602
read(1,'(A242)')ar__2
read(ar__2,3)
+ SDSS_2(i__),RAdeg(i__),DEdeg(i__),NUV(i__),e_NUV(i__),
+ FUV(i__),e_FUV(i__),umag(i__),e_umag(i__),gmag(i__),
+ e_gmag(i__),rmag(i__),e_rmag(i__),imag(i__),e_imag(i__),
+ zmag(i__),e_zmag(i__),Ymag(i__),e_Ymag(i__),Jmag(i__),
+ e_Jmag(i__),Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),
+ Type(i__)
c ..............Just test output...........
write(6,3)
+ SDSS_2(i__),RAdeg(i__),DEdeg(i__),NUV(i__),e_NUV(i__),
+ FUV(i__),e_FUV(i__),umag(i__),e_umag(i__),gmag(i__),
+ e_gmag(i__),rmag(i__),e_rmag(i__),imag(i__),e_imag(i__),
+ zmag(i__),e_zmag(i__),Ymag(i__),e_Ymag(i__),Jmag(i__),
+ e_Jmag(i__),Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),
+ Type(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! *White-dwarf masses, effective temperatures,
* surface gravities, spectral types and distances
* of the 1602 WDMS binaries in our catalogue, as
* determined from spectral modelling
C Format for file interpretation
4 format(
+ I1,5X,A19,1X,A6,1X,I5,1X,I4,1X,I3,1X,I5,1X,I5,1X,F4.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,I4,1X,I4,1X,I2,1X,I4,1X,I4,1X,A9,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,3592
read(1,'(A118)')ar__3
read(ar__3,4)
+ Sol(i__),SDSS_3(i__),Type_1(i__),MJD(i__),Plate(i__),
+ Fiber(i__),Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),
+ Mwd(i__),e_Mwd(i__),dwd(i__),e_dwd(i__),MType(i__),dsec(i__),
+ e_dsec(i__),Flag(i__),Notes(i__)
if(ar__3(1:1) .EQ. '') Sol(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ Sol(i__),SDSS_3(i__),Type_1(i__),MJD(i__),Plate(i__),
+ Fiber(i__),Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),
+ Mwd(i__),e_Mwd(i__),dwd(i__),e_dwd(i__),MType(i__),dsec(i__),
+ e_dsec(i__),Flag(i__),Notes(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Radial velocities measured from the NaI
* {lambda}{lambda}8183.27,8194.81 doublet and the
* H{alpha} emission for 1068 systems in our catalogue
C Format for file interpretation
5 format(
+ 4X,A19,1X,I1,1X,F12.4,1X,F7.2,1X,F5.2,1X,F7.2,1X,F5.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,1225
read(1,'(A68)')ar__4
read(ar__4,5)
+ SDSS_4(i__),m_SDSS(i__),HJD(i__),RVNa(i__),e_RVNa(i__),
+ RVHa(i__),e_RVHa(i__),Rel(i__)
if(ar__4(25:25) .EQ. '') m_SDSS(i__) = iNULL__
c ..............Just test output...........
write(6,5)
+ SDSS_4(i__),m_SDSS(i__),HJD(i__),RVNa(i__),e_RVNa(i__),
+ RVHa(i__),e_RVHa(i__),Rel(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end