FORTRAN Generation
(/./ftp/cats/J/ApJS/234/16)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/234/16 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-Sep-09
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/234/16  H{alpha} & H{beta} spectral regions of low-z QSOs  (Chen+, 2018)
*================================================================================
*A catalog of quasar properties from the Baryon Oscillation Spectroscopic Survey.
*    Chen Z.-F., Pan D.-S., Pang T.-T., Huang Y.
*   <Astrophys. J. Suppl. Ser., 234, 16 (2018)>
*   =2018ApJS..234...16C
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! The catalog including the measurements around
                                H{alpha} spectral region

      integer*4 nr__
      parameter (nr__=2860)	! Number of records
      character*243 ar__   	! Full-size record

      character*18  SDSS       (nr__) ! SDSS DR12 designation
*                                     (HHMMSS.ss+DDMMSS.s; J2000)
      character*15  SpID       (nr__) ! Plate-MJD-fiber ID (G1)
      real*4        zem        (nr__) ! [0/0.5] DR12Q redshift from visual inspection
      real*4        zOIII      (nr__) ! [0/0.5] Redshift from [OIII]5007{AA} (G2)
      real*4        L5100      (nr__) ! ([10-7W]) [42.2/49.1] Log continuum luminosity
*                                     at 5100{AA}; erg/s (G3)
      real*4        logM2      (nr__) ! ([Msun]) [0/9.5] Log black hole mass
*                                     from Equation 2 (1)
      real*4        logM1      (nr__) ! ([Msun]) [0/10.4] Log black hole mass
*                                     from Equation 1 (2)
      integer*4     FWHMB      (nr__) ! (km/s) [0/19938] Broad H{alpha} FWHM with
*                                     single Gaussian (3)
      integer*4     FWHMB3     (nr__) ! (km/s) [0/19326] Broad H{alpha} FWHM with
*                                     3 Gaussians (3)
      real*4        logLHaB    (nr__) ! ([10-7W]) [0/44.2] Log broad H{alpha} line luminosity;
*                                     erg/s (3)
      real*4        e_logLHaB  (nr__) ! ([10-7W]) [0/0.4] The 1{sigma} uncertainty in logLHaB
      real*4        EWHaB      (nr__) ! (0.1nm) [0/797] Equivalent width of broad H{alpha};
*                                     Angstroms (3)
      real*4        e_EWHaB    (nr__) ! (0.1nm) [0/114] The 1{sigma} uncertainty in EWHaB
      real*4        logLHa_    (nr__) ! ([10-7W]) [0/44.2] Log broad + narrow H{alpha} line
*                                     luminosity; erg/s (3)
      real*4        logLHaN    (nr__) ! ([10-7W]) [0/43.2] Log narrow H{alpha} line luminosity;
*                                     erg/s (4)
      real*4        e_logLHaN  (nr__) ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logLHaN
      real*4        EWHaN      (nr__) ! (0.1nm) [0/401] Equivalent width of narrow H{alpha};
*                                     Angstroms (4)
      real*4        e_EWHaN    (nr__) ! (0.1nm) [0/36] The 1{sigma} uncertainty in EWHaN
      integer*4     FWHMN      (nr__) ! (km/s) [0/1191] FWHM of [NII]6585{AA} line (5)
      real*4        logL6585   (nr__) ! ([10-7W]) [0/43] Log [NII]6585{AA} line luminosity (5)
      real*4        e_logL6585 (nr__) ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logL6585
      real*4        EW6585     (nr__) ! (0.1nm) [0/173] Equivalent width of [NII]6585{AA} (5)
      real*4        e_EW6585   (nr__) ! (0.1nm) [0/26] The 1{sigma} uncertainty in EW6585
      real*4        logL6549   (nr__) ! ([10-7W]) [0/43] Log [NII]6549{AA} line luminosity (6)
      real*4        e_logL6549 (nr__) ! ([10-7W]) [0/4] The 1{sigma} uncertainty in logL6549
      real*4        EW6549     (nr__) ! (0.1nm) [0/59] Equivalent width of [NII]6549{AA} (6)
      real*4        e_EW6549   (nr__) ! (0.1nm) [0/38] The 1{sigma} uncertainty in EW6549
      real*4        logL6732   (nr__) ! ([10-7W]) [0/42] Log [SII]6732{AA} line luminosity (7)
      real*4        e_logL6732 (nr__) ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logL6732
      real*4        EW6732     (nr__) ! (0.1nm) [0/85] Equivalent width of [SII]6732{AA} (7)
      real*4        e_EW6732   (nr__) ! (0.1nm) [0/13] The 1{sigma} uncertainty in EW6732
      real*4        logL6718   (nr__) ! ([10-7W]) [0/42.1] Log [SII]6718{AA} line
*                                    luminosity (8)
      real*4        e_logL6718 (nr__) ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logL6718
      real*4        EW6718     (nr__) ! (0.1nm) [0/86] Equivalent width of [SII]6718{AA} (8)
      real*4        e_EW6718   (nr__) ! (0.1nm) [0/12] The 1{sigma} uncertainty in EW6718
      real*4        chi3       (nr__) ! [0.4/25.1] Reduced {chi}^2^ of 3 Gaussians
*                                     to broad H{alpha}
      real*4        chi        (nr__) ! [0.5/78] Reduced {chi}^2^ of a Gaussian
*                                     to broad H{alpha}
      integer*4     ctpix      (nr__) ! [208/264] Number of good pixels (9)
      real*4        S_N        (nr__) ! [2/90] Signal-to-Noise (9)
*Note (1): With the FWHM of the broad H{alpha} modeled by a combination of
*          three Gaussian functions. 0.0 indicates that the broad H{alpha}
*          has EW_Ha_<3{sigma}EW_Ha_.
*Note (2): With the (a,b)=(7.00,0.50) and the FWHM of the broad H{beta}
*          modeled by a combination of three Gaussian functions. 0.0 indicates
*          that the broad H{beta} has EW_H{beta}_<3{sigma}EW_H{beta}_.
*Note (3): 0.0 indicates that the broad H{alpha} component has
*          EW_H{alpha}_<3{sigma}EW_H{alpha}_.
*Note (4): 0.0 indicates that the narrow H{alpha} component has
*          EW_H{alpha}_<2{sigma}EW_H{alpha}_.
*Note (5): 0.0 indicates that the [NII]{lambda}6585 has
*           EW_[NII]{lambda}6585_<2{sigma}_EW[NII]{lambda}6585_
*Note (6): 0.0 indicates that the EW6549<2{sigma}EW6549.
*Note (7): 0.0 indicates that the EW6732<2{sigma}EW6732.
*Note (8): 0.0 indicates that the EW6718<2{sigma}EW6718.
*Note (9): In the H{alpha} spectral region; 6400-6800{AA}.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! The catalog including the measurements around
                                H{beta} spectral region

      integer*4 nr__1
      parameter (nr__1=21107)	! Number of records
      character*202 ar__1  	! Full-size record

      character*18  SDSS_1     (nr__1) ! SDSS DR12 designation
*                                     (HHMMSS.ss+DDMMSS.s; J2000)
      character*15  SpID_1     (nr__1) ! Plate-MJD-fiber ID (G1)
      real*4        zem_1      (nr__1) ! [0.04/0.9] DR12Q redshift from visual
*                                    inspection
      real*4        zOIII_1    (nr__1) ! [0/0.9] Redshift from [OIII]5007{AA} (G2)
      real*4        L5100_1    (nr__1) ! ([10-7W]) [42.5/46] Log Continuum luminosity
*                                    at 5100{AA}; erg/s (G3)
      real*4        logM2_1    (nr__1) ! ([Msun]) [0/9.7] Log black hole mass
*                                    from Equation 2 (1)
      real*4        logM1_1    (nr__1) ! ([Msun]) [0/9.7] Log black hole mass
*                                    from Equation 1 (2)
      integer*4     FWHMB_1    (nr__1) ! (km/s) [0/12737] Broad H{beta} FWHM with
*                                    single Gaussian (3)
      integer*4     FWHMB3_1   (nr__1) ! (km/s) [0/12156] Broad H{beta} FWHM with
*                                    3 Gaussians (3)
      real*4        logLHbB    (nr__1) ! ([10-7W]) [0/45] Log broad H{beta} line luminosity;
*                                    erg/s (3)
      real*4        e_logLHbB  (nr__1) ! ([10-7W]) [0.02/1] The 1{sigma} uncertainty in logLHbB
      real*4        EWHbB      (nr__1) ! (0.1nm) [0/165] Equivalent width of broad H{beta};
*                                    Angstroms (3)
      real*4        e_EWHbB    (nr__1) ! (0.1nm) [0/49] The 1{sigma} uncertainty in EWHbB
      real*4        logLHb_    (nr__1) ! ([10-7W]) [0/44.2] Log broad + narrow H{beta} line
*                                    luminosity; erg/s (3)
      real*4        logLHbN    (nr__1) ! ([10-7W]) [0/43.1] Log narrow H{beta} line luminosity;
*                                    erg/s
      real*4        e_logLHbN  (nr__1) ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logLHbN
      real*4        EWHbN      (nr__1) ! (0.1nm) [0/98] Equivalent width of narrow H{beta};
*                                    Angstroms (4)
      real*4        e_EWHbN    (nr__1) ! (0.1nm) [0/20] The 1{sigma} uncertainty in EWHbN
      integer*4     FWHMN_1    (nr__1) ! (km/s) [0/1195] FWHM of [NII]6585{AA} line (5)
      real*4        logL5007   (nr__1) ! ([10-7W]) [0/44] Log [OIII]5007{AA} line luminosity
      real*4        e_logL5007 (nr__1) ! ([10-7W]) [0/10] The 1{sigma} uncertainty in logL5007
      real*4        EW5007     (nr__1) ! (0.1nm) [0/1844] Equivalent width of
*                                     [OIII]5007{AA} (5)
      real*4        e_EW5007   (nr__1) ! (0.1nm) [0/103] The 1{sigma} uncertainty in EW5007
      real*4        logL4959   (nr__1) ! ([10-7W]) [0/44] Log [OIII]4959{AA} line luminosity
      real*4        e_logL4959 (nr__1) ! ([10-7W]) [0/28] The 1{sigma} uncertainty in logL4959
      real*4        EW4959     (nr__1) ! (0.1nm) [0/611] Equivalent width of
*                                    [OIII]4959{AA} (5)
      real*4        e_EW4959   (nr__1) ! (0.1nm) [0/88] The 1{sigma} uncertainty in EW4959
      real*4        chi3_1     (nr__1) ! [0.5/20] Reduced {chi}^2^ of 3 Gaussians
*                                    to broad H{beta}
      real*4        chi_1      (nr__1) ! [0.5/61] Reduced {chi}^2^ of a Gaussian
*                                    to broad H{beta}
      integer*4     ctpix_1    (nr__1) ! [267/355] Number of good pixels (6)
      real*4        S_N_1      (nr__1) ! [2/96] Signal-to-Noise (6)
*Note (1): With the  (a,b)=(6.91,0.50) and the FWHM of the broad H{beta}
*          modeled by a combination of three Gaussian functions. 0.0 indicates
*          that the broad H{beta} has EW_Hb_< 3{sigma}EW_Hb_.
*Note (2): With the (a,b)=(6.96,0.50) and the FWHM of the broad H{beta}
*          modeled by a combination of three Gaussian functions. 0.0 indicates
*          that the broad H{beta} has EW_H{beta}_< 3{sigma}EW_H{beta}_.
*Note (3): 0.0 indicates that the broad H{beta} component has
*          EW_H{beta}_<3{sigma}EW_H{beta}_.
*Note (4): 0.0 indicates that the broad H{beta} component has
*          EW_H{beta}_<2{sigma}EW_H{beta}_.
*Note (5): Which are the sum of the core and wing components. 0.0 indicates
*          that the EW5007<2{sigma}EW5007.
*Note (6): In the H{beta} spectral region; 4700-5100{AA}.

C=============================================================================

C  Loading file 'table1.dat'	! The catalog including the measurements around
*                                H{alpha} spectral region

C  Format for file interpretation

    1 format(
     +  A18,1X,A15,1X,F6.4,1X,F6.4,1X,F5.2,1X,F4.2,1X,F5.2,1X,I5,1X,
     +  I5,1X,F5.2,1X,F4.2,1X,F5.1,1X,F5.1,1X,F5.2,1X,F5.2,1X,F4.2,1X,
     +  F5.1,1X,F4.1,1X,I4,1X,F5.2,1X,F4.2,1X,F5.1,1X,F4.1,1X,F5.2,1X,
     +  F4.2,1X,F4.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F4.1,1X,F4.1,1X,F5.2,
     +  1X,F4.2,1X,F4.1,1X,F4.1,1X,F6.3,1X,F6.3,1X,I3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,2860
        read(1,'(A243)')ar__
        read(ar__,1)
     +  SDSS(i__),SpID(i__),zem(i__),zOIII(i__),L5100(i__),logM2(i__),
     +  logM1(i__),FWHMB(i__),FWHMB3(i__),logLHaB(i__),e_logLHaB(i__),
     +  EWHaB(i__),e_EWHaB(i__),logLHa_(i__),logLHaN(i__),
     +  e_logLHaN(i__),EWHaN(i__),e_EWHaN(i__),FWHMN(i__),
     +  logL6585(i__),e_logL6585(i__),EW6585(i__),e_EW6585(i__),
     +  logL6549(i__),e_logL6549(i__),EW6549(i__),e_EW6549(i__),
     +  logL6732(i__),e_logL6732(i__),EW6732(i__),e_EW6732(i__),
     +  logL6718(i__),e_logL6718(i__),EW6718(i__),e_EW6718(i__),
     +  chi3(i__),chi(i__),ctpix(i__),S_N(i__)
c    ..............Just test output...........
        write(6,1)
     +  SDSS(i__),SpID(i__),zem(i__),zOIII(i__),L5100(i__),logM2(i__),
     +  logM1(i__),FWHMB(i__),FWHMB3(i__),logLHaB(i__),e_logLHaB(i__),
     +  EWHaB(i__),e_EWHaB(i__),logLHa_(i__),logLHaN(i__),
     +  e_logLHaN(i__),EWHaN(i__),e_EWHaN(i__),FWHMN(i__),
     +  logL6585(i__),e_logL6585(i__),EW6585(i__),e_EW6585(i__),
     +  logL6549(i__),e_logL6549(i__),EW6549(i__),e_EW6549(i__),
     +  logL6732(i__),e_logL6732(i__),EW6732(i__),e_EW6732(i__),
     +  logL6718(i__),e_logL6718(i__),EW6718(i__),e_EW6718(i__),
     +  chi3(i__),chi(i__),ctpix(i__),S_N(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table2.dat'	! The catalog including the measurements around
*                                H{beta} spectral region

C  Format for file interpretation

    2 format(
     +  A18,1X,A15,1X,F6.4,1X,F6.4,1X,F5.2,1X,F4.2,1X,F4.2,1X,I5,1X,
     +  I5,1X,F5.2,1X,F4.2,1X,F5.1,1X,F4.1,1X,F5.2,1X,F5.2,1X,F4.2,1X,
     +  F4.1,1X,F4.1,1X,I4,1X,F5.2,1X,F4.2,1X,F6.1,1X,F5.1,1X,F5.2,1X,
     +  F5.2,1X,F5.1,1X,F4.1,1X,F6.3,1X,F6.3,1X,I3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,21107
        read(1,'(A202)')ar__1
        read(ar__1,2)
     +  SDSS_1(i__),SpID_1(i__),zem_1(i__),zOIII_1(i__),L5100_1(i__),
     +  logM2_1(i__),logM1_1(i__),FWHMB_1(i__),FWHMB3_1(i__),
     +  logLHbB(i__),e_logLHbB(i__),EWHbB(i__),e_EWHbB(i__),
     +  logLHb_(i__),logLHbN(i__),e_logLHbN(i__),EWHbN(i__),
     +  e_EWHbN(i__),FWHMN_1(i__),logL5007(i__),e_logL5007(i__),
     +  EW5007(i__),e_EW5007(i__),logL4959(i__),e_logL4959(i__),
     +  EW4959(i__),e_EW4959(i__),chi3_1(i__),chi_1(i__),ctpix_1(i__),
     +  S_N_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  SDSS_1(i__),SpID_1(i__),zem_1(i__),zOIII_1(i__),L5100_1(i__),
     +  logM2_1(i__),logM1_1(i__),FWHMB_1(i__),FWHMB3_1(i__),
     +  logLHbB(i__),e_logLHbB(i__),EWHbB(i__),e_EWHbB(i__),
     +  logLHb_(i__),logLHbN(i__),e_logLHbN(i__),EWHbN(i__),
     +  e_EWHbN(i__),FWHMN_1(i__),logL5007(i__),e_logL5007(i__),
     +  EW5007(i__),e_EW5007(i__),logL4959(i__),e_logL4959(i__),
     +  EW4959(i__),e_EW4959(i__),chi3_1(i__),chi_1(i__),ctpix_1(i__),
     +  S_N_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end