FORTRAN Generation
(/./ftp/cats/J/A_A/427/107)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/427/107 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 2024-Apr-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__=--2147483648.)  	! NULL real number
      parameter  (iNULL__=(-2147483647-1))	! NULL int  number
      integer    idig			! testing NULL number

C=============================================================================
Cat. J/A+A/427/107       Polarization in BAL QSOs                 (Lamy+, 2004)
*================================================================================
*Polarization properties of broad absorption line QSOs : new statistical clues.
*    Lamy H., Hutsemekers D.
*   <Astron. Astrophys. 427, 107 (2004)>
*   =2004A&A...427..107L
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Optical measurements for the BAL QSO sample

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*10  Name        ! Quasar name
      integer*4     RAh         ! (h) Right ascension (J2000)
      integer*4     RAm         ! (min) Right ascension (J2000)
      real*4        RAs         ! (s) Right ascension (J2000)
      character*1   DE_         ! Declination sign (J2000)
      integer*4     DEd         ! (deg) Declination (J2000)
      integer*4     DEm         ! (arcmin) Declination (J2000)
      real*4        DEs         ! (arcsec) ? Declination (J2000)
      real*4        z           ! Redshift
      integer*4     Type        ! BAL QSO type (1)
      real*4        Pol         ! (%) ? Polarization degree
      real*4        e_Pol       ! (%) ? Uncertainty of polarization degree
      real*4        Pol0        ! (%) ? Debiased polarization degree
      integer*4     r_Pol       ! ? References for polarization data (2)
      integer*4     BI          ! (km/s) ? Balnicity index (3)
      real*4        DI          ! ? Detachment index (4)
      integer*4     HWHM_CIV    ! (km/s) ? CIV Half Width at Half Maximum
      integer*4     HWHM__CIII_ ! (km/s) ? [CIII] Half Width at Half Maximum
      real*4        HREW_CIV    ! (0.1nm) ? CIV Rest Equivalent Width (5)
      real*4        REW__CIII_  ! (0.1nm) ? [CIII] Rest Equivalent Width
      real*4        REW_FeII2400 ! (0.1nm) ? FeII 2400 Rest Equivalent Width
      real*4        REW_FeII2070 ! (0.1nm) ? FeII 2070 Rest Equivalent Width
      real*4        alphaB      ! ? Continuum slope B
      integer*4     Vmax        ! (km/s) ? Maximum velocity in the CIV trough
      integer*4     Ref2        ! References for optical indices (6)
      real*4        BMAG        ! (mag) ? B absolute Magnitude
*Note (1): Object Type:
*     2 = HIBAL QSOs
*     3 = Strong LIBAL QSOs
*     4 = Weak LIBAL QSOs
*     5 = Marginal LIBAL QSOs
*     6 = Unclassified BAL QSO.
*Note (2): Reference for polarization data:
*     1 = Hutsemekers et al., 1998, Cat. <J/A+A/340/371>, Paper I
*     2 = Lamy & Hutsemekers, 2000, Cat. <J/A+AS/142/451>
*     3 = Schmidt & Hines, 1999ApJ...512..125S
*     4 = Ogle et al., 1999ApJS..125....1O (Keck observations)
*     5 = Ogle et al., 1999ApJS..125....1O (Palomar observations)
*     6 = Brotherton et al., 2002AJ....124.2575B
*     7 = Sluse et al., 2004, A&A, in prep.
*     8 = Hutsemekers et al., 2004, in prep.
*Note (3): The Balnicity index is a modified velocity equivalent width
*     of the C IV BAL.
*Note (4): The Detachment index measure the onset velocity of the
*    strongest C IV BAL trough in units of the adjacent emission line
*    half-width, that is the degree of detachment of the absorption line
*    relative to the emission one.
*Note (5): Only the red half of the CIV absorption line is considered
*    in the measurement of the Equivalent Width.
*Note (6): Reference for the spectra from which indices are measured.
*    They are encoded with 3 digits: each non-zero digit designates
*    a specific reference:
*      1: Hutsemekers et al., 1998, Cat. <J/A+A/340/371>, Paper I
*      2: Hutsemekers & Lamy, 2000A&A...358..835H, Paper II
*      3: Weymann et al., 1991ApJ...373...23W
*      4: Becker et al. 2000ApJ...538...72B, 2001, Cat. <J/ApJS/135/227>
*      5: Menou et al., 2001ApJ...561..645M
*      6: Korista et al., 1993ApJS...88..357K
*      7: Ogle et al., 1999ApJS..125....1O
*      8: spectra scanned from Large Bright Quasar Survey (LBQS)
*      9: Brotherton et al., 2002AJ....124.2575B

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

C  Loading file 'table1.dat'	! Optical measurements for the BAL QSO sample

C  Format for file interpretation

    1 format(
     +  A10,2X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,3X,F4.2,4X,I1,
     +  4X,F4.2,4X,F4.2,4X,F4.2,4X,I1,3X,I5,3X,F5.2,4X,I4,4X,I4,4X,
     +  F4.1,4X,F4.1,4X,F5.2,3X,F5.2,3X,F5.2,3X,I5,4X,I3,4X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,139
        read(1,'(A177)')ar__
        read(ar__,1)
     +  Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,Type,Pol,e_Pol,Pol0,r_Pol,
     +  BI,DI,HWHM_CIV,HWHM__CIII_,HREW_CIV,REW__CIII_,REW_FeII2400,
     +  REW_FeII2070,alphaB,Vmax,Ref2,BMAG
        if(ar__(33:36) .EQ. '') DEs = rNULL__
        if(ar__(53:56) .EQ. '') Pol = rNULL__
        if(ar__(61:64) .EQ. '') e_Pol = rNULL__
        if(ar__(69:72) .EQ. '') Pol0 = rNULL__
        if(ar__(77:77) .EQ. '') r_Pol = iNULL__
        if(ar__(81:85) .EQ. '') BI = iNULL__
        if(ar__(89:93) .EQ. '') DI = rNULL__
        if(ar__(98:101) .EQ. '') HWHM_CIV = iNULL__
        if(ar__(106:109) .EQ. '') HWHM__CIII_ = iNULL__
        if(ar__(114:117) .EQ. '') HREW_CIV = rNULL__
        if(ar__(122:125) .EQ. '') REW__CIII_ = rNULL__
        if(ar__(130:134) .EQ. '') REW_FeII2400 = rNULL__
        if(ar__(138:142) .EQ. '') REW_FeII2070 = rNULL__
        if(ar__(146:150) .EQ. '') alphaB = rNULL__
        if(ar__(154:158) .EQ. '') Vmax = iNULL__
        if(ar__(170:175) .EQ. '') BMAG = rNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,Type,Pol,e_Pol,Pol0,r_Pol,
     +  BI,DI,HWHM_CIV,HWHM__CIII_,HREW_CIV,REW__CIII_,REW_FeII2400,
     +  REW_FeII2070,alphaB,Vmax,Ref2,BMAG
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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