FORTRAN Generation
(/./ftp/cats/J/MNRAS/458/2221)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/458/2221 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-Mar-28
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/458/2221    ATLAS3D Project. XXXI                    (Nyland+, 2016)
*================================================================================
*The ATLAS^3D^ Project.
*XXXI. Nuclear radio emission in nearby early-type galaxies.
*    Nyland K., Young L.M., Wrobel J.M., Sarzi M., Morganti R., Alatalo K.,
*    Blitz L., Bournaud F., Bureau M., Cappellari M., Crocker A.F., Davies R.L.,
*    Davis T.A., De Zeeuw P.T., Duc P.-A., Emsellem E., Khochfar S.,
*    Krajnovic D., Kuntschner H., McDermid R.M., Naab T., Oosterloo T.,
*    Scott N., Serra P., Weijmans A.-M.
*   <Mon. Not. R. Astron. Soc., 458, 2221-2268 (2016)>
*   =2016MNRAS.458.2221N    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'tablea1.dat'	! VLA 5GHz sample and flux density measurements

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

      character*1   f_Name      ! [a] Flag on Name (1)
      character*9   Name        ! Galaxy name
      character*1   n_Name      ! [b] Note on Name (1)
      real*4        Dist        ! (Mpc) Atlas3D distance (Paper I,
*                                  Cappellari et al. 2011, Cat. J/MNRAS/413/813)
      integer*4     Virgo       ! [0/1]  Virgo membership
      character*1   F_S         ! kinematic class (Paper III, Emsellem et al.,
*                                   2011, Cat. J/MNRAS/414/888) (G1)
      real*4        logMJAM     ! ([Msun]) ?=- Dynamically modelled stellar mass
*                                   (Paper XV, Cappellari et al.,
*                                   2013MNRAS.432.1709C)
      character*1   l_logMH2    ! Limit flag on logMH2
      real*4        logMH2      ! ([Msun]) Molecular hydrogen mass (Paper IV,
*                                   Young et al., 2011, Cat. J/MNRAS/414/940)
      character*1   l_SW91      ! Limit flag on SW91
      real*4        SW91        ! (mJy) ?=- Integrated flux density from the 5GHz
*                                   study of nearby ETGs carried-out by
*                                   Wrobel & Heeschen (1991AJ....101..148W)
*                                   at 5 arcsec resolution
      real*4        e_SW91      ! (mJy) ? rms uncertainty on SW91
      integer*4     rms         ! (uJy/beam) ?=- Average rms noise (3)
      character*1   l_Speak     ! Limit flag on Speak (3)
      real*4        Speak       ! (mJy/beam) ?=- Peak flux density (3)
      real*4        e_Speak     ! (mJy/beam) ? rms uncertainty on Speak (3)
      real*4        Sint        ! (mJy) ?=- Integrated flux density (3)(4)
      real*4        e_Sint      ! (mJy) ? rms uncertainty on Sint (3)
      character*1   l_logL      ! Limit flag on logL (3)
      real*4        logL        ! ([W/Hz]) ?=- 5 GHz radio luminosity (3)(5)
*Note (1): Notes as follows:
*   a = Extended source not well represented by a single two-dimensional Gaussian
*        model. The peak and integrated flux densities were calculated by drawing
*        an aperture at the 3xrms level around the source in the casa Viewer
*        and then using the IMSTAT task to determine the flux parameters.
*   b = Multicomponent source. The integrated flux density refers to the sum of
*        all components. See Table A3 for information on the properties of
*        individual components.
*Note (3): refer to data from our new high-resolution 5GHz VLA observations.
*Note (4): Note that measurements of the integrated flux density are only given
*   for sources that are formally resolved by JMFIT.
*Note (5): When an integrated flux density is given, L is based on the integrated
*   flux density. If only a peak flux density is given (either a measurement or
*   an upper limit), then L is based on the peak flux density.

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

C  Declarations for 'tablea2.dat'	! 5GHz spatial parameters of detections

      integer*4 nr__1
      parameter (nr__1=53)	! Number of records
      character*108 ar__1  	! 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*1   f_Name_1    ! [a] Flag on Name (1)
      character*9   Name_1      ! galaxy name
      character*1   n_Name_1    ! [b] Note on Name (1)
      character*1   Morph       ! [RU] Radio morphology based on the output of
*                                   the JMFIT task in aips (G2)
      integer*4     RAh         ! (h) Right ascension (J2000) of the emission at
*                                   the peak flux density (G3)
      integer*4     RAm         ! (min) Right ascension (J2000) of the emission
*                                   at the peak flux density (G3)
      real*4        RAs         ! (s) Right ascension (J2000) of the emission
*                                   at the peak flux density (G3)
      character*1   DE_         ! Declination sign (J2000) of the central
*                                       position of the emission (G3)
      integer*4     DEd         ! (deg) Declination of the central position of
*                                   the emission (G3)
      integer*4     DEm         ! (arcmin) Declination of the central position of
*                                   the emission (G3)
      real*4        DEs         ! (arcsec) Declination of the central position of
*                                   the emission (G3)
      real*4        Beam1       ! (arcsec) Angular dimension of the synthesized beam
*                                   (major axis)
      real*4        Beam2       ! (arcsec) Angular dimension of the synthesized beam
*                                   (minor axis)
      real*4        BPA         ! (deg) Beam position angle, measured anti clockwise
*                                   from North
      character*1   l_thetaM    ! Limit flag on thetaM
      real*4        thetaM      ! (arcsec) Radius upper limit or major axis of the
*                                   emission (4)
      real*4        e_thetaM    ! (arcsec) ? rms uncertainty on thetaM
      real*4        thetam_1    ! (arcsec) ? Minor axis of the emission (4)
      real*4        e_thetam_1  ! (arcsec) ? rms uncertainty on thetam
      real*4        PA          ! (deg) ?=- Position angle of the emission from JMFIT
      real*4        e_PA        ! (deg) ? rms uncertainty on PA
      character*1   l_Maj       ! Limit flag on Maj
      real*8        Maj         ! (pc) Linear radius upper limit or major axis of
*                                   the emission in physical units
      real*8        Min         ! (pc) ? Linear minor axis of the emission
*                                   in physical units
*Note (1): Notes as follows:
*  a = Extended source not well-represented by a two-dimensional Gaussian model.
*       The source dimensions were measured using the casa Viewer.
*  b = Multicomponent source. The spatial dimensions refer to the source closest
*       to the centre of the galaxy based on ground-based optical measurements.
*       See Table A4 for information on individual components.
*Note (4): Angular dimensions of the emission (major x minor axis). If JMFIT
*  was only able to deconvolve the major axis of the source, then the minor
*  axis extent is given as 0.00. The errors are from JMFIT and are only given
*  if the emission was successfully deconvolved in at least one dimension and
*  categorized as resolved. For non-Gaussian sources, source dimensions were
*  determined using the casa Viewer and no error is reported.

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

C  Declarations for 'tablea3.dat'	! 5GHz image properties of sources with multiple
                               components

      integer*4 nr__2
      parameter (nr__2=16)	! Number of records
      character*82 ar__2  	! Full-size record

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

      character*8   Name_2      ! Galaxy name
      character*19  Comp        ! Radio component name. The brightest
*                                     component in each galaxy is listed first
      integer*4     RAh_1       ! (h) Right ascension (J2000) at the location of
*                                     the peak flux density
      integer*4     RAm_1       ! (min) Right ascension (J2000) at the location of
*                                     the peak flux density
      real*4        RAs_1       ! (s) Right ascension (J2000) at the location of
*                                     the peak flux density
      character*1   DE__1       ! Declination sign (J2000) at the location of
*                                     the peak flux density
      integer*4     DEd_1       ! (deg) Declination (J2000) at the location of the
*                                     peak flux density
      integer*4     DEm_1       ! (arcmin) Declination (J2000) at the location of the
*                                     peak flux density
      real*4        DEs_1       ! (arcsec) Declination (J2000) at the location of the
*                                     peak flux density
      real*4        Speak_1     ! (mJy/beam) Peak flux density
      real*4        e_Speak_1   ! (mJy/beam) rms uncertainty on Speak
      real*4        Sint_1      ! (mJy) ?=-  Integrated flux density
      real*4        e_Sint_1    ! (mJy) ? rms uncertainty on Sint
      real*4        logL_1      ! ([W/Hz]) Radio luminosity (1)
*Note (1): When an integrated flux density is given, L is based on the integrated
*  flux density. If only a peak flux density is given (either a measurement or
*  an upper limit), then L is based on the peak flux density.

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

C  Declarations for 'tablea4.dat'	! 5GHz spatial properties of sources with
                               multiple components

      integer*4 nr__3
      parameter (nr__3=16)	! Number of records
      character*80 ar__3  	! Full-size record

      character*8   Name_3      ! Galaxy name
      character*19  Comp_1      ! Radio component name. The brightest component
*                                  in each galaxy is listed first
      character*1   Morph_1     ! [RU] Radio morphology as reported by JMFIT (G2)
      character*1   l_thetaM_1  ! Limit flag on thetaM
      real*4        thetaM_2    ! (arcsec) Radius upper limit or major axis of the
*                                  emission (G3)
      real*4        e_thetaM_2  ! (arcsec) ? rms uncertainty on thetaM
      real*4        thetam_3    ! (arcsec) ? Minor axis of the emission (G3)
      real*4        e_thetam_3  ! (arcsec) ? rms uncertainty on thetam
      real*4        PA_1        ! (deg) ?=- Position angle of the emission from JMFIT
      real*4        e_PA_1      ! (deg) ? rms uncertainty on PA
      character*1   l_Maj_1     ! Limit flag on Maj
      real*8        Maj_1       ! (pc) Linear radius upper limit or major axis of
*                                  the emission in physical units
      real*4        Min_1       ! (pc) ? Linear minor axis of the emission in
*                                  physical units

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

C  Declarations for 'tablea5.dat'	! *Archival high-resolution radio continuum data

      integer*4 nr__4
      parameter (nr__4=29)	! Number of records
      character*76 ar__4  	! Full-size record

      character*8   Name_4      ! Galaxy name
      character*1   n_Name_2    ! [abc] Note on Name (1)
      real*4        Dist_1      ! (Mpc) Atlas3D distance (Paper I,
*                                  Cappellari et al. 2011, Cat. J/MNRAS/413/813))
      integer*4     Virgo_1     ! [0/1]  Virgo membership
      character*1   F_S_1       ! Kinematic class (Paper III, Emsellem et al.,
*                                  2011, Cat. J/MNRAS/414/888) (G1)
      real*4        logMJAM_1   ! ([Msun]) ?=- Dynamically modelled stellar mass
*                                  (Paper XV, Cappellari et al.,
*                                  2013MNRAS.432.1709C))
      character*1   l_logMH2_1  ! Limit flag on logMH2
      real*4        logMH2_1    ! ([Msun]) Molecular hydrogen mass  (Paper IV,
*                                  Young et al., 2011, Cat. J/MNRAS/414/940)
      real*4        SW91_1      ! (mJy) ?=- Integrated flux density from the 5GHz study
*                                  of nearby ETGs carried-out by Wrobel &
*                                  Heeschen (1991AJ....101..148W)
*                                  at 5 arcsec resolution
      real*4        e_SW91_1    ! (mJy) ? rms uncertainty on SW91
      real*4        Freq        ! (GHz) Original observing frequency
      real*4        Res         ! (arcsec) Approximate synthesized beam major axis
      character*1   l_S5GHz     ! Limit flag on S5GHz
      real*8        S5GHz       ! (mJy) ?=- Integrated nuclear flux density at 5GHz (2)
      real*4        e_S5GHz     ! (mJy/beam) ? rms uncertainty on S5GHz
      character*1   l_logL_1    ! Limit flag on logL
      real*4        logL_2      ! ([W/Hz]) ?=- Nuclear radio luminosity at 5GHz
      character*1   Ref         ! References (3)
*Note (1): Notes as follows:
*   a = For NGC 4261, a milliarcsecond-scale spatial resolution VLBA radio flux
*        density is reported.
*   b = Although NGC 4697 and NGC 4477 were included in our new 5GHz VLA
*        observations, the sensitivity and quality of the final images were not
*        sufficient to detect them. We include a literature detection of NGC 4697
*        (Wrobel et al., 2008ApJ...686..838W) as well as a detection of NGC 4477
*        from project 12B-191 based on our own independent analysis.
*   c = For NGC 5846, the flux density has been extrapolated from to 5GHz using
*        a spectral index of -0.03. This spectral index was calculated based on
*        high-resolution measurements at 8.4 and 15GHz in Filho et al.
*        (2004A&A...418..429F) and Nagar et al. (2005, Cat. J/A+A/435/521),
*        respectively.
*Note (2): For measurements originally at frequencies other than 5GHz, the flux
*   densities listed here have been scaled to 5GHz using the source spectral
*   index reported in the literature. If no such spectral index information is
*   available, the flux density is scaled to 5GHz assuming a flat spectral index
*   of {alpha}=-0.1, where S~{nu}^{alpha}. Uncertainties are reported as given in
*   the literature when possible. If errors were not reported in the literature,
*   we estimate the flux density uncertainty as described in Section 5.
*Note (3): References as follows:
*   1 = Filho, Barthel & Ho (2006A&A...451...71F)
*   2 = Filho et al. (2004A&A...418..429F)
*   3 = Kharb et al. (2012AJ....143...78K)
*   4 = Nagar et al. (2001ApJ...559L..87N)
*   5 = Nagar et al. (2005, Cat. J/A+A/435/521)
*   6 = Project 12B-191
*   7 = Wrobel et al. (2008ApJ...686..838W)

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

C  Declarations for 'tablea6.dat'	! Additional galaxy properties

      integer*4 nr__5
      parameter (nr__5=148)	! Number of records
      character*59 ar__5  	! Full-size record

      character*9   Name_5      ! Galaxy name
      real*4        logMBH      ! ([Msun]) Black hole mass
      integer*4     r_logMBH    ! ? Black hole mass reference (1)
      character*1   l_logLX     ! Limit flag on logLX
      real*4        logLX       ! ([10-7W]) ?=- 2-10keV X-ray luminosity from
*                                          the literature
      integer*4     r_logLX     ! ?=- X-ray luminosity reference (2)
      character*1   l_log_LBol_LEdd ! Limit flag on log(LBol/LEdd)
      real*4        log_LBol_LEdd ! ?=- Eddington ratio (3)
      character*1   l_logRX     ! Limit flag on logRX
      real*4        logRX       ! ([-]) ?=- Radio-X-ray ratio as defined in
*                                          Terashima & Wilson
*                                          (2003ApJ...583..145T)
      real*4        logL_OIII_  ! ([10-7J]) ?=- Nuclear [OIII] luminosity
      real*4        EW_OIII_    ! (0.1nm) ?=- Nuclear [OIII] equivalent width
      character*5   Class       ! Nuclear emission line classification
*Note (1): for galaxies with dynamical mass estimates available in the
* literature. If no reference is given, the black hole mass is estimated from
* the MBH-{sigma}* relation (McConnell & Ma, 2013ApJ...764..184M).
*   References as follows:
*     1 = Krajnovic et al. (2002A&A...390..423K)
*     2 = Schulze & Gebhardt (2011ApJ...729...21S)
*     3 = Bower et al. (2001ApJ...550...75B)
*     4 = Barth et al. (2001ApJ...555..685B)
*     5 = van den Bosch & de Zeeuw (2010MNRAS.401.1770V)
*     6 = Nowak et al. (2010MNRAS.403..646N)
*     7 = Gultekin et al. (2009ApJ...695.1577G)
*     8 = Walsh et al. (2012ApJ...753...79W)
*     9 = Ferrarese, Ford & Jaffe (1996ApJ...470..444F)
*    10 = Walsh, Barth & Sarzi (2010ApJ...721..762W)
*    11 = Sarzi et al. (2001ApJ...550...65S)
*    12 = Rusli et al. (2013AJ....146..160R)
*    13 = Gebhardt et al. (2011ApJ...729..119G)
*    14 = Davis et al. (2014MNRAS.444.3427D)
*Note (2): Details on the computation of this parameter are provided in
*  Section 6.3.2.
*Note (3): References as follows:
*    15 = Pellegrini (2010ApJ...717..640P)
*    16 = Miller et al. (2012ApJ...747...57M)
*    17 = Kharb et al. (2012AJ....143...78K)

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

C  Declarations for 'tablec1.dat'	! Relative contour levels in the 5GHz continuum maps

      integer*4 nr__6
      parameter (nr__6=53)	! Number of records
      character*55 ar__6  	! Full-size record

      character*9   Name_6      ! Galaxy name
      integer*4     rms_1       ! (uJy/beam) rms noise
      character*42  RC          ! Relative contours

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

C  Loading file 'tablea1.dat'	! VLA 5GHz sample and flux density measurements

C  Format for file interpretation

    1 format(
     +  A1,A9,A1,1X,F4.1,1X,I1,1X,A1,1X,F5.2,1X,A1,F4.2,1X,A1,F4.1,1X,
     +  F3.1,1X,I2,1X,A1,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,1X,A1,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea1.dat')
      write(6,*) '....Loading file: tablea1.dat'
      do i__=1,125
        read(1,'(A76)')ar__
        read(ar__,1)
     +  f_Name,Name,n_Name,Dist,Virgo,F_S,logMJAM,l_logMH2,logMH2,
     +  l_SW91,SW91,e_SW91,rms,l_Speak,Speak,e_Speak,Sint,e_Sint,
     +  l_logL,logL
        if (idig(ar__(22:26)).EQ.0) logMJAM =  rNULL__
        if (idig(ar__(35:38)).EQ.0) SW91 =  rNULL__
        if(ar__(40:42) .EQ. '') e_SW91 = rNULL__
        if (rms .EQ. 45) rms =  iNULL__
        if (idig(ar__(48:52)).EQ.0) Speak =  rNULL__
        if(ar__(54:57) .EQ. '') e_Speak = rNULL__
        if (idig(ar__(60:64)).EQ.0) Sint =  rNULL__
        if(ar__(66:69) .EQ. '') e_Sint = rNULL__
        if (idig(ar__(72:76)).EQ.0) logL =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  f_Name,Name,n_Name,Dist,Virgo,F_S,logMJAM,l_logMH2,logMH2,
     +  l_SW91,SW91,e_SW91,rms,l_Speak,Speak,e_Speak,Sint,e_Sint,
     +  l_logL,logL
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea2.dat'	! 5GHz spatial parameters of detections

C  Format for file interpretation

    2 format(
     +  A1,A9,A1,1X,A1,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,
     +  F4.2,1X,F4.2,1X,F6.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
     +  F6.2,1X,F5.2,1X,A1,F7.2,1X,F7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea2.dat')
      write(6,*) '....Loading file: tablea2.dat'
      do i__=1,53
        read(1,'(A108)')ar__1
        read(ar__1,2)
     +  f_Name_1,Name_1,n_Name_1,Morph,RAh,RAm,RAs,DE_,DEd,DEm,DEs,
     +  Beam1,Beam2,BPA,l_thetaM,thetaM,e_thetaM,thetam_1,e_thetam_1,
     +  PA,e_PA,l_Maj,Maj,Min
        if(ar__1(65:68) .EQ. '') e_thetaM = rNULL__
        if(ar__1(70:73) .EQ. '') thetam_1 = rNULL__
        if(ar__1(75:78) .EQ. '') e_thetam_1 = rNULL__
        if (idig(ar__1(80:85)).EQ.0) PA =  rNULL__
        if(ar__1(87:91) .EQ. '') e_PA = rNULL__
        if(ar__1(102:108) .EQ. '') Min = 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,2)
     +  f_Name_1,Name_1,n_Name_1,Morph,RAh,RAm,RAs,DE_,DEd,DEm,DEs,
     +  Beam1,Beam2,BPA,l_thetaM,thetaM,e_thetaM,thetam_1,e_thetam_1,
     +  PA,e_PA,l_Maj,Maj,Min
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea3.dat'	! 5GHz image properties of sources with multiple
*                               components

C  Format for file interpretation

    3 format(
     +  A8,1X,A19,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2,
     +  1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea3.dat')
      write(6,*) '....Loading file: tablea3.dat'
      do i__=1,16
        read(1,'(A82)')ar__2
        read(ar__2,3)
     +  Name_2,Comp,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Speak_1,
     +  e_Speak_1,Sint_1,e_Sint_1,logL_1
        if (idig(ar__2(67:71)).EQ.0) Sint_1 =  rNULL__
        if(ar__2(73:76) .EQ. '') e_Sint_1 = rNULL__
        RAdeg_1 = rNULL__
        DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh_1*15.
        if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
        if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
        if(DEd_1 .GE. 0) DEdeg_1=DEd_1
        if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
        if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
        if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c    ..............Just test output...........
        write(6,3)
     +  Name_2,Comp,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Speak_1,
     +  e_Speak_1,Sint_1,e_Sint_1,logL_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea4.dat'	! 5GHz spatial properties of sources with
*                               multiple components

C  Format for file interpretation

    4 format(
     +  A8,1X,A19,1X,A1,1X,A1,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.2,1X,
     +  F5.2,1X,A1,F7.2,1X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea4.dat')
      write(6,*) '....Loading file: tablea4.dat'
      do i__=1,16
        read(1,'(A80)')ar__3
        read(ar__3,4)
     +  Name_3,Comp_1,Morph_1,l_thetaM_1,thetaM_2,e_thetaM_2,thetam_3,
     +  e_thetam_3,PA_1,e_PA_1,l_Maj_1,Maj_1,Min_1
        if(ar__3(38:41) .EQ. '') e_thetaM_2 = rNULL__
        if(ar__3(43:46) .EQ. '') thetam_3 = rNULL__
        if(ar__3(48:51) .EQ. '') e_thetam_3 = rNULL__
        if (idig(ar__3(53:58)).EQ.0) PA_1 =  rNULL__
        if(ar__3(60:64) .EQ. '') e_PA_1 = rNULL__
        if(ar__3(75:80) .EQ. '') Min_1 = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Name_3,Comp_1,Morph_1,l_thetaM_1,thetaM_2,e_thetaM_2,thetam_3,
     +  e_thetam_3,PA_1,e_PA_1,l_Maj_1,Maj_1,Min_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea5.dat'	! *Archival high-resolution radio continuum data

C  Format for file interpretation

    5 format(
     +  A8,A1,1X,F4.1,1X,I1,1X,A1,1X,F5.2,1X,A1,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F4.1,1X,F5.3,1X,A1,F7.2,1X,F5.2,1X,A1,F5.2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea5.dat')
      write(6,*) '....Loading file: tablea5.dat'
      do i__=1,29
        read(1,'(A76)')ar__4
        read(ar__4,5)
     +  Name_4,n_Name_2,Dist_1,Virgo_1,F_S_1,logMJAM_1,l_logMH2_1,
     +  logMH2_1,SW91_1,e_SW91_1,Freq,Res,l_S5GHz,S5GHz,e_S5GHz,
     +  l_logL_1,logL_2,Ref
        if (idig(ar__4(20:24)).EQ.0) logMJAM_1 =  rNULL__
        if (idig(ar__4(32:36)).EQ.0) SW91_1 =  rNULL__
        if(ar__4(38:41) .EQ. '') e_SW91_1 = rNULL__
        if (idig(ar__4(55:61)).EQ.0) S5GHz =  rNULL__
        if(ar__4(63:67) .EQ. '') e_S5GHz = rNULL__
        if (idig(ar__4(70:74)).EQ.0) logL_2 =  rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Name_4,n_Name_2,Dist_1,Virgo_1,F_S_1,logMJAM_1,l_logMH2_1,
     +  logMH2_1,SW91_1,e_SW91_1,Freq,Res,l_S5GHz,S5GHz,e_S5GHz,
     +  l_logL_1,logL_2,Ref
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea6.dat'	! Additional galaxy properties

C  Format for file interpretation

    6 format(
     +  A9,1X,F4.2,1X,I2,1X,A1,F5.2,1X,I2,1X,A1,F5.2,1X,A1,F5.2,1X,
     +  F5.2,1X,F5.2,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea6.dat')
      write(6,*) '....Loading file: tablea6.dat'
      do i__=1,148
        read(1,'(A59)')ar__5
        read(ar__5,6)
     +  Name_5,logMBH,r_logMBH,l_logLX,logLX,r_logLX,l_log_LBol_LEdd,
     +  log_LBol_LEdd,l_logRX,logRX,logL_OIII_,EW_OIII_,Class
        if(ar__5(16:17) .EQ. '') r_logMBH = iNULL__
        if (idig(ar__5(20:24)).EQ.0) logLX =  rNULL__
        if (r_logLX .EQ. 45) r_logLX =  iNULL__
        if (idig(ar__5(30:34)).EQ.0) log_LBol_LEdd =  rNULL__
        if (idig(ar__5(37:41)).EQ.0) logRX =  rNULL__
        if (idig(ar__5(43:47)).EQ.0) logL_OIII_ =  rNULL__
        if (idig(ar__5(49:53)).EQ.0) EW_OIII_ =  rNULL__
c    ..............Just test output...........
        write(6,6)
     +  Name_5,logMBH,r_logMBH,l_logLX,logLX,r_logLX,l_log_LBol_LEdd,
     +  log_LBol_LEdd,l_logRX,logRX,logL_OIII_,EW_OIII_,Class
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablec1.dat'	! Relative contour levels in the 5GHz continuum maps

C  Format for file interpretation

    7 format(A9,1X,I2,1X,A42)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablec1.dat')
      write(6,*) '....Loading file: tablec1.dat'
      do i__=1,53
        read(1,'(A55)')ar__6
        read(ar__6,7)Name_6,rms_1,RC
c    ..............Just test output...........
        write(6,7)Name_6,rms_1,RC
c    .......End.of.Just test output...........
      end do
      close(1)

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

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end