FORTRAN Generation
(/./ftp/cats/IX/1)

Conversion of standardized ReadMe file for file /./ftp/cats/IX/1 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 2018-Sep-25
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. IX/1    Wisconsin soft X-ray diffuse background all-sky Survey (McCammon+ 1983)
*================================================================================
*Wisconsin Soft X-Ray Diffuse Background All-Sky Survey
*    McCammon D., Burrows D.N., Sanders W.T., Kraushaar W.L.
*   <Astrophys. J. 269, 107 (1983)>
*   =1983ApJ...269..107M
*   =1997yCat.9001....0M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'effaprod.dat'	! Effective Area Product as a function of energy
                                for the different X-ray bands

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

      real*4        E          (nr__) ! (keV) [0.0025/10.6] Photon energy
      real*4        Area_B     (nr__) ! (cm2.sr) Effective area in the B(130-188eV) band
      real*4        Area_C     (nr__) ! (cm2.sr) Effective area in the C(160-284eV)-band
      real*4        Area_M1    (nr__) ! (cm2.sr) Effective area in the M1(440-930eV)-band
      real*4        Area_M2    (nr__) ! (cm2.sr) Effective area in the M2(600-1100eV)-band
      real*4        Area_I     (nr__) ! (cm2.sr) Effective area in the I(770-1500eV)-band
      real*4        Area_J     (nr__) ! (cm2.sr) Effective area in the J(1100-2200eV)-band
      real*4        Area_26    (nr__) ! (cm2.sr) Effective area in the 1.8-6.3keV band

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

C  Declarations for 'p13049.dat'	! Observations from flight 13.049 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight     (nr__1) ! Flight number
      integer*4     time_m     (nr__1) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s     (nr__1) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag       (nr__1) ! [0-7] Octal flag word (1)
      real*4        GLON       (nr__1) ! (deg) [0/360] Galactic longitude
      real*4        GLAT       (nr__1) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg      (nr__1) ! (deg) Right Ascension (B1950)
      real*4        DEdeg      (nr__1) ! (deg) Declination (B1950)
      real*4        Az         (nr__1) ! (deg) [0/360] Azimuth angle
      real*4        Ze         (nr__1) ! (deg) [0/180] Zenith angle
      real*4        Alt        (nr__1) ! (km) [95/220] Rocket altitude
      real*4        Ncol       (nr__1) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol       (nr__1) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol      (nr__1) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B       (nr__1) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B       (nr__1) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate      (nr__1) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate    (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm       (nr__1) ! [0/1] Atmospheric transmission for band B
      real*4        Crate      (nr__1) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate    (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm       (nr__1) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate      (nr__1) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate    (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm       (nr__1) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate     (nr__1) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate   (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm      (nr__1) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate     (nr__1) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate   (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm      (nr__1) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate      (nr__1) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate    (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm       (nr__1) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate      (nr__1) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate    (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm       (nr__1) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate    (nr__1) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate   (nr__1) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm     (nr__1) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p13083.dat'	! Observations from flight 13.083 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_1   (nr__2) ! Flight number
      integer*4     time_m_1   (nr__2) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_1   (nr__2) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_1     (nr__2) ! [0-7] Octal flag word (1)
      real*4        GLON_1     (nr__2) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_1     (nr__2) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_1    (nr__2) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_1    (nr__2) ! (deg) Declination (B1950)
      real*4        Az_1       (nr__2) ! (deg) [0/360] Azimuth angle
      real*4        Ze_1       (nr__2) ! (deg) [0/180] Zenith angle
      real*4        Alt_1      (nr__2) ! (km) [95/220] Rocket altitude
      real*4        Ncol_1     (nr__2) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_1     (nr__2) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_1    (nr__2) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_1     (nr__2) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_1     (nr__2) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_1    (nr__2) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_1  (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_1     (nr__2) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_1    (nr__2) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_1  (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_1     (nr__2) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_1    (nr__2) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_1  (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_1     (nr__2) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_1   (nr__2) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_1 (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_1    (nr__2) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_1   (nr__2) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_1 (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_1    (nr__2) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_1    (nr__2) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_1  (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_1     (nr__2) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_1    (nr__2) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_1  (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_1     (nr__2) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_1  (nr__2) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_1 (nr__2) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_1   (nr__2) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p13084.dat'	! Observations from flight 13.084 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_2   (nr__3) ! Flight number
      integer*4     time_m_2   (nr__3) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_2   (nr__3) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_2     (nr__3) ! [0-7] Octal flag word (1)
      real*4        GLON_2     (nr__3) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_2     (nr__3) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_2    (nr__3) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_2    (nr__3) ! (deg) Declination (B1950)
      real*4        Az_2       (nr__3) ! (deg) [0/360] Azimuth angle
      real*4        Ze_2       (nr__3) ! (deg) [0/180] Zenith angle
      real*4        Alt_2      (nr__3) ! (km) [95/220] Rocket altitude
      real*4        Ncol_2     (nr__3) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_2     (nr__3) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_2    (nr__3) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_2     (nr__3) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_2     (nr__3) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_2    (nr__3) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_2  (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_2     (nr__3) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_2    (nr__3) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_2  (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_2     (nr__3) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_2    (nr__3) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_2  (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_2     (nr__3) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_2   (nr__3) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_2 (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_2    (nr__3) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_2   (nr__3) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_2 (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_2    (nr__3) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_2    (nr__3) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_2  (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_2     (nr__3) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_2    (nr__3) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_2  (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_2     (nr__3) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_2  (nr__3) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_2 (nr__3) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_2   (nr__3) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p13102.dat'	! Observations from flight 13.102 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_3   (nr__4) ! Flight number
      integer*4     time_m_3   (nr__4) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_3   (nr__4) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_3     (nr__4) ! [0-7] Octal flag word (1)
      real*4        GLON_3     (nr__4) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_3     (nr__4) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_3    (nr__4) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_3    (nr__4) ! (deg) Declination (B1950)
      real*4        Az_3       (nr__4) ! (deg) [0/360] Azimuth angle
      real*4        Ze_3       (nr__4) ! (deg) [0/180] Zenith angle
      real*4        Alt_3      (nr__4) ! (km) [95/220] Rocket altitude
      real*4        Ncol_3     (nr__4) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_3     (nr__4) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_3    (nr__4) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_3     (nr__4) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_3     (nr__4) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_3    (nr__4) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_3  (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_3     (nr__4) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_3    (nr__4) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_3  (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_3     (nr__4) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_3    (nr__4) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_3  (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_3     (nr__4) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_3   (nr__4) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_3 (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_3    (nr__4) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_3   (nr__4) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_3 (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_3    (nr__4) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_3    (nr__4) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_3  (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_3     (nr__4) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_3    (nr__4) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_3  (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_3     (nr__4) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_3  (nr__4) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_3 (nr__4) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_3   (nr__4) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p13103.dat'	! Observations from flight 13.103 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_4   (nr__5) ! Flight number
      integer*4     time_m_4   (nr__5) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_4   (nr__5) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_4     (nr__5) ! [0-7] Octal flag word (1)
      real*4        GLON_4     (nr__5) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_4     (nr__5) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_4    (nr__5) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_4    (nr__5) ! (deg) Declination (B1950)
      real*4        Az_4       (nr__5) ! (deg) [0/360] Azimuth angle
      real*4        Ze_4       (nr__5) ! (deg) [0/180] Zenith angle
      real*4        Alt_4      (nr__5) ! (km) [95/220] Rocket altitude
      real*4        Ncol_4     (nr__5) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_4     (nr__5) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_4    (nr__5) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_4     (nr__5) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_4     (nr__5) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_4    (nr__5) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_4  (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_4     (nr__5) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_4    (nr__5) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_4  (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_4     (nr__5) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_4    (nr__5) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_4  (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_4     (nr__5) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_4   (nr__5) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_4 (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_4    (nr__5) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_4   (nr__5) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_4 (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_4    (nr__5) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_4    (nr__5) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_4  (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_4     (nr__5) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_4    (nr__5) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_4  (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_4     (nr__5) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_4  (nr__5) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_4 (nr__5) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_4   (nr__5) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p13122.dat'	! Observations from flight 13.122 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_5   (nr__6) ! Flight number
      integer*4     time_m_5   (nr__6) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_5   (nr__6) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_5     (nr__6) ! [0-7] Octal flag word (1)
      real*4        GLON_5     (nr__6) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_5     (nr__6) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_5    (nr__6) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_5    (nr__6) ! (deg) Declination (B1950)
      real*4        Az_5       (nr__6) ! (deg) [0/360] Azimuth angle
      real*4        Ze_5       (nr__6) ! (deg) [0/180] Zenith angle
      real*4        Alt_5      (nr__6) ! (km) [95/220] Rocket altitude
      real*4        Ncol_5     (nr__6) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_5     (nr__6) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_5    (nr__6) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_5     (nr__6) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_5     (nr__6) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_5    (nr__6) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_5  (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_5     (nr__6) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_5    (nr__6) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_5  (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_5     (nr__6) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_5    (nr__6) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_5  (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_5     (nr__6) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_5   (nr__6) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_5 (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_5    (nr__6) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_5   (nr__6) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_5 (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_5    (nr__6) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_5    (nr__6) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_5  (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_5     (nr__6) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_5    (nr__6) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_5  (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_5     (nr__6) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_5  (nr__6) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_5 (nr__6) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_5   (nr__6) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p13137.dat'	! Observations from flight 13.137 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_6   (nr__7) ! Flight number
      integer*4     time_m_6   (nr__7) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_6   (nr__7) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_6     (nr__7) ! [0-7] Octal flag word (1)
      real*4        GLON_6     (nr__7) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_6     (nr__7) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_6    (nr__7) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_6    (nr__7) ! (deg) Declination (B1950)
      real*4        Az_6       (nr__7) ! (deg) [0/360] Azimuth angle
      real*4        Ze_6       (nr__7) ! (deg) [0/180] Zenith angle
      real*4        Alt_6      (nr__7) ! (km) [95/220] Rocket altitude
      real*4        Ncol_6     (nr__7) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_6     (nr__7) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_6    (nr__7) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_6     (nr__7) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_6     (nr__7) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_6    (nr__7) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_6  (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_6     (nr__7) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_6    (nr__7) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_6  (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_6     (nr__7) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_6    (nr__7) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_6  (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_6     (nr__7) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_6   (nr__7) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_6 (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_6    (nr__7) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_6   (nr__7) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_6 (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_6    (nr__7) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_6    (nr__7) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_6  (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_6     (nr__7) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_6    (nr__7) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_6  (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_6     (nr__7) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_6  (nr__7) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_6 (nr__7) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_6   (nr__7) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p25045.dat'	! Observations from flight 25.045 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_7   (nr__8) ! Flight number
      integer*4     time_m_7   (nr__8) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_7   (nr__8) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_7     (nr__8) ! [0-7] Octal flag word (1)
      real*4        GLON_7     (nr__8) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_7     (nr__8) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_7    (nr__8) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_7    (nr__8) ! (deg) Declination (B1950)
      real*4        Az_7       (nr__8) ! (deg) [0/360] Azimuth angle
      real*4        Ze_7       (nr__8) ! (deg) [0/180] Zenith angle
      real*4        Alt_7      (nr__8) ! (km) [95/220] Rocket altitude
      real*4        Ncol_7     (nr__8) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_7     (nr__8) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_7    (nr__8) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_7     (nr__8) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_7     (nr__8) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_7    (nr__8) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_7  (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_7     (nr__8) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_7    (nr__8) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_7  (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_7     (nr__8) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_7    (nr__8) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_7  (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_7     (nr__8) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_7   (nr__8) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_7 (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_7    (nr__8) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_7   (nr__8) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_7 (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_7    (nr__8) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_7    (nr__8) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_7  (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_7     (nr__8) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_7    (nr__8) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_7  (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_7     (nr__8) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_7  (nr__8) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_7 (nr__8) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_7   (nr__8) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p25051.dat'	! Observations from flight 25.051 (PASS3)

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

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_8   (nr__9) ! Flight number
      integer*4     time_m_8   (nr__9) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_8   (nr__9) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_8     (nr__9) ! [0-7] Octal flag word (1)
      real*4        GLON_8     (nr__9) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_8     (nr__9) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_8    (nr__9) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_8    (nr__9) ! (deg) Declination (B1950)
      real*4        Az_8       (nr__9) ! (deg) [0/360] Azimuth angle
      real*4        Ze_8       (nr__9) ! (deg) [0/180] Zenith angle
      real*4        Alt_8      (nr__9) ! (km) [95/220] Rocket altitude
      real*4        Ncol_8     (nr__9) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_8     (nr__9) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_8    (nr__9) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_8     (nr__9) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_8     (nr__9) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_8    (nr__9) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_8  (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_8     (nr__9) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_8    (nr__9) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_8  (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_8     (nr__9) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_8    (nr__9) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_8  (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_8     (nr__9) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_8   (nr__9) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_8 (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_8    (nr__9) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_8   (nr__9) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_8 (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_8    (nr__9) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_8    (nr__9) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_8  (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_8     (nr__9) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_8    (nr__9) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_8  (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_8     (nr__9) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_8  (nr__9) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_8 (nr__9) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_8   (nr__9) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'p26061.dat'	! Observations from flight 26.061 (PASS3)

      integer*4 nr__10
      parameter (nr__10=1245)	! Number of records
      character*400 ar__10 	! Full-size record

C  B1950 position composed of: RAdeg DEdeg
      real*4        Flight_9   (nr__10) ! Flight number
      integer*4     time_m_9   (nr__10) ! (min) Minutes past the hour for the center of the
*                                    0.2s observation
      real*4        time_s_9   (nr__10) ! (s) Seconds past the minute for the center of
*                                    the 0.2s observation
      character*6   Flag_9     (nr__10) ! [0-7] Octal flag word (1)
      real*4        GLON_9     (nr__10) ! (deg) [0/360] Galactic longitude
      real*4        GLAT_9     (nr__10) ! (deg) [-90/90] Galactic latitude
      real*4        RAdeg_9    (nr__10) ! (deg) Right Ascension (B1950)
      real*4        DEdeg_9    (nr__10) ! (deg) Declination (B1950)
      real*4        Az_9       (nr__10) ! (deg) [0/360] Azimuth angle
      real*4        Ze_9       (nr__10) ! (deg) [0/180] Zenith angle
      real*4        Alt_9      (nr__10) ! (km) [95/220] Rocket altitude
      real*4        Ncol_9     (nr__10) ! (cm-2) Column density of nitrogen atoms (2)
      real*4        Ocol_9     (nr__10) ! (cm-2) Column density of oxygen atoms (2)
      real*4        Arcol_9    (nr__10) ! (cm-2) Column density of argon atoms (3)
      real*4        Az_B_9     (nr__10) ! (deg) [0/360] Azimuth angle, relative to the Earth's
*                                    magnetic field
      real*4        Ze_B_9     (nr__10) ! (deg) [0/180] Zenith angle, relative to the Earth's
*                                    magnetic field
      real*4        Brate_9    (nr__10) ! (ct/s) ?=-99 Count rate in band B(130-188eV)
      real*4        e_Brate_9  (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for Brate
      real*4        Btrm_9     (nr__10) ! [0/1] Atmospheric transmission for band B
      real*4        Crate_9    (nr__10) ! (ct/s) ?=-99 Count rate in band C(160-284eV)
      real*4        e_Crate_9  (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for Crate
      real*4        Ctrm_9     (nr__10) ! [0/1] Atmospheric transmission for band C
      real*4        Mrate_9    (nr__10) ! (ct/s) ?=-99 Count rate in band M
      real*4        e_Mrate_9  (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for Mrate
      real*4        Mtrm_9     (nr__10) ! [0/1] Atmospheric transmission for band M
      real*4        M1rate_9   (nr__10) ! (ct/s) ?=-99 Count rate in band M1(440-930eV)
      real*4        e_M1rate_9 (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for M1rate
      real*4        M1trm_9    (nr__10) ! [0/1] Atmospheric transmission for band M1
      real*4        M2rate_9   (nr__10) ! (ct/s) ?=-99 Count rate in band M2(600-1100eV)
      real*4        e_M2rate_9 (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for M2rate
      real*4        M2trm_9    (nr__10) ! [0/1] Atmospheric transmission for band M2
      real*4        Irate_9    (nr__10) ! (ct/s) ?=-99 Count rate in band I(770-1500eV)
      real*4        e_Irate_9  (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for Irate
      real*4        Itrm_9     (nr__10) ! [0/1] Atmospheric transmission for band I
      real*4        Jrate_9    (nr__10) ! (ct/s) ?=-99 Count rate in band J(1100-2200eV)
      real*4        e_Jrate_9  (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for Jrate
      real*4        Jtrm_9     (nr__10) ! [0/1] Atmospheric transmission for band J
      real*4        v26rate_9  (nr__10) ! (ct/s) ?=-99 Count rate in band 1.8-6.3keV
      real*4        e_26rate_9 (nr__10) ! (ct/s) ?=-99 1-{sigma} uncertainty for 26rate
      real*4        v26trm_9   (nr__10) ! [0/1] Atmospheric transmission for 1.8-6.3keV
*Note (1): these 6 characters should be interpreted as an 18-bit octal word,
*   with only the least significant (rightmost) 16 bits being used.
*   They contain the following information (bits numbered 1-16 from
*   left to right - most significant to least significant):
*   1 = set to 1 if the observation contains bad telemetry.
*   2 = set to 1 if there are fewer than 50 useful telemetry subframes in the
*       observation (each observation can contain a maximum of 125 subframes
*       for flights 13.083 and 13.102, or 250 for the other flights).
*   3 = set to 1 if the livetime correction is uncertain.
*   5 = set to 1 if the B and/or C-band data are definitely contaminated.
*   6 = set to 1 if the B and/or C-band data are possibly contaminated.
*  10 = set to 1 if the payload was pointed at the Earth.
*  11 = set to 1 if a B and/or C-band source was in the field of view. These are!
*       HZ43, Am Her, the Cygnus Loop, Vela-Puppis SNR's, Sco X-1, and U Gem.
*  12 = set to 1 if the M and/or I-band data show evidence for the presence
*       of a source in the field of view.
*  13 = set to 1 if the J-and/or 2-6keV band data show evidence for the
*       presence of a source in the field of view.
*  14 = set to 1 if the following sources are in the field of view:
*       sources listed in the 4-th Uhuru catalog (4U) with >=20ct/s;
*       Vela-Puppis SNR's; Cygnus Loop.
*
*   Example:  The flag word of the first record of file 31 (13.083) is written
*   as:    004000, indicating that bit 5 is set because the B and C-band data
*   for that observation are known to be contaminated.
*
*   Note that records were not cast onto the maps if the following flags
*   were set:
*   * For the B and C-bands: 1, 2, 3, 5, 10, 11
*   * For the M1-through 2-6keV bands: 1, 2, 3, 10, 14
*Note (2): number of atoms along the line of sight in the atmosphere,
*   equals twice the number of molecules per cm^2^ calculated from
*   CIRA72 model atmosphere, as described in the paper
*Note (3): number of atoms along the line of sight in the atmosphere,
*   calculated from CIRA72 model atmosphere, as described in the paper

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

C  Declarations for 'maps.dat'	! Summary of FITS maps

      integer*4 nr__11
      parameter (nr__11=36)	! Number of records
      character*83 ar__11 	! Full-size record

      integer*4     Nx         (nr__11) ! [92/182] Number of pixels along X-axis
      integer*4     Ny         (nr__11) ! [92/92] Number of pixels along Y-axis
      integer*4     Nz         (nr__11) ! [2/2] Number of frames (map and error map)
      character*12  FileName   (nr__11) ! Name of FITS file
      character*61  Text       (nr__11) ! Explanations of file

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

C  Loading file 'effaprod.dat'	! Effective Area Product as a function of energy
*                                for the different X-ray bands

C  Format for file interpretation

    1 format(
     +  1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,
     +  E9.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'effaprod.dat')
      write(6,*) '....Loading file: effaprod.dat'
      do i__=1,560
        read(1,'(A80)')ar__
        read(ar__,1)
     +  E(i__),Area_B(i__),Area_C(i__),Area_M1(i__),Area_M2(i__),
     +  Area_I(i__),Area_J(i__),Area_26(i__)
c    ..............Just test output...........
        write(6,1)
     +  E(i__),Area_B(i__),Area_C(i__),Area_M1(i__),Area_M2(i__),
     +  Area_I(i__),Area_J(i__),Area_26(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p13049.dat'	! Observations from flight 13.049 (PASS3)

C  Format for file interpretation

    2 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p13049.dat')
      write(6,*) '....Loading file: p13049.dat'
      do i__=1,1176
        read(1,'(A400)')ar__1
        read(ar__1,2)
     +  Flight(i__),time_m(i__),time_s(i__),Flag(i__),GLON(i__),
     +  GLAT(i__),RAdeg(i__),DEdeg(i__),Az(i__),Ze(i__),Alt(i__),
     +  Ncol(i__),Ocol(i__),Arcol(i__),Az_B(i__),Ze_B(i__),Brate(i__),
     +  e_Brate(i__),Btrm(i__),Crate(i__),e_Crate(i__),Ctrm(i__),
     +  Mrate(i__),e_Mrate(i__),Mtrm(i__),M1rate(i__),e_M1rate(i__),
     +  M1trm(i__),M2rate(i__),e_M2rate(i__),M2trm(i__),Irate(i__),
     +  e_Irate(i__),Itrm(i__),Jrate(i__),e_Jrate(i__),Jtrm(i__),
     +  v26rate(i__),e_26rate(i__),v26trm(i__)
c    ..............Just test output...........
        write(6,2)
     +  Flight(i__),time_m(i__),time_s(i__),Flag(i__),GLON(i__),
     +  GLAT(i__),RAdeg(i__),DEdeg(i__),Az(i__),Ze(i__),Alt(i__),
     +  Ncol(i__),Ocol(i__),Arcol(i__),Az_B(i__),Ze_B(i__),Brate(i__),
     +  e_Brate(i__),Btrm(i__),Crate(i__),e_Crate(i__),Ctrm(i__),
     +  Mrate(i__),e_Mrate(i__),Mtrm(i__),M1rate(i__),e_M1rate(i__),
     +  M1trm(i__),M2rate(i__),e_M2rate(i__),M2trm(i__),Irate(i__),
     +  e_Irate(i__),Itrm(i__),Jrate(i__),e_Jrate(i__),Jtrm(i__),
     +  v26rate(i__),e_26rate(i__),v26trm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p13083.dat'	! Observations from flight 13.083 (PASS3)

C  Format for file interpretation

    3 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p13083.dat')
      write(6,*) '....Loading file: p13083.dat'
      do i__=1,1186
        read(1,'(A400)')ar__2
        read(ar__2,3)
     +  Flight_1(i__),time_m_1(i__),time_s_1(i__),Flag_1(i__),
     +  GLON_1(i__),GLAT_1(i__),RAdeg_1(i__),DEdeg_1(i__),Az_1(i__),
     +  Ze_1(i__),Alt_1(i__),Ncol_1(i__),Ocol_1(i__),Arcol_1(i__),
     +  Az_B_1(i__),Ze_B_1(i__),Brate_1(i__),e_Brate_1(i__),
     +  Btrm_1(i__),Crate_1(i__),e_Crate_1(i__),Ctrm_1(i__),
     +  Mrate_1(i__),e_Mrate_1(i__),Mtrm_1(i__),M1rate_1(i__),
     +  e_M1rate_1(i__),M1trm_1(i__),M2rate_1(i__),e_M2rate_1(i__),
     +  M2trm_1(i__),Irate_1(i__),e_Irate_1(i__),Itrm_1(i__),
     +  Jrate_1(i__),e_Jrate_1(i__),Jtrm_1(i__),v26rate_1(i__),
     +  e_26rate_1(i__),v26trm_1(i__)
c    ..............Just test output...........
        write(6,3)
     +  Flight_1(i__),time_m_1(i__),time_s_1(i__),Flag_1(i__),
     +  GLON_1(i__),GLAT_1(i__),RAdeg_1(i__),DEdeg_1(i__),Az_1(i__),
     +  Ze_1(i__),Alt_1(i__),Ncol_1(i__),Ocol_1(i__),Arcol_1(i__),
     +  Az_B_1(i__),Ze_B_1(i__),Brate_1(i__),e_Brate_1(i__),
     +  Btrm_1(i__),Crate_1(i__),e_Crate_1(i__),Ctrm_1(i__),
     +  Mrate_1(i__),e_Mrate_1(i__),Mtrm_1(i__),M1rate_1(i__),
     +  e_M1rate_1(i__),M1trm_1(i__),M2rate_1(i__),e_M2rate_1(i__),
     +  M2trm_1(i__),Irate_1(i__),e_Irate_1(i__),Itrm_1(i__),
     +  Jrate_1(i__),e_Jrate_1(i__),Jtrm_1(i__),v26rate_1(i__),
     +  e_26rate_1(i__),v26trm_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p13084.dat'	! Observations from flight 13.084 (PASS3)

C  Format for file interpretation

    4 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p13084.dat')
      write(6,*) '....Loading file: p13084.dat'
      do i__=1,1356
        read(1,'(A400)')ar__3
        read(ar__3,4)
     +  Flight_2(i__),time_m_2(i__),time_s_2(i__),Flag_2(i__),
     +  GLON_2(i__),GLAT_2(i__),RAdeg_2(i__),DEdeg_2(i__),Az_2(i__),
     +  Ze_2(i__),Alt_2(i__),Ncol_2(i__),Ocol_2(i__),Arcol_2(i__),
     +  Az_B_2(i__),Ze_B_2(i__),Brate_2(i__),e_Brate_2(i__),
     +  Btrm_2(i__),Crate_2(i__),e_Crate_2(i__),Ctrm_2(i__),
     +  Mrate_2(i__),e_Mrate_2(i__),Mtrm_2(i__),M1rate_2(i__),
     +  e_M1rate_2(i__),M1trm_2(i__),M2rate_2(i__),e_M2rate_2(i__),
     +  M2trm_2(i__),Irate_2(i__),e_Irate_2(i__),Itrm_2(i__),
     +  Jrate_2(i__),e_Jrate_2(i__),Jtrm_2(i__),v26rate_2(i__),
     +  e_26rate_2(i__),v26trm_2(i__)
c    ..............Just test output...........
        write(6,4)
     +  Flight_2(i__),time_m_2(i__),time_s_2(i__),Flag_2(i__),
     +  GLON_2(i__),GLAT_2(i__),RAdeg_2(i__),DEdeg_2(i__),Az_2(i__),
     +  Ze_2(i__),Alt_2(i__),Ncol_2(i__),Ocol_2(i__),Arcol_2(i__),
     +  Az_B_2(i__),Ze_B_2(i__),Brate_2(i__),e_Brate_2(i__),
     +  Btrm_2(i__),Crate_2(i__),e_Crate_2(i__),Ctrm_2(i__),
     +  Mrate_2(i__),e_Mrate_2(i__),Mtrm_2(i__),M1rate_2(i__),
     +  e_M1rate_2(i__),M1trm_2(i__),M2rate_2(i__),e_M2rate_2(i__),
     +  M2trm_2(i__),Irate_2(i__),e_Irate_2(i__),Itrm_2(i__),
     +  Jrate_2(i__),e_Jrate_2(i__),Jtrm_2(i__),v26rate_2(i__),
     +  e_26rate_2(i__),v26trm_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p13102.dat'	! Observations from flight 13.102 (PASS3)

C  Format for file interpretation

    5 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p13102.dat')
      write(6,*) '....Loading file: p13102.dat'
      do i__=1,1211
        read(1,'(A400)')ar__4
        read(ar__4,5)
     +  Flight_3(i__),time_m_3(i__),time_s_3(i__),Flag_3(i__),
     +  GLON_3(i__),GLAT_3(i__),RAdeg_3(i__),DEdeg_3(i__),Az_3(i__),
     +  Ze_3(i__),Alt_3(i__),Ncol_3(i__),Ocol_3(i__),Arcol_3(i__),
     +  Az_B_3(i__),Ze_B_3(i__),Brate_3(i__),e_Brate_3(i__),
     +  Btrm_3(i__),Crate_3(i__),e_Crate_3(i__),Ctrm_3(i__),
     +  Mrate_3(i__),e_Mrate_3(i__),Mtrm_3(i__),M1rate_3(i__),
     +  e_M1rate_3(i__),M1trm_3(i__),M2rate_3(i__),e_M2rate_3(i__),
     +  M2trm_3(i__),Irate_3(i__),e_Irate_3(i__),Itrm_3(i__),
     +  Jrate_3(i__),e_Jrate_3(i__),Jtrm_3(i__),v26rate_3(i__),
     +  e_26rate_3(i__),v26trm_3(i__)
c    ..............Just test output...........
        write(6,5)
     +  Flight_3(i__),time_m_3(i__),time_s_3(i__),Flag_3(i__),
     +  GLON_3(i__),GLAT_3(i__),RAdeg_3(i__),DEdeg_3(i__),Az_3(i__),
     +  Ze_3(i__),Alt_3(i__),Ncol_3(i__),Ocol_3(i__),Arcol_3(i__),
     +  Az_B_3(i__),Ze_B_3(i__),Brate_3(i__),e_Brate_3(i__),
     +  Btrm_3(i__),Crate_3(i__),e_Crate_3(i__),Ctrm_3(i__),
     +  Mrate_3(i__),e_Mrate_3(i__),Mtrm_3(i__),M1rate_3(i__),
     +  e_M1rate_3(i__),M1trm_3(i__),M2rate_3(i__),e_M2rate_3(i__),
     +  M2trm_3(i__),Irate_3(i__),e_Irate_3(i__),Itrm_3(i__),
     +  Jrate_3(i__),e_Jrate_3(i__),Jtrm_3(i__),v26rate_3(i__),
     +  e_26rate_3(i__),v26trm_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p13103.dat'	! Observations from flight 13.103 (PASS3)

C  Format for file interpretation

    6 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p13103.dat')
      write(6,*) '....Loading file: p13103.dat'
      do i__=1,1026
        read(1,'(A400)')ar__5
        read(ar__5,6)
     +  Flight_4(i__),time_m_4(i__),time_s_4(i__),Flag_4(i__),
     +  GLON_4(i__),GLAT_4(i__),RAdeg_4(i__),DEdeg_4(i__),Az_4(i__),
     +  Ze_4(i__),Alt_4(i__),Ncol_4(i__),Ocol_4(i__),Arcol_4(i__),
     +  Az_B_4(i__),Ze_B_4(i__),Brate_4(i__),e_Brate_4(i__),
     +  Btrm_4(i__),Crate_4(i__),e_Crate_4(i__),Ctrm_4(i__),
     +  Mrate_4(i__),e_Mrate_4(i__),Mtrm_4(i__),M1rate_4(i__),
     +  e_M1rate_4(i__),M1trm_4(i__),M2rate_4(i__),e_M2rate_4(i__),
     +  M2trm_4(i__),Irate_4(i__),e_Irate_4(i__),Itrm_4(i__),
     +  Jrate_4(i__),e_Jrate_4(i__),Jtrm_4(i__),v26rate_4(i__),
     +  e_26rate_4(i__),v26trm_4(i__)
c    ..............Just test output...........
        write(6,6)
     +  Flight_4(i__),time_m_4(i__),time_s_4(i__),Flag_4(i__),
     +  GLON_4(i__),GLAT_4(i__),RAdeg_4(i__),DEdeg_4(i__),Az_4(i__),
     +  Ze_4(i__),Alt_4(i__),Ncol_4(i__),Ocol_4(i__),Arcol_4(i__),
     +  Az_B_4(i__),Ze_B_4(i__),Brate_4(i__),e_Brate_4(i__),
     +  Btrm_4(i__),Crate_4(i__),e_Crate_4(i__),Ctrm_4(i__),
     +  Mrate_4(i__),e_Mrate_4(i__),Mtrm_4(i__),M1rate_4(i__),
     +  e_M1rate_4(i__),M1trm_4(i__),M2rate_4(i__),e_M2rate_4(i__),
     +  M2trm_4(i__),Irate_4(i__),e_Irate_4(i__),Itrm_4(i__),
     +  Jrate_4(i__),e_Jrate_4(i__),Jtrm_4(i__),v26rate_4(i__),
     +  e_26rate_4(i__),v26trm_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p13122.dat'	! Observations from flight 13.122 (PASS3)

C  Format for file interpretation

    7 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p13122.dat')
      write(6,*) '....Loading file: p13122.dat'
      do i__=1,1296
        read(1,'(A400)')ar__6
        read(ar__6,7)
     +  Flight_5(i__),time_m_5(i__),time_s_5(i__),Flag_5(i__),
     +  GLON_5(i__),GLAT_5(i__),RAdeg_5(i__),DEdeg_5(i__),Az_5(i__),
     +  Ze_5(i__),Alt_5(i__),Ncol_5(i__),Ocol_5(i__),Arcol_5(i__),
     +  Az_B_5(i__),Ze_B_5(i__),Brate_5(i__),e_Brate_5(i__),
     +  Btrm_5(i__),Crate_5(i__),e_Crate_5(i__),Ctrm_5(i__),
     +  Mrate_5(i__),e_Mrate_5(i__),Mtrm_5(i__),M1rate_5(i__),
     +  e_M1rate_5(i__),M1trm_5(i__),M2rate_5(i__),e_M2rate_5(i__),
     +  M2trm_5(i__),Irate_5(i__),e_Irate_5(i__),Itrm_5(i__),
     +  Jrate_5(i__),e_Jrate_5(i__),Jtrm_5(i__),v26rate_5(i__),
     +  e_26rate_5(i__),v26trm_5(i__)
c    ..............Just test output...........
        write(6,7)
     +  Flight_5(i__),time_m_5(i__),time_s_5(i__),Flag_5(i__),
     +  GLON_5(i__),GLAT_5(i__),RAdeg_5(i__),DEdeg_5(i__),Az_5(i__),
     +  Ze_5(i__),Alt_5(i__),Ncol_5(i__),Ocol_5(i__),Arcol_5(i__),
     +  Az_B_5(i__),Ze_B_5(i__),Brate_5(i__),e_Brate_5(i__),
     +  Btrm_5(i__),Crate_5(i__),e_Crate_5(i__),Ctrm_5(i__),
     +  Mrate_5(i__),e_Mrate_5(i__),Mtrm_5(i__),M1rate_5(i__),
     +  e_M1rate_5(i__),M1trm_5(i__),M2rate_5(i__),e_M2rate_5(i__),
     +  M2trm_5(i__),Irate_5(i__),e_Irate_5(i__),Itrm_5(i__),
     +  Jrate_5(i__),e_Jrate_5(i__),Jtrm_5(i__),v26rate_5(i__),
     +  e_26rate_5(i__),v26trm_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p13137.dat'	! Observations from flight 13.137 (PASS3)

C  Format for file interpretation

    8 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p13137.dat')
      write(6,*) '....Loading file: p13137.dat'
      do i__=1,1326
        read(1,'(A400)')ar__7
        read(ar__7,8)
     +  Flight_6(i__),time_m_6(i__),time_s_6(i__),Flag_6(i__),
     +  GLON_6(i__),GLAT_6(i__),RAdeg_6(i__),DEdeg_6(i__),Az_6(i__),
     +  Ze_6(i__),Alt_6(i__),Ncol_6(i__),Ocol_6(i__),Arcol_6(i__),
     +  Az_B_6(i__),Ze_B_6(i__),Brate_6(i__),e_Brate_6(i__),
     +  Btrm_6(i__),Crate_6(i__),e_Crate_6(i__),Ctrm_6(i__),
     +  Mrate_6(i__),e_Mrate_6(i__),Mtrm_6(i__),M1rate_6(i__),
     +  e_M1rate_6(i__),M1trm_6(i__),M2rate_6(i__),e_M2rate_6(i__),
     +  M2trm_6(i__),Irate_6(i__),e_Irate_6(i__),Itrm_6(i__),
     +  Jrate_6(i__),e_Jrate_6(i__),Jtrm_6(i__),v26rate_6(i__),
     +  e_26rate_6(i__),v26trm_6(i__)
c    ..............Just test output...........
        write(6,8)
     +  Flight_6(i__),time_m_6(i__),time_s_6(i__),Flag_6(i__),
     +  GLON_6(i__),GLAT_6(i__),RAdeg_6(i__),DEdeg_6(i__),Az_6(i__),
     +  Ze_6(i__),Alt_6(i__),Ncol_6(i__),Ocol_6(i__),Arcol_6(i__),
     +  Az_B_6(i__),Ze_B_6(i__),Brate_6(i__),e_Brate_6(i__),
     +  Btrm_6(i__),Crate_6(i__),e_Crate_6(i__),Ctrm_6(i__),
     +  Mrate_6(i__),e_Mrate_6(i__),Mtrm_6(i__),M1rate_6(i__),
     +  e_M1rate_6(i__),M1trm_6(i__),M2rate_6(i__),e_M2rate_6(i__),
     +  M2trm_6(i__),Irate_6(i__),e_Irate_6(i__),Itrm_6(i__),
     +  Jrate_6(i__),e_Jrate_6(i__),Jtrm_6(i__),v26rate_6(i__),
     +  e_26rate_6(i__),v26trm_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p25045.dat'	! Observations from flight 25.045 (PASS3)

C  Format for file interpretation

    9 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p25045.dat')
      write(6,*) '....Loading file: p25045.dat'
      do i__=1,1401
        read(1,'(A400)')ar__8
        read(ar__8,9)
     +  Flight_7(i__),time_m_7(i__),time_s_7(i__),Flag_7(i__),
     +  GLON_7(i__),GLAT_7(i__),RAdeg_7(i__),DEdeg_7(i__),Az_7(i__),
     +  Ze_7(i__),Alt_7(i__),Ncol_7(i__),Ocol_7(i__),Arcol_7(i__),
     +  Az_B_7(i__),Ze_B_7(i__),Brate_7(i__),e_Brate_7(i__),
     +  Btrm_7(i__),Crate_7(i__),e_Crate_7(i__),Ctrm_7(i__),
     +  Mrate_7(i__),e_Mrate_7(i__),Mtrm_7(i__),M1rate_7(i__),
     +  e_M1rate_7(i__),M1trm_7(i__),M2rate_7(i__),e_M2rate_7(i__),
     +  M2trm_7(i__),Irate_7(i__),e_Irate_7(i__),Itrm_7(i__),
     +  Jrate_7(i__),e_Jrate_7(i__),Jtrm_7(i__),v26rate_7(i__),
     +  e_26rate_7(i__),v26trm_7(i__)
c    ..............Just test output...........
        write(6,9)
     +  Flight_7(i__),time_m_7(i__),time_s_7(i__),Flag_7(i__),
     +  GLON_7(i__),GLAT_7(i__),RAdeg_7(i__),DEdeg_7(i__),Az_7(i__),
     +  Ze_7(i__),Alt_7(i__),Ncol_7(i__),Ocol_7(i__),Arcol_7(i__),
     +  Az_B_7(i__),Ze_B_7(i__),Brate_7(i__),e_Brate_7(i__),
     +  Btrm_7(i__),Crate_7(i__),e_Crate_7(i__),Ctrm_7(i__),
     +  Mrate_7(i__),e_Mrate_7(i__),Mtrm_7(i__),M1rate_7(i__),
     +  e_M1rate_7(i__),M1trm_7(i__),M2rate_7(i__),e_M2rate_7(i__),
     +  M2trm_7(i__),Irate_7(i__),e_Irate_7(i__),Itrm_7(i__),
     +  Jrate_7(i__),e_Jrate_7(i__),Jtrm_7(i__),v26rate_7(i__),
     +  e_26rate_7(i__),v26trm_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p25051.dat'	! Observations from flight 25.051 (PASS3)

C  Format for file interpretation

   10 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p25051.dat')
      write(6,*) '....Loading file: p25051.dat'
      do i__=1,1451
        read(1,'(A400)')ar__9
        read(ar__9,10)
     +  Flight_8(i__),time_m_8(i__),time_s_8(i__),Flag_8(i__),
     +  GLON_8(i__),GLAT_8(i__),RAdeg_8(i__),DEdeg_8(i__),Az_8(i__),
     +  Ze_8(i__),Alt_8(i__),Ncol_8(i__),Ocol_8(i__),Arcol_8(i__),
     +  Az_B_8(i__),Ze_B_8(i__),Brate_8(i__),e_Brate_8(i__),
     +  Btrm_8(i__),Crate_8(i__),e_Crate_8(i__),Ctrm_8(i__),
     +  Mrate_8(i__),e_Mrate_8(i__),Mtrm_8(i__),M1rate_8(i__),
     +  e_M1rate_8(i__),M1trm_8(i__),M2rate_8(i__),e_M2rate_8(i__),
     +  M2trm_8(i__),Irate_8(i__),e_Irate_8(i__),Itrm_8(i__),
     +  Jrate_8(i__),e_Jrate_8(i__),Jtrm_8(i__),v26rate_8(i__),
     +  e_26rate_8(i__),v26trm_8(i__)
c    ..............Just test output...........
        write(6,10)
     +  Flight_8(i__),time_m_8(i__),time_s_8(i__),Flag_8(i__),
     +  GLON_8(i__),GLAT_8(i__),RAdeg_8(i__),DEdeg_8(i__),Az_8(i__),
     +  Ze_8(i__),Alt_8(i__),Ncol_8(i__),Ocol_8(i__),Arcol_8(i__),
     +  Az_B_8(i__),Ze_B_8(i__),Brate_8(i__),e_Brate_8(i__),
     +  Btrm_8(i__),Crate_8(i__),e_Crate_8(i__),Ctrm_8(i__),
     +  Mrate_8(i__),e_Mrate_8(i__),Mtrm_8(i__),M1rate_8(i__),
     +  e_M1rate_8(i__),M1trm_8(i__),M2rate_8(i__),e_M2rate_8(i__),
     +  M2trm_8(i__),Irate_8(i__),e_Irate_8(i__),Itrm_8(i__),
     +  Jrate_8(i__),e_Jrate_8(i__),Jtrm_8(i__),v26rate_8(i__),
     +  e_26rate_8(i__),v26trm_8(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'p26061.dat'	! Observations from flight 26.061 (PASS3)

C  Format for file interpretation

   11 format(
     +  4X,F6.3,8X,I2,5X,F5.2,4X,A6,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4,
     +  E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'p26061.dat')
      write(6,*) '....Loading file: p26061.dat'
      do i__=1,1245
        read(1,'(A400)')ar__10
        read(ar__10,11)
     +  Flight_9(i__),time_m_9(i__),time_s_9(i__),Flag_9(i__),
     +  GLON_9(i__),GLAT_9(i__),RAdeg_9(i__),DEdeg_9(i__),Az_9(i__),
     +  Ze_9(i__),Alt_9(i__),Ncol_9(i__),Ocol_9(i__),Arcol_9(i__),
     +  Az_B_9(i__),Ze_B_9(i__),Brate_9(i__),e_Brate_9(i__),
     +  Btrm_9(i__),Crate_9(i__),e_Crate_9(i__),Ctrm_9(i__),
     +  Mrate_9(i__),e_Mrate_9(i__),Mtrm_9(i__),M1rate_9(i__),
     +  e_M1rate_9(i__),M1trm_9(i__),M2rate_9(i__),e_M2rate_9(i__),
     +  M2trm_9(i__),Irate_9(i__),e_Irate_9(i__),Itrm_9(i__),
     +  Jrate_9(i__),e_Jrate_9(i__),Jtrm_9(i__),v26rate_9(i__),
     +  e_26rate_9(i__),v26trm_9(i__)
c    ..............Just test output...........
        write(6,11)
     +  Flight_9(i__),time_m_9(i__),time_s_9(i__),Flag_9(i__),
     +  GLON_9(i__),GLAT_9(i__),RAdeg_9(i__),DEdeg_9(i__),Az_9(i__),
     +  Ze_9(i__),Alt_9(i__),Ncol_9(i__),Ocol_9(i__),Arcol_9(i__),
     +  Az_B_9(i__),Ze_B_9(i__),Brate_9(i__),e_Brate_9(i__),
     +  Btrm_9(i__),Crate_9(i__),e_Crate_9(i__),Ctrm_9(i__),
     +  Mrate_9(i__),e_Mrate_9(i__),Mtrm_9(i__),M1rate_9(i__),
     +  e_M1rate_9(i__),M1trm_9(i__),M2rate_9(i__),e_M2rate_9(i__),
     +  M2trm_9(i__),Irate_9(i__),e_Irate_9(i__),Itrm_9(i__),
     +  Jrate_9(i__),e_Jrate_9(i__),Jtrm_9(i__),v26rate_9(i__),
     +  e_26rate_9(i__),v26trm_9(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'maps.dat'	! Summary of FITS maps

C  Format for file interpretation

   12 format(I3,1X,I2,1X,I1,1X,A12,1X,A61)

C  Effective file loading

      open(unit=1,status='old',file=
     +'maps.dat')
      write(6,*) '....Loading file: maps.dat'
      do i__=1,36
        read(1,'(A83)')ar__11
        read(ar__11,12)Nx(i__),Ny(i__),Nz(i__),FileName(i__),Text(i__)
c    ..............Just test output...........
        write(6,12)Nx(i__),Ny(i__),Nz(i__),FileName(i__),Text(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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