Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/216/4 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-Apr-21
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/ApJS/216/4 SDSS-RM project: technical overview (Shen+, 2015)
*================================================================================
*The Sloan Digital Sky Survey Reverberation Mapping project: technical overview.
* Shen Y., Brandt W.N., Dawson K.S., Hall P.B., McGreer I.D., Anderson S.F.,
* Chen Y., Denney K.D., Eftekharzadeh S., Fan X., Gao Y., Green P.J.,
* Greene J.E., Ho L.C., Horne K., Jiang L., Kelly B.C., Kinemuchi K.,
* Kochanek C.S., Paris I., Peters C.M., Peterson B.M., Petitjean P.,
* Ponder K., Richards G.T., Schneider D.P., Seth A., Smith R.N.,
* Strauss M.A., Tao C., Trump J.R., Wood-Vasey W.M., Zu Y., Eisenstein D.J.,
* Pan K., Bizyaev D., Malanushenko V., Malanushenko E., Oravetz D.
* <Astrophys. J. Suppl. Ser., 216, 4 (2015)>
* =2015ApJS..216....4S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! The tiled sample
integer*4 nr__
parameter (nr__=1000) ! Number of records
character*758 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__) ! [0/999] Index of objects in this catalog (RMID)
real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000) (1)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000) (1)
real*4 z (nr__) ! [0.1/4.4]?=0 Redshift (2)
character*8 Type (nr__) ! Source type (3)
real*4 umag (nr__) ! (mag) [16.8/26]?=0 u-band PSF magnitude (asinh) in
* SDSS filter; undereddened
real*4 gmag (nr__) ! (mag) [16.2/23.7]?=0 g-band PSF magnitude (asinh) in
* SDSS filter; undereddened
real*4 rmag (nr__) ! (mag) [16/22.8]?=0 r-band PSF magnitude (asinh) in
* SDSS filter; undereddened
real*4 imag (nr__) ! (mag) [15.8/21.7]?=0 i-band PSF magnitude (asinh) in
* SDSS filter; undereddened
real*4 zmag (nr__) ! (mag) [15.8/22.9]?=0 z-band PSF magnitude (asinh) in
* SDSS filter; undereddened
integer*4 OT (nr__) ! [0/6] Object type (0=unknown, 3=extended source or
* 6=point source)
character*4 Sample (nr__) ! Sample ID (4)
integer*4 Plate (32,nr__) ! [7338/7340] SDSS plate(s) for 32 epochs (5)
integer*4 Fiber (32,nr__) ! [1/1000] SDSS fiber(s) for 32 epochs (5)
integer*4 MJD (32,nr__) ! [56660/56837] SDSS Modified Julian Date(s) for
* 28 epochs (5)
real*4 S_N (32,nr__) ! [-1.5/81.7] Median S/N/pixel of each epoch
* (over full wavelength range)
*Note (1): Object positions are from the corrected SDSS astrometry discussed in
* DR10 (Ahn et al. 2014ApJS..211...17A).
*Note (2): Redshifts are measured by running the SDSS redshift pipeline (Bolton
* et al. 2012AJ....144..144B) on the coadded spectra from all epochs
* and visually verified.
*Note (3): Types are:
* RM_TILE1 = RM targets tile 1 (230 QSOs)
* RM_TILE2 = RM targets tile 2 (619 QSOs)
* LRG = luminous red galaxy target (1 source)
* STD = Spectrophotometric standard F-stars (70 sources)
* SKY = Sky fibers (80 sources)
*Note (4): Sample are:
* boss = BOSS DR9 quasar sample (Paris et al. 2012, VII/269);
* dr7 = SDSS Data Release 7 (DR7) quasars (Shen et al. 2011, J/ApJS/194/45);
* mmt = Follow-up spectroscopy with MMT/Hectospec
* (P. J. Green et al., in preparation)
* egs = quasars in the EGS field discovered by the DEEP2 survey
* (Newman et al. 2013ApJS..208....5N)
*Note (5): The Plate, Fiber and MJD columns record the spectroscopic
* identifications for all epochs of each target.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Log of BOSS RM observations
integer*4 nr__1
parameter (nr__1=32) ! Number of records
character*56 ar__1 ! Full-size record
character*7 Date (nr__1) ! UT months of observation (in year 2014)
integer*4 Epoch (nr__1) ! [1/32] Epoch number (1)
character*1 f_Epoch (nr__1) ! [ab] low quality spectra (2)
integer*4 MJDs (nr__1) ! (d) [56659/56837] Start Modified Julian Date
integer*4 MJDe (nr__1) ! (d) [56660/56837] End Modified Julian Date
real*8 v_MJD_ (nr__1) ! (d) Mean Modified Julian Date
character*10 File (nr__1) ! Plate-MJDe file identifier
real*4 S_N_1 (nr__1) ! [0.5/34.2] Plate S/N
integer*4 N (nr__1) ! [2/11] Number of coadded 15 minute
* sub-exposures in each epoch
real*4 Seeing (nr__1) ! (arcsec) [1.3/3] Median seeing during the observation
* estimated from the guide camera images
*Note (1): Except for Epochs 1 and 6, all epochs were observed on the same night.
*Note (2): Flag as follows:
* a = Epochs that have a plate S/N<5 due to poor observing conditions. Most of
* the targets on these plates should be discarded in the cross-correlation
* analysis due to the poor spectral quality.
* b = Epoch 56825 suffered from a software error with no simultaneous
* calibration exposures taken for science exposures. Afternoon calibration
* exposures were used in reduction.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sp.dat' ! List of spectra for the 32 epochs by object
integer*4 nr__2
parameter (nr__2=32000) ! Number of records
character*24 ar__2 ! Full-size record
integer*4 Seq_1 (nr__2) ! [0/999] Index of objects in this catalog (RMID)
integer*4 Plate_1 (nr__2) ! SDSS plate
integer*4 Fiber_1 (nr__2) ! SDSS fiber
integer*4 MJD_1 (nr__2) ! SDSS Modified Julian Date
real*4 S_N_2 (nr__2) ! Median S/N per pixel (over full wavelength range)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Mask of sky pixels
integer*4 nr__3
parameter (nr__3=589) ! Number of records
character*16 ar__3 ! Full-size record
real*8 lambda (nr__3) ! (0.1nm) [5458.8/10296.8] Wavelength in angstroms (1)
real*4 loglambda (nr__3) ! ([0.1nm]) [3.7371/4.0127] Logarithmic wavelength (1)
*Note (1): Central wavelengths of sky pixels within [5000,10300]{AA},
* generated using the sky spectra from the SDSS-RM program and
* following the approach in Wild & Hewett (2005MNRAS.358.1083W).
C=============================================================================
C Loading file 'table1.dat' ! The tiled sample
C Format for file interpretation
1 format(
+ I3,1X,F10.6,1X,F9.6,1X,F6.4,1X,A8,1X,F6.3,1X,F6.3,1X,F6.3,1X,
+ F6.3,1X,F6.3,1X,I1,1X,A4,1X,32I5,1X,32I5,1X,32I6,1X,32F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,1000
read(1,'(A758)')ar__
read(ar__,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),z(i__),Type(i__),umag(i__),
+ gmag(i__),rmag(i__),imag(i__),zmag(i__),OT(i__),Sample(i__),
+ (Plate(j__,i__),j__=1,32),(Fiber(j__,i__),j__=1,32),(MJD(j__,
+ i__),j__=1,32),(S_N(j__,i__),j__=1,32)
c ..............Just test output...........
write(6,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),z(i__),Type(i__),umag(i__),
+ gmag(i__),rmag(i__),imag(i__),zmag(i__),OT(i__),Sample(i__),
+ (Plate(j__,i__),j__=1,32),(Fiber(j__,i__),j__=1,32),(MJD(j__,
+ i__),j__=1,32),(S_N(j__,i__),j__=1,32)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Log of BOSS RM observations
C Format for file interpretation
2 format(
+ A7,1X,I2,A1,1X,I5,1X,I5,1X,F9.3,1X,A10,1X,F4.1,1X,I2,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,32
read(1,'(A56)')ar__1
read(ar__1,2)
+ Date(i__),Epoch(i__),f_Epoch(i__),MJDs(i__),MJDe(i__),
+ v_MJD_(i__),File(i__),S_N_1(i__),N(i__),Seeing(i__)
c ..............Just test output...........
write(6,2)
+ Date(i__),Epoch(i__),f_Epoch(i__),MJDs(i__),MJDe(i__),
+ v_MJD_(i__),File(i__),S_N_1(i__),N(i__),Seeing(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'sp.dat' ! List of spectra for the 32 epochs by object
C Format for file interpretation
3 format(I3,1X,I4,1X,I4,1X,I5,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'sp.dat')
write(6,*) '....Loading file: sp.dat'
do i__=1,32000
read(1,'(A24)')ar__2
read(ar__2,3)
+ Seq_1(i__),Plate_1(i__),Fiber_1(i__),MJD_1(i__),S_N_2(i__)
c ..............Just test output...........
write(6,3)
+ Seq_1(i__),Plate_1(i__),Fiber_1(i__),MJD_1(i__),S_N_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Mask of sky pixels
C Format for file interpretation
4 format(F9.3,1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,589
read(1,'(A16)')ar__3
read(ar__3,4)lambda(i__),loglambda(i__)
c ..............Just test output...........
write(6,4)lambda(i__),loglambda(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end