Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/666/806 into FORTRAN code for reading data files line by line.
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.5, on 2013-Jun-19
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJ/666/806 SED of Spitzer quasars (QUEST) (Netzer+, 2007)
*================================================================================
*Spitzer Quasar and ULIRG Evolution Study (QUEST).
*II. The spectral energy distributions of Palomar-Green quasars.
* Netzer H., Lutz D., Schweitzer M., Contursi A., Sturm E., Tacconi L.J.,
* Veilleux S., Kim D.-C., Rupke D., Baker A.J., Dasyra K., Mazzarella J.,
* Lord S.
* <Astrophys. J., 666, 806-816 (2007)>
* =2007ApJ...666..806N
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! QSO sample and supplementary data
integer*4 nr__
parameter (nr__=29) ! Number of records
character*77 ar__ ! Full-size record
character*11 Name ! QSO name
character*8 OName ! Other name
real*4 z ! Redshift
real*4 log_L5100 ! ([10-7W]) Continuum luminosity {lambda}L{lambda} at
* 5100{AA} rest wavelength (1)
real*4 SJ ! (mJy) Near-infrared flux in J band
real*4 SH ! (mJy) Near-infrared flux in H band
real*4 SK ! (mJy) Near-infrared flux in K band
real*4 SL ! (mJy) ?=- Near-infrared flux in L band
character*7 Ref ! Near-infrared fluxes references (2)
character*10 Class ! Classification as far-infrared strong/weak
* (relative to the mid-infrared), or undetected
*Note (1): from spectra by T. Boroson 2006, private communication.
*Note (2): References as follows:
* 2M = 2MASS magnitudes, extended source catalog (Jarrett et al.
* 2000AJ....119.2498J, cat. VII/233); K20 isophotal magnitudes
* for slightly extended sources, point source catalog (II/246)
* magnitudes otherwise
* H82 = Hyland & Allen (1982MNRAS.199..943H)
* M83 = McAlary et al. (1983ApJS...52..341M)
* N87 = Neugebauer et al. (1987ApJS...63..615N)
* S89 = Sanders et al. (1989ApJ...347...29S).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Average QSO SEDs
integer*4 nr__1
parameter (nr__1=354) ! Number of records
character*18 ar__1 ! Full-size record
real*4 lambda ! (um) Rest wavelength
real*4 ObsF ! ? Observed QSO SED, {lambda}F_{lambda}_ (1)
real*4 IntF ! Intrinsic AGN SED for FIR-weak QSOs,
* {lambda}F_{lambda}_ (2)
*Note (1): Average observed SED of 28 PG QSOs (see Fig. 3, top, smoothed),
* in arbitrary units
*Note (2): Average intrinsic AGN SED of 8 'FIR-weak' PG QSOs, after
* subtraction of a starburst component (see Fig. 6, right, smoothed),
* in arbitrary units
C=============================================================================
C Loading file 'table1.dat' ! QSO sample and supplementary data
C Format for file interpretation
1 format(
+ A11,1X,A8,1X,F6.4,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F6.2,1X,
+ A7,1X,A10)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,29
read(1,'(A77)')ar__
read(ar__,1)Name,OName,z,log_L5100,SJ,SH,SK,SL,Ref,Class
if (SL .EQ. 2.2233e-322) SL = rNULL__
c ..............Just test output...........
write(6,1)Name,OName,z,log_L5100,SJ,SH,SK,SL,Ref,Class
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Average QSO SEDs
C Format for file interpretation
2 format(F6.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,354
read(1,'(A18)')ar__1
read(ar__1,2)lambda,ObsF,IntF
if(ar__1(8:12) .EQ. '') ObsF = rNULL__
c ..............Just test output...........
write(6,2)lambda,ObsF,IntF
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end