FORTRAN Generation
(/./ftp/cats/J/MNRAS/459/220)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/459/220 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-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/MNRAS/459/220 Optical redshifts of OH- and HI-emitting galaxies (Suess+, 2016)
*================================================================================
*Identifying OH imposters in the ALFALFA Neutral Hydrogen Survey.
*    Suess K.A., Darling J., Haynes M.P., Giovanelli R.
*   <Mon. Not. R. Astron. Soc., 459, 220-231 (2016)>
*   =2016MNRAS.459..220S    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Optical counterpart redshifts of HI-emitting
                                 galaxies

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

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

      character*17  Name       (nr__) ! Object name
      character*1   n_Name     (nr__) ! [2-6] Note on Name (1)
      integer*4     RAh        (nr__) ! (h) Optical counterpart right ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Optical counterpart right ascension (J2000)
      real*4        RAs        (nr__) ! (s) Optical counterpart right ascension (J2000)
      character*1   DE_        (nr__) ! Optical counterpart declination sign (J2000)
      integer*4     DEd        (nr__) ! (deg) Optical counterpart declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Optical counterpart declination (J2000)
      real*4        DEs        (nr__) ! (arcsec) Optical counterpart declination (J2000)
      character*6   n_Pos      (nr__) ! Note on Pos (2)
      integer*4     Vopt       (nr__) ! (km/s) Optical velocity (uncertainty is 130km/s)
      integer*4     VHI        (nr__) ! (km/s) ? HI velocity from Haynes et al.
*                                 (2011, Cat. J/AJ/142/170)
      integer*4     e_VHI      (nr__) ! (km/s) ? rms uncertainty on VHI
*Note (1): 2 Notes as follows:
*  2 = This spectrum showed both emission and absorption lines. The emission
*       lines were more prevalent, but the offset absorption lines clearly
*       showed in the blue half of the spectrum. Absorption was especially
*       evident in H{beta} and blueward Balmer series lines, where broad emission
*       lines were nearly divided in two by the offset absorption. The emission
*       velocity occurred at 5035km/s, a match for the ALFALFA HI line; the
*       absorption lines were offset to 5770km/s. The velocity listed in Table 3
*       is the emission velocity, which matches the ALFALFA HI velocity. The
*       galaxy exhibits disturbed morphology in SDSS, and the absorption and
*       emission lines may be offset because they come from different nuclei.
*  3 = This object was identified as a clear HI detection; in SDSS it appears
*       to be a blue compact dwarf galaxy.
*  4 = The optical counterpart for this object was confirmed to be the blue
*       extended object to the west of the central ALFALFA coordinates.
*  5 = Two objects were observed within the ALFALFA beam radius. One observed
*       galaxy was to the east of the ALFALFA target coordinates, and one to the
*       south. Both galaxies matched the HI velocity, so the coordinates and
*       measured velocities for both are included in this table.
*  6 = Most observed objects have at least five identifying spectral lines. For
*       this object, only H{alpha} was identified. Despite the scarcity of
*       observable optical lines in this spectra, we are confident of the HI
*       confirmation. The singular line is bright, shows extension along the
*       spatial axis of the spectrum, and it is very close to the expected HI
*       redshift. In visible light, the object appears to be a low surface
*       brightness galaxy.
*Note (2): For 249263, East and North optical counterparts.

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

C  Declarations for 'table4.dat'	! Ambiguous optical counterparts

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

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

      character*17  Name_1     (nr__1) ! Object name
      character*2   m_Name     (nr__1) ! [AB: ] Multiplicity index on Name
      character*1   n_Name_1   (nr__1) ! [2-9] Note on Name (1)
      integer*4     RAh_1      (nr__1) ! (h) Observed right ascension (J2000)
      integer*4     RAm_1      (nr__1) ! (min) Observed right ascension (J2000)
      real*4        RAs_1      (nr__1) ! (s) Observed right ascension (J2000)
      character*1   DE__1      (nr__1) ! Observed declination sign (J2000)
      integer*4     DEd_1      (nr__1) ! (deg) Observed declination (J2000)
      integer*4     DEm_1      (nr__1) ! (arcmin) Observed declination (J2000)
      real*4        DEs_1      (nr__1) ! (arcsec) Observed declination (J2000)
      character*2   n_Pos_1    (nr__1) ! Note on Pos (2)
      integer*4     Vopt_1     (nr__1) ! (km/s) ? Optical velocity (uncertainty is 130km/s)
      integer*4     VHI_1      (nr__1) ! (km/s) ? HI velocity from Haynes et al.
*                                 (2011, Cat. J/AJ/142/170)
      integer*4     e_VHI_1    (nr__1) ! (km/s) ? rms uncertainty on VHI
      integer*4     VOH        (nr__1) ! (km/s) ? OH velocity
      integer*4     e_VOH      (nr__1) ! (km/s) ? rms uncertainty on VOH
      character*17  Com        (nr__1) ! Comments (3)
*Note (1): Notes as follows:
*  2 = Possible pointing error during observations; this object may not be the
*       WISE bright source or the ALFALFA detection.
*  3 = The {alpha}.40 data release of the ALFALFA catalogue incorrectly states
*       that this object is an OHM (Haynes et al., 2011, Cat. J/AJ/142/170).
*       While the measured optical velocity is much higher than the ALFALFA
*       velocity, it does not match the OH velocity and the object's identity
*       remains unknown.
*  4 = The velocity determination for this object was measured from only one
*       line, presumed to be H{alpha}.
*  5 = Bleed-in from a nearby star obscured optical lines for this object.
*  6 = Two objects were observed within beam uncertainty of the ALFALFA
*       detection. The first had no visible optical lines, and the second (the
*       velocity listed in Table 4) showed broad line emission that matched
*       neither the OH nor the HI velocity.
*  7 = Due to high redshift, this object is likely an AGN.
*  8 = Guiding errors during observing rendered these frames unusable; no
*       further observations were made.
*  9 = Two possible optical counterparts were observed. One matched neither OH
*       nor HI velocities, and the other had no clear optical lines.
*Note (2): A: or B: for HI011145+290458
*Note (3): Comments as follows:
*   neither = observed velocity matches neither HI nor OH
*   no lines = objects for which no optical lines were observable and no velocity
*               determination could be made
*   ML = ALFALFA HI detection is a marginal line, likely not a real detections
*   UL = ALFALFA line is uncertain, and could be HI or another line

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

C  Loading file 'table3.dat'	! Optical counterpart redshifts of HI-emitting
*                                 galaxies

C  Format for file interpretation

    1 format(
     +  A17,A1,1X,I2,I2,F5.2,A1,I2,I2,F4.1,1X,A6,1X,I5,1X,I5,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,130
        read(1,'(A60)')ar__
        read(ar__,1)
     +  Name(i__),n_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),n_Pos(i__),Vopt(i__),VHI(i__),
     +  e_VHI(i__)
        if(ar__(52:56) .EQ. '') VHI(i__) = iNULL__
        if(ar__(58:60) .EQ. '') e_VHI(i__) = iNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),n_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),n_Pos(i__),Vopt(i__),VHI(i__),
     +  e_VHI(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Ambiguous optical counterparts

C  Format for file interpretation

    2 format(
     +  A17,A2,A1,1X,I2,I2,F5.2,A1,I2,I2,F4.1,1X,A2,1X,I6,2X,I5,1X,I2,
     +  1X,I6,1X,I2,1X,A17)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,62
        read(1,'(A87)')ar__1
        read(ar__1,2)
     +  Name_1(i__),m_Name(i__),n_Name_1(i__),RAh_1(i__),RAm_1(i__),
     +  RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
     +  n_Pos_1(i__),Vopt_1(i__),VHI_1(i__),e_VHI_1(i__),VOH(i__),
     +  e_VOH(i__),Com(i__)
        if(ar__1(44:49) .EQ. '') Vopt_1(i__) = iNULL__
        if(ar__1(52:56) .EQ. '') VHI_1(i__) = iNULL__
        if(ar__1(58:59) .EQ. '') e_VHI_1(i__) = iNULL__
        if(ar__1(61:66) .EQ. '') VOH(i__) = iNULL__
        if(ar__1(68:69) .EQ. '') e_VOH(i__) = iNULL__
        RAdeg_1(i__) = rNULL__
        DEdeg_1(i__) = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_1 from input data
c  (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
        if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
        if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Name_1(i__),m_Name(i__),n_Name_1(i__),RAh_1(i__),RAm_1(i__),
     +  RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
     +  n_Pos_1(i__),Vopt_1(i__),VHI_1(i__),e_VHI_1(i__),VOH(i__),
     +  e_VOH(i__),Com(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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