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 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 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      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*1   Type       (nr__) ! Sample type (G1)
      integer*4     Seq        (nr__) ! [1/170] Index number within Type
      character*15  Name       (nr__) ! Galaxy Name
      integer*4     RAh        (nr__) ! (h) 2MASS Hour of Right Ascension (J2000)
      integer*4     RAm        (nr__) ! (min) 2MASS Minute of Right Ascension (J2000)
      real*8        RAs        (nr__) ! (s) 2MASS Second of Right Ascension (J2000)
      character*1   DE_        (nr__) ! 2MASS Sign of the Declination (J2000)
      integer*4     DEd        (nr__) ! (deg) 2MASS Degree of Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) 2MASS Arcminute of Declination (J2000)
      real*8        DEs        (nr__) ! (arcsec) 2MASS Arcsecond of Declination (J2000)
      real*8        z          (nr__) ! [-0.0008/0.5] Redshift from NED
      real*8        Dist       (nr__) ! (Mpc) [0.6/1360]? Redshift-independent distance
      integer*4     r_Dist     (nr__) ! [1/36]? Reference for Dist
*                                     (see refs.dat file)
      character*5   Sp         (nr__) ! Spectral Classification
      character*1   l_logLx    (nr__) ! Limit flag on logLx
      real*4        logLx      (nr__) ! ([10-7W]) [35.9/45.7]? Log absorption-corrected
*                                     2-10keV luminosity; erg/s
      character*1   CT         (nr__) ! Compton Thick? Y: yes; N: no
      integer*4     r_logLx    (nr__) ! [101/171]? Reference for logLx
*                                     (see refs.dat file)
      character*1   l_logLnev1 (nr__) ! Limit flag on logLnev1
      real*4        logLnev1   (nr__) ! ([10-7W]) [36.4/43.1]? Log Spitzer/IRS NeV 14.3um
*                                     luminosity; erg/s
      character*1   l_logLnev2 (nr__) ! Limit flag on logLnev2
      real*4        logLnev2   (nr__) ! ([10-7W]) [37.1/43.1]? Log Spitzer/IRS NeV 24.3um
*                                     luminosity; erg/s
      character*1   l_logLoiv  (nr__) ! Limit flag on logLoiv
      real*4        logLoiv    (nr__) ! ([10-7W]) [36.6/43.8]? Log Spitzer/IRS OIV 25.9um
*                                     luminosity; erg/s
      real*4        Z_1        (nr__) ! (Sun) [0.02/2.8]? Metallicity derived from optical
*                                     spectra
      real*4        e_Z        (nr__) ! (Sun) [0.001/8]? Uncertainty in Z
      integer*4     r_Z        (nr__) ! ? Reference for Z (see refs.dat file)
      character*18  AName      (nr__) ! 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     (nr__1) ! Sample type (G1)
      integer*4     Seq_1      (nr__1) ! [1/170] Index number within Type
      character*15  Name_1     (nr__1) ! Galaxy Name
      character*1   l_Foiii52  (nr__1) ! Limit flag on Foiii52
      real*8        Foiii52    (nr__1) ! (10-17W/m2) [1/1363]? Herschel/PACS [OIII]52um line flux
      real*4        e_Foiii52  (nr__1) ! (10-17W/m2) [0.6/184]? Uncertainty in Foiii52
      character*1   l_Fniii57  (nr__1) ! Limit flag on Fniii57
      real*8        Fniii57    (nr__1) ! (10-17W/m2) [0.9/3288]? Herschel/PACS [NIII]57um
*                                     line flux
      real*4        e_Fniii57  (nr__1) ! (10-17W/m2) [0.4/179]? Uncertainty in Fniii57
      character*1   l_Foi63    (nr__1) ! Limit flag on Foi63
      real*8        Foi63      (nr__1) ! (10-17W/m2) [0.5/13352]? Herschel/PACS [OI]63um
*                                     line flux
      real*4        e_Foi63    (nr__1) ! (10-17W/m2) [0.2/404]? Uncertainty in Foi63
      character*1   l_Foiii88  (nr__1) ! Limit flag on Foiii88
      real*8        Foiii88    (nr__1) ! (10-17W/m2) [0.4/7653]? Herschel/PACS [OIII]88um
*                                     line flux
      real*4        e_Foiii88  (nr__1) ! (10-17W/m2) [0.2/236]? Uncertainty in Foiii88
      character*1   l_Fnii122  (nr__1) ! Limit flag on Fnii1222
      real*8        Fnii122    (nr__1) ! (10-17W/m2) [0.2/2002]? Herschel/PACS [NII]122um
*                                     line flux
      real*4        e_Fnii122  (nr__1) ! (10-17W/m2) [0.06/50]? Uncertainty in Fnii122
      character*1   l_Foi145   (nr__1) ! Limit flag on Foi145
      real*8        Foi145     (nr__1) ! (10-17W/m2) [0.1/1607]? Herschel/PACS [OI]145um
*                                     line flux
      real*4        e_Foi145   (nr__1) ! (10-17W/m2) [0.05/328]? Uncertainty in Foi145
      character*1   l_Fcii158  (nr__1) ! Limit flag on Fcii158
      real*8        Fcii158    (nr__1) ! (10-17W/m2) [0.5/12875]? Herschel/PACS [CII]158um
*                                     line flux
      real*4        e_Fcii158  (nr__1) ! (10-17W/m2) [0.08/345]? Uncertainty in Fcii158
      character*3   Aper       (nr__1) ! 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     (nr__2) ! Sample type (G1)
      integer*4     Seq_2      (nr__2) ! [1/170] Index number within Type
      character*15  Name_2     (nr__2) ! Galaxy Name
      character*1   l_Woiii52  (nr__2) ! Limit flag on Woiii52
      real*4        Woiii52    (nr__2) ! (um) [0.005/0.6]? Herschel/PACS [OIII]52um line
*                               equivalent width
      real*4        e_Woiii52  (nr__2) ! (um) [0.003/0.6]? Uncertainty in Woiii52
      character*1   l_Wniii57  (nr__2) ! Limit flag on Wniii57
      real*4        Wniii57    (nr__2) ! (um) [0.003/0.1]? Herschel/PACS [NIII]57um line
*                               equivalent width
      real*4        e_Wniii57  (nr__2) ! (um) [0.001/0.02]? Uncertainty in Wniii57
      character*1   l_Woi63    (nr__2) ! Limit flag on Woi63
      real*4        Woi63      (nr__2) ! (um) [0.002/1.5]? Herschel/PACS [OI]63um line
*                               equivalent width
      real*4        e_Woi63    (nr__2) ! (um) [0.001/1.4]? Uncertainty in Woi63
      character*1   l_Woiii88  (nr__2) ! Limit flag on Woiii88
      real*4        Woiii88    (nr__2) ! (um) [0.004/2]? Herschel/PACS [OIII]88um line
*                               equivalent width
      real*4        e_Woiii88  (nr__2) ! (um) [0.001/2]? Uncertainty in Woiii88
      character*1   l_Wnii122  (nr__2) ! Limit flag on Wnii122
      real*4        Wnii122    (nr__2) ! (um) [0.004/1.2]? Herschel/PACS [NII]122um line
*                               equivalent width
      real*4        e_Wnii122  (nr__2) ! (um) [0.001/0.08]? Uncertainty in Wnii122
      character*1   l_Woi145   (nr__2) ! Limit flag on Woi145
      real*4        Woi145     (nr__2) ! (um) [0.005/1.1]? Herschel/PACS [OI]145um line
*                               equivalent width
      real*4        e_Woi145   (nr__2) ! (um) [0.001/0.07]? Uncertainty in Woi145
      character*1   l_Wcii158  (nr__2) ! Limit flag on Wcii158
      real*4        Wcii158    (nr__2) ! (um) [0.09/5.3]? Herschel/PACS [CII]158um line
*                               equivalent width
      real*4        e_Wcii158  (nr__2) ! (um) [0.003/2]? Uncertainty in Wcii158
      character*3   Aper_1     (nr__2) ! 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     (nr__3) ! Sample type (G1)
      integer*4     Seq_3      (nr__3) ! [1/170] Index number within Type
      character*15  Name_3     (nr__3) ! Galaxy Name
      real*4        Coiii52    (nr__3) ! (Jy) [0.3/451]? Herschel/PACS [OIII]52um line
*                                   adjacent continuum flux
      real*4        e_Coiii52  (nr__3) ! (Jy) [0.04/15]? Uncertainty in Coiii52
      real*8        Cniii57    (nr__3) ! (Jy) [0.3/1039]? Herschel/PACS [NIII]57um line
*                                   adjacent continuum flux
      real*4        e_Cniii57  (nr__3) ! (Jy) [0.04/9]? Uncertainty in Cniii57
      character*1   l_Coi63    (nr__3) ! Limit flag on Coi63
      real*8        Coi63      (nr__3) ! (Jy) [0.04/1132]? Herschel/PACS [OI]63um line
*                                   adjacent continuum flux
      real*4        e_Coi63    (nr__3) ! (Jy) [0.04/16]? Uncertainty in Coi63
      character*1   l_Coiii88  (nr__3) ! Limit flag on Coiii88
      real*8        Coiii88    (nr__3) ! (Jy) [0.07/1207]? Herschel/PACS [OIII]88um line
*                                   adjacent continuum flux
      real*4        e_Coiii88  (nr__3) ! (Jy) [0.02/8]? Uncertainty Coiii88
      character*1   l_Cnii122  (nr__3) ! Limit flag on Cnii122
      real*8        Cnii122    (nr__3) ! (Jy) [0.06/1074]? Herschel/PACS [NII]122um line
*                                   adjacent continuum flux
      real*4        e_Cnii122  (nr__3) ! (Jy) [0.01/4]? Uncertainty in Cnii122
      character*1   l_Coi145   (nr__3) ! Limit flag on Coi145
      real*4        Coi145     (nr__3) ! (Jy) [0.05/934]? Herschel/PACS [OI]145um line
*                                   adjacent continuum flux
      real*4        e_Coi145   (nr__3) ! (Jy) [0.01/5]? Uncertainty in Coi145
      character*1   l_Ccii158  (nr__3) ! Limit flag on Ccii158
      real*4        Ccii158    (nr__3) ! (Jy) [0.07/798]? Herschel/PACS [CII]158um line
*                                   adjacent continuum flux
      real*4        e_Ccii158  (nr__3) ! (Jy) [0.02/36]? Uncertainty in Ccii158
      character*1   r_Ccii158  (nr__3) ! Reference for Ccii158 (see refs.dat file)
      character*3   Aper_2     (nr__3) ! 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     (nr__4) ! Sample type (G1)
      integer*4     Seq_4      (nr__4) ! [1/170] Index number within Type
      character*15  Name_4     (nr__4) ! Galaxy Name
      real*4        Soiii52    (nr__4) ! (km/s) [124/775]? Herschel/PACS [OIII]52um line
*                                 Gaussian sigma
      real*4        e_Soiii52  (nr__4) ! (km/s) [7/182]? Uncertainty in error for [OIII]52
      real*4        Sniii57    (nr__4) ! (km/s) [81/390]? Herschel/PACS [NIII]57um line
*                                 Gaussian sigma
      real*4        e_Sniii57  (nr__4) ! (km/s) [0.2/153]? Uncertainty in error for [NIII]57
      real*4        Soi63      (nr__4) ! (km/s) [60/633]? Herschel/PACS [OI]63um line
*                                 Gaussian sigma
      real*4        e_Soi63    (nr__4) ! (km/s) [0.1/66]? Uncertainty in error for [OI]63
      real*4        Soiii88    (nr__4) ! (km/s) [73/421]? Herschel/PACS [OIII]88um line
*                                 Gaussian sigma
      real*4        e_Soiii88  (nr__4) ! (km/s) [0.2/112]? Uncertainty in error for [OIII]88
      real*4        Snii122    (nr__4) ! (km/s) [102/911]? Herschel/PACS [NII]122um line
*                                 Gaussian sigma
      real*4        e_Snii122  (nr__4) ! (km/s) [0.4/135]? Uncertainty in error for [NII]122
      real*4        Soi145     (nr__4) ! (km/s) [108/344]? Herschel/PACS [OI]145um line
*                                 Gaussian sigma
      real*4        e_Soi145   (nr__4) ! (km/s) [0.5/614]? Uncertainty in error for [OI]145
      real*4        Scii158    (nr__4) ! (km/s) [89/340]? Herschel/PACS [CII]158um line
*                                 Gaussian sigma
      real*4        e_Scii158  (nr__4) ! (km/s) [0.2/100]? Uncertainty in error for [CII]158
      character*3   Aper_3     (nr__4) ! 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     (nr__5) ! Sample type (G1)
      integer*4     Seq_5      (nr__5) ! [1/170] Index number within Type
      character*15  Name_5     (nr__5) ! Galaxy Name
      character*1   l_Fsiv10   (nr__5) ! Limit flag on Fsiv10
      real*4        Fsiv10     (nr__5) ! (10-17W/m2) [0.07/537]? Spitzer/IRS [SIV]10.5um
*                                       line flux
      real*4        e_Fsiv10   (nr__5) ! (10-17W/m2) [0/54]? Uncertainty in Fsiv10
      character*1   l_Fneii12  (nr__5) ! Limit flag on Fneii12
      real*8        Fneii12    (nr__5) ! (10-17W/m2) [0.2/2833]? Spitzer/IRS [NeII]12.8um
*                                       line flux
      real*4        e_Fneii12  (nr__5) ! (10-17W/m2) [0.02/68]? Uncertainty in Fneii12
      character*1   l_Fnev14   (nr__5) ! Limit flag on Fnev14
      real*4        Fnev14     (nr__5) ! (10-17W/m2) [0.03/899]? Spitzer/IRS [NeV]14.3um
*                                       line flux
      real*4        e_Fnev14   (nr__5) ! (10-17W/m2) [0.01/23]? Uncertainty in Fnev14
      character*1   l_Fneiii15 (nr__5) ! Limit flag on Fneiii15
      real*8        Fneiii15   (nr__5) ! (10-17W/m2) [0.08/1371]? Spitzer/IRS [NeIII]15.6um
*                                       line flux
      real*4        e_Fneiii15 (nr__5) ! (10-17W/m2) [0.01/30]? Uncertainty in Fneiii15
      character*1   l_Fsiii18  (nr__5) ! Limit flag on Fsiii18
      real*4        Fsiii18    (nr__5) ! (10-17W/m2) [0.1/667]? Spitzer/IRS [SIII]18.7um
*                                       line flux
      real*4        e_Fsiii18  (nr__5) ! (10-17W/m2) [0/27]? Uncertainty in Fsiii18
      character*1   l_Fnev24   (nr__5) ! Limit flag on Fnev24
      real*4        Fnev24     (nr__5) ! (10-17W/m2) [0.07/816]? Spitzer/IRS [NeV]24.3um
*                                       line flux
      real*4        e_Fnev24   (nr__5) ! (10-17W/m2) [0.04/59]? Uncertainty in Fnev24
      character*1   l_Foiv25   (nr__5) ! Limit flag on Foiv25
      real*8        Foiv25     (nr__5) ! (10-17W/m2) [0.07/2030]? Spitzer/IRS [OIV]25.9um
*                                       line flux
      real*4        e_Foiv25   (nr__5) ! (10-17W/m2) [0.01/49]? Uncertainty in Foiv25
      character*1   l_Fsiii33  (nr__5) ! Limit flag on Fsiii33
      real*8        Fsiii33    (nr__5) ! (10-17W/m2) [0.09/1813]? Spitzer/IRS [SIII]33.5um
*                                       line flux
      real*4        e_Fsiii33  (nr__5) ! (10-17W/m2) [0.08/102]? Uncertainty in Fsiii33
      character*1   l_Fsii34   (nr__5) ! Limit flag on Fsii34
      real*8        Fsii34     (nr__5) ! (10-17W/m2) [1/2412]? Spitzer/IRS [SiII]34.8um
*                                       line flux
      real*4        e_Fsii34   (nr__5) ! (10-17W/m2) [0.3/81]? Uncertainty in Fsii34
      character*5   Mod        (nr__5) ! IRS Mode (1)
      character*23  Ref        (nr__5) ! Reference for the Spitzer spectra
*                                       (see refs.dat file)
      character*1   f_Name     (nr__5) ! 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    (nr__6) ! (deg) Right Ascension J2000
      real*8        DEdeg_1    (nr__6) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*14  Name_6     (nr__6) ! Galaxy Name
      integer*4     RAh_1      (nr__6) ! (h) 2MASS Hour of Right Ascension (J2000)
      integer*4     RAm_1      (nr__6) ! (min) 2MASS Minute of Right Ascension (J2000)
      real*8        RAs_1      (nr__6) ! (s) 2MASS Second of Right Ascension (J2000)
      character*1   DE__1      (nr__6) ! 2MASS Sign of the Declination (J2000)
      integer*4     DEd_1      (nr__6) ! (deg) 2MASS Degree of Declination (J2000)
      integer*4     DEm_1      (nr__6) ! (arcmin) 2MASS Arcminute of Declination (J2000)
      real*8        DEs_1      (nr__6) ! (arcsec) 2MASS Arcsecond of Declination (J2000)
      real*8        z_2        (nr__6) ! [-0.0008/0.2] Redshift from NED
      character*5   Sp_1       (nr__6) ! Spectral Classification
      character*1   l_Fnii205  (nr__6) ! Limit flag for Fnii205
      real*4        Fnii205    (nr__6) ! (10-17W/m2) [0.4/438]? Flux of the [NII]205um line
      real*4        e_Fnii205  (nr__6) ! (10-17W/m2) [0.1/12]? Uncertainty in Fnii205
      character*1   l_Fci371   (nr__6) ! Limit flag for Fci371
      real*4        Fci371     (nr__6) ! (10-17W/m2) [0.3/108]? Flux of the [CI]371um line
      real*4        e_Fci371   (nr__6) ! (10-17W/m2) [0/3]? Uncertainty in Fci371
      character*1   l_Fci609   (nr__6) ! Limit flag for Fci609
      real*4        Fci609     (nr__6) ! (10-17W/m2) [0.3/43]? Flux of the [CI]609um line
      real*4        e_Fci609   (nr__6) ! (10-17W/m2) [0/3]? Uncertainty in Fci609
      character*5   Ref_1      (nr__6) ! 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     (nr__7) ! Galaxy Name
      character*5   Sp_2       (nr__7) ! Spectral Classification
      real*4        Rnii       (nr__7) ! [0.2/3]? Line ratio [NII] 205/122
      real*4        e_Rnii     (nr__7) ! [0.01/2]? Uncertainty in Rnii
      character*1   l_N1       (nr__7) ! Limit flag on N1
      real*4        N1         (nr__7) ! ([cm-3]) [0.1/3]? Log density derived
*                                 from [NII] 205/122 ratio
      real*4        e_N1       (nr__7) ! ([cm-3]) [0/3]? Lower error boundary in N1
      real*4        E_N1_1     (nr__7) ! ([cm-3]) [1/4]? Upper error boundary in N1
      real*4        Rsiii      (nr__7) ! [0.6/15]? Line ratio [SIII] 33/18
      real*4        e_Rsiii    (nr__7) ! [0.04/6]? Uncertainy in Rsiii
      character*1   l_N2       (nr__7) ! Limit flag on N2
      real*4        N2         (nr__7) ! ([cm-3]) [0.3/4]? Log density derived
*                                 from [SIII] 33/18 ratio
      real*4        e_N2       (nr__7) ! ([cm-3]) [0/4]? Lower error boundary in N2
      real*4        E_N2_1     (nr__7) ! ([cm-3]) [2/4]? Upper error boundary in N2
      real*4        Roiii      (nr__7) ! [0.1/2.2]? Line ratio [OIII]88/52
      real*4        e_Roiii    (nr__7) ! [0.09/0.7]? Error on the [OIII]88/52 line ratio
      character*1   l_N3       (nr__7) ! Limit flag on N3
      real*4        N3         (nr__7) ! ([cm-3]) [0.7/5]? Log density derived
*                                 from [OIII] 88/52 ratio
      real*4        e_N3       (nr__7) ! ([cm-3]) [0/4]? Lower error boundary in N3
      real*4        E_N3_1     (nr__7) ! ([cm-3]) [1/6]? Upper error boundary in N3
      real*4        Rnev       (nr__7) ! [0.4/6]? Line ratio [NeV] 24/14
      real*4        e_Rnev     (nr__7) ! [0.02/2]? Error in Rnev
      character*1   l_N4       (nr__7) ! Limit flag on N4
      real*4        N4         (nr__7) ! ([cm-3]) [1/4]? Log density derived
*                                 from [NeV] 24/14 ratio
      real*4        e_N4       (nr__7) ! ([cm-3]) [0/4]? Lower error boundary in N4
      real*4        E_N4_1     (nr__7) ! ([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     (nr__8) ! Sample type (G1)
      integer*4     Seq_6      (nr__8) ! [1/170] Index number within Type
      character*15  Name_8     (nr__8) ! Galaxy Name
      character*1   l_Foiii52_1(nr__8) ! Limit flag on Foiii52
      real*8        Foiii52_1  (nr__8) ! (10-17W/m2) [1/1363]? Herschel/PACS [OIII]52um line flux
      real*4        e_Foiii52_1(nr__8) ! (10-17W/m2) [0.6/184]? Uncertainty in Foiii52
      character*1   l_Fniii57_1(nr__8) ! Limit flag on Fniii57
      real*8        Fniii57_1  (nr__8) ! (10-17W/m2) [0.9/3288]? Herschel/PACS [NIII]57um
*                                     line flux
      real*4        e_Fniii57_1(nr__8) ! (10-17W/m2) [0.4/179]? Uncertainty in Fniii57
      character*1   l_Foi63_1  (nr__8) ! Limit flag on Foi63
      real*8        Foi63_1    (nr__8) ! (10-17W/m2) [0.5/13352]? Herschel/PACS [OI]63um
*                                     line flux
      real*4        e_Foi63_1  (nr__8) ! (10-17W/m2) [0.2/404]? Uncertainty in Foi63
      character*1   l_Foiii88_1(nr__8) ! Limit flag on Foiii88
      real*8        Foiii88_1  (nr__8) ! (10-17W/m2) [0.4/7653]? Herschel/PACS [OIII]88um
*                                     line flux
      real*4        e_Foiii88_1(nr__8) ! (10-17W/m2) [0.2/236]? Uncertainty in Foiii88
      character*1   l_Fnii122_1(nr__8) ! Limit flag on Fnii1222
      real*8        Fnii122_1  (nr__8) ! (10-17W/m2) [0.2/2002]? Herschel/PACS [NII]122um
*                                     line flux
      real*4        e_Fnii122_1(nr__8) ! (10-17W/m2) [0.06/50]? Uncertainty in Fnii122
      character*1   l_Foi145_1 (nr__8) ! Limit flag on Foi145
      real*8        Foi145_1   (nr__8) ! (10-17W/m2) [0.1/1607]? Herschel/PACS [OI]145um
*                                     line flux
      real*4        e_Foi145_1 (nr__8) ! (10-17W/m2) [0.05/328]? Uncertainty in Foi145
      character*1   l_Fcii158_1(nr__8) ! Limit flag on Fcii158
      real*8        Fcii158_1  (nr__8) ! (10-17W/m2) [0.5/12875]? Herschel/PACS [CII]158um
*                                     line flux
      real*4        e_Fcii158_1(nr__8) ! (10-17W/m2) [0.08/345]? Uncertainty in Fcii158
      character*3   Aper_4     (nr__8) ! 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      (nr__9) ! Reference code
      character*19  BibCode    (nr__9) ! Bibcode
      character*28  Auth       (nr__9) ! First author's name(s)
      character*52  Comm       (nr__9) ! 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(i__),Seq(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),z(i__),Dist(i__),
     +  r_Dist(i__),Sp(i__),l_logLx(i__),logLx(i__),CT(i__),
     +  r_logLx(i__),l_logLnev1(i__),logLnev1(i__),l_logLnev2(i__),
     +  logLnev2(i__),l_logLoiv(i__),logLoiv(i__),Z_1(i__),e_Z(i__),
     +  r_Z(i__),AName(i__)
        if(ar__(63:70) .EQ. '') Dist(i__) = rNULL__
        if(ar__(72:73) .EQ. '') r_Dist(i__) = iNULL__
        if(ar__(83:87) .EQ. '') logLx(i__) = rNULL__
        if(ar__(91:93) .EQ. '') r_logLx(i__) = iNULL__
        if(ar__(97:101) .EQ. '') logLnev1(i__) = rNULL__
        if(ar__(105:109) .EQ. '') logLnev2(i__) = rNULL__
        if(ar__(113:117) .EQ. '') logLoiv(i__) = rNULL__
        if(ar__(119:123) .EQ. '') Z_1(i__) = rNULL__
        if(ar__(125:129) .EQ. '') e_Z(i__) = rNULL__
        if(ar__(131:133) .EQ. '') r_Z(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)
     +  Type(i__),Seq(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),z(i__),Dist(i__),
     +  r_Dist(i__),Sp(i__),l_logLx(i__),logLx(i__),CT(i__),
     +  r_logLx(i__),l_logLnev1(i__),logLnev1(i__),l_logLnev2(i__),
     +  logLnev2(i__),l_logLoiv(i__),logLoiv(i__),Z_1(i__),e_Z(i__),
     +  r_Z(i__),AName(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'	! 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(i__),Seq_1(i__),Name_1(i__),l_Foiii52(i__),
     +  Foiii52(i__),e_Foiii52(i__),l_Fniii57(i__),Fniii57(i__),
     +  e_Fniii57(i__),l_Foi63(i__),Foi63(i__),e_Foi63(i__),
     +  l_Foiii88(i__),Foiii88(i__),e_Foiii88(i__),l_Fnii122(i__),
     +  Fnii122(i__),e_Fnii122(i__),l_Foi145(i__),Foi145(i__),
     +  e_Foi145(i__),l_Fcii158(i__),Fcii158(i__),e_Fcii158(i__),
     +  Aper(i__)
        if(ar__1(25:31) .EQ. '') Foiii52(i__) = rNULL__
        if(ar__1(33:38) .EQ. '') e_Foiii52(i__) = rNULL__
        if(ar__1(42:48) .EQ. '') Fniii57(i__) = rNULL__
        if(ar__1(50:55) .EQ. '') e_Fniii57(i__) = rNULL__
        if(ar__1(59:66) .EQ. '') Foi63(i__) = rNULL__
        if(ar__1(68:73) .EQ. '') e_Foi63(i__) = rNULL__
        if(ar__1(77:83) .EQ. '') Foiii88(i__) = rNULL__
        if(ar__1(85:90) .EQ. '') e_Foiii88(i__) = rNULL__
        if(ar__1(94:100) .EQ. '') Fnii122(i__) = rNULL__
        if(ar__1(102:106) .EQ. '') e_Fnii122(i__) = rNULL__
        if(ar__1(110:116) .EQ. '') Foi145(i__) = rNULL__
        if(ar__1(118:123) .EQ. '') e_Foi145(i__) = rNULL__
        if(ar__1(127:134) .EQ. '') Fcii158(i__) = rNULL__
        if(ar__1(136:141) .EQ. '') e_Fcii158(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Type_1(i__),Seq_1(i__),Name_1(i__),l_Foiii52(i__),
     +  Foiii52(i__),e_Foiii52(i__),l_Fniii57(i__),Fniii57(i__),
     +  e_Fniii57(i__),l_Foi63(i__),Foi63(i__),e_Foi63(i__),
     +  l_Foiii88(i__),Foiii88(i__),e_Foiii88(i__),l_Fnii122(i__),
     +  Fnii122(i__),e_Fnii122(i__),l_Foi145(i__),Foi145(i__),
     +  e_Foi145(i__),l_Fcii158(i__),Fcii158(i__),e_Fcii158(i__),
     +  Aper(i__)
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(i__),Seq_2(i__),Name_2(i__),l_Woiii52(i__),
     +  Woiii52(i__),e_Woiii52(i__),l_Wniii57(i__),Wniii57(i__),
     +  e_Wniii57(i__),l_Woi63(i__),Woi63(i__),e_Woi63(i__),
     +  l_Woiii88(i__),Woiii88(i__),e_Woiii88(i__),l_Wnii122(i__),
     +  Wnii122(i__),e_Wnii122(i__),l_Woi145(i__),Woi145(i__),
     +  e_Woi145(i__),l_Wcii158(i__),Wcii158(i__),e_Wcii158(i__),
     +  Aper_1(i__)
        if(ar__2(25:29) .EQ. '') Woiii52(i__) = rNULL__
        if(ar__2(31:35) .EQ. '') e_Woiii52(i__) = rNULL__
        if(ar__2(39:43) .EQ. '') Wniii57(i__) = rNULL__
        if(ar__2(45:49) .EQ. '') e_Wniii57(i__) = rNULL__
        if(ar__2(53:57) .EQ. '') Woi63(i__) = rNULL__
        if(ar__2(59:63) .EQ. '') e_Woi63(i__) = rNULL__
        if(ar__2(67:71) .EQ. '') Woiii88(i__) = rNULL__
        if(ar__2(73:77) .EQ. '') e_Woiii88(i__) = rNULL__
        if(ar__2(81:85) .EQ. '') Wnii122(i__) = rNULL__
        if(ar__2(87:91) .EQ. '') e_Wnii122(i__) = rNULL__
        if(ar__2(95:99) .EQ. '') Woi145(i__) = rNULL__
        if(ar__2(101:105) .EQ. '') e_Woi145(i__) = rNULL__
        if(ar__2(109:113) .EQ. '') Wcii158(i__) = rNULL__
        if(ar__2(115:119) .EQ. '') e_Wcii158(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Type_2(i__),Seq_2(i__),Name_2(i__),l_Woiii52(i__),
     +  Woiii52(i__),e_Woiii52(i__),l_Wniii57(i__),Wniii57(i__),
     +  e_Wniii57(i__),l_Woi63(i__),Woi63(i__),e_Woi63(i__),
     +  l_Woiii88(i__),Woiii88(i__),e_Woiii88(i__),l_Wnii122(i__),
     +  Wnii122(i__),e_Wnii122(i__),l_Woi145(i__),Woi145(i__),
     +  e_Woi145(i__),l_Wcii158(i__),Wcii158(i__),e_Wcii158(i__),
     +  Aper_1(i__)
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(i__),Seq_3(i__),Name_3(i__),Coiii52(i__),
     +  e_Coiii52(i__),Cniii57(i__),e_Cniii57(i__),l_Coi63(i__),
     +  Coi63(i__),e_Coi63(i__),l_Coiii88(i__),Coiii88(i__),
     +  e_Coiii88(i__),l_Cnii122(i__),Cnii122(i__),e_Cnii122(i__),
     +  l_Coi145(i__),Coi145(i__),e_Coi145(i__),l_Ccii158(i__),
     +  Ccii158(i__),e_Ccii158(i__),r_Ccii158(i__),Aper_2(i__)
        if(ar__3(23:28) .EQ. '') Coiii52(i__) = rNULL__
        if(ar__3(30:34) .EQ. '') e_Coiii52(i__) = rNULL__
        if(ar__3(36:42) .EQ. '') Cniii57(i__) = rNULL__
        if(ar__3(44:47) .EQ. '') e_Cniii57(i__) = rNULL__
        if(ar__3(51:57) .EQ. '') Coi63(i__) = rNULL__
        if(ar__3(59:63) .EQ. '') e_Coi63(i__) = rNULL__
        if(ar__3(67:73) .EQ. '') Coiii88(i__) = rNULL__
        if(ar__3(75:78) .EQ. '') e_Coiii88(i__) = rNULL__
        if(ar__3(82:88) .EQ. '') Cnii122(i__) = rNULL__
        if(ar__3(90:93) .EQ. '') e_Cnii122(i__) = rNULL__
        if(ar__3(97:102) .EQ. '') Coi145(i__) = rNULL__
        if(ar__3(104:107) .EQ. '') e_Coi145(i__) = rNULL__
        if(ar__3(111:116) .EQ. '') Ccii158(i__) = rNULL__
        if(ar__3(118:122) .EQ. '') e_Ccii158(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Type_3(i__),Seq_3(i__),Name_3(i__),Coiii52(i__),
     +  e_Coiii52(i__),Cniii57(i__),e_Cniii57(i__),l_Coi63(i__),
     +  Coi63(i__),e_Coi63(i__),l_Coiii88(i__),Coiii88(i__),
     +  e_Coiii88(i__),l_Cnii122(i__),Cnii122(i__),e_Cnii122(i__),
     +  l_Coi145(i__),Coi145(i__),e_Coi145(i__),l_Ccii158(i__),
     +  Ccii158(i__),e_Ccii158(i__),r_Ccii158(i__),Aper_2(i__)
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(i__),Seq_4(i__),Name_4(i__),Soiii52(i__),
     +  e_Soiii52(i__),Sniii57(i__),e_Sniii57(i__),Soi63(i__),
     +  e_Soi63(i__),Soiii88(i__),e_Soiii88(i__),Snii122(i__),
     +  e_Snii122(i__),Soi145(i__),e_Soi145(i__),Scii158(i__),
     +  e_Scii158(i__),Aper_3(i__)
        if(ar__4(23:27) .EQ. '') Soiii52(i__) = rNULL__
        if(ar__4(29:33) .EQ. '') e_Soiii52(i__) = rNULL__
        if(ar__4(35:39) .EQ. '') Sniii57(i__) = rNULL__
        if(ar__4(41:45) .EQ. '') e_Sniii57(i__) = rNULL__
        if(ar__4(47:51) .EQ. '') Soi63(i__) = rNULL__
        if(ar__4(53:56) .EQ. '') e_Soi63(i__) = rNULL__
        if(ar__4(58:62) .EQ. '') Soiii88(i__) = rNULL__
        if(ar__4(64:68) .EQ. '') e_Soiii88(i__) = rNULL__
        if(ar__4(70:74) .EQ. '') Snii122(i__) = rNULL__
        if(ar__4(76:80) .EQ. '') e_Snii122(i__) = rNULL__
        if(ar__4(82:86) .EQ. '') Soi145(i__) = rNULL__
        if(ar__4(88:92) .EQ. '') e_Soi145(i__) = rNULL__
        if(ar__4(94:98) .EQ. '') Scii158(i__) = rNULL__
        if(ar__4(100:104) .EQ. '') e_Scii158(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Type_4(i__),Seq_4(i__),Name_4(i__),Soiii52(i__),
     +  e_Soiii52(i__),Sniii57(i__),e_Sniii57(i__),Soi63(i__),
     +  e_Soi63(i__),Soiii88(i__),e_Soiii88(i__),Snii122(i__),
     +  e_Snii122(i__),Soi145(i__),e_Soi145(i__),Scii158(i__),
     +  e_Scii158(i__),Aper_3(i__)
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(i__),Seq_5(i__),Name_5(i__),l_Fsiv10(i__),Fsiv10(i__),
     +  e_Fsiv10(i__),l_Fneii12(i__),Fneii12(i__),e_Fneii12(i__),
     +  l_Fnev14(i__),Fnev14(i__),e_Fnev14(i__),l_Fneiii15(i__),
     +  Fneiii15(i__),e_Fneiii15(i__),l_Fsiii18(i__),Fsiii18(i__),
     +  e_Fsiii18(i__),l_Fnev24(i__),Fnev24(i__),e_Fnev24(i__),
     +  l_Foiv25(i__),Foiv25(i__),e_Foiv25(i__),l_Fsiii33(i__),
     +  Fsiii33(i__),e_Fsiii33(i__),l_Fsii34(i__),Fsii34(i__),
     +  e_Fsii34(i__),Mod(i__),Ref(i__),f_Name(i__)
        if(ar__5(25:30) .EQ. '') Fsiv10(i__) = rNULL__
        if(ar__5(32:36) .EQ. '') e_Fsiv10(i__) = rNULL__
        if(ar__5(40:46) .EQ. '') Fneii12(i__) = rNULL__
        if(ar__5(48:52) .EQ. '') e_Fneii12(i__) = rNULL__
        if(ar__5(56:61) .EQ. '') Fnev14(i__) = rNULL__
        if(ar__5(63:67) .EQ. '') e_Fnev14(i__) = rNULL__
        if(ar__5(71:77) .EQ. '') Fneiii15(i__) = rNULL__
        if(ar__5(79:83) .EQ. '') e_Fneiii15(i__) = rNULL__
        if(ar__5(87:92) .EQ. '') Fsiii18(i__) = rNULL__
        if(ar__5(94:98) .EQ. '') e_Fsiii18(i__) = rNULL__
        if(ar__5(102:107) .EQ. '') Fnev24(i__) = rNULL__
        if(ar__5(109:113) .EQ. '') e_Fnev24(i__) = rNULL__
        if(ar__5(117:123) .EQ. '') Foiv25(i__) = rNULL__
        if(ar__5(125:129) .EQ. '') e_Foiv25(i__) = rNULL__
        if(ar__5(133:139) .EQ. '') Fsiii33(i__) = rNULL__
        if(ar__5(141:146) .EQ. '') e_Fsiii33(i__) = rNULL__
        if(ar__5(150:156) .EQ. '') Fsii34(i__) = rNULL__
        if(ar__5(158:162) .EQ. '') e_Fsii34(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  Type_5(i__),Seq_5(i__),Name_5(i__),l_Fsiv10(i__),Fsiv10(i__),
     +  e_Fsiv10(i__),l_Fneii12(i__),Fneii12(i__),e_Fneii12(i__),
     +  l_Fnev14(i__),Fnev14(i__),e_Fnev14(i__),l_Fneiii15(i__),
     +  Fneiii15(i__),e_Fneiii15(i__),l_Fsiii18(i__),Fsiii18(i__),
     +  e_Fsiii18(i__),l_Fnev24(i__),Fnev24(i__),e_Fnev24(i__),
     +  l_Foiv25(i__),Foiv25(i__),e_Foiv25(i__),l_Fsiii33(i__),
     +  Fsiii33(i__),e_Fsiii33(i__),l_Fsii34(i__),Fsii34(i__),
     +  e_Fsii34(i__),Mod(i__),Ref(i__),f_Name(i__)
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(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
     +  DEd_1(i__),DEm_1(i__),DEs_1(i__),z_2(i__),Sp_1(i__),
     +  l_Fnii205(i__),Fnii205(i__),e_Fnii205(i__),l_Fci371(i__),
     +  Fci371(i__),e_Fci371(i__),l_Fci609(i__),Fci609(i__),
     +  e_Fci609(i__),Ref_1(i__)
        if(ar__6(65:70) .EQ. '') Fnii205(i__) = rNULL__
        if(ar__6(72:76) .EQ. '') e_Fnii205(i__) = rNULL__
        if(ar__6(80:85) .EQ. '') Fci371(i__) = rNULL__
        if(ar__6(87:90) .EQ. '') e_Fci371(i__) = rNULL__
        if(ar__6(94:98) .EQ. '') Fci609(i__) = rNULL__
        if(ar__6(100:103) .EQ. '') e_Fci609(i__) = rNULL__
        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,7)
     +  Name_6(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
     +  DEd_1(i__),DEm_1(i__),DEs_1(i__),z_2(i__),Sp_1(i__),
     +  l_Fnii205(i__),Fnii205(i__),e_Fnii205(i__),l_Fci371(i__),
     +  Fci371(i__),e_Fci371(i__),l_Fci609(i__),Fci609(i__),
     +  e_Fci609(i__),Ref_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
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(i__),Sp_2(i__),Rnii(i__),e_Rnii(i__),l_N1(i__),N1(i__),
     +  e_N1(i__),E_N1_1(i__),Rsiii(i__),e_Rsiii(i__),l_N2(i__),
     +  N2(i__),e_N2(i__),E_N2_1(i__),Roiii(i__),e_Roiii(i__),
     +  l_N3(i__),N3(i__),e_N3(i__),E_N3_1(i__),Rnev(i__),e_Rnev(i__),
     +  l_N4(i__),N4(i__),e_N4(i__),E_N4_1(i__)
        if(ar__7(23:26) .EQ. '') Rnii(i__) = rNULL__
        if(ar__7(28:31) .EQ. '') e_Rnii(i__) = rNULL__
        if(ar__7(35:38) .EQ. '') N1(i__) = rNULL__
        if(ar__7(40:43) .EQ. '') e_N1(i__) = rNULL__
        if(ar__7(45:48) .EQ. '') E_N1_1(i__) = rNULL__
        if(ar__7(50:54) .EQ. '') Rsiii(i__) = rNULL__
        if(ar__7(56:59) .EQ. '') e_Rsiii(i__) = rNULL__
        if(ar__7(63:66) .EQ. '') N2(i__) = rNULL__
        if(ar__7(68:71) .EQ. '') e_N2(i__) = rNULL__
        if(ar__7(73:76) .EQ. '') E_N2_1(i__) = rNULL__
        if(ar__7(78:81) .EQ. '') Roiii(i__) = rNULL__
        if(ar__7(83:86) .EQ. '') e_Roiii(i__) = rNULL__
        if(ar__7(90:93) .EQ. '') N3(i__) = rNULL__
        if(ar__7(95:98) .EQ. '') e_N3(i__) = rNULL__
        if(ar__7(100:103) .EQ. '') E_N3_1(i__) = rNULL__
        if(ar__7(105:108) .EQ. '') Rnev(i__) = rNULL__
        if(ar__7(110:113) .EQ. '') e_Rnev(i__) = rNULL__
        if(ar__7(117:120) .EQ. '') N4(i__) = rNULL__
        if(ar__7(122:125) .EQ. '') e_N4(i__) = rNULL__
        if(ar__7(127:130) .EQ. '') E_N4_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  Name_7(i__),Sp_2(i__),Rnii(i__),e_Rnii(i__),l_N1(i__),N1(i__),
     +  e_N1(i__),E_N1_1(i__),Rsiii(i__),e_Rsiii(i__),l_N2(i__),
     +  N2(i__),e_N2(i__),E_N2_1(i__),Roiii(i__),e_Roiii(i__),
     +  l_N3(i__),N3(i__),e_N3(i__),E_N3_1(i__),Rnev(i__),e_Rnev(i__),
     +  l_N4(i__),N4(i__),e_N4(i__),E_N4_1(i__)
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(i__),Seq_6(i__),Name_8(i__),l_Foiii52_1(i__),
     +  Foiii52_1(i__),e_Foiii52_1(i__),l_Fniii57_1(i__),
     +  Fniii57_1(i__),e_Fniii57_1(i__),l_Foi63_1(i__),Foi63_1(i__),
     +  e_Foi63_1(i__),l_Foiii88_1(i__),Foiii88_1(i__),
     +  e_Foiii88_1(i__),l_Fnii122_1(i__),Fnii122_1(i__),
     +  e_Fnii122_1(i__),l_Foi145_1(i__),Foi145_1(i__),
     +  e_Foi145_1(i__),l_Fcii158_1(i__),Fcii158_1(i__),
     +  e_Fcii158_1(i__),Aper_4(i__)
        if(ar__8(25:31) .EQ. '') Foiii52_1(i__) = rNULL__
        if(ar__8(33:38) .EQ. '') e_Foiii52_1(i__) = rNULL__
        if(ar__8(42:48) .EQ. '') Fniii57_1(i__) = rNULL__
        if(ar__8(50:55) .EQ. '') e_Fniii57_1(i__) = rNULL__
        if(ar__8(59:66) .EQ. '') Foi63_1(i__) = rNULL__
        if(ar__8(68:73) .EQ. '') e_Foi63_1(i__) = rNULL__
        if(ar__8(77:83) .EQ. '') Foiii88_1(i__) = rNULL__
        if(ar__8(85:90) .EQ. '') e_Foiii88_1(i__) = rNULL__
        if(ar__8(94:100) .EQ. '') Fnii122_1(i__) = rNULL__
        if(ar__8(102:106) .EQ. '') e_Fnii122_1(i__) = rNULL__
        if(ar__8(110:116) .EQ. '') Foi145_1(i__) = rNULL__
        if(ar__8(118:123) .EQ. '') e_Foi145_1(i__) = rNULL__
        if(ar__8(127:134) .EQ. '') Fcii158_1(i__) = rNULL__
        if(ar__8(136:141) .EQ. '') e_Fcii158_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  Type_6(i__),Seq_6(i__),Name_8(i__),l_Foiii52_1(i__),
     +  Foiii52_1(i__),e_Foiii52_1(i__),l_Fniii57_1(i__),
     +  Fniii57_1(i__),e_Fniii57_1(i__),l_Foi63_1(i__),Foi63_1(i__),
     +  e_Foi63_1(i__),l_Foiii88_1(i__),Foiii88_1(i__),
     +  e_Foiii88_1(i__),l_Fnii122_1(i__),Fnii122_1(i__),
     +  e_Fnii122_1(i__),l_Foi145_1(i__),Foi145_1(i__),
     +  e_Foi145_1(i__),l_Fcii158_1(i__),Fcii158_1(i__),
     +  e_Fcii158_1(i__),Aper_4(i__)
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(i__),BibCode(i__),Auth(i__),Comm(i__)
c    ..............Just test output...........
        write(6,10)Ref_2(i__),BibCode(i__),Auth(i__),Comm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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