Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/832/67 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-18
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/832/67 Double-peaked narrow lines in AGN. II. z<0.1 (Nevin+, 2016)
*================================================================================
*The origin of double-peaked narrow lines in active galactic nuclei.
*II. Kinematic classifications for the population at z < 0.1.
* Nevin R., Comerford J., Muller-Sanchez F., Barrows R., Cooper M.
* <Astrophys. J., 832, 67-67 (2016)>
* =2016ApJ...832...67N (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Summary of longslit observations
integer*4 nr__
parameter (nr__=71) ! Number of records
character*64 ar__ ! Full-size record
character*19 SDSS (nr__) ! SDSS ID (JHHMMSS.ss+DDMMSS.s) (1)
character*10 Name (nr__) ! Galaxy Name (JHHMM+DDMM)
real*8 z (nr__) ! [0.008/0.098] Redshift determined from the
* velocity of the stellar absorption lines
real*4 e_z (nr__) ! [0.000006/0.00004] Uncertainty in z
character*16 Tel (nr__) ! Observatory and spectrograph (2)
*Note (1): J101836.77+512753.0 is very likely a misprint for J101835.77+512753.0
* (see Liu+, 2010, J/ApJ/708/427); corrected at CDS.
*Note (2): Observatory and spectrograph used:
* APO/DIS = APO Dual Imaging Spectrograph, 0.42"/pixel in the
* blue channel, 0.4"/pixel in the red channel;
* Keck/DEIMOS = Keck DEep Imaging Multi-Object Spectrograph, 0.12"/pixel;
* Lick/Kast = Lick Kast Spectrograph, 0.78"/pixel;
* MMT/Blue Channel = MMT Blue Channel Spectrograph, 0.29"/pixel;
* Palomar/DBSP = Palomar Double Spectrograph, 0.39"/pixel.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Measured luminosities
integer*4 nr__1
parameter (nr__1=71) ! Number of records
character*91 ar__1 ! Full-size record
character*10 Name_1 (nr__1) ! Galaxy Name
real*4 LOIII (nr__1) ! (10-7W) Observed [OIII] luminosity
real*4 e_LOIII (nr__1) ! (10-7W) Uncertainty in LOIII
real*4 LOIIIc (nr__1) ! (10-7W) Dereddened [OIII] luminosity
real*4 e_LOIIIc (nr__1) ! (10-7W) Uncertainty in LOIIIc
real*4 Lbol (nr__1) ! (10-7W) ? Bolometric luminosity
real*4 e_Lbol (nr__1) ! (10-7W) ? Uncertainty in Lbol
real*4 scale (nr__1) ! (kpc/arcsec) [0.1/1.9] Conversion from arcsec to kpc
* for Galaxy
integer*4 ak (nr__1) ! (pix) [4/55] Akaike width of the narrow-line
* region (2)
integer*4 e_ak (nr__1) ! (pix) [1/8] Uncertainty in ak
integer*4 Rad (nr__1) ! (pc) [950/10650] Radius of the narrow-line region
integer*4 e_Rad (nr__1) ! (pc) [80/1850] Uncertainty in Rad
*Note (2): The Akaike Information Criterion (AIC) is a least squares statistic
* that introduces a penalty for additional parameters, defined by
* Akaike (1974ITAC...19..716A). Numerically, AIC={chi}^2^+2k,
* where k is the number of parameters and {chi}^2^ is the
* chi-square statistic. See section 2.3 for further explanations.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Values for the parameters for quantitative
classification for both observed position angles
integer*4 nr__2
parameter (nr__2=142) ! Number of records
character*79 ar__2 ! Full-size record
character*10 Name_2 (nr__2) ! Galaxy Name (1)
integer*4 PAobs (nr__2) ! (deg) [0/252] Observation position angle
character*5 Nrow (nr__2) ! Number of Rows with >2 Gaussian components (2)
real*4 Vr (nr__2) ! (km/s) [16.6/664.4] Line of sight velocity from the
* single Gaussian centroid
real*4 e_Vr (nr__2) ! (km/s) [0.1/266] Uncertainty in Vr
real*4 sigma (nr__2) ! (km/s) [114/1039] Dispersion of Single Gauss fit
real*4 e_sigma (nr__2) ! (km/s) [0.8/507] Uncertainty in sigma
real*4 sigma1 (nr__2) ! (km/s) [122/1202] Dispersion of Double Gauss fit,
* component 1
real*4 e_sigma1 (nr__2) ! (km/s) [0.9/564] Uncertainty in sigma1
real*4 sigma2 (nr__2) ! (km/s) [61/1051] Dispersion of Double Gauss fit,
* component 2
real*4 e_sigma2 (nr__2) ! (km/s) [0.5/450] Uncertainty in sigma2
integer*4 PAgal (nr__2) ! (deg) [3/173] Galaxy photometric major axis position
* angle from SDSS r band photometry
integer*4 PAOIII (nr__2) ! (deg) [5/176] Galaxy [OIII] major axis position angle
*Note (1): J0731+4258 is very likely a misprint for J0731+4528 and
* J1350+3211 is very likely a misprint for J1350+2311; corrected at CDS.
*Note (2): The number of rows that are best fit by >2 Gaussians within the
* resolved center of emission (FWHM_cont +/-1 row). For instance, in
* row one, 5/9 indicates that out of the nine rows of the resolved
* center of emission, 5 rows are better fit by 3 Gaussians.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Official kinematic classification for each galaxy
integer*4 nr__3
parameter (nr__3=71) ! Number of records
character*76 ar__3 ! Full-size record
character*10 Name_3 (nr__3) ! Galaxy Name (1)
character*1 l_NG (nr__3) ! Limit flag on NG
integer*4 NG (nr__3) ! [2] Number of Gaussian components (2)
character*2 l_Vr (nr__3) ! [ <=>] Limit flag on Vr
integer*4 Vr_1 (nr__3) ! (km/s) [400] Maximum line of sight radial velocity
* measured from the single Gaussian fit
* for both PA
character*2 l_sigma (nr__3) ! [ <=>] Limit flag on sigma
integer*4 sigma_1 (nr__3) ! (km/s) [500] Largest dispersion of the single Gaussian
* fit from both observed PA
character*2 l_sigma12 (nr__3) ! [ <=>] Limit flag on sigma12
integer*4 sigma12 (nr__3) ! (km/s) [500] Largest dispersion of either of the two
* Gaussian fit from both PA
character*3 Align (nr__3) ! Alignment flag: yes=PAOIII equal PAgal (3)
character*32 Class (nr__3) ! Classification (includes 35 "Outflow" and
* 26 "outflow Composite")
*Note (1): The 1{sigma} error on the measured value straddles the classification
* cutoff for some objects, e.g. J0135+1435. However, as the measured
* value of J0135+1435 is greater than 500 km/s we classify this galaxy
* accordingly as Outflow Composite. For all other galaxies in our sample
* in which the measured value straddles a given classification cutoff
* within error are classified according to the measured value.
*Note (2): The number of components is >2 if at either position angle, more
* than half the rows within the spatial center of the galaxy have >2
* components as the best fit.
*Note (3): Alignment is determined from PA[OIII] and PAgal. If these two
* measurements are within 20{deg} of one another, the galaxy is aligned.
C=============================================================================
C Loading file 'table1.dat' ! Summary of longslit observations
C Format for file interpretation
1 format(A19,1X,A10,1X,F8.6,1X,E7.1,1X,A16)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,71
read(1,'(A64)')ar__
read(ar__,1)SDSS(i__),Name(i__),z(i__),e_z(i__),Tel(i__)
c ..............Just test output...........
write(6,1)SDSS(i__),Name(i__),z(i__),e_z(i__),Tel(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Measured luminosities
C Format for file interpretation
2 format(
+ A10,1X,E10.4,1X,E8.2,1X,E9.3,1X,E8.2,1X,E10.4,1X,E9.3,1X,F4.2,
+ 1X,I2,1X,I1,1X,I5,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,71
read(1,'(A91)')ar__1
read(ar__1,2)
+ Name_1(i__),LOIII(i__),e_LOIII(i__),LOIIIc(i__),e_LOIIIc(i__),
+ Lbol(i__),e_Lbol(i__),scale(i__),ak(i__),e_ak(i__),Rad(i__),
+ e_Rad(i__)
if(ar__1(51:60) .EQ. '') Lbol(i__) = rNULL__
if(ar__1(62:70) .EQ. '') e_Lbol(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),LOIII(i__),e_LOIII(i__),LOIIIc(i__),e_LOIIIc(i__),
+ Lbol(i__),e_Lbol(i__),scale(i__),ak(i__),e_ak(i__),Rad(i__),
+ e_Rad(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Values for the parameters for quantitative
* classification for both observed position angles
C Format for file interpretation
3 format(
+ A10,1X,I3,1X,A5,1X,F5.1,1X,F5.1,1X,F6.1,1X,F5.1,1X,F6.1,1X,
+ F5.1,1X,F6.1,1X,F5.1,1X,I3,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,142
read(1,'(A79)')ar__2
read(ar__2,3)
+ Name_2(i__),PAobs(i__),Nrow(i__),Vr(i__),e_Vr(i__),sigma(i__),
+ e_sigma(i__),sigma1(i__),e_sigma1(i__),sigma2(i__),
+ e_sigma2(i__),PAgal(i__),PAOIII(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),PAobs(i__),Nrow(i__),Vr(i__),e_Vr(i__),sigma(i__),
+ e_sigma(i__),sigma1(i__),e_sigma1(i__),sigma2(i__),
+ e_sigma2(i__),PAgal(i__),PAOIII(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Official kinematic classification for each galaxy
C Format for file interpretation
4 format(
+ A10,1X,A1,1X,I1,1X,A2,1X,I3,1X,A2,1X,I3,3X,A2,1X,I3,3X,A3,1X,
+ A32)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,71
read(1,'(A76)')ar__3
read(ar__3,4)
+ Name_3(i__),l_NG(i__),NG(i__),l_Vr(i__),Vr_1(i__),
+ l_sigma(i__),sigma_1(i__),l_sigma12(i__),sigma12(i__),
+ Align(i__),Class(i__)
c ..............Just test output...........
write(6,4)
+ Name_3(i__),l_NG(i__),NG(i__),l_Vr(i__),Vr_1(i__),
+ l_sigma(i__),sigma_1(i__),l_sigma12(i__),sigma12(i__),
+ Align(i__),Class(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end