Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/702/441 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-Apr-20
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/702/441 Obscured active galaxies (Greene+, 2009)
*================================================================================
*The growth of black holes: insights from obscured active galaxies.
* Greene J.E., Zakamska N.L., Liu X., Barth A.J., Ho L.C.
* <Astrophys. J., 702, 441-459 (2009)>
* =2009ApJ...702..441G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Sample and observations
integer*4 nr__
parameter (nr__=15) ! Number of records
character*70 ar__ ! Full-size record
character*19 SDSS (nr__) ! SDSS identification (JHHMMSS.ss+DDMMSS.s)
real*4 z (nr__) ! Redshift
real*4 LOIII (nr__) ! ([10-7W]) L_[OIII]_ (in erg/s) (1)
character*10 Date (nr__) ! ("YYYY/MM/DD") Date of LDSS-3 observation
integer*4 Ts (nr__) ! (s) Total duration of spectroscopic exposure
integer*4 N (nr__) ! Number of slit positions
integer*4 Ti (nr__) ! (s) ? Total integrated time
integer*4 PA1 (nr__) ! (deg) Position angle of slit 1 position (2)
integer*4 PA2 (nr__) ! (deg) ? Position angle of slit 2 position (2)
integer*4 PA3 (nr__) ! (deg) ? Position angle of slit 3 position (2)
*Note (1): L_[OIII]_ as measured from the SDSS spectrum of Reyes et al.
* 2008, Cat. J/AJ/136/2373; used in the object selection.
*Note (2): measured E from N.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Galaxy properties
integer*4 nr__1
parameter (nr__1=15) ! Number of records
character*80 ar__1 ! Full-size record
character*19 SDSS_1 (nr__1) ! SDSS identification (JHHMMSS.ss+DDMMSS.s)
character*1 F (nr__1) ! Filter used for structural measurements
integer*4 n_1 (nr__1) ! ? best-fit Sersic index for the bulge-like
* component from the GALFIT model
real*4 rMAG (nr__1) ! (mag) Absolute magnitude of the bulge component (3)
real*4 e_rMAG (nr__1) ! (mag) rMAG uncertainty
real*4 k (nr__1) ! (mag) k-correction used to correct the absolute
* magnitudes in F
real*4 Af (nr__1) ! (mag) Galactic reddening correction in F (4)
real*4 OIIIm (nr__1) ! (mag) Magnitude of the [OIII] line (5)
real*4 re (nr__1) ! (arcsec) ? Effective radius of the bulge-like
* component of our best-fit model
real*4 e_re (nr__1) ! (arcsec) ? re uncertainty
real*4 S (nr__1) ! (kpc/arcsec) Scale-factor to convert re to a physical size
real*4 B_T (nr__1) ! ? Bulge-to-Total ratio for objects with
* sensible GALFIT decompositions (6)
integer*4 Sig (nr__1) ! (km/s) ? Stellar velocity dispersion within 2.25"
integer*4 e_Sig (nr__1) ! (km/s) ? Sig uncertainty
*Note (3): Absolute magnitude as derived from our best-fit GALFIT model,
* with an [OIII] correction (Column OIIIm), correction for Galactic
* reddening (Column Af) and a k-correction (Column k) derived using the
* code of Blanton & Roweis 2007AJ....133..734B. Cases without sensible
* multi-component models list total Petrosian magnitudes and are not
* assigned a Sersic index, effective radius, or B/T.
*Note (4): Galactic reddening correction in the filter shown in Column F.
* Derived from the dust maps of Schlegel et al. 1998 and extinction law
* of Cardelli et al. 1989.
*Note (5): Magnitude of the [OIII] line used to correct the total r-band
* (z<~0.38) or i-band magnitude.
*Note (6): The B/T for SDSS J125334.49-034158.1 is very low because there is
* also a bar component. If included as bulge light, the B/T increases to
* 0.05.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Additional SDSS sample (version corrected by the
author)
integer*4 nr__2
parameter (nr__2=111) ! Number of records
character*83 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*18 SDSS_2 (nr__2) ! Official SDSS name (HHMMSS.ss+DDMMSS.s; J2000)
real*8 RAdeg (nr__2) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__2) ! (deg) Declination in decimal degrees (J2000)
integer*4 Plate (nr__2) ! SDSS plate
integer*4 Fiber (nr__2) ! SDSS fiber
integer*4 MJD (nr__2) ! SDSS Modified Julian Date
real*4 z_1 (nr__2) ! Redshift as measured by the SDSS pipeline
real*4 LOIII_1 (nr__2) ! ([10-7W]) Log [OIII] luminosity in erg/s
integer*4 Sig_1 (nr__2) ! (km/s) Stellar velocity dispersion measured from
* the SDSS spectra
integer*4 e_Sig_1 (nr__2) ! (km/s) Uncertainty in Sig
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Line emission and blue light
integer*4 nr__3
parameter (nr__3=15) ! Number of records
character*79 ar__3 ! Full-size record
character*19 SDSS_3 (nr__3) ! SDSS identification (JHHMMSS.ss+DDMMSS.s)
real*4 logLOIII (nr__3) ! ([10-7W]) Luminosity in erg/s in the [OIII] line (7)
integer*4 FWHM (nr__3) ! (km/s) The [OIII] line FWHM, corrected for
* instrumental resolution (7)
real*4 OIII_Hb (nr__3) ! Ratio of [OIII] to H{beta} emission (7)
real*4 OII_OIII (nr__3) ! Ratio of [OII] to [OIII] emission from the
* SDSS spectra
real*4 logL5100d (nr__3) ! ([10-7W]) ? L5100{AA} (in erg/s) when detected in our
* three-component continuum model
real*4 logL5100H (nr__3) ! ([10-7W]) ? L5100{AA} (erg/s) inferred from broad
* H{beta} when detected
real*4 logL5100l (nr__3) ! ([10-7W]) Final adopted upper limit on contribution to
* L5100{AA} (erg/s) from scattered light (8)
character*1 n_logL5100l(nr__3) ! [*] Possible star forming populations (9)
real*4 logDn (nr__3) ! ([-]) ? The D_n_(4000) index (10)
real*4 HdA (nr__3) ! ? The H{delta}A index (11)
real*4 logM300 (nr__3) ! ([-]) ? The fraction of total galaxy mass formed in
* the last 300Myr, based on stellar population
* synthesis
*Note (7): Within the 2.25" aperture.
*Note (8): Final adopted upper limit on contribution to L_5100{AA}_ from
* scattered light, as derived from limits on a scattered broad H{beta}
* component.
*Note (9): * indicates systems where the blue light is substantially more
* than can be accounted for with scattered light. These systems may have
* a substantial star-forming population.
*Note (10): The Dn(4000) index, measured from our best-fit continuum model
* to the SDSS spectra following Balogh et al. (1999ApJ...527...54B).
*Note (11): The H{delta}A index, measured from our best-fit continuum model
* to the SDSS spectra following Worthey & Ottaviani (1997, Cat.
* J/ApJS/111/377). Positive values of H{delta}A indicate absorption.
C=============================================================================
C Loading file 'table1.dat' ! Sample and observations
C Format for file interpretation
1 format(
+ 5X,A19,1X,F5.3,1X,F5.2,1X,A10,1X,I4,1X,I1,1X,I3,1X,I3,1X,I3,
+ 1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,15
read(1,'(A70)')ar__
read(ar__,1)
+ SDSS(i__),z(i__),LOIII(i__),Date(i__),Ts(i__),N(i__),Ti(i__),
+ PA1(i__),PA2(i__),PA3(i__)
if(ar__(56:58) .EQ. '') Ti(i__) = iNULL__
if(ar__(64:66) .EQ. '') PA2(i__) = iNULL__
if(ar__(68:70) .EQ. '') PA3(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ SDSS(i__),z(i__),LOIII(i__),Date(i__),Ts(i__),N(i__),Ti(i__),
+ PA1(i__),PA2(i__),PA3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Galaxy properties
C Format for file interpretation
2 format(
+ 5X,A19,1X,A1,1X,I1,1X,F6.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.1,1X,
+ F3.1,1X,F3.1,1X,F3.1,1X,F5.3,1X,I3,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,15
read(1,'(A80)')ar__1
read(ar__1,2)
+ SDSS_1(i__),F(i__),n_1(i__),rMAG(i__),e_rMAG(i__),k(i__),
+ Af(i__),OIIIm(i__),re(i__),e_re(i__),S(i__),B_T(i__),Sig(i__),
+ e_Sig(i__)
if(ar__1(28:28) .EQ. '') n_1(i__) = iNULL__
if(ar__1(57:59) .EQ. '') re(i__) = rNULL__
if(ar__1(61:63) .EQ. '') e_re(i__) = rNULL__
if(ar__1(69:73) .EQ. '') B_T(i__) = rNULL__
if(ar__1(75:77) .EQ. '') Sig(i__) = iNULL__
if(ar__1(79:80) .EQ. '') e_Sig(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ SDSS_1(i__),F(i__),n_1(i__),rMAG(i__),e_rMAG(i__),k(i__),
+ Af(i__),OIIIm(i__),re(i__),e_re(i__),S(i__),B_T(i__),Sig(i__),
+ e_Sig(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Additional SDSS sample (version corrected by the
* author)
C Format for file interpretation
3 format(
+ A18,2X,F10.6,3X,F9.6,3X,I3,2X,I4,1X,I5,1X,F6.4,1X,F5.2,2X,I3,
+ 2X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,111
read(1,'(A83)')ar__2
read(ar__2,3)
+ SDSS_2(i__),RAdeg(i__),DEdeg(i__),Plate(i__),Fiber(i__),
+ MJD(i__),z_1(i__),LOIII_1(i__),Sig_1(i__),e_Sig_1(i__)
c ..............Just test output...........
write(6,3)
+ SDSS_2(i__),RAdeg(i__),DEdeg(i__),Plate(i__),Fiber(i__),
+ MJD(i__),z_1(i__),LOIII_1(i__),Sig_1(i__),e_Sig_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Line emission and blue light
C Format for file interpretation
4 format(
+ 5X,A19,1X,F5.2,1X,I4,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.2,1X,F4.1,
+ A1,1X,F3.1,1X,F4.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,15
read(1,'(A79)')ar__3
read(ar__3,4)
+ SDSS_3(i__),logLOIII(i__),FWHM(i__),OIII_Hb(i__),
+ OII_OIII(i__),logL5100d(i__),logL5100H(i__),logL5100l(i__),
+ n_logL5100l(i__),logDn(i__),HdA(i__),logM300(i__)
if(ar__3(49:53) .EQ. '') logL5100d(i__) = rNULL__
if(ar__3(55:59) .EQ. '') logL5100H(i__) = rNULL__
if(ar__3(67:69) .EQ. '') logDn(i__) = rNULL__
if(ar__3(71:74) .EQ. '') HdA(i__) = rNULL__
if(ar__3(76:79) .EQ. '') logM300(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ SDSS_3(i__),logLOIII(i__),FWHM(i__),OIII_Hb(i__),
+ OII_OIII(i__),logL5100d(i__),logL5100H(i__),logL5100l(i__),
+ n_logL5100l(i__),logDn(i__),HdA(i__),logM300(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end