FORTRAN Generation
(/./ftp/cats/J/ApJS/154/673)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/154/673 into FORTRAN code for loading all data files into arrays.

Note that special values are assigned to unknown or unspecified numbers (also called NULL numbers); when necessary, the coordinate components making up the right ascension and declination are converted into floating-point numbers representing these angles in degrees.



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-19
C=============================================================================
*  This code was generated from the ReadMe file documenting a catalogue
*  according to the "Standard for Documentation of Astronomical Catalogues"
*  currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
*  (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
*  Please report problems or questions to   
C=============================================================================

      implicit none
*  Unspecified or NULL values, generally corresponding to blank columns,
*  are assigned one of the following special values:
*     rNULL__    for unknown or NULL floating-point values
*     iNULL__    for unknown or NULL   integer      values
      real*4     rNULL__
      integer*4  iNULL__
      parameter  (rNULL__=--2147483648.)  	! NULL real number
      parameter  (iNULL__=(-2147483647-1))	! NULL int  number
      integer    idig			! testing NULL number

C=============================================================================
Cat. J/ApJS/154/673      COBE DIRBE Point Source Catalog          (Smith+, 2004)
*================================================================================
*The COBE DIRBE Point Source Catalog.
*    Smith B.J., Price S.D., Baker R.I.
*   <Astrophys. J. Suppl. Ser., 154, 673 (2004)>
*   =2004ApJS..154..673S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! DIRBE Photometry and Variability

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

      character*17  DIRBE      (nr__) ! DIRBE name (G1)
      character*16  v2MASS     (nr__) ! 2MASS name (2)
      character*23  IRAS       (nr__) ! IRAS/MSX name (3)
      real*8        F1_25      (nr__) ! (Jy) ?=-99.9 1.25{mu}m weighted mean flux density (4)
      real*8        e_F1_25    (nr__) ! (Jy) ?=-99.9 RMS error of F1.25
      real*8        eF1_25     (nr__) ! (Jy) ?=-99.9 Average 1.25{mu}m light curve errorbar
      real*4        D1_25      (nr__) ! (mag) ?=-99.9 1.25 light curve var. amplitude (5)
      real*4        e_D1_25    (nr__) ! (mag) ?=-99.90 Uncertainty in the D1.25
      integer*4     o_F1_25    (nr__) ! Number of 1.25{mu}m unfiltered datapoints
      real*8        F2_2       (nr__) ! (Jy) 2.2{mu}m weighted mean flux density (4)
      real*8        e_F2_2     (nr__) ! (Jy) ?=-99.9 RMS error of F2.2
      real*8        eF2_2      (nr__) ! (Jy) ?=-99.9 Average 2.2{mu}m light curve errorbar
      real*4        D2_2       (nr__) ! (mag) ?=-99.9 2.2um light curve variation amplitude (5)
      real*4        e_D2_2     (nr__) ! (mag) ?=-99.90 Uncertainty in the D2.2
      integer*4     o_F2_2     (nr__) ! Number of 2.2{mu}m unfiltered datapoints
      real*8        F3_5       (nr__) ! (Jy) ?=-99.9 3.5{mu}m weighted mean flux density (4)
      real*8        e_F3_5     (nr__) ! (Jy) ?=-99.9 RMS error of F3.5
      real*8        eF3_5      (nr__) ! (Jy) ?=-99.9 Average 3.5{mu}m light curve errorbar
      real*4        D3_5       (nr__) ! (mag) ?=-99.9 3.5um light curve variation amplitude (5)
      real*4        e_D3_5     (nr__) ! (mag) ?=-99.90 Uncertainty in the D3.5
      integer*4     o_F3_5     (nr__) ! Number of 3.5 {mu}m unfiltered datapoints
      real*8        F4_9       (nr__) ! (Jy) ?=-99.9 4.9{mu}m weighted mean flux density (4)
      real*8        e_F4_9     (nr__) ! (Jy) ?=-99.9 RMS error of F4.9
      real*8        eF4_9      (nr__) ! (Jy) ?=-99.9 Average 4.9{mu}m light curve errorbar
      real*4        D4_9       (nr__) ! (mag) ?=-99.9 4.9um light curve variation amplitude (5)
      real*4        e_D4_9     (nr__) ! (mag) ?=-99.90 Uncertainty in the D4.9
      integer*4     o_F4_9     (nr__) ! Number of 4.9{mu}m unfiltered datapoints
      real*8        F12        (nr__) ! (Jy) ?=-99.9 12{mu}m weighted mean flux density (4)
      real*8        e_F12      (nr__) ! (Jy) ?=-99.9 RMS error of F12
      real*8        eF12       (nr__) ! (Jy) ?=-99.9 Average 12{mu}m light curve errorbar
      real*4        D12        (nr__) ! (mag) ?=-99.9 12um light curve variation amplitude (5)
      real*4        e_D12      (nr__) ! (mag) ?=-99.90 Uncertainty in the D12
      integer*4     o_F12      (nr__) ! Number of 12{mu}m unfiltered datapoints
      real*8        F25        (nr__) ! (Jy) ?=-99.9 25{mu}m weighted mean flux density (4)
      real*8        e_F25      (nr__) ! (Jy) ?=-99.9 RMS error of F25
      real*8        eF25       (nr__) ! (Jy) ?=-99.9 Average 25{mu}m light curve errorbar
      real*4        D25        (nr__) ! (mag) ?=-99.9 25um light curve variation amplitude (5)
      real*4        e_D25      (nr__) ! (mag) ?=-99.90 Uncertainty in the D25
      integer*4     o_F25      (nr__) ! Number of 25{mu}m unfiltered datapoints
      real*8        F60        (nr__) ! (Jy) ?=-99.9 60{mu}m weighted mean flux density (4)
      real*8        e_F60      (nr__) ! (Jy) ?=-99.9 RMS error of F60
      real*8        eF60       (nr__) ! (Jy) ?=-99.9 Average 60{mu}m light curve errorbar
      real*4        D60        (nr__) ! (mag) ?=-99.9 60um light curve variation amplitude (5)
      real*4        e_D60      (nr__) ! (mag) ?=-99.90 Uncertainty in the D60
      integer*4     o_F60      (nr__) ! Number of 60{mu}m unfiltered datapoints
      real*8        F100       (nr__) ! (Jy) ?=-99.9 100{mu}m weighted mean flux density (4)
      real*8        e_F100     (nr__) ! (Jy) ?=-99.9 RMS error of F100
      real*8        eF100      (nr__) ! (Jy) ?=-99.9 Average 100{mu}m light curve errorbar
      real*4        D100       (nr__) ! (mag) ?=-99.9 100um light curve variation amplitude (5)
      real*4        e_D100     (nr__) ! (mag) ?=-99.90 Uncertainty in the D100
      integer*4     o_F100     (nr__) ! Number of 100{mu}m unfiltered datapoints
      real*8        F140       (nr__) ! (Jy) ?=-99.9 140{mu}m weighted mean flux density (4)
      real*8        e_F140     (nr__) ! (Jy) ?=-99.9 RMS error of F140
      real*8        eF140      (nr__) ! (Jy) ?=-99.9 Average 140{mu}m light curve errorbar
      real*4        D140       (nr__) ! (mag) ?=-99.9 140um light curve variation amplitude (5)
      real*4        e_D140     (nr__) ! (mag) ?=-99.90 Uncertainty in the D140
      integer*4     o_F140     (nr__) ! Number of 140{mu}m unfiltered datapoints
      real*8        F240       (nr__) ! (Jy) ?=-99.9 240{mu}m weighted mean flux density (4)
      real*8        e_F240     (nr__) ! (Jy) ?=-99.9 RMS error of F240
      real*8        eF240      (nr__) ! (Jy) ?=-99.9 Average 240{mu}m light curve errorbar
      real*4        D240       (nr__) ! (mag) ?=-99.9 240um light curve variation amplitude (5)
      real*4        e_D240     (nr__) ! (mag) ?=-99.90 Uncertainty in the D240
      integer*4     o_F240     (nr__) ! Number of 240{mu}m unfiltered datapoints
      integer*4     f_F1_25    (nr__) ! 1.25{mu}m confusion flags (6)
      integer*4     f_F2_2     (nr__) ! 2.2{mu}m confusion flags (6)
      integer*4     f_F3_5     (nr__) ! 3.5{mu}m confusion flags (6)
      integer*4     f_F4_9     (nr__) ! 4.9{mu}m confusion flags (6)
      integer*4     f_F12      (nr__) ! 12{mu}m confusion flags (6)
      integer*4     f_F25      (nr__) ! 25{mu}m confusion flags (6)
      integer*4     f_F60      (nr__) ! 60{mu}m confusion flags (6)
      integer*4     f_F100     (nr__) ! 100{mu}m confusion flags (6)
      integer*4     f_F140     (nr__) ! 140{mu}m confusion flags (6)
      integer*4     f_F240     (nr__) ! 240{mu}m confusion flags (6)
*Note (2): Because of the method that the Catalog was constructed,
*    all DIRBE Point Source Catalog sources have either an IRAS/MSX or a
*    2MASS counterpart (with the exception of the four special sources
*    noted in the text). A 2MASS name is not given if there was no 2MASS
*    counterpart to the IRAS/MSX source with J<=4.51 and/or K<=3.81.
*Note (3): Note that an IRAS/MSX name is only given if there is an IRAS/MSX
*    counterpart to the 2MASS source with F(12um)>=15Jy and/or F(25um)>=27.5Jy.
*Note (4): -99.9 = no information is available.
*Note (5): Calculated from the light curve after it has been
*    block-averaged in weekly blocks using a running average. Note that
*    this may not be the full amplitude of variation for the star, since
*    the light curve may be incomplete or not cover a full pulsation
*    period. This parameter is set to -99.9 if the weekly-averaged light
*    curve at minimum is detected at less than a 3{sigma} level.
*Note (6): Confusion flags are as follows:
*   Flag1 = 1: IRAS/2MASS/MSX/CIO/(Egan & Price, Cat. <V/98>)
*      companion greater than the DIRBE noise level within a 0.5 degree
*      radius. This flag was not set at 140 or 240{mu}m
*   Flag2 = 1: Greater than 3{sigma} discrepancy between the
*      time-averaged DIRBE flux density and the IRAS or 2MASS flux density.
*      This flag is not set for the Catalog of Infrared Observation
*      comparison or the MSX catalog comparison, because of possible large
*      uncertainties in the photometry. It was also not set at 140 or 240
*      microns, due to the lack of an all-sky survey at those wavelengths.
*      This flag may be an indicator of extended emission around a point
*      source rather than a discrete companion. Alternatively, it may
*      indicate strong variability.
*   Flag3 = 1: This flag was set if the rms/<errorbar> in the light
*    curve was M=3 for any 2-week period containing at least 10 datapoints.
*    This is an indication that the light curve may be affected by a nearby
*    companion that was missed by the filtering process.

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

C  Declarations for 'table4.dat'	! 2MASS/IRAS/MSX Photometry

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

      character*17  DIRBE_1    (nr__1) ! DIRBE name (G1) (G2)
      character*16  v2MASS_1   (nr__1) ! 2MASS name
      character*23  IRAS_1     (nr__1) ! IRAS name
      real*8        Jmag       (nr__1) ! (mag) ?=-99.99 2MASS J Magnitude
      real*8        e_Jmag     (nr__1) ! (mag) ?=-99.990 Uncertainty in Jmag
      real*8        Hmag       (nr__1) ! (mag) ?=-99.99 2MASS H Magnitude
      real*8        e_Hmag     (nr__1) ! (mag) ?=-99.990 Uncertainty in Hmag
      real*8        Kmag       (nr__1) ! (mag) ?=-99.99 2MASS Ks Magnitude
      real*8        e_Kmag     (nr__1) ! (mag) ?=-99.990 Uncertainty in Kmag
      real*8        F12i       (nr__1) ! (Jy) ?=-99.99 The IRAS 12{mu}m flux density
      real*8        F25i       (nr__1) ! (Jy) ?=-99.99 The IRAS 25{mu}m flux density
      real*8        F60i       (nr__1) ! (Jy) ?=-99.99 The IRAS 60{mu}m flux density
      real*8        F100i      (nr__1) ! (Jy) ?=-99.99 IRAS 100{mu}m flux density
      integer*4     q_F12i     (nr__1) ! ?=-99.9 IRAS 12{mu}m quality flag
      integer*4     q_F25i     (nr__1) ! ?=-99.9 IRAS 25{mu}m quality flag
      integer*4     q_F60i     (nr__1) ! ?=-99.9 IRAS 60{mu}m quality flag
      integer*4     q_F100i    (nr__1) ! ?=-99.9 IRAS 100{mu}m quality flag
      integer*4     Var        (nr__1) ! ?=-99.9 IRAS VAR parameter
      character*1   Conf       (nr__1) ! IRAS confusion flag
      integer*4     CIRR1      (nr__1) ! ?=-99.9 IRAS CIRR1 flag
      integer*4     CIRR2      (nr__1) ! ?=-99.9 IRAS CIRR2 flag
      integer*4     CIRR3      (nr__1) ! ?=-99.9 IRAS CIRR3 flag

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

C  Declarations for 'table5.dat'	! Companion Information

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

      character*17  DIRBE_2    (nr__2) ! DIRBE name (G1) (1)
      real*4        Wave       (nr__2) ! (um) DIRBE band central wavelength (2)
      character*23  Name       (nr__2) ! Name of the companion (3)
      real*8        Flux       (nr__2) ! (Jy) Catalog companion flux density
      real*8        Off        (nr__2) ! (deg) Companion offset from the DIRBE Catalog position
*Note (1): This table contains information about companion sources
*    within the DIRBE beam. DIRBE Catalog sources are only included in this
*    table if they have been flagged because they have a companion in the
*    DIRBE beam above 25, 20, 20, 10, 30, or 55 Jy at 1.25, 2.2, 3.5, 4.9,
*    12, and 25{mu}m, respectively.
*Note (2): Note that these are the DIRBE wavelengths, not the wavelength
*    of the original observation, which may be slightly different.
*Note (3): the companions are from:
*    - At 1.25 and 2.2{mu}m, companions are from the 2MASS Catalog.
*    - At 12 or 25{mu}m, the companions are from the IRAS Point Source
*      Catalog or the IRAS Small Scale Structure Catalog.
*    - At 3.5{mu}m, companions are from the Catalog of Infrared Observations.
*    - At 4.9{mu}m, companions are either from the Catalog of Infrared
*      Observations, or from the synthetic Egan and Price (1996, Cat. <V/98>) 
*      4.2{mu}m list.
*    Names of sources from the Egan and Price (1996, Cat. <V/98>) list are
*    preceded by `E+P' to distinguish them from those in the Catalog of
*    Infrared Observations.

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

C  Declarations for 'table8.dat'	! DIRBE Catalog Associations

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

C  Position composed of: GLON GLAT
      character*17  DIRBE_3    (nr__3) ! The DIRBE name (G1) (G2)
      character*16  v2MASS_2   (nr__3) ! The 2MASS name
      character*23  IRAS_2     (nr__3) ! The IRAS name
      character*26  SName      (nr__3) ! The SIMBAD Name
      character*3   SObj       (nr__3) ! The SIMBAD object type
      character*11  SType      (nr__3) ! The SIMBAD spectral type
      real*4        Off_1      (nr__3) ! (arcsec) Offset between SIMBAD and DIRBE positions
      character*26  SName2     (nr__3) ! The second SIMBAD name (2)
      character*3   SObj2      (nr__3) ! The second SIMBAD source object type
      character*11  SType2     (nr__3) ! The second SIMBAD source spectral type
      real*4        SOff       (nr__3) ! (arcsec) ? Offset between the two SIMBAD positions
      character*1   LRS        (nr__3) ! The IRAS LRS type (3)
      character*12  Name_1     (nr__3) ! The GCVS, NSV, or NSVS name
      character*10  VarType    (nr__3) ! The GCVS variability type
      character*2   OH         (nr__3) ! OH/IR star compilation (4)
      character*1   l_Per      (nr__3) ! Limit flag on Per
      real*8        Per        (nr__3) ! (d) ?=-99.9 The GCVS/NSV/NSVS period
      real*4        EVel       (nr__3) ! (km/s) ?=-99.9 The OH expansion velocity (5)
      character*1   Cat        (nr__3) ! [TU0] The 2MASS associated optical catalog (6)
      real*4        Bmag       (nr__3) ! (mag) ?=-99.99 2MASS optical B band magnitude
      real*4        Vmag       (nr__3) ! (mag) ?=-99.99 2MASS optical V band magnitude
      real*4        GLON       (nr__3) ! (deg) Galactic longitude
      real*4        GLAT       (nr__3) ! (deg) Galactic latitude
*Note (2): If within 5 arcsec of first position.
*Note (3): From Kwok et al. (1997, Cat. <J/ApJS/112/557>; see Cat. <III/197>)
*Note (4): Set to `OH' if in Chen et al. (2001, Cat. <J/A+A/368/1006>)
*           OH/IR star compilation
*Note (5): From Chen et al. (2001, Cat. <J/A+A/368/1006>) OH/IR star compilation.
*Note (6): T = Tycho-2 (Cat. I/259), U = USNO-A2.0 (Cat. I/252).
*     See 2MASS Explanatory Supplement for more information.

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

C  Declarations for 'filters.dat'	! Transmission of DIRBE filters

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

      real*8        lambda     (nr__4) ! (um) Wavelength
      real*4        T1         (nr__4) ! [0,1] Transmission for filter#1 (1.25um)
      real*4        T2         (nr__4) ! [0,1] Transmission for filter#2 (2.2um)
      real*4        T3         (nr__4) ! [0,1] Transmission for filter#3 (3.5um)
      real*4        T4         (nr__4) ! [0,1] Transmission for filter#4 (4.9um)
      real*4        T5         (nr__4) ! [0,1] Transmission for filter#5 (12um)
      real*4        T6         (nr__4) ! [0,1] Transmission for filter#6 (25um)
      real*4        T7         (nr__4) ! [0,1] Transmission for filter#7 (60um)
      real*4        T8         (nr__4) ! [0,1] Transmission for filter#8 (100um)
      real*4        T9         (nr__4) ! [0,1] Transmission for filter#9 (140um)
      real*4        T10        (nr__4) ! [0,1] Transmission for filter#10 (240um)

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

C  Loading file 'table1.dat'	! DIRBE Photometry and Variability

C  Format for file interpretation

    1 format(
     +  A17,1X,A16,1X,A23,1X,F9.1,1X,F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,
     +  I4,1X,F9.1,1X,F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,I4,1X,F9.1,1X,
     +  F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,I4,1X,F9.1,1X,F9.1,1X,F9.1,1X,
     +  F6.2,1X,F6.2,1X,I4,1X,F9.1,1X,F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,
     +  I4,1X,F9.1,1X,F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,I4,1X,F9.1,1X,
     +  F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,I4,1X,F9.1,1X,F9.1,1X,F9.1,1X,
     +  F6.2,1X,F6.2,1X,I4,1X,F9.1,1X,F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,
     +  I4,1X,F9.1,1X,F9.1,1X,F9.1,1X,F6.2,1X,F6.2,1X,I4,1X,I3,1X,I3,
     +  1X,I3,1X,I3,1X,I3,1X,I3,1X,I3,1X,I3,1X,I3,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,11788
        read(1,'(A592)')ar__
        read(ar__,1)
     +  DIRBE(i__),v2MASS(i__),IRAS(i__),F1_25(i__),e_F1_25(i__),
     +  eF1_25(i__),D1_25(i__),e_D1_25(i__),o_F1_25(i__),F2_2(i__),
     +  e_F2_2(i__),eF2_2(i__),D2_2(i__),e_D2_2(i__),o_F2_2(i__),
     +  F3_5(i__),e_F3_5(i__),eF3_5(i__),D3_5(i__),e_D3_5(i__),
     +  o_F3_5(i__),F4_9(i__),e_F4_9(i__),eF4_9(i__),D4_9(i__),
     +  e_D4_9(i__),o_F4_9(i__),F12(i__),e_F12(i__),eF12(i__),
     +  D12(i__),e_D12(i__),o_F12(i__),F25(i__),e_F25(i__),eF25(i__),
     +  D25(i__),e_D25(i__),o_F25(i__),F60(i__),e_F60(i__),eF60(i__),
     +  D60(i__),e_D60(i__),o_F60(i__),F100(i__),e_F100(i__),
     +  eF100(i__),D100(i__),e_D100(i__),o_F100(i__),F140(i__),
     +  e_F140(i__),eF140(i__),D140(i__),e_D140(i__),o_F140(i__),
     +  F240(i__),e_F240(i__),eF240(i__),D240(i__),e_D240(i__),
     +  o_F240(i__),f_F1_25(i__),f_F2_2(i__),f_F3_5(i__),f_F4_9(i__),
     +  f_F12(i__),f_F25(i__),f_F60(i__),f_F100(i__),f_F140(i__),
     +  f_F240(i__)
c    ..............Just test output...........
        write(6,1)
     +  DIRBE(i__),v2MASS(i__),IRAS(i__),F1_25(i__),e_F1_25(i__),
     +  eF1_25(i__),D1_25(i__),e_D1_25(i__),o_F1_25(i__),F2_2(i__),
     +  e_F2_2(i__),eF2_2(i__),D2_2(i__),e_D2_2(i__),o_F2_2(i__),
     +  F3_5(i__),e_F3_5(i__),eF3_5(i__),D3_5(i__),e_D3_5(i__),
     +  o_F3_5(i__),F4_9(i__),e_F4_9(i__),eF4_9(i__),D4_9(i__),
     +  e_D4_9(i__),o_F4_9(i__),F12(i__),e_F12(i__),eF12(i__),
     +  D12(i__),e_D12(i__),o_F12(i__),F25(i__),e_F25(i__),eF25(i__),
     +  D25(i__),e_D25(i__),o_F25(i__),F60(i__),e_F60(i__),eF60(i__),
     +  D60(i__),e_D60(i__),o_F60(i__),F100(i__),e_F100(i__),
     +  eF100(i__),D100(i__),e_D100(i__),o_F100(i__),F140(i__),
     +  e_F140(i__),eF140(i__),D140(i__),e_D140(i__),o_F140(i__),
     +  F240(i__),e_F240(i__),eF240(i__),D240(i__),e_D240(i__),
     +  o_F240(i__),f_F1_25(i__),f_F2_2(i__),f_F3_5(i__),f_F4_9(i__),
     +  f_F12(i__),f_F25(i__),f_F60(i__),f_F100(i__),f_F140(i__),
     +  f_F240(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! 2MASS/IRAS/MSX Photometry

C  Format for file interpretation

    2 format(
     +  A17,1X,A16,1X,A23,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,
     +  F7.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,F9.3,1X,I3,1X,I3,1X,I3,1X,I4,
     +  1X,I3,1X,A1,1X,I3,1X,I3,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,11788
        read(1,'(A181)')ar__1
        read(ar__1,2)
     +  DIRBE_1(i__),v2MASS_1(i__),IRAS_1(i__),Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),F12i(i__),
     +  F25i(i__),F60i(i__),F100i(i__),q_F12i(i__),q_F25i(i__),
     +  q_F60i(i__),q_F100i(i__),Var(i__),Conf(i__),CIRR1(i__),
     +  CIRR2(i__),CIRR3(i__)
c    ..............Just test output...........
        write(6,2)
     +  DIRBE_1(i__),v2MASS_1(i__),IRAS_1(i__),Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),F12i(i__),
     +  F25i(i__),F60i(i__),F100i(i__),q_F12i(i__),q_F25i(i__),
     +  q_F60i(i__),q_F100i(i__),Var(i__),Conf(i__),CIRR1(i__),
     +  CIRR2(i__),CIRR3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Companion Information

C  Format for file interpretation

    3 format(A17,1X,F6.2,1X,A23,1X,F10.1,1X,F9.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,94658
        read(1,'(A69)')ar__2
        read(ar__2,3)
     +  DIRBE_2(i__),Wave(i__),Name(i__),Flux(i__),Off(i__)
c    ..............Just test output...........
        write(6,3)
     +  DIRBE_2(i__),Wave(i__),Name(i__),Flux(i__),Off(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! DIRBE Catalog Associations

C  Format for file interpretation

    4 format(
     +  A17,1X,A16,1X,A23,1X,A26,1X,A3,1X,A11,1X,F4.1,1X,A26,1X,A3,1X,
     +  A11,1X,F4.1,1X,A1,1X,A12,1X,A10,1X,A2,1X,A1,F8.2,1X,F5.1,1X,
     +  A1,4X,F6.2,4X,F6.2,1X,F6.2,1X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,11788
        read(1,'(A235)')ar__3
        read(ar__3,4)
     +  DIRBE_3(i__),v2MASS_2(i__),IRAS_2(i__),SName(i__),SObj(i__),
     +  SType(i__),Off_1(i__),SName2(i__),SObj2(i__),SType2(i__),
     +  SOff(i__),LRS(i__),Name_1(i__),VarType(i__),OH(i__),
     +  l_Per(i__),Per(i__),EVel(i__),Cat(i__),Bmag(i__),Vmag(i__),
     +  GLON(i__),GLAT(i__)
        if(ar__3(151:154) .EQ. '') SOff(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  DIRBE_3(i__),v2MASS_2(i__),IRAS_2(i__),SName(i__),SObj(i__),
     +  SType(i__),Off_1(i__),SName2(i__),SObj2(i__),SType2(i__),
     +  SOff(i__),LRS(i__),Name_1(i__),VarType(i__),OH(i__),
     +  l_Per(i__),Per(i__),EVel(i__),Cat(i__),Bmag(i__),Vmag(i__),
     +  GLON(i__),GLAT(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'filters.dat'	! Transmission of DIRBE filters

C  Format for file interpretation

    5 format(
     +  F7.3,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'filters.dat')
      write(6,*) '....Loading file: filters.dat'
      do i__=1,800
        read(1,'(A58)')ar__4
        read(ar__4,5)
     +  lambda(i__),T1(i__),T2(i__),T3(i__),T4(i__),T5(i__),T6(i__),
     +  T7(i__),T8(i__),T9(i__),T10(i__)
c    ..............Just test output...........
        write(6,5)
     +  lambda(i__),T1(i__),T2(i__),T3(i__),T4(i__),T5(i__),T6(i__),
     +  T7(i__),T8(i__),T9(i__),T10(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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