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 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.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        ! SDSS DR12 designation
*                                     (HHMMSS.ss+DDMMSS.s; J2000)
      character*15  SpID        ! Plate-MJD-fiber ID (G1)
      real*4        zem         ! [0/0.5] DR12Q redshift from visual inspection
      real*4        zOIII       ! [0/0.5] Redshift from [OIII]5007{AA} (G2)
      real*4        L5100       ! ([10-7W]) [42.2/49.1] Log continuum luminosity
*                                     at 5100{AA}; erg/s (G3)
      real*4        logM2       ! ([Msun]) [0/9.5] Log black hole mass
*                                     from Equation 2 (1)
      real*4        logM1       ! ([Msun]) [0/10.4] Log black hole mass
*                                     from Equation 1 (2)
      integer*4     FWHMB       ! (km/s) [0/19938] Broad H{alpha} FWHM with
*                                     single Gaussian (3)
      integer*4     FWHMB3      ! (km/s) [0/19326] Broad H{alpha} FWHM with
*                                     3 Gaussians (3)
      real*4        logLHaB     ! ([10-7W]) [0/44.2] Log broad H{alpha} line luminosity;
*                                     erg/s (3)
      real*4        e_logLHaB   ! ([10-7W]) [0/0.4] The 1{sigma} uncertainty in logLHaB
      real*4        EWHaB       ! (0.1nm) [0/797] Equivalent width of broad H{alpha};
*                                     Angstroms (3)
      real*4        e_EWHaB     ! (0.1nm) [0/114] The 1{sigma} uncertainty in EWHaB
      real*4        logLHa_     ! ([10-7W]) [0/44.2] Log broad + narrow H{alpha} line
*                                     luminosity; erg/s (3)
      real*4        logLHaN     ! ([10-7W]) [0/43.2] Log narrow H{alpha} line luminosity;
*                                     erg/s (4)
      real*4        e_logLHaN   ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logLHaN
      real*4        EWHaN       ! (0.1nm) [0/401] Equivalent width of narrow H{alpha};
*                                     Angstroms (4)
      real*4        e_EWHaN     ! (0.1nm) [0/36] The 1{sigma} uncertainty in EWHaN
      integer*4     FWHMN       ! (km/s) [0/1191] FWHM of [NII]6585{AA} line (5)
      real*4        logL6585    ! ([10-7W]) [0/43] Log [NII]6585{AA} line luminosity (5)
      real*4        e_logL6585  ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logL6585
      real*4        EW6585      ! (0.1nm) [0/173] Equivalent width of [NII]6585{AA} (5)
      real*4        e_EW6585    ! (0.1nm) [0/26] The 1{sigma} uncertainty in EW6585
      real*4        logL6549    ! ([10-7W]) [0/43] Log [NII]6549{AA} line luminosity (6)
      real*4        e_logL6549  ! ([10-7W]) [0/4] The 1{sigma} uncertainty in logL6549
      real*4        EW6549      ! (0.1nm) [0/59] Equivalent width of [NII]6549{AA} (6)
      real*4        e_EW6549    ! (0.1nm) [0/38] The 1{sigma} uncertainty in EW6549
      real*4        logL6732    ! ([10-7W]) [0/42] Log [SII]6732{AA} line luminosity (7)
      real*4        e_logL6732  ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logL6732
      real*4        EW6732      ! (0.1nm) [0/85] Equivalent width of [SII]6732{AA} (7)
      real*4        e_EW6732    ! (0.1nm) [0/13] The 1{sigma} uncertainty in EW6732
      real*4        logL6718    ! ([10-7W]) [0/42.1] Log [SII]6718{AA} line
*                                    luminosity (8)
      real*4        e_logL6718  ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logL6718
      real*4        EW6718      ! (0.1nm) [0/86] Equivalent width of [SII]6718{AA} (8)
      real*4        e_EW6718    ! (0.1nm) [0/12] The 1{sigma} uncertainty in EW6718
      real*4        chi3        ! [0.4/25.1] Reduced {chi}^2^ of 3 Gaussians
*                                     to broad H{alpha}
      real*4        chi         ! [0.5/78] Reduced {chi}^2^ of a Gaussian
*                                     to broad H{alpha}
      integer*4     ctpix       ! [208/264] Number of good pixels (9)
      real*4        S_N         ! [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      ! SDSS DR12 designation
*                                     (HHMMSS.ss+DDMMSS.s; J2000)
      character*15  SpID_1      ! Plate-MJD-fiber ID (G1)
      real*4        zem_1       ! [0.04/0.9] DR12Q redshift from visual
*                                    inspection
      real*4        zOIII_1     ! [0/0.9] Redshift from [OIII]5007{AA} (G2)
      real*4        L5100_1     ! ([10-7W]) [42.5/46] Log Continuum luminosity
*                                    at 5100{AA}; erg/s (G3)
      real*4        logM2_1     ! ([Msun]) [0/9.7] Log black hole mass
*                                    from Equation 2 (1)
      real*4        logM1_1     ! ([Msun]) [0/9.7] Log black hole mass
*                                    from Equation 1 (2)
      integer*4     FWHMB_1     ! (km/s) [0/12737] Broad H{beta} FWHM with
*                                    single Gaussian (3)
      integer*4     FWHMB3_1    ! (km/s) [0/12156] Broad H{beta} FWHM with
*                                    3 Gaussians (3)
      real*4        logLHbB     ! ([10-7W]) [0/45] Log broad H{beta} line luminosity;
*                                    erg/s (3)
      real*4        e_logLHbB   ! ([10-7W]) [0.02/1] The 1{sigma} uncertainty in logLHbB
      real*4        EWHbB       ! (0.1nm) [0/165] Equivalent width of broad H{beta};
*                                    Angstroms (3)
      real*4        e_EWHbB     ! (0.1nm) [0/49] The 1{sigma} uncertainty in EWHbB
      real*4        logLHb_     ! ([10-7W]) [0/44.2] Log broad + narrow H{beta} line
*                                    luminosity; erg/s (3)
      real*4        logLHbN     ! ([10-7W]) [0/43.1] Log narrow H{beta} line luminosity;
*                                    erg/s
      real*4        e_logLHbN   ! ([10-7W]) [0/0.6] The 1{sigma} uncertainty in logLHbN
      real*4        EWHbN       ! (0.1nm) [0/98] Equivalent width of narrow H{beta};
*                                    Angstroms (4)
      real*4        e_EWHbN     ! (0.1nm) [0/20] The 1{sigma} uncertainty in EWHbN
      integer*4     FWHMN_1     ! (km/s) [0/1195] FWHM of [NII]6585{AA} line (5)
      real*4        logL5007    ! ([10-7W]) [0/44] Log [OIII]5007{AA} line luminosity
      real*4        e_logL5007  ! ([10-7W]) [0/10] The 1{sigma} uncertainty in logL5007
      real*4        EW5007      ! (0.1nm) [0/1844] Equivalent width of
*                                     [OIII]5007{AA} (5)
      real*4        e_EW5007    ! (0.1nm) [0/103] The 1{sigma} uncertainty in EW5007
      real*4        logL4959    ! ([10-7W]) [0/44] Log [OIII]4959{AA} line luminosity
      real*4        e_logL4959  ! ([10-7W]) [0/28] The 1{sigma} uncertainty in logL4959
      real*4        EW4959      ! (0.1nm) [0/611] Equivalent width of
*                                    [OIII]4959{AA} (5)
      real*4        e_EW4959    ! (0.1nm) [0/88] The 1{sigma} uncertainty in EW4959
      real*4        chi3_1      ! [0.5/20] Reduced {chi}^2^ of 3 Gaussians
*                                    to broad H{beta}
      real*4        chi_1       ! [0.5/61] Reduced {chi}^2^ of a Gaussian
*                                    to broad H{beta}
      integer*4     ctpix_1     ! [267/355] Number of good pixels (6)
      real*4        S_N_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,SpID,zem,zOIII,L5100,logM2,logM1,FWHMB,FWHMB3,logLHaB,
     +  e_logLHaB,EWHaB,e_EWHaB,logLHa_,logLHaN,e_logLHaN,EWHaN,
     +  e_EWHaN,FWHMN,logL6585,e_logL6585,EW6585,e_EW6585,logL6549,
     +  e_logL6549,EW6549,e_EW6549,logL6732,e_logL6732,EW6732,
     +  e_EW6732,logL6718,e_logL6718,EW6718,e_EW6718,chi3,chi,ctpix,
     +  S_N
c    ..............Just test output...........
        write(6,1)
     +  SDSS,SpID,zem,zOIII,L5100,logM2,logM1,FWHMB,FWHMB3,logLHaB,
     +  e_logLHaB,EWHaB,e_EWHaB,logLHa_,logLHaN,e_logLHaN,EWHaN,
     +  e_EWHaN,FWHMN,logL6585,e_logL6585,EW6585,e_EW6585,logL6549,
     +  e_logL6549,EW6549,e_EW6549,logL6732,e_logL6732,EW6732,
     +  e_EW6732,logL6718,e_logL6718,EW6718,e_EW6718,chi3,chi,ctpix,
     +  S_N
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,SpID_1,zem_1,zOIII_1,L5100_1,logM2_1,logM1_1,FWHMB_1,
     +  FWHMB3_1,logLHbB,e_logLHbB,EWHbB,e_EWHbB,logLHb_,logLHbN,
     +  e_logLHbN,EWHbN,e_EWHbN,FWHMN_1,logL5007,e_logL5007,EW5007,
     +  e_EW5007,logL4959,e_logL4959,EW4959,e_EW4959,chi3_1,chi_1,
     +  ctpix_1,S_N_1
c    ..............Just test output...........
        write(6,2)
     +  SDSS_1,SpID_1,zem_1,zOIII_1,L5100_1,logM2_1,logM1_1,FWHMB_1,
     +  FWHMB3_1,logLHbB,e_logLHbB,EWHbB,e_EWHbB,logLHb_,logLHbN,
     +  e_logLHbN,EWHbN,e_EWHbN,FWHMN_1,logL5007,e_logL5007,EW5007,
     +  e_EW5007,logL4959,e_logL4959,EW4959,e_EW4959,chi3_1,chi_1,
     +  ctpix_1,S_N_1
c    .......End.of.Just test output...........
      end do
      close(1)

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