FORTRAN Generation
(/./ftp/cats/J/ApJS/226/19)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/226/19 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-Oct-13
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/226/19   FIR spectra of AGNs from Herschel   (Fernandez-Ontiveros+, 2016)
*================================================================================
*Far-infrared line spectra of active galaxies from the Herschel/PACS
*spectrometer: the complete database.
*    Fernandez-Ontiveros J.A., Spinoglio L., Pereira-Santaella M., Malkan M.A.,
*    Andreani P., Dasyra K.M.
*   <Astrophys. J. Suppl. Ser., 226, 19-19 (2016)>
*   =2016ApJS..226...19F    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Sample of active, starburst, and dwarf galaxies
                               observed with the Herschel/PACS spectrometer

      integer*4 nr__
      parameter (nr__=241)	! Number of records
      character*152 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*1   Type        ! Sample type (G1)
      integer*4     Seq         ! [1/170] Index number within Type
      character*15  Name        ! Galaxy Name
      integer*4     RAh         ! (h) 2MASS Hour of Right Ascension (J2000)
      integer*4     RAm         ! (min) 2MASS Minute of Right Ascension (J2000)
      real*8        RAs         ! (s) 2MASS Second of Right Ascension (J2000)
      character*1   DE_         ! 2MASS Sign of the Declination (J2000)
      integer*4     DEd         ! (deg) 2MASS Degree of Declination (J2000)
      integer*4     DEm         ! (arcmin) 2MASS Arcminute of Declination (J2000)
      real*8        DEs         ! (arcsec) 2MASS Arcsecond of Declination (J2000)
      real*8        z           ! [-0.0008/0.5] Redshift from NED
      real*8        Dist        ! (Mpc) [0.6/1360]? Redshift-independent distance
      integer*4     r_Dist      ! [1/36]? Reference for Dist
*                                     (see refs.dat file)
      character*5   Sp          ! Spectral Classification
      character*1   l_logLx     ! Limit flag on logLx
      real*4        logLx       ! ([10-7W]) [35.9/45.7]? Log absorption-corrected
*                                     2-10keV luminosity; erg/s
      character*1   CT          ! Compton Thick? Y: yes; N: no
      integer*4     r_logLx     ! [101/171]? Reference for logLx
*                                     (see refs.dat file)
      character*1   l_logLnev1  ! Limit flag on logLnev1
      real*4        logLnev1    ! ([10-7W]) [36.4/43.1]? Log Spitzer/IRS NeV 14.3um
*                                     luminosity; erg/s
      character*1   l_logLnev2  ! Limit flag on logLnev2
      real*4        logLnev2    ! ([10-7W]) [37.1/43.1]? Log Spitzer/IRS NeV 24.3um
*                                     luminosity; erg/s
      character*1   l_logLoiv   ! Limit flag on logLoiv
      real*4        logLoiv     ! ([10-7W]) [36.6/43.8]? Log Spitzer/IRS OIV 25.9um
*                                     luminosity; erg/s
      real*4        Z_1         ! (Sun) [0.02/2.8]? Metallicity derived from optical
*                                     spectra
      real*4        e_Z         ! (Sun) [0.001/8]? Uncertainty in Z
      integer*4     r_Z         ! ? Reference for Z (see refs.dat file)
      character*18  AName       ! Alternative Galaxy Name

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

C  Declarations for 'table4.dat'	! Far-IR fine structure line flux fluxes

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

      character*1   Type_1      ! Sample type (G1)
      integer*4     Seq_1       ! [1/170] Index number within Type
      character*15  Name_1      ! Galaxy Name
      character*1   l_Foiii52   ! Limit flag on Foiii52
      real*8        Foiii52     ! (10-17W/m2) [1/1363]? Herschel/PACS [OIII]52um line flux
      real*4        e_Foiii52   ! (10-17W/m2) [0.6/184]? Uncertainty in Foiii52
      character*1   l_Fniii57   ! Limit flag on Fniii57
      real*8        Fniii57     ! (10-17W/m2) [0.9/3288]? Herschel/PACS [NIII]57um
*                                     line flux
      real*4        e_Fniii57   ! (10-17W/m2) [0.4/179]? Uncertainty in Fniii57
      character*1   l_Foi63     ! Limit flag on Foi63
      real*8        Foi63       ! (10-17W/m2) [0.5/13352]? Herschel/PACS [OI]63um
*                                     line flux
      real*4        e_Foi63     ! (10-17W/m2) [0.2/404]? Uncertainty in Foi63
      character*1   l_Foiii88   ! Limit flag on Foiii88
      real*8        Foiii88     ! (10-17W/m2) [0.4/7653]? Herschel/PACS [OIII]88um
*                                     line flux
      real*4        e_Foiii88   ! (10-17W/m2) [0.2/236]? Uncertainty in Foiii88
      character*1   l_Fnii122   ! Limit flag on Fnii1222
      real*8        Fnii122     ! (10-17W/m2) [0.2/2002]? Herschel/PACS [NII]122um
*                                     line flux
      real*4        e_Fnii122   ! (10-17W/m2) [0.06/50]? Uncertainty in Fnii122
      character*1   l_Foi145    ! Limit flag on Foi145
      real*8        Foi145      ! (10-17W/m2) [0.1/1607]? Herschel/PACS [OI]145um
*                                     line flux
      real*4        e_Foi145    ! (10-17W/m2) [0.05/328]? Uncertainty in Foi145
      character*1   l_Fcii158   ! Limit flag on Fcii158
      real*8        Fcii158     ! (10-17W/m2) [0.5/12875]? Herschel/PACS [CII]158um
*                                     line flux
      real*4        e_Fcii158   ! (10-17W/m2) [0.08/345]? Uncertainty in Fcii158
      character*3   Aper        ! Aperture used (G2)

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

C  Declarations for 'table5.dat'	! Equivalent width for far-IR line equivalent widths

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

      character*1   Type_2      ! Sample type (G1)
      integer*4     Seq_2       ! [1/170] Index number within Type
      character*15  Name_2      ! Galaxy Name
      character*1   l_Woiii52   ! Limit flag on Woiii52
      real*4        Woiii52     ! (um) [0.005/0.6]? Herschel/PACS [OIII]52um line
*                               equivalent width
      real*4        e_Woiii52   ! (um) [0.003/0.6]? Uncertainty in Woiii52
      character*1   l_Wniii57   ! Limit flag on Wniii57
      real*4        Wniii57     ! (um) [0.003/0.1]? Herschel/PACS [NIII]57um line
*                               equivalent width
      real*4        e_Wniii57   ! (um) [0.001/0.02]? Uncertainty in Wniii57
      character*1   l_Woi63     ! Limit flag on Woi63
      real*4        Woi63       ! (um) [0.002/1.5]? Herschel/PACS [OI]63um line
*                               equivalent width
      real*4        e_Woi63     ! (um) [0.001/1.4]? Uncertainty in Woi63
      character*1   l_Woiii88   ! Limit flag on Woiii88
      real*4        Woiii88     ! (um) [0.004/2]? Herschel/PACS [OIII]88um line
*                               equivalent width
      real*4        e_Woiii88   ! (um) [0.001/2]? Uncertainty in Woiii88
      character*1   l_Wnii122   ! Limit flag on Wnii122
      real*4        Wnii122     ! (um) [0.004/1.2]? Herschel/PACS [NII]122um line
*                               equivalent width
      real*4        e_Wnii122   ! (um) [0.001/0.08]? Uncertainty in Wnii122
      character*1   l_Woi145    ! Limit flag on Woi145
      real*4        Woi145      ! (um) [0.005/1.1]? Herschel/PACS [OI]145um line
*                               equivalent width
      real*4        e_Woi145    ! (um) [0.001/0.07]? Uncertainty in Woi145
      character*1   l_Wcii158   ! Limit flag on Wcii158
      real*4        Wcii158     ! (um) [0.09/5.3]? Herschel/PACS [CII]158um line
*                               equivalent width
      real*4        e_Wcii158   ! (um) [0.003/2]? Uncertainty in Wcii158
      character*3   Aper_1      ! Aperture used (G2)

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

C  Declarations for 'table6.dat'	! Far-IR continuum fluxes

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

      character*1   Type_3      ! Sample type (G1)
      integer*4     Seq_3       ! [1/170] Index number within Type
      character*15  Name_3      ! Galaxy Name
      real*4        Coiii52     ! (Jy) [0.3/451]? Herschel/PACS [OIII]52um line
*                                   adjacent continuum flux
      real*4        e_Coiii52   ! (Jy) [0.04/15]? Uncertainty in Coiii52
      real*8        Cniii57     ! (Jy) [0.3/1039]? Herschel/PACS [NIII]57um line
*                                   adjacent continuum flux
      real*4        e_Cniii57   ! (Jy) [0.04/9]? Uncertainty in Cniii57
      character*1   l_Coi63     ! Limit flag on Coi63
      real*8        Coi63       ! (Jy) [0.04/1132]? Herschel/PACS [OI]63um line
*                                   adjacent continuum flux
      real*4        e_Coi63     ! (Jy) [0.04/16]? Uncertainty in Coi63
      character*1   l_Coiii88   ! Limit flag on Coiii88
      real*8        Coiii88     ! (Jy) [0.07/1207]? Herschel/PACS [OIII]88um line
*                                   adjacent continuum flux
      real*4        e_Coiii88   ! (Jy) [0.02/8]? Uncertainty Coiii88
      character*1   l_Cnii122   ! Limit flag on Cnii122
      real*8        Cnii122     ! (Jy) [0.06/1074]? Herschel/PACS [NII]122um line
*                                   adjacent continuum flux
      real*4        e_Cnii122   ! (Jy) [0.01/4]? Uncertainty in Cnii122
      character*1   l_Coi145    ! Limit flag on Coi145
      real*4        Coi145      ! (Jy) [0.05/934]? Herschel/PACS [OI]145um line
*                                   adjacent continuum flux
      real*4        e_Coi145    ! (Jy) [0.01/5]? Uncertainty in Coi145
      character*1   l_Ccii158   ! Limit flag on Ccii158
      real*4        Ccii158     ! (Jy) [0.07/798]? Herschel/PACS [CII]158um line
*                                   adjacent continuum flux
      real*4        e_Ccii158   ! (Jy) [0.02/36]? Uncertainty in Ccii158
      character*1   r_Ccii158   ! Reference for Ccii158 (see refs.dat file)
      character*3   Aper_2      ! Aperture used (G2)

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

C  Declarations for 'table7.dat'	! Gaussian sigma for far-IR fine structure lines

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

      character*1   Type_4      ! Sample type (G1)
      integer*4     Seq_4       ! [1/170] Index number within Type
      character*15  Name_4      ! Galaxy Name
      real*4        Soiii52     ! (km/s) [124/775]? Herschel/PACS [OIII]52um line
*                                 Gaussian sigma
      real*4        e_Soiii52   ! (km/s) [7/182]? Uncertainty in error for [OIII]52
      real*4        Sniii57     ! (km/s) [81/390]? Herschel/PACS [NIII]57um line
*                                 Gaussian sigma
      real*4        e_Sniii57   ! (km/s) [0.2/153]? Uncertainty in error for [NIII]57
      real*4        Soi63       ! (km/s) [60/633]? Herschel/PACS [OI]63um line
*                                 Gaussian sigma
      real*4        e_Soi63     ! (km/s) [0.1/66]? Uncertainty in error for [OI]63
      real*4        Soiii88     ! (km/s) [73/421]? Herschel/PACS [OIII]88um line
*                                 Gaussian sigma
      real*4        e_Soiii88   ! (km/s) [0.2/112]? Uncertainty in error for [OIII]88
      real*4        Snii122     ! (km/s) [102/911]? Herschel/PACS [NII]122um line
*                                 Gaussian sigma
      real*4        e_Snii122   ! (km/s) [0.4/135]? Uncertainty in error for [NII]122
      real*4        Soi145      ! (km/s) [108/344]? Herschel/PACS [OI]145um line
*                                 Gaussian sigma
      real*4        e_Soi145    ! (km/s) [0.5/614]? Uncertainty in error for [OI]145
      real*4        Scii158     ! (km/s) [89/340]? Herschel/PACS [CII]158um line
*                                 Gaussian sigma
      real*4        e_Scii158   ! (km/s) [0.2/100]? Uncertainty in error for [CII]158
      character*3   Aper_3      ! Aperture used (G2)

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

C  Declarations for 'table8.dat'	! Mid-IR fine structure line fluxes

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

      character*1   Type_5      ! Sample type (G1)
      integer*4     Seq_5       ! [1/170] Index number within Type
      character*15  Name_5      ! Galaxy Name
      character*1   l_Fsiv10    ! Limit flag on Fsiv10
      real*4        Fsiv10      ! (10-17W/m2) [0.07/537]? Spitzer/IRS [SIV]10.5um
*                                       line flux
      real*4        e_Fsiv10    ! (10-17W/m2) [0/54]? Uncertainty in Fsiv10
      character*1   l_Fneii12   ! Limit flag on Fneii12
      real*8        Fneii12     ! (10-17W/m2) [0.2/2833]? Spitzer/IRS [NeII]12.8um
*                                       line flux
      real*4        e_Fneii12   ! (10-17W/m2) [0.02/68]? Uncertainty in Fneii12
      character*1   l_Fnev14    ! Limit flag on Fnev14
      real*4        Fnev14      ! (10-17W/m2) [0.03/899]? Spitzer/IRS [NeV]14.3um
*                                       line flux
      real*4        e_Fnev14    ! (10-17W/m2) [0.01/23]? Uncertainty in Fnev14
      character*1   l_Fneiii15  ! Limit flag on Fneiii15
      real*8        Fneiii15    ! (10-17W/m2) [0.08/1371]? Spitzer/IRS [NeIII]15.6um
*                                       line flux
      real*4        e_Fneiii15  ! (10-17W/m2) [0.01/30]? Uncertainty in Fneiii15
      character*1   l_Fsiii18   ! Limit flag on Fsiii18
      real*4        Fsiii18     ! (10-17W/m2) [0.1/667]? Spitzer/IRS [SIII]18.7um
*                                       line flux
      real*4        e_Fsiii18   ! (10-17W/m2) [0/27]? Uncertainty in Fsiii18
      character*1   l_Fnev24    ! Limit flag on Fnev24
      real*4        Fnev24      ! (10-17W/m2) [0.07/816]? Spitzer/IRS [NeV]24.3um
*                                       line flux
      real*4        e_Fnev24    ! (10-17W/m2) [0.04/59]? Uncertainty in Fnev24
      character*1   l_Foiv25    ! Limit flag on Foiv25
      real*8        Foiv25      ! (10-17W/m2) [0.07/2030]? Spitzer/IRS [OIV]25.9um
*                                       line flux
      real*4        e_Foiv25    ! (10-17W/m2) [0.01/49]? Uncertainty in Foiv25
      character*1   l_Fsiii33   ! Limit flag on Fsiii33
      real*8        Fsiii33     ! (10-17W/m2) [0.09/1813]? Spitzer/IRS [SIII]33.5um
*                                       line flux
      real*4        e_Fsiii33   ! (10-17W/m2) [0.08/102]? Uncertainty in Fsiii33
      character*1   l_Fsii34    ! Limit flag on Fsii34
      real*8        Fsii34      ! (10-17W/m2) [1/2412]? Spitzer/IRS [SiII]34.8um
*                                       line flux
      real*4        e_Fsii34    ! (10-17W/m2) [0.3/81]? Uncertainty in Fsii34
      character*5   Mod         ! IRS Mode (1)
      character*23  Ref         ! Reference for the Spitzer spectra
*                                       (see refs.dat file)
      character*1   f_Name      ! Flag on Name (2)
*Note (1):
*    HR = High spectral resolution;
*    LR = Low spectral resolution.
*Note (2): Flag as follows:
*     b = Low S/N spectrum;
*     c = Only continuum, no line detected.

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

C  Declarations for 'table9.dat'	! [NII]205um, [CI]371um, and [CI]609um line fluxes
                               from Herschel/SPIRE

      integer*4 nr__6
      parameter (nr__6=91)	! Number of records
      character*109 ar__6  	! 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*14  Name_6      ! Galaxy Name
      integer*4     RAh_1       ! (h) 2MASS Hour of Right Ascension (J2000)
      integer*4     RAm_1       ! (min) 2MASS Minute of Right Ascension (J2000)
      real*8        RAs_1       ! (s) 2MASS Second of Right Ascension (J2000)
      character*1   DE__1       ! 2MASS Sign of the Declination (J2000)
      integer*4     DEd_1       ! (deg) 2MASS Degree of Declination (J2000)
      integer*4     DEm_1       ! (arcmin) 2MASS Arcminute of Declination (J2000)
      real*8        DEs_1       ! (arcsec) 2MASS Arcsecond of Declination (J2000)
      real*8        z_2         ! [-0.0008/0.2] Redshift from NED
      character*5   Sp_1        ! Spectral Classification
      character*1   l_Fnii205   ! Limit flag for Fnii205
      real*4        Fnii205     ! (10-17W/m2) [0.4/438]? Flux of the [NII]205um line
      real*4        e_Fnii205   ! (10-17W/m2) [0.1/12]? Uncertainty in Fnii205
      character*1   l_Fci371    ! Limit flag for Fci371
      real*4        Fci371      ! (10-17W/m2) [0.3/108]? Flux of the [CI]371um line
      real*4        e_Fci371    ! (10-17W/m2) [0/3]? Uncertainty in Fci371
      character*1   l_Fci609    ! Limit flag for Fci609
      real*4        Fci609      ! (10-17W/m2) [0.3/43]? Flux of the [CI]609um line
      real*4        e_Fci609    ! (10-17W/m2) [0/3]? Uncertainty in Fci609
      character*5   Ref_1       ! Line flux references (1)
*Note (1): References as follows:
*    1 = Gerin & Phillips (2000ApJ...537..644G);
*    2 = Israel & Baas (2002A&A...383...82I);
*    3 = Israel (2009A&A...493..525I);
*    4 = Israel et al. (2015A&A...578A..95I);
*    5 = Kamenetzky et al. (2014ApJ...795..174K);
*    6 = Kamenetzky et al. (2016ApJ...829...93K);
*    7 = Rigopoulou et al. (2013MNRAS.434.2051R);

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

C  Declarations for 'table10.dat'	! Density determinations

      integer*4 nr__7
      parameter (nr__7=132)	! Number of records
      character*130 ar__7  	! Full-size record

      character*15  Name_7      ! Galaxy Name
      character*5   Sp_2        ! Spectral Classification
      real*4        Rnii        ! [0.2/3]? Line ratio [NII] 205/122
      real*4        e_Rnii      ! [0.01/2]? Uncertainty in Rnii
      character*1   l_N1        ! Limit flag on N1
      real*4        N1          ! ([cm-3]) [0.1/3]? Log density derived
*                                 from [NII] 205/122 ratio
      real*4        e_N1        ! ([cm-3]) [0/3]? Lower error boundary in N1
      real*4        E_N1_1      ! ([cm-3]) [1/4]? Upper error boundary in N1
      real*4        Rsiii       ! [0.6/15]? Line ratio [SIII] 33/18
      real*4        e_Rsiii     ! [0.04/6]? Uncertainy in Rsiii
      character*1   l_N2        ! Limit flag on N2
      real*4        N2          ! ([cm-3]) [0.3/4]? Log density derived
*                                 from [SIII] 33/18 ratio
      real*4        e_N2        ! ([cm-3]) [0/4]? Lower error boundary in N2
      real*4        E_N2_1      ! ([cm-3]) [2/4]? Upper error boundary in N2
      real*4        Roiii       ! [0.1/2.2]? Line ratio [OIII]88/52
      real*4        e_Roiii     ! [0.09/0.7]? Error on the [OIII]88/52 line ratio
      character*1   l_N3        ! Limit flag on N3
      real*4        N3          ! ([cm-3]) [0.7/5]? Log density derived
*                                 from [OIII] 88/52 ratio
      real*4        e_N3        ! ([cm-3]) [0/4]? Lower error boundary in N3
      real*4        E_N3_1      ! ([cm-3]) [1/6]? Upper error boundary in N3
      real*4        Rnev        ! [0.4/6]? Line ratio [NeV] 24/14
      real*4        e_Rnev      ! [0.02/2]? Error in Rnev
      character*1   l_N4        ! Limit flag on N4
      real*4        N4          ! ([cm-3]) [1/4]? Log density derived
*                                 from [NeV] 24/14 ratio
      real*4        e_N4        ! ([cm-3]) [0/4]? Lower error boundary in N4
      real*4        E_N4_1      ! ([cm-3]) [2/5]? Upper error boundary in N4

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

C  Declarations for 'table11.dat'	! Far-IR fine structure line fluxes from
                               Herschel/PACS for the central spaxel, 3x3 spaxel
                               array, and 5x5 spaxel array

      integer*4 nr__8
      parameter (nr__8=594)	! Number of records
      character*145 ar__8  	! Full-size record

      character*1   Type_6      ! Sample type (G1)
      integer*4     Seq_6       ! [1/170] Index number within Type
      character*15  Name_8      ! Galaxy Name
      character*1   l_Foiii52_1 ! Limit flag on Foiii52
      real*8        Foiii52_1   ! (10-17W/m2) [1/1363]? Herschel/PACS [OIII]52um line flux
      real*4        e_Foiii52_1 ! (10-17W/m2) [0.6/184]? Uncertainty in Foiii52
      character*1   l_Fniii57_1 ! Limit flag on Fniii57
      real*8        Fniii57_1   ! (10-17W/m2) [0.9/3288]? Herschel/PACS [NIII]57um
*                                     line flux
      real*4        e_Fniii57_1 ! (10-17W/m2) [0.4/179]? Uncertainty in Fniii57
      character*1   l_Foi63_1   ! Limit flag on Foi63
      real*8        Foi63_1     ! (10-17W/m2) [0.5/13352]? Herschel/PACS [OI]63um
*                                     line flux
      real*4        e_Foi63_1   ! (10-17W/m2) [0.2/404]? Uncertainty in Foi63
      character*1   l_Foiii88_1 ! Limit flag on Foiii88
      real*8        Foiii88_1   ! (10-17W/m2) [0.4/7653]? Herschel/PACS [OIII]88um
*                                     line flux
      real*4        e_Foiii88_1 ! (10-17W/m2) [0.2/236]? Uncertainty in Foiii88
      character*1   l_Fnii122_1 ! Limit flag on Fnii1222
      real*8        Fnii122_1   ! (10-17W/m2) [0.2/2002]? Herschel/PACS [NII]122um
*                                     line flux
      real*4        e_Fnii122_1 ! (10-17W/m2) [0.06/50]? Uncertainty in Fnii122
      character*1   l_Foi145_1  ! Limit flag on Foi145
      real*8        Foi145_1    ! (10-17W/m2) [0.1/1607]? Herschel/PACS [OI]145um
*                                     line flux
      real*4        e_Foi145_1  ! (10-17W/m2) [0.05/328]? Uncertainty in Foi145
      character*1   l_Fcii158_1 ! Limit flag on Fcii158
      real*8        Fcii158_1   ! (10-17W/m2) [0.5/12875]? Herschel/PACS [CII]158um
*                                     line flux
      real*4        e_Fcii158_1 ! (10-17W/m2) [0.08/345]? Uncertainty in Fcii158
      character*3   Aper_4      ! Aperture used (G2)

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

C  Declarations for 'refs.dat'	! References

      integer*4 nr__9
      parameter (nr__9=179)	! Number of records
      character*105 ar__9  	! Full-size record

      character*3   Ref_2       ! Reference code
      character*19  BibCode     ! Bibcode
      character*28  Auth        ! First author's name(s)
      character*52  Comm        ! Comment

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

C  Loading file 'table3.dat'	! Sample of active, starburst, and dwarf galaxies
*                               observed with the Herschel/PACS spectrometer

C  Format for file interpretation

    1 format(
     +  A1,1X,I3,1X,A15,1X,I2,1X,I2,1X,F7.4,1X,A1,I2,1X,I2,1X,F7.4,1X,
     +  F10.7,1X,F8.3,1X,I2,1X,A5,1X,A1,1X,F5.2,1X,A1,1X,I3,1X,A1,1X,
     +  F5.2,1X,A1,1X,F5.2,1X,A1,1X,F5.2,1X,F5.3,1X,F5.3,1X,I3,1X,A18)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,241
        read(1,'(A152)')ar__
        read(ar__,1)
     +  Type,Seq,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,Dist,r_Dist,Sp,
     +  l_logLx,logLx,CT,r_logLx,l_logLnev1,logLnev1,l_logLnev2,
     +  logLnev2,l_logLoiv,logLoiv,Z_1,e_Z,r_Z,AName
        if(ar__(63:70) .EQ. '') Dist = rNULL__
        if(ar__(72:73) .EQ. '') r_Dist = iNULL__
        if(ar__(83:87) .EQ. '') logLx = rNULL__
        if(ar__(91:93) .EQ. '') r_logLx = iNULL__
        if(ar__(97:101) .EQ. '') logLnev1 = rNULL__
        if(ar__(105:109) .EQ. '') logLnev2 = rNULL__
        if(ar__(113:117) .EQ. '') logLoiv = rNULL__
        if(ar__(119:123) .EQ. '') Z_1 = rNULL__
        if(ar__(125:129) .EQ. '') e_Z = rNULL__
        if(ar__(131:133) .EQ. '') r_Z = iNULL__
        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)
     +  Type,Seq,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,Dist,r_Dist,Sp,
     +  l_logLx,logLx,CT,r_logLx,l_logLnev1,logLnev1,l_logLnev2,
     +  logLnev2,l_logLoiv,logLoiv,Z_1,e_Z,r_Z,AName
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Far-IR fine structure line flux fluxes

C  Format for file interpretation

    2 format(
     +  A1,1X,I3,1X,A15,1X,A1,1X,F7.2,1X,F6.2,1X,A1,1X,F7.2,1X,F6.2,
     +  1X,A1,1X,F8.2,1X,F6.2,1X,A1,1X,F7.2,1X,F6.2,1X,A1,1X,F7.2,1X,
     +  F5.2,1X,A1,1X,F7.2,1X,F6.2,1X,A1,1X,F8.2,1X,F6.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,198
        read(1,'(A145)')ar__1
        read(ar__1,2)
     +  Type_1,Seq_1,Name_1,l_Foiii52,Foiii52,e_Foiii52,l_Fniii57,
     +  Fniii57,e_Fniii57,l_Foi63,Foi63,e_Foi63,l_Foiii88,Foiii88,
     +  e_Foiii88,l_Fnii122,Fnii122,e_Fnii122,l_Foi145,Foi145,
     +  e_Foi145,l_Fcii158,Fcii158,e_Fcii158,Aper
        if(ar__1(25:31) .EQ. '') Foiii52 = rNULL__
        if(ar__1(33:38) .EQ. '') e_Foiii52 = rNULL__
        if(ar__1(42:48) .EQ. '') Fniii57 = rNULL__
        if(ar__1(50:55) .EQ. '') e_Fniii57 = rNULL__
        if(ar__1(59:66) .EQ. '') Foi63 = rNULL__
        if(ar__1(68:73) .EQ. '') e_Foi63 = rNULL__
        if(ar__1(77:83) .EQ. '') Foiii88 = rNULL__
        if(ar__1(85:90) .EQ. '') e_Foiii88 = rNULL__
        if(ar__1(94:100) .EQ. '') Fnii122 = rNULL__
        if(ar__1(102:106) .EQ. '') e_Fnii122 = rNULL__
        if(ar__1(110:116) .EQ. '') Foi145 = rNULL__
        if(ar__1(118:123) .EQ. '') e_Foi145 = rNULL__
        if(ar__1(127:134) .EQ. '') Fcii158 = rNULL__
        if(ar__1(136:141) .EQ. '') e_Fcii158 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Type_1,Seq_1,Name_1,l_Foiii52,Foiii52,e_Foiii52,l_Fniii57,
     +  Fniii57,e_Fniii57,l_Foi63,Foi63,e_Foi63,l_Foiii88,Foiii88,
     +  e_Foiii88,l_Fnii122,Fnii122,e_Fnii122,l_Foi145,Foi145,
     +  e_Foi145,l_Fcii158,Fcii158,e_Fcii158,Aper
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Equivalent width for far-IR line equivalent widths

C  Format for file interpretation

    3 format(
     +  A1,1X,I3,1X,A15,1X,A1,1X,F5.3,1X,F5.3,1X,A1,1X,F5.3,1X,F5.3,
     +  1X,A1,1X,F5.3,1X,F5.3,1X,A1,1X,F5.3,1X,F5.3,1X,A1,1X,F5.3,1X,
     +  F5.3,1X,A1,1X,F5.3,1X,F5.3,1X,A1,1X,F5.3,1X,F5.3,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,198
        read(1,'(A123)')ar__2
        read(ar__2,3)
     +  Type_2,Seq_2,Name_2,l_Woiii52,Woiii52,e_Woiii52,l_Wniii57,
     +  Wniii57,e_Wniii57,l_Woi63,Woi63,e_Woi63,l_Woiii88,Woiii88,
     +  e_Woiii88,l_Wnii122,Wnii122,e_Wnii122,l_Woi145,Woi145,
     +  e_Woi145,l_Wcii158,Wcii158,e_Wcii158,Aper_1
        if(ar__2(25:29) .EQ. '') Woiii52 = rNULL__
        if(ar__2(31:35) .EQ. '') e_Woiii52 = rNULL__
        if(ar__2(39:43) .EQ. '') Wniii57 = rNULL__
        if(ar__2(45:49) .EQ. '') e_Wniii57 = rNULL__
        if(ar__2(53:57) .EQ. '') Woi63 = rNULL__
        if(ar__2(59:63) .EQ. '') e_Woi63 = rNULL__
        if(ar__2(67:71) .EQ. '') Woiii88 = rNULL__
        if(ar__2(73:77) .EQ. '') e_Woiii88 = rNULL__
        if(ar__2(81:85) .EQ. '') Wnii122 = rNULL__
        if(ar__2(87:91) .EQ. '') e_Wnii122 = rNULL__
        if(ar__2(95:99) .EQ. '') Woi145 = rNULL__
        if(ar__2(101:105) .EQ. '') e_Woi145 = rNULL__
        if(ar__2(109:113) .EQ. '') Wcii158 = rNULL__
        if(ar__2(115:119) .EQ. '') e_Wcii158 = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Type_2,Seq_2,Name_2,l_Woiii52,Woiii52,e_Woiii52,l_Wniii57,
     +  Wniii57,e_Wniii57,l_Woi63,Woi63,e_Woi63,l_Woiii88,Woiii88,
     +  e_Woiii88,l_Wnii122,Wnii122,e_Wnii122,l_Woi145,Woi145,
     +  e_Woi145,l_Wcii158,Wcii158,e_Wcii158,Aper_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Far-IR continuum fluxes

C  Format for file interpretation

    4 format(
     +  A1,1X,I3,1X,A15,1X,F6.2,1X,F5.2,1X,F7.2,1X,F4.2,1X,A1,1X,F7.2,
     +  1X,F5.2,1X,A1,1X,F7.2,1X,F4.2,1X,A1,1X,F7.2,1X,F4.2,1X,A1,1X,
     +  F6.2,1X,F4.2,1X,A1,1X,F6.2,1X,F5.2,1X,A1,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,198
        read(1,'(A128)')ar__3
        read(ar__3,4)
     +  Type_3,Seq_3,Name_3,Coiii52,e_Coiii52,Cniii57,e_Cniii57,
     +  l_Coi63,Coi63,e_Coi63,l_Coiii88,Coiii88,e_Coiii88,l_Cnii122,
     +  Cnii122,e_Cnii122,l_Coi145,Coi145,e_Coi145,l_Ccii158,Ccii158,
     +  e_Ccii158,r_Ccii158,Aper_2
        if(ar__3(23:28) .EQ. '') Coiii52 = rNULL__
        if(ar__3(30:34) .EQ. '') e_Coiii52 = rNULL__
        if(ar__3(36:42) .EQ. '') Cniii57 = rNULL__
        if(ar__3(44:47) .EQ. '') e_Cniii57 = rNULL__
        if(ar__3(51:57) .EQ. '') Coi63 = rNULL__
        if(ar__3(59:63) .EQ. '') e_Coi63 = rNULL__
        if(ar__3(67:73) .EQ. '') Coiii88 = rNULL__
        if(ar__3(75:78) .EQ. '') e_Coiii88 = rNULL__
        if(ar__3(82:88) .EQ. '') Cnii122 = rNULL__
        if(ar__3(90:93) .EQ. '') e_Cnii122 = rNULL__
        if(ar__3(97:102) .EQ. '') Coi145 = rNULL__
        if(ar__3(104:107) .EQ. '') e_Coi145 = rNULL__
        if(ar__3(111:116) .EQ. '') Ccii158 = rNULL__
        if(ar__3(118:122) .EQ. '') e_Ccii158 = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Type_3,Seq_3,Name_3,Coiii52,e_Coiii52,Cniii57,e_Cniii57,
     +  l_Coi63,Coi63,e_Coi63,l_Coiii88,Coiii88,e_Coiii88,l_Cnii122,
     +  Cnii122,e_Cnii122,l_Coi145,Coi145,e_Coi145,l_Ccii158,Ccii158,
     +  e_Ccii158,r_Ccii158,Aper_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Gaussian sigma for far-IR fine structure lines

C  Format for file interpretation

    5 format(
     +  A1,1X,I3,1X,A15,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,
     +  F4.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,
     +  1X,F5.1,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,190
        read(1,'(A108)')ar__4
        read(ar__4,5)
     +  Type_4,Seq_4,Name_4,Soiii52,e_Soiii52,Sniii57,e_Sniii57,Soi63,
     +  e_Soi63,Soiii88,e_Soiii88,Snii122,e_Snii122,Soi145,e_Soi145,
     +  Scii158,e_Scii158,Aper_3
        if(ar__4(23:27) .EQ. '') Soiii52 = rNULL__
        if(ar__4(29:33) .EQ. '') e_Soiii52 = rNULL__
        if(ar__4(35:39) .EQ. '') Sniii57 = rNULL__
        if(ar__4(41:45) .EQ. '') e_Sniii57 = rNULL__
        if(ar__4(47:51) .EQ. '') Soi63 = rNULL__
        if(ar__4(53:56) .EQ. '') e_Soi63 = rNULL__
        if(ar__4(58:62) .EQ. '') Soiii88 = rNULL__
        if(ar__4(64:68) .EQ. '') e_Soiii88 = rNULL__
        if(ar__4(70:74) .EQ. '') Snii122 = rNULL__
        if(ar__4(76:80) .EQ. '') e_Snii122 = rNULL__
        if(ar__4(82:86) .EQ. '') Soi145 = rNULL__
        if(ar__4(88:92) .EQ. '') e_Soi145 = rNULL__
        if(ar__4(94:98) .EQ. '') Scii158 = rNULL__
        if(ar__4(100:104) .EQ. '') e_Scii158 = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Type_4,Seq_4,Name_4,Soiii52,e_Soiii52,Sniii57,e_Sniii57,Soi63,
     +  e_Soi63,Soiii88,e_Soiii88,Snii122,e_Snii122,Soi145,e_Soi145,
     +  Scii158,e_Scii158,Aper_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Mid-IR fine structure line fluxes

C  Format for file interpretation

    6 format(
     +  A1,1X,I3,1X,A15,1X,A1,1X,F6.2,1X,F5.2,1X,A1,1X,F7.2,1X,F5.2,
     +  1X,A1,1X,F6.2,1X,F5.2,1X,A1,1X,F7.2,1X,F5.2,1X,A1,1X,F6.2,1X,
     +  F5.2,1X,A1,1X,F6.2,1X,F5.2,1X,A1,1X,F7.2,1X,F5.2,1X,A1,1X,
     +  F7.2,1X,F6.2,1X,A1,1X,F7.2,1X,F5.2,1X,A5,1X,A23,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,198
        read(1,'(A194)')ar__5
        read(ar__5,6)
     +  Type_5,Seq_5,Name_5,l_Fsiv10,Fsiv10,e_Fsiv10,l_Fneii12,
     +  Fneii12,e_Fneii12,l_Fnev14,Fnev14,e_Fnev14,l_Fneiii15,
     +  Fneiii15,e_Fneiii15,l_Fsiii18,Fsiii18,e_Fsiii18,l_Fnev24,
     +  Fnev24,e_Fnev24,l_Foiv25,Foiv25,e_Foiv25,l_Fsiii33,Fsiii33,
     +  e_Fsiii33,l_Fsii34,Fsii34,e_Fsii34,Mod,Ref,f_Name
        if(ar__5(25:30) .EQ. '') Fsiv10 = rNULL__
        if(ar__5(32:36) .EQ. '') e_Fsiv10 = rNULL__
        if(ar__5(40:46) .EQ. '') Fneii12 = rNULL__
        if(ar__5(48:52) .EQ. '') e_Fneii12 = rNULL__
        if(ar__5(56:61) .EQ. '') Fnev14 = rNULL__
        if(ar__5(63:67) .EQ. '') e_Fnev14 = rNULL__
        if(ar__5(71:77) .EQ. '') Fneiii15 = rNULL__
        if(ar__5(79:83) .EQ. '') e_Fneiii15 = rNULL__
        if(ar__5(87:92) .EQ. '') Fsiii18 = rNULL__
        if(ar__5(94:98) .EQ. '') e_Fsiii18 = rNULL__
        if(ar__5(102:107) .EQ. '') Fnev24 = rNULL__
        if(ar__5(109:113) .EQ. '') e_Fnev24 = rNULL__
        if(ar__5(117:123) .EQ. '') Foiv25 = rNULL__
        if(ar__5(125:129) .EQ. '') e_Foiv25 = rNULL__
        if(ar__5(133:139) .EQ. '') Fsiii33 = rNULL__
        if(ar__5(141:146) .EQ. '') e_Fsiii33 = rNULL__
        if(ar__5(150:156) .EQ. '') Fsii34 = rNULL__
        if(ar__5(158:162) .EQ. '') e_Fsii34 = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  Type_5,Seq_5,Name_5,l_Fsiv10,Fsiv10,e_Fsiv10,l_Fneii12,
     +  Fneii12,e_Fneii12,l_Fnev14,Fnev14,e_Fnev14,l_Fneiii15,
     +  Fneiii15,e_Fneiii15,l_Fsiii18,Fsiii18,e_Fsiii18,l_Fnev24,
     +  Fnev24,e_Fnev24,l_Foiv25,Foiv25,e_Foiv25,l_Fsiii33,Fsiii33,
     +  e_Fsiii33,l_Fsii34,Fsii34,e_Fsii34,Mod,Ref,f_Name
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! [NII]205um, [CI]371um, and [CI]609um line fluxes
*                               from Herschel/SPIRE

C  Format for file interpretation

    7 format(
     +  A14,1X,I2,1X,I2,1X,F7.4,1X,A1,1X,I2,1X,I2,1X,F7.4,1X,F10.7,1X,
     +  A5,1X,A1,1X,F6.2,1X,F5.2,1X,A1,1X,F6.2,1X,F4.2,1X,A1,1X,F5.2,
     +  1X,F4.2,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,91
        read(1,'(A109)')ar__6
        read(ar__6,7)
     +  Name_6,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,z_2,Sp_1,
     +  l_Fnii205,Fnii205,e_Fnii205,l_Fci371,Fci371,e_Fci371,l_Fci609,
     +  Fci609,e_Fci609,Ref_1
        if(ar__6(65:70) .EQ. '') Fnii205 = rNULL__
        if(ar__6(72:76) .EQ. '') e_Fnii205 = rNULL__
        if(ar__6(80:85) .EQ. '') Fci371 = rNULL__
        if(ar__6(87:90) .EQ. '') e_Fci371 = rNULL__
        if(ar__6(94:98) .EQ. '') Fci609 = rNULL__
        if(ar__6(100:103) .EQ. '') e_Fci609 = 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,7)
     +  Name_6,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,z_2,Sp_1,
     +  l_Fnii205,Fnii205,e_Fnii205,l_Fci371,Fci371,e_Fci371,l_Fci609,
     +  Fci609,e_Fci609,Ref_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 'table10.dat'	! Density determinations

C  Format for file interpretation

    8 format(
     +  A15,1X,A5,1X,F4.2,1X,F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
     +  A1,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1,1X,F4.2,1X,
     +  F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,132
        read(1,'(A130)')ar__7
        read(ar__7,8)
     +  Name_7,Sp_2,Rnii,e_Rnii,l_N1,N1,e_N1,E_N1_1,Rsiii,e_Rsiii,
     +  l_N2,N2,e_N2,E_N2_1,Roiii,e_Roiii,l_N3,N3,e_N3,E_N3_1,Rnev,
     +  e_Rnev,l_N4,N4,e_N4,E_N4_1
        if(ar__7(23:26) .EQ. '') Rnii = rNULL__
        if(ar__7(28:31) .EQ. '') e_Rnii = rNULL__
        if(ar__7(35:38) .EQ. '') N1 = rNULL__
        if(ar__7(40:43) .EQ. '') e_N1 = rNULL__
        if(ar__7(45:48) .EQ. '') E_N1_1 = rNULL__
        if(ar__7(50:54) .EQ. '') Rsiii = rNULL__
        if(ar__7(56:59) .EQ. '') e_Rsiii = rNULL__
        if(ar__7(63:66) .EQ. '') N2 = rNULL__
        if(ar__7(68:71) .EQ. '') e_N2 = rNULL__
        if(ar__7(73:76) .EQ. '') E_N2_1 = rNULL__
        if(ar__7(78:81) .EQ. '') Roiii = rNULL__
        if(ar__7(83:86) .EQ. '') e_Roiii = rNULL__
        if(ar__7(90:93) .EQ. '') N3 = rNULL__
        if(ar__7(95:98) .EQ. '') e_N3 = rNULL__
        if(ar__7(100:103) .EQ. '') E_N3_1 = rNULL__
        if(ar__7(105:108) .EQ. '') Rnev = rNULL__
        if(ar__7(110:113) .EQ. '') e_Rnev = rNULL__
        if(ar__7(117:120) .EQ. '') N4 = rNULL__
        if(ar__7(122:125) .EQ. '') e_N4 = rNULL__
        if(ar__7(127:130) .EQ. '') E_N4_1 = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  Name_7,Sp_2,Rnii,e_Rnii,l_N1,N1,e_N1,E_N1_1,Rsiii,e_Rsiii,
     +  l_N2,N2,e_N2,E_N2_1,Roiii,e_Roiii,l_N3,N3,e_N3,E_N3_1,Rnev,
     +  e_Rnev,l_N4,N4,e_N4,E_N4_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table11.dat'	! Far-IR fine structure line fluxes from
*                               Herschel/PACS for the central spaxel, 3x3 spaxel
*                               array, and 5x5 spaxel array

C  Format for file interpretation

    9 format(
     +  A1,1X,I3,1X,A15,1X,A1,1X,F7.2,1X,F6.2,1X,A1,1X,F7.2,1X,F6.2,
     +  1X,A1,1X,F8.2,1X,F6.2,1X,A1,1X,F7.2,1X,F6.2,1X,A1,1X,F7.2,1X,
     +  F5.2,1X,A1,1X,F7.2,1X,F6.2,1X,A1,1X,F8.2,1X,F6.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table11.dat')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,594
        read(1,'(A145)')ar__8
        read(ar__8,9)
     +  Type_6,Seq_6,Name_8,l_Foiii52_1,Foiii52_1,e_Foiii52_1,
     +  l_Fniii57_1,Fniii57_1,e_Fniii57_1,l_Foi63_1,Foi63_1,e_Foi63_1,
     +  l_Foiii88_1,Foiii88_1,e_Foiii88_1,l_Fnii122_1,Fnii122_1,
     +  e_Fnii122_1,l_Foi145_1,Foi145_1,e_Foi145_1,l_Fcii158_1,
     +  Fcii158_1,e_Fcii158_1,Aper_4
        if(ar__8(25:31) .EQ. '') Foiii52_1 = rNULL__
        if(ar__8(33:38) .EQ. '') e_Foiii52_1 = rNULL__
        if(ar__8(42:48) .EQ. '') Fniii57_1 = rNULL__
        if(ar__8(50:55) .EQ. '') e_Fniii57_1 = rNULL__
        if(ar__8(59:66) .EQ. '') Foi63_1 = rNULL__
        if(ar__8(68:73) .EQ. '') e_Foi63_1 = rNULL__
        if(ar__8(77:83) .EQ. '') Foiii88_1 = rNULL__
        if(ar__8(85:90) .EQ. '') e_Foiii88_1 = rNULL__
        if(ar__8(94:100) .EQ. '') Fnii122_1 = rNULL__
        if(ar__8(102:106) .EQ. '') e_Fnii122_1 = rNULL__
        if(ar__8(110:116) .EQ. '') Foi145_1 = rNULL__
        if(ar__8(118:123) .EQ. '') e_Foi145_1 = rNULL__
        if(ar__8(127:134) .EQ. '') Fcii158_1 = rNULL__
        if(ar__8(136:141) .EQ. '') e_Fcii158_1 = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  Type_6,Seq_6,Name_8,l_Foiii52_1,Foiii52_1,e_Foiii52_1,
     +  l_Fniii57_1,Fniii57_1,e_Fniii57_1,l_Foi63_1,Foi63_1,e_Foi63_1,
     +  l_Foiii88_1,Foiii88_1,e_Foiii88_1,l_Fnii122_1,Fnii122_1,
     +  e_Fnii122_1,l_Foi145_1,Foi145_1,e_Foi145_1,l_Fcii158_1,
     +  Fcii158_1,e_Fcii158_1,Aper_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

   10 format(A3,1X,A19,1X,A28,1X,A52)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,179
        read(1,'(A105)')ar__9
        read(ar__9,10)Ref_2,BibCode,Auth,Comm
c    ..............Just test output...........
        write(6,10)Ref_2,BibCode,Auth,Comm
c    .......End.of.Just test output...........
      end do
      close(1)

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