Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/719/1672 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-Mar-28 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/719/1672 SDSS binary quasars at high redshift. I. (Hennawi+, 2010) *================================================================================ *Binary quasars at high redshift. *I. 24 new quasar pairs at z ~ 3-4. * Hennawi J.F., Myers A.D., Shen Y., Strauss M.A., Djorgovski S.G., Fan X., * Glikman E., Mahabal A., Martin C.L., Richards G.T., Schneider D.P., * Shankar F. * <Astrophys. J., 719, 1672-1692 (2010)> * =2010ApJ...719.1672H C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! High-redshift binary quasars in the SDSS imaging footprint integer*4 nr__ parameter (nr__=54) ! Number of records character*111 ar__ ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*16 v_HMS2010_ (nr__) ! Brighter quasar designation * (SDSS JHHMM+DDMMA or PSS JHHMM+DDMMA) (1) character*1 f__HMS2010_(nr__) ! [a] PSS quasar (2) integer*4 RAh (nr__) ! (h) Hour of right ascension (J2000) integer*4 RAm (nr__) ! (min) Minute of right ascension (J2000) real*4 RAs (nr__) ! (s) Second of right ascension (J2000) character*1 DE_ (nr__) ! Sign of declination (J2000) integer*4 DEd (nr__) ! (deg) Degree of declination (J2000) integer*4 DEm (nr__) ! (arcmin) Arcminute of declination (J2000) real*4 DEs (nr__) ! (arcsec) Arcsecond of declination (J2000) real*4 umag (nr__) ! (mag) ? SDSS u-band magnitude (3) real*4 gmag (nr__) ! (mag) ? SDSS g-band magnitude (3) character*1 u_rmag (nr__) ! [~] Uncertainty flag on rmag real*4 rmag (nr__) ! (mag) ? SDSS r-band magnitude (3) real*4 imag (nr__) ! (mag) ? SDSS i-band magnitude (3) real*4 zmag (nr__) ! (mag) ? SDSS z-band magnitude (3) real*4 chi2s (nr__) ! ? The minimum stellar locus distance * ({chi}^2^_star_) real*4 chi2p (nr__) ! ? Quasar locus distance ({chi}^2^_phot_) real*4 z (nr__) ! Redshift integer*4 e_z (nr__) ! (km/s) Estimation of the redshift uncertainty integer*4 DelV (nr__) ! (km/s) ? Velocity difference {Delta}V between the * two members of a pair (primary only) real*4 DelT (nr__) ! (arcsec) ? Angular separation {Delta}{theta} integer*4 RT (nr__) ! (kpc) ? Transverse proper separation *Note (1): The brightest member of the pair in the i band is always labeled * "A" and the fainter member is labeled "B."; see Note (2) for the PSS * designation. *Note (2): * a = This quasar pair was selected via a different method during a pilot * search for proto-clusters associated with z>4 quasars (see Section 4.1). * The redshifts quoted are approximate and we did not measure a velocity * difference. The fainter quasar PSS 1315+2924B is not detected in the * SDSS imaging, its quoted r-band magnitude is an estimate from Keck * LRIS R-band imaging. *Note (3): Extinction-corrected SDSS five-band point-spread function (PSF) * photometry is given in the columns umag, gmag, rmag, imag, and zmag. C============================================================================= C Loading file 'table3.dat' ! High-redshift binary quasars in the SDSS * imaging footprint C Format for file interpretation 1 format( + A16,1X,A1,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2, + 1X,F5.2,1X,A1,F5.2,1X,F5.2,1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.3,1X, + I4,1X,I4,1X,F4.1,1X,I3) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,54 read(1,'(A111)')ar__ read(ar__,1) + v_HMS2010_(i__),f__HMS2010_(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),umag(i__),gmag(i__), + u_rmag(i__),rmag(i__),imag(i__),zmag(i__),chi2s(i__), + chi2p(i__),z(i__),e_z(i__),DelV(i__),DelT(i__),RT(i__) if(ar__(46:50) .EQ. '') umag(i__) = rNULL__ if(ar__(52:56) .EQ. '') gmag(i__) = rNULL__ if(ar__(59:63) .EQ. '') rmag(i__) = rNULL__ if(ar__(65:69) .EQ. '') imag(i__) = rNULL__ if(ar__(71:75) .EQ. '') zmag(i__) = rNULL__ if(ar__(77:81) .EQ. '') chi2s(i__) = rNULL__ if(ar__(83:86) .EQ. '') chi2p(i__) = rNULL__ if(ar__(99:102) .EQ. '') DelV(i__) = iNULL__ if(ar__(104:107) .EQ. '') DelT(i__) = rNULL__ if(ar__(109:111) .EQ. '') RT(i__) = iNULL__ RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,1) + v_HMS2010_(i__),f__HMS2010_(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),umag(i__),gmag(i__), + u_rmag(i__),rmag(i__),imag(i__),zmag(i__),chi2s(i__), + chi2p(i__),z(i__),e_z(i__),DelV(i__),DelT(i__),RT(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end