FORTRAN Generation
(/./ftp/cats/I/312)

Conversion of standardized ReadMe file for file /./ftp/cats/I/312 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 2017-Nov-20
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. I/312             PPMX Catalog of positions and proper motions   (Roeser+ 2008)
*================================================================================
*PPM-Extended, a catalogue of positions and proper motions
*     Roeser S., Schilbach E., Schwan H., Kharchenko N.V.,  Piskunov A.E.,
*     Scholz R.-D.
*    <Astron. Astrophys. 488, 401 (2008)>
*    =2008A&A...488..401R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'n0000.dat'	! Declination from +00 to +07.5{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX       (nr__) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX     (nr__) ! [pf] if name is identical (346 cases)
      integer*8     RAmas      (nr__) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas      (nr__) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA       (nr__) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE       (nr__) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA       (nr__) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE       (nr__) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas    (nr__) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas    (nr__) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA     (nr__) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE     (nr__) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag       (nr__) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag       (nr__) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag       (nr__) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag     (nr__) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag       (nr__) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag     (nr__) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag       (nr__) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag     (nr__) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag       (nr__) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag     (nr__) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag       (nr__) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag     (nr__) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs       (nr__) ! ? Number of observation in LSQ fit
      character*1   P          (nr__) ! [P] P if LSQ fit is bad
      character*1   sub        (nr__) ! [SHO] subset flag (1)
      character*1   r_ID       (nr__) ! [HTPSGA] source catalogue flag (2)
      character*12  ID         (nr__) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n0730.dat'	! Declination from +07.5 to +15{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_1     (nr__1) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_1   (nr__1) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_1    (nr__1) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_1    (nr__1) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_1     (nr__1) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_1     (nr__1) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_1     (nr__1) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_1     (nr__1) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_1  (nr__1) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_1  (nr__1) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_1   (nr__1) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_1   (nr__1) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_1     (nr__1) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_1     (nr__1) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_1     (nr__1) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_1   (nr__1) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_1     (nr__1) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_1   (nr__1) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_1     (nr__1) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_1   (nr__1) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_1     (nr__1) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_1   (nr__1) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_1     (nr__1) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_1   (nr__1) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_1     (nr__1) ! ? Number of observation in LSQ fit
      character*1   P_1        (nr__1) ! [P] P if LSQ fit is bad
      character*1   sub_1      (nr__1) ! [SHO] subset flag (1)
      character*1   r_ID_1     (nr__1) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_1       (nr__1) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n1500.dat'	! Declination from +15 to +22.5{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_2     (nr__2) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_2   (nr__2) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_2    (nr__2) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_2    (nr__2) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_2     (nr__2) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_2     (nr__2) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_2     (nr__2) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_2     (nr__2) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_2  (nr__2) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_2  (nr__2) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_2   (nr__2) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_2   (nr__2) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_2     (nr__2) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_2     (nr__2) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_2     (nr__2) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_2   (nr__2) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_2     (nr__2) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_2   (nr__2) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_2     (nr__2) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_2   (nr__2) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_2     (nr__2) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_2   (nr__2) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_2     (nr__2) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_2   (nr__2) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_2     (nr__2) ! ? Number of observation in LSQ fit
      character*1   P_2        (nr__2) ! [P] P if LSQ fit is bad
      character*1   sub_2      (nr__2) ! [SHO] subset flag (1)
      character*1   r_ID_2     (nr__2) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_2       (nr__2) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n2230.dat'	! Declination from +22.5 to +30{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_3     (nr__3) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_3   (nr__3) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_3    (nr__3) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_3    (nr__3) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_3     (nr__3) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_3     (nr__3) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_3     (nr__3) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_3     (nr__3) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_3  (nr__3) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_3  (nr__3) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_3   (nr__3) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_3   (nr__3) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_3     (nr__3) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_3     (nr__3) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_3     (nr__3) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_3   (nr__3) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_3     (nr__3) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_3   (nr__3) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_3     (nr__3) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_3   (nr__3) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_3     (nr__3) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_3   (nr__3) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_3     (nr__3) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_3   (nr__3) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_3     (nr__3) ! ? Number of observation in LSQ fit
      character*1   P_3        (nr__3) ! [P] P if LSQ fit is bad
      character*1   sub_3      (nr__3) ! [SHO] subset flag (1)
      character*1   r_ID_3     (nr__3) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_3       (nr__3) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n3000.dat'	! Declination from +30 to +37.5{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_4     (nr__4) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_4   (nr__4) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_4    (nr__4) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_4    (nr__4) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_4     (nr__4) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_4     (nr__4) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_4     (nr__4) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_4     (nr__4) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_4  (nr__4) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_4  (nr__4) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_4   (nr__4) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_4   (nr__4) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_4     (nr__4) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_4     (nr__4) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_4     (nr__4) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_4   (nr__4) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_4     (nr__4) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_4   (nr__4) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_4     (nr__4) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_4   (nr__4) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_4     (nr__4) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_4   (nr__4) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_4     (nr__4) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_4   (nr__4) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_4     (nr__4) ! ? Number of observation in LSQ fit
      character*1   P_4        (nr__4) ! [P] P if LSQ fit is bad
      character*1   sub_4      (nr__4) ! [SHO] subset flag (1)
      character*1   r_ID_4     (nr__4) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_4       (nr__4) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n3730.dat'	! Declination from +37.5 to +45{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_5     (nr__5) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_5   (nr__5) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_5    (nr__5) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_5    (nr__5) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_5     (nr__5) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_5     (nr__5) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_5     (nr__5) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_5     (nr__5) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_5  (nr__5) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_5  (nr__5) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_5   (nr__5) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_5   (nr__5) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_5     (nr__5) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_5     (nr__5) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_5     (nr__5) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_5   (nr__5) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_5     (nr__5) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_5   (nr__5) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_5     (nr__5) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_5   (nr__5) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_5     (nr__5) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_5   (nr__5) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_5     (nr__5) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_5   (nr__5) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_5     (nr__5) ! ? Number of observation in LSQ fit
      character*1   P_5        (nr__5) ! [P] P if LSQ fit is bad
      character*1   sub_5      (nr__5) ! [SHO] subset flag (1)
      character*1   r_ID_5     (nr__5) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_5       (nr__5) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n4500.dat'	! Declination from +45 to +52.5{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_6     (nr__6) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_6   (nr__6) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_6    (nr__6) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_6    (nr__6) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_6     (nr__6) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_6     (nr__6) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_6     (nr__6) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_6     (nr__6) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_6  (nr__6) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_6  (nr__6) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_6   (nr__6) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_6   (nr__6) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_6     (nr__6) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_6     (nr__6) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_6     (nr__6) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_6   (nr__6) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_6     (nr__6) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_6   (nr__6) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_6     (nr__6) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_6   (nr__6) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_6     (nr__6) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_6   (nr__6) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_6     (nr__6) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_6   (nr__6) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_6     (nr__6) ! ? Number of observation in LSQ fit
      character*1   P_6        (nr__6) ! [P] P if LSQ fit is bad
      character*1   sub_6      (nr__6) ! [SHO] subset flag (1)
      character*1   r_ID_6     (nr__6) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_6       (nr__6) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n5230.dat'	! Declination from +52.5 to +60{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_7     (nr__7) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_7   (nr__7) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_7    (nr__7) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_7    (nr__7) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_7     (nr__7) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_7     (nr__7) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_7     (nr__7) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_7     (nr__7) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_7  (nr__7) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_7  (nr__7) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_7   (nr__7) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_7   (nr__7) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_7     (nr__7) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_7     (nr__7) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_7     (nr__7) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_7   (nr__7) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_7     (nr__7) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_7   (nr__7) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_7     (nr__7) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_7   (nr__7) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_7     (nr__7) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_7   (nr__7) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_7     (nr__7) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_7   (nr__7) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_7     (nr__7) ! ? Number of observation in LSQ fit
      character*1   P_7        (nr__7) ! [P] P if LSQ fit is bad
      character*1   sub_7      (nr__7) ! [SHO] subset flag (1)
      character*1   r_ID_7     (nr__7) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_7       (nr__7) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n6000.dat'	! Declination from +60 to +67.5{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_8     (nr__8) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_8   (nr__8) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_8    (nr__8) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_8    (nr__8) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_8     (nr__8) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_8     (nr__8) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_8     (nr__8) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_8     (nr__8) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_8  (nr__8) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_8  (nr__8) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_8   (nr__8) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_8   (nr__8) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_8     (nr__8) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_8     (nr__8) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_8     (nr__8) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_8   (nr__8) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_8     (nr__8) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_8   (nr__8) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_8     (nr__8) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_8   (nr__8) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_8     (nr__8) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_8   (nr__8) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_8     (nr__8) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_8   (nr__8) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_8     (nr__8) ! ? Number of observation in LSQ fit
      character*1   P_8        (nr__8) ! [P] P if LSQ fit is bad
      character*1   sub_8      (nr__8) ! [SHO] subset flag (1)
      character*1   r_ID_8     (nr__8) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_8       (nr__8) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n6730.dat'	! Declination from +67.5 to +75{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_9     (nr__9) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_9   (nr__9) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_9    (nr__9) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_9    (nr__9) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_9     (nr__9) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_9     (nr__9) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_9     (nr__9) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_9     (nr__9) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_9  (nr__9) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_9  (nr__9) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_9   (nr__9) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_9   (nr__9) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_9     (nr__9) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_9     (nr__9) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_9     (nr__9) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_9   (nr__9) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_9     (nr__9) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_9   (nr__9) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_9     (nr__9) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_9   (nr__9) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_9     (nr__9) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_9   (nr__9) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_9     (nr__9) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_9   (nr__9) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_9     (nr__9) ! ? Number of observation in LSQ fit
      character*1   P_9        (nr__9) ! [P] P if LSQ fit is bad
      character*1   sub_9      (nr__9) ! [SHO] subset flag (1)
      character*1   r_ID_9     (nr__9) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_9       (nr__9) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n7500.dat'	! Declination from +75 to +82.5{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_10    (nr__10) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_10  (nr__10) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_10   (nr__10) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_10   (nr__10) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_10    (nr__10) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_10    (nr__10) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_10    (nr__10) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_10    (nr__10) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_10 (nr__10) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_10 (nr__10) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_10  (nr__10) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_10  (nr__10) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_10    (nr__10) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_10    (nr__10) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_10    (nr__10) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_10  (nr__10) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_10    (nr__10) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_10  (nr__10) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_10    (nr__10) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_10  (nr__10) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_10    (nr__10) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_10  (nr__10) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_10    (nr__10) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_10  (nr__10) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_10    (nr__10) ! ? Number of observation in LSQ fit
      character*1   P_10       (nr__10) ! [P] P if LSQ fit is bad
      character*1   sub_10     (nr__10) ! [SHO] subset flag (1)
      character*1   r_ID_10    (nr__10) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_10      (nr__10) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 'n8230.dat'	! Declination from +82.5 to +90{deg}

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

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_11    (nr__11) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_11  (nr__11) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_11   (nr__11) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_11   (nr__11) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_11    (nr__11) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_11    (nr__11) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_11    (nr__11) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_11    (nr__11) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_11 (nr__11) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_11 (nr__11) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_11  (nr__11) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_11  (nr__11) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_11    (nr__11) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_11    (nr__11) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_11    (nr__11) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_11  (nr__11) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_11    (nr__11) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_11  (nr__11) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_11    (nr__11) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_11  (nr__11) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_11    (nr__11) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_11  (nr__11) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_11    (nr__11) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_11  (nr__11) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_11    (nr__11) ! ? Number of observation in LSQ fit
      character*1   P_11       (nr__11) ! [P] P if LSQ fit is bad
      character*1   sub_11     (nr__11) ! [SHO] subset flag (1)
      character*1   r_ID_11    (nr__11) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_11      (nr__11) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's0000.dat'	! Declination from -00 to -07.5{deg}

      integer*4 nr__12
      parameter (nr__12=970673)	! Number of records
      character*166 ar__12 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_12    (nr__12) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_12  (nr__12) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_12   (nr__12) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_12   (nr__12) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_12    (nr__12) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_12    (nr__12) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_12    (nr__12) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_12    (nr__12) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_12 (nr__12) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_12 (nr__12) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_12  (nr__12) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_12  (nr__12) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_12    (nr__12) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_12    (nr__12) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_12    (nr__12) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_12  (nr__12) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_12    (nr__12) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_12  (nr__12) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_12    (nr__12) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_12  (nr__12) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_12    (nr__12) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_12  (nr__12) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_12    (nr__12) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_12  (nr__12) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_12    (nr__12) ! ? Number of observation in LSQ fit
      character*1   P_12       (nr__12) ! [P] P if LSQ fit is bad
      character*1   sub_12     (nr__12) ! [SHO] subset flag (1)
      character*1   r_ID_12    (nr__12) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_12      (nr__12) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's0730.dat'	! Declination from -07.5 to -15{deg}

      integer*4 nr__13
      parameter (nr__13=970094)	! Number of records
      character*166 ar__13 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_13    (nr__13) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_13  (nr__13) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_13   (nr__13) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_13   (nr__13) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_13    (nr__13) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_13    (nr__13) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_13    (nr__13) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_13    (nr__13) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_13 (nr__13) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_13 (nr__13) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_13  (nr__13) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_13  (nr__13) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_13    (nr__13) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_13    (nr__13) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_13    (nr__13) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_13  (nr__13) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_13    (nr__13) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_13  (nr__13) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_13    (nr__13) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_13  (nr__13) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_13    (nr__13) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_13  (nr__13) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_13    (nr__13) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_13  (nr__13) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_13    (nr__13) ! ? Number of observation in LSQ fit
      character*1   P_13       (nr__13) ! [P] P if LSQ fit is bad
      character*1   sub_13     (nr__13) ! [SHO] subset flag (1)
      character*1   r_ID_13    (nr__13) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_13      (nr__13) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's1500.dat'	! Declination from -15 to -22.5{deg}

      integer*4 nr__14
      parameter (nr__14=1010372)	! Number of records
      character*166 ar__14 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_14    (nr__14) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_14  (nr__14) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_14   (nr__14) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_14   (nr__14) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_14    (nr__14) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_14    (nr__14) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_14    (nr__14) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_14    (nr__14) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_14 (nr__14) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_14 (nr__14) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_14  (nr__14) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_14  (nr__14) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_14    (nr__14) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_14    (nr__14) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_14    (nr__14) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_14  (nr__14) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_14    (nr__14) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_14  (nr__14) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_14    (nr__14) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_14  (nr__14) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_14    (nr__14) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_14  (nr__14) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_14    (nr__14) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_14  (nr__14) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_14    (nr__14) ! ? Number of observation in LSQ fit
      character*1   P_14       (nr__14) ! [P] P if LSQ fit is bad
      character*1   sub_14     (nr__14) ! [SHO] subset flag (1)
      character*1   r_ID_14    (nr__14) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_14      (nr__14) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's2230.dat'	! Declination from -22.5 to -30{deg}

      integer*4 nr__15
      parameter (nr__15=1019597)	! Number of records
      character*166 ar__15 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_15    (nr__15) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_15  (nr__15) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_15   (nr__15) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_15   (nr__15) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_15    (nr__15) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_15    (nr__15) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_15    (nr__15) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_15    (nr__15) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_15 (nr__15) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_15 (nr__15) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_15  (nr__15) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_15  (nr__15) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_15    (nr__15) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_15    (nr__15) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_15    (nr__15) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_15  (nr__15) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_15    (nr__15) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_15  (nr__15) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_15    (nr__15) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_15  (nr__15) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_15    (nr__15) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_15  (nr__15) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_15    (nr__15) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_15  (nr__15) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_15    (nr__15) ! ? Number of observation in LSQ fit
      character*1   P_15       (nr__15) ! [P] P if LSQ fit is bad
      character*1   sub_15     (nr__15) ! [SHO] subset flag (1)
      character*1   r_ID_15    (nr__15) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_15      (nr__15) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's3000.dat'	! Declination from -30 to -37.5{deg}

      integer*4 nr__16
      parameter (nr__16=975900)	! Number of records
      character*166 ar__16 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_16    (nr__16) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_16  (nr__16) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_16   (nr__16) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_16   (nr__16) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_16    (nr__16) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_16    (nr__16) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_16    (nr__16) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_16    (nr__16) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_16 (nr__16) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_16 (nr__16) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_16  (nr__16) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_16  (nr__16) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_16    (nr__16) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_16    (nr__16) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_16    (nr__16) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_16  (nr__16) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_16    (nr__16) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_16  (nr__16) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_16    (nr__16) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_16  (nr__16) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_16    (nr__16) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_16  (nr__16) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_16    (nr__16) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_16  (nr__16) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_16    (nr__16) ! ? Number of observation in LSQ fit
      character*1   P_16       (nr__16) ! [P] P if LSQ fit is bad
      character*1   sub_16     (nr__16) ! [SHO] subset flag (1)
      character*1   r_ID_16    (nr__16) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_16      (nr__16) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's3730.dat'	! Declination from -37.5 to -45{deg}

      integer*4 nr__17
      parameter (nr__17=991495)	! Number of records
      character*166 ar__17 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_17    (nr__17) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_17  (nr__17) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_17   (nr__17) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_17   (nr__17) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_17    (nr__17) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_17    (nr__17) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_17    (nr__17) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_17    (nr__17) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_17 (nr__17) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_17 (nr__17) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_17  (nr__17) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_17  (nr__17) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_17    (nr__17) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_17    (nr__17) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_17    (nr__17) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_17  (nr__17) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_17    (nr__17) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_17  (nr__17) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_17    (nr__17) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_17  (nr__17) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_17    (nr__17) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_17  (nr__17) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_17    (nr__17) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_17  (nr__17) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_17    (nr__17) ! ? Number of observation in LSQ fit
      character*1   P_17       (nr__17) ! [P] P if LSQ fit is bad
      character*1   sub_17     (nr__17) ! [SHO] subset flag (1)
      character*1   r_ID_17    (nr__17) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_17      (nr__17) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's4500.dat'	! Declination from -45 to -52.5{deg}

      integer*4 nr__18
      parameter (nr__18=951351)	! Number of records
      character*166 ar__18 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_18    (nr__18) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_18  (nr__18) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_18   (nr__18) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_18   (nr__18) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_18    (nr__18) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_18    (nr__18) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_18    (nr__18) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_18    (nr__18) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_18 (nr__18) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_18 (nr__18) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_18  (nr__18) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_18  (nr__18) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_18    (nr__18) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_18    (nr__18) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_18    (nr__18) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_18  (nr__18) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_18    (nr__18) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_18  (nr__18) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_18    (nr__18) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_18  (nr__18) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_18    (nr__18) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_18  (nr__18) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_18    (nr__18) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_18  (nr__18) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_18    (nr__18) ! ? Number of observation in LSQ fit
      character*1   P_18       (nr__18) ! [P] P if LSQ fit is bad
      character*1   sub_18     (nr__18) ! [SHO] subset flag (1)
      character*1   r_ID_18    (nr__18) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_18      (nr__18) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's5230.dat'	! Declination from -52.5 to -60{deg}

      integer*4 nr__19
      parameter (nr__19=782446)	! Number of records
      character*166 ar__19 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_19    (nr__19) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_19  (nr__19) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_19   (nr__19) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_19   (nr__19) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_19    (nr__19) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_19    (nr__19) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_19    (nr__19) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_19    (nr__19) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_19 (nr__19) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_19 (nr__19) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_19  (nr__19) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_19  (nr__19) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_19    (nr__19) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_19    (nr__19) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_19    (nr__19) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_19  (nr__19) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_19    (nr__19) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_19  (nr__19) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_19    (nr__19) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_19  (nr__19) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_19    (nr__19) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_19  (nr__19) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_19    (nr__19) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_19  (nr__19) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_19    (nr__19) ! ? Number of observation in LSQ fit
      character*1   P_19       (nr__19) ! [P] P if LSQ fit is bad
      character*1   sub_19     (nr__19) ! [SHO] subset flag (1)
      character*1   r_ID_19    (nr__19) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_19      (nr__19) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's6000.dat'	! Declination from -60 to -67.5{deg}

      integer*4 nr__20
      parameter (nr__20=640425)	! Number of records
      character*166 ar__20 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_20    (nr__20) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_20  (nr__20) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_20   (nr__20) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_20   (nr__20) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_20    (nr__20) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_20    (nr__20) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_20    (nr__20) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_20    (nr__20) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_20 (nr__20) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_20 (nr__20) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_20  (nr__20) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_20  (nr__20) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_20    (nr__20) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_20    (nr__20) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_20    (nr__20) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_20  (nr__20) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_20    (nr__20) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_20  (nr__20) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_20    (nr__20) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_20  (nr__20) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_20    (nr__20) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_20  (nr__20) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_20    (nr__20) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_20  (nr__20) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_20    (nr__20) ! ? Number of observation in LSQ fit
      character*1   P_20       (nr__20) ! [P] P if LSQ fit is bad
      character*1   sub_20     (nr__20) ! [SHO] subset flag (1)
      character*1   r_ID_20    (nr__20) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_20      (nr__20) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's6730.dat'	! Declination from -67.5 to -75{deg}

      integer*4 nr__21
      parameter (nr__21=452628)	! Number of records
      character*166 ar__21 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_21    (nr__21) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_21  (nr__21) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_21   (nr__21) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_21   (nr__21) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_21    (nr__21) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_21    (nr__21) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_21    (nr__21) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_21    (nr__21) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_21 (nr__21) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_21 (nr__21) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_21  (nr__21) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_21  (nr__21) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_21    (nr__21) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_21    (nr__21) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_21    (nr__21) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_21  (nr__21) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_21    (nr__21) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_21  (nr__21) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_21    (nr__21) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_21  (nr__21) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_21    (nr__21) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_21  (nr__21) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_21    (nr__21) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_21  (nr__21) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_21    (nr__21) ! ? Number of observation in LSQ fit
      character*1   P_21       (nr__21) ! [P] P if LSQ fit is bad
      character*1   sub_21     (nr__21) ! [SHO] subset flag (1)
      character*1   r_ID_21    (nr__21) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_21      (nr__21) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's7500.dat'	! Declination from -75 to -82.5{deg}

      integer*4 nr__22
      parameter (nr__22=269692)	! Number of records
      character*166 ar__22 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_22    (nr__22) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_22  (nr__22) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_22   (nr__22) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_22   (nr__22) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_22    (nr__22) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_22    (nr__22) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_22    (nr__22) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_22    (nr__22) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_22 (nr__22) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_22 (nr__22) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_22  (nr__22) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_22  (nr__22) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_22    (nr__22) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_22    (nr__22) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_22    (nr__22) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_22  (nr__22) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_22    (nr__22) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_22  (nr__22) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_22    (nr__22) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_22  (nr__22) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_22    (nr__22) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_22  (nr__22) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_22    (nr__22) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_22  (nr__22) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_22    (nr__22) ! ? Number of observation in LSQ fit
      character*1   P_22       (nr__22) ! [P] P if LSQ fit is bad
      character*1   sub_22     (nr__22) ! [SHO] subset flag (1)
      character*1   r_ID_22    (nr__22) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_22      (nr__22) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Declarations for 's8230.dat'	! Declination from -82.5 to -90{deg}

      integer*4 nr__23
      parameter (nr__23=92360)	! Number of records
      character*166 ar__23 	! Full-size record

C  J2000.0 position composed of: RAmas DEmas (Epoch=2000.0)
      character*15  PPMX_23    (nr__23) ! Name (IAU convention HHMMSS.S+DDMMSS)
      character*1   m_PPMX_23  (nr__23) ! [pf] if name is identical (346 cases)
      integer*8     RAmas_23   (nr__23) ! (mas) Right Ascension J2000.0, epoch 2000.0
      integer*8     DEmas_23   (nr__23) ! (mas) Declination J2000.0, epoch 2000.0
      integer*4     pmRA_23    (nr__23) ! (0.01mas/yr) Proper Motion in RA*cos(DEmas) (mas/100yr)
      integer*4     pmDE_23    (nr__23) ! (0.01mas/yr) Proper Motion in DE (mas/100yr)
      integer*4     epRA_23    (nr__23) ! (0.01yr) Mean Epoch (RA) - 2000.00
      integer*4     epDE_23    (nr__23) ! (0.01yr) Mean Epoch (DE) - 2000.00
      integer*4     e_RAmas_23 (nr__23) ! (mas) Mean error in RA*cos(DEmas) at epRA
      integer*4     e_DEmas_23 (nr__23) ! (mas) Mean error in DE at epDE
      integer*4     e_pmRA_23  (nr__23) ! (0.1mas/yr) Mean error in pmRA*cos(DEmas) (mas/10yr)
      integer*4     e_pmDE_23  (nr__23) ! (0.1mas/yr) Mean error in pmDE (mas/10yr)
      integer*4     Cmag_23    (nr__23) ! (mmag) ?=-9999 Catalogue magnitude from source
      integer*4     Rmag_23    (nr__23) ! (mmag) ?=-9999 Calculated Ru magnitude (579-642nm)
*                                            from source
      integer*4     Bmag_23    (nr__23) ! (mmag) ?=-9999 B magnitude in Johnson system
      integer*4     e_Bmag_23  (nr__23) ! (mmag) ?=-99 Standard error on B magnitude
      integer*4     Vmag_23    (nr__23) ! (mmag) ?=-9999 V magnitude in Johnson system
      integer*4     e_Vmag_23  (nr__23) ! (mmag) ?=-99 Standard error on V magnitude
      integer*4     Jmag_23    (nr__23) ! (mmag) ?=-9999 J magnitude from 2MASS
      integer*4     e_Jmag_23  (nr__23) ! (mmag) ?=-99 Standard error on J magnitude
      integer*4     Hmag_23    (nr__23) ! (mmag) ?=-9999 H magnitude from 2MASS
      integer*4     e_Hmag_23  (nr__23) ! (mmag) ?=-99 Standard error on H magnitude
      integer*4     Kmag_23    (nr__23) ! (mmag) ?=-9999 K magnitude from 2MASS
      integer*4     e_Kmag_23  (nr__23) ! (mmag) ?=-99 Standard error on K magnitude
      integer*4     Nobs_23    (nr__23) ! ? Number of observation in LSQ fit
      character*1   P_23       (nr__23) ! [P] P if LSQ fit is bad
      character*1   sub_23     (nr__23) ! [SHO] subset flag (1)
      character*1   r_ID_23    (nr__23) ! [HTPSGA] source catalogue flag (2)
      character*12  ID_23      (nr__23) ! star identifier in the source catalogue
*Note (1): PPMX consist of 3 parts designated as follows:
*    S = a survey S of stars brighter than R=12.8, consisting of 5,620,114 stars
*    H = 874,934 stars fainter than the survey limit have measurements in the AC
*        and therefore highly accurate proper motions
*    O = all other stars (11,593,871).
*Note (2): These columns give the star identifiers in the following catalogues:
*      H = the ARIHIP catalogue  (Wielen et al., 2001VeARI..40....1W)
*      T = the Tycho-2 catalogue (Hoeg et al. 2000, Cat. I/259)
*      P = the PPM catalogue (Roeser & Bastian, 1991, Cat. I/146, I/193)
*      S = STARNET catalogue (Roeser 1996, IAU Symp. 172, 481)
*      G = the Guide Star Catalog 1.2 (Morrison et al., 2001, Cat. I/254)
*      A = the AC2000.0 catalogue  (Urban et al., 2001, Cat. I/275)

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

C  Loading file 'n0000.dat'	! Declination from +00 to +07.5{deg}

C  Format for file interpretation

    1 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n0000.dat')
      write(6,*) '....Loading file: n0000.dat'
      do i__=1,1131366
        read(1,'(A166)')ar__
        read(ar__,1)
     +  PPMX(i__),m_PPMX(i__),RAmas(i__),DEmas(i__),pmRA(i__),
     +  pmDE(i__),epRA(i__),epDE(i__),e_RAmas(i__),e_DEmas(i__),
     +  e_pmRA(i__),e_pmDE(i__),Cmag(i__),Rmag(i__),Bmag(i__),
     +  e_Bmag(i__),Vmag(i__),e_Vmag(i__),Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),Nobs(i__),P(i__),
     +  sub(i__),r_ID(i__),ID(i__)
        if(ar__(146:147) .EQ. '') Nobs(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  PPMX(i__),m_PPMX(i__),RAmas(i__),DEmas(i__),pmRA(i__),
     +  pmDE(i__),epRA(i__),epDE(i__),e_RAmas(i__),e_DEmas(i__),
     +  e_pmRA(i__),e_pmDE(i__),Cmag(i__),Rmag(i__),Bmag(i__),
     +  e_Bmag(i__),Vmag(i__),e_Vmag(i__),Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),Nobs(i__),P(i__),
     +  sub(i__),r_ID(i__),ID(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n0730.dat'	! Declination from +07.5 to +15{deg}

C  Format for file interpretation

    2 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n0730.dat')
      write(6,*) '....Loading file: n0730.dat'
      do i__=1,1040650
        read(1,'(A166)')ar__1
        read(ar__1,2)
     +  PPMX_1(i__),m_PPMX_1(i__),RAmas_1(i__),DEmas_1(i__),
     +  pmRA_1(i__),pmDE_1(i__),epRA_1(i__),epDE_1(i__),
     +  e_RAmas_1(i__),e_DEmas_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),
     +  Cmag_1(i__),Rmag_1(i__),Bmag_1(i__),e_Bmag_1(i__),Vmag_1(i__),
     +  e_Vmag_1(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),
     +  e_Hmag_1(i__),Kmag_1(i__),e_Kmag_1(i__),Nobs_1(i__),P_1(i__),
     +  sub_1(i__),r_ID_1(i__),ID_1(i__)
        if(ar__1(146:147) .EQ. '') Nobs_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  PPMX_1(i__),m_PPMX_1(i__),RAmas_1(i__),DEmas_1(i__),
     +  pmRA_1(i__),pmDE_1(i__),epRA_1(i__),epDE_1(i__),
     +  e_RAmas_1(i__),e_DEmas_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),
     +  Cmag_1(i__),Rmag_1(i__),Bmag_1(i__),e_Bmag_1(i__),Vmag_1(i__),
     +  e_Vmag_1(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),
     +  e_Hmag_1(i__),Kmag_1(i__),e_Kmag_1(i__),Nobs_1(i__),P_1(i__),
     +  sub_1(i__),r_ID_1(i__),ID_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n1500.dat'	! Declination from +15 to +22.5{deg}

C  Format for file interpretation

    3 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n1500.dat')
      write(6,*) '....Loading file: n1500.dat'
      do i__=1,1005212
        read(1,'(A166)')ar__2
        read(ar__2,3)
     +  PPMX_2(i__),m_PPMX_2(i__),RAmas_2(i__),DEmas_2(i__),
     +  pmRA_2(i__),pmDE_2(i__),epRA_2(i__),epDE_2(i__),
     +  e_RAmas_2(i__),e_DEmas_2(i__),e_pmRA_2(i__),e_pmDE_2(i__),
     +  Cmag_2(i__),Rmag_2(i__),Bmag_2(i__),e_Bmag_2(i__),Vmag_2(i__),
     +  e_Vmag_2(i__),Jmag_2(i__),e_Jmag_2(i__),Hmag_2(i__),
     +  e_Hmag_2(i__),Kmag_2(i__),e_Kmag_2(i__),Nobs_2(i__),P_2(i__),
     +  sub_2(i__),r_ID_2(i__),ID_2(i__)
        if(ar__2(146:147) .EQ. '') Nobs_2(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  PPMX_2(i__),m_PPMX_2(i__),RAmas_2(i__),DEmas_2(i__),
     +  pmRA_2(i__),pmDE_2(i__),epRA_2(i__),epDE_2(i__),
     +  e_RAmas_2(i__),e_DEmas_2(i__),e_pmRA_2(i__),e_pmDE_2(i__),
     +  Cmag_2(i__),Rmag_2(i__),Bmag_2(i__),e_Bmag_2(i__),Vmag_2(i__),
     +  e_Vmag_2(i__),Jmag_2(i__),e_Jmag_2(i__),Hmag_2(i__),
     +  e_Hmag_2(i__),Kmag_2(i__),e_Kmag_2(i__),Nobs_2(i__),P_2(i__),
     +  sub_2(i__),r_ID_2(i__),ID_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n2230.dat'	! Declination from +22.5 to +30{deg}

C  Format for file interpretation

    4 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n2230.dat')
      write(6,*) '....Loading file: n2230.dat'
      do i__=1,1040381
        read(1,'(A166)')ar__3
        read(ar__3,4)
     +  PPMX_3(i__),m_PPMX_3(i__),RAmas_3(i__),DEmas_3(i__),
     +  pmRA_3(i__),pmDE_3(i__),epRA_3(i__),epDE_3(i__),
     +  e_RAmas_3(i__),e_DEmas_3(i__),e_pmRA_3(i__),e_pmDE_3(i__),
     +  Cmag_3(i__),Rmag_3(i__),Bmag_3(i__),e_Bmag_3(i__),Vmag_3(i__),
     +  e_Vmag_3(i__),Jmag_3(i__),e_Jmag_3(i__),Hmag_3(i__),
     +  e_Hmag_3(i__),Kmag_3(i__),e_Kmag_3(i__),Nobs_3(i__),P_3(i__),
     +  sub_3(i__),r_ID_3(i__),ID_3(i__)
        if(ar__3(146:147) .EQ. '') Nobs_3(i__) = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  PPMX_3(i__),m_PPMX_3(i__),RAmas_3(i__),DEmas_3(i__),
     +  pmRA_3(i__),pmDE_3(i__),epRA_3(i__),epDE_3(i__),
     +  e_RAmas_3(i__),e_DEmas_3(i__),e_pmRA_3(i__),e_pmDE_3(i__),
     +  Cmag_3(i__),Rmag_3(i__),Bmag_3(i__),e_Bmag_3(i__),Vmag_3(i__),
     +  e_Vmag_3(i__),Jmag_3(i__),e_Jmag_3(i__),Hmag_3(i__),
     +  e_Hmag_3(i__),Kmag_3(i__),e_Kmag_3(i__),Nobs_3(i__),P_3(i__),
     +  sub_3(i__),r_ID_3(i__),ID_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n3000.dat'	! Declination from +30 to +37.5{deg}

C  Format for file interpretation

    5 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n3000.dat')
      write(6,*) '....Loading file: n3000.dat'
      do i__=1,1020072
        read(1,'(A166)')ar__4
        read(ar__4,5)
     +  PPMX_4(i__),m_PPMX_4(i__),RAmas_4(i__),DEmas_4(i__),
     +  pmRA_4(i__),pmDE_4(i__),epRA_4(i__),epDE_4(i__),
     +  e_RAmas_4(i__),e_DEmas_4(i__),e_pmRA_4(i__),e_pmDE_4(i__),
     +  Cmag_4(i__),Rmag_4(i__),Bmag_4(i__),e_Bmag_4(i__),Vmag_4(i__),
     +  e_Vmag_4(i__),Jmag_4(i__),e_Jmag_4(i__),Hmag_4(i__),
     +  e_Hmag_4(i__),Kmag_4(i__),e_Kmag_4(i__),Nobs_4(i__),P_4(i__),
     +  sub_4(i__),r_ID_4(i__),ID_4(i__)
        if(ar__4(146:147) .EQ. '') Nobs_4(i__) = iNULL__
c    ..............Just test output...........
        write(6,5)
     +  PPMX_4(i__),m_PPMX_4(i__),RAmas_4(i__),DEmas_4(i__),
     +  pmRA_4(i__),pmDE_4(i__),epRA_4(i__),epDE_4(i__),
     +  e_RAmas_4(i__),e_DEmas_4(i__),e_pmRA_4(i__),e_pmDE_4(i__),
     +  Cmag_4(i__),Rmag_4(i__),Bmag_4(i__),e_Bmag_4(i__),Vmag_4(i__),
     +  e_Vmag_4(i__),Jmag_4(i__),e_Jmag_4(i__),Hmag_4(i__),
     +  e_Hmag_4(i__),Kmag_4(i__),e_Kmag_4(i__),Nobs_4(i__),P_4(i__),
     +  sub_4(i__),r_ID_4(i__),ID_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n3730.dat'	! Declination from +37.5 to +45{deg}

C  Format for file interpretation

    6 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n3730.dat')
      write(6,*) '....Loading file: n3730.dat'
      do i__=1,982895
        read(1,'(A166)')ar__5
        read(ar__5,6)
     +  PPMX_5(i__),m_PPMX_5(i__),RAmas_5(i__),DEmas_5(i__),
     +  pmRA_5(i__),pmDE_5(i__),epRA_5(i__),epDE_5(i__),
     +  e_RAmas_5(i__),e_DEmas_5(i__),e_pmRA_5(i__),e_pmDE_5(i__),
     +  Cmag_5(i__),Rmag_5(i__),Bmag_5(i__),e_Bmag_5(i__),Vmag_5(i__),
     +  e_Vmag_5(i__),Jmag_5(i__),e_Jmag_5(i__),Hmag_5(i__),
     +  e_Hmag_5(i__),Kmag_5(i__),e_Kmag_5(i__),Nobs_5(i__),P_5(i__),
     +  sub_5(i__),r_ID_5(i__),ID_5(i__)
        if(ar__5(146:147) .EQ. '') Nobs_5(i__) = iNULL__
c    ..............Just test output...........
        write(6,6)
     +  PPMX_5(i__),m_PPMX_5(i__),RAmas_5(i__),DEmas_5(i__),
     +  pmRA_5(i__),pmDE_5(i__),epRA_5(i__),epDE_5(i__),
     +  e_RAmas_5(i__),e_DEmas_5(i__),e_pmRA_5(i__),e_pmDE_5(i__),
     +  Cmag_5(i__),Rmag_5(i__),Bmag_5(i__),e_Bmag_5(i__),Vmag_5(i__),
     +  e_Vmag_5(i__),Jmag_5(i__),e_Jmag_5(i__),Hmag_5(i__),
     +  e_Hmag_5(i__),Kmag_5(i__),e_Kmag_5(i__),Nobs_5(i__),P_5(i__),
     +  sub_5(i__),r_ID_5(i__),ID_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n4500.dat'	! Declination from +45 to +52.5{deg}

C  Format for file interpretation

    7 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n4500.dat')
      write(6,*) '....Loading file: n4500.dat'
      do i__=1,863340
        read(1,'(A166)')ar__6
        read(ar__6,7)
     +  PPMX_6(i__),m_PPMX_6(i__),RAmas_6(i__),DEmas_6(i__),
     +  pmRA_6(i__),pmDE_6(i__),epRA_6(i__),epDE_6(i__),
     +  e_RAmas_6(i__),e_DEmas_6(i__),e_pmRA_6(i__),e_pmDE_6(i__),
     +  Cmag_6(i__),Rmag_6(i__),Bmag_6(i__),e_Bmag_6(i__),Vmag_6(i__),
     +  e_Vmag_6(i__),Jmag_6(i__),e_Jmag_6(i__),Hmag_6(i__),
     +  e_Hmag_6(i__),Kmag_6(i__),e_Kmag_6(i__),Nobs_6(i__),P_6(i__),
     +  sub_6(i__),r_ID_6(i__),ID_6(i__)
        if(ar__6(146:147) .EQ. '') Nobs_6(i__) = iNULL__
c    ..............Just test output...........
        write(6,7)
     +  PPMX_6(i__),m_PPMX_6(i__),RAmas_6(i__),DEmas_6(i__),
     +  pmRA_6(i__),pmDE_6(i__),epRA_6(i__),epDE_6(i__),
     +  e_RAmas_6(i__),e_DEmas_6(i__),e_pmRA_6(i__),e_pmDE_6(i__),
     +  Cmag_6(i__),Rmag_6(i__),Bmag_6(i__),e_Bmag_6(i__),Vmag_6(i__),
     +  e_Vmag_6(i__),Jmag_6(i__),e_Jmag_6(i__),Hmag_6(i__),
     +  e_Hmag_6(i__),Kmag_6(i__),e_Kmag_6(i__),Nobs_6(i__),P_6(i__),
     +  sub_6(i__),r_ID_6(i__),ID_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n5230.dat'	! Declination from +52.5 to +60{deg}

C  Format for file interpretation

    8 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n5230.dat')
      write(6,*) '....Loading file: n5230.dat'
      do i__=1,663079
        read(1,'(A166)')ar__7
        read(ar__7,8)
     +  PPMX_7(i__),m_PPMX_7(i__),RAmas_7(i__),DEmas_7(i__),
     +  pmRA_7(i__),pmDE_7(i__),epRA_7(i__),epDE_7(i__),
     +  e_RAmas_7(i__),e_DEmas_7(i__),e_pmRA_7(i__),e_pmDE_7(i__),
     +  Cmag_7(i__),Rmag_7(i__),Bmag_7(i__),e_Bmag_7(i__),Vmag_7(i__),
     +  e_Vmag_7(i__),Jmag_7(i__),e_Jmag_7(i__),Hmag_7(i__),
     +  e_Hmag_7(i__),Kmag_7(i__),e_Kmag_7(i__),Nobs_7(i__),P_7(i__),
     +  sub_7(i__),r_ID_7(i__),ID_7(i__)
        if(ar__7(146:147) .EQ. '') Nobs_7(i__) = iNULL__
c    ..............Just test output...........
        write(6,8)
     +  PPMX_7(i__),m_PPMX_7(i__),RAmas_7(i__),DEmas_7(i__),
     +  pmRA_7(i__),pmDE_7(i__),epRA_7(i__),epDE_7(i__),
     +  e_RAmas_7(i__),e_DEmas_7(i__),e_pmRA_7(i__),e_pmDE_7(i__),
     +  Cmag_7(i__),Rmag_7(i__),Bmag_7(i__),e_Bmag_7(i__),Vmag_7(i__),
     +  e_Vmag_7(i__),Jmag_7(i__),e_Jmag_7(i__),Hmag_7(i__),
     +  e_Hmag_7(i__),Kmag_7(i__),e_Kmag_7(i__),Nobs_7(i__),P_7(i__),
     +  sub_7(i__),r_ID_7(i__),ID_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n6000.dat'	! Declination from +60 to +67.5{deg}

C  Format for file interpretation

    9 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n6000.dat')
      write(6,*) '....Loading file: n6000.dat'
      do i__=1,512841
        read(1,'(A166)')ar__8
        read(ar__8,9)
     +  PPMX_8(i__),m_PPMX_8(i__),RAmas_8(i__),DEmas_8(i__),
     +  pmRA_8(i__),pmDE_8(i__),epRA_8(i__),epDE_8(i__),
     +  e_RAmas_8(i__),e_DEmas_8(i__),e_pmRA_8(i__),e_pmDE_8(i__),
     +  Cmag_8(i__),Rmag_8(i__),Bmag_8(i__),e_Bmag_8(i__),Vmag_8(i__),
     +  e_Vmag_8(i__),Jmag_8(i__),e_Jmag_8(i__),Hmag_8(i__),
     +  e_Hmag_8(i__),Kmag_8(i__),e_Kmag_8(i__),Nobs_8(i__),P_8(i__),
     +  sub_8(i__),r_ID_8(i__),ID_8(i__)
        if(ar__8(146:147) .EQ. '') Nobs_8(i__) = iNULL__
c    ..............Just test output...........
        write(6,9)
     +  PPMX_8(i__),m_PPMX_8(i__),RAmas_8(i__),DEmas_8(i__),
     +  pmRA_8(i__),pmDE_8(i__),epRA_8(i__),epDE_8(i__),
     +  e_RAmas_8(i__),e_DEmas_8(i__),e_pmRA_8(i__),e_pmDE_8(i__),
     +  Cmag_8(i__),Rmag_8(i__),Bmag_8(i__),e_Bmag_8(i__),Vmag_8(i__),
     +  e_Vmag_8(i__),Jmag_8(i__),e_Jmag_8(i__),Hmag_8(i__),
     +  e_Hmag_8(i__),Kmag_8(i__),e_Kmag_8(i__),Nobs_8(i__),P_8(i__),
     +  sub_8(i__),r_ID_8(i__),ID_8(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n6730.dat'	! Declination from +67.5 to +75{deg}

C  Format for file interpretation

   10 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n6730.dat')
      write(6,*) '....Loading file: n6730.dat'
      do i__=1,369595
        read(1,'(A166)')ar__9
        read(ar__9,10)
     +  PPMX_9(i__),m_PPMX_9(i__),RAmas_9(i__),DEmas_9(i__),
     +  pmRA_9(i__),pmDE_9(i__),epRA_9(i__),epDE_9(i__),
     +  e_RAmas_9(i__),e_DEmas_9(i__),e_pmRA_9(i__),e_pmDE_9(i__),
     +  Cmag_9(i__),Rmag_9(i__),Bmag_9(i__),e_Bmag_9(i__),Vmag_9(i__),
     +  e_Vmag_9(i__),Jmag_9(i__),e_Jmag_9(i__),Hmag_9(i__),
     +  e_Hmag_9(i__),Kmag_9(i__),e_Kmag_9(i__),Nobs_9(i__),P_9(i__),
     +  sub_9(i__),r_ID_9(i__),ID_9(i__)
        if(ar__9(146:147) .EQ. '') Nobs_9(i__) = iNULL__
c    ..............Just test output...........
        write(6,10)
     +  PPMX_9(i__),m_PPMX_9(i__),RAmas_9(i__),DEmas_9(i__),
     +  pmRA_9(i__),pmDE_9(i__),epRA_9(i__),epDE_9(i__),
     +  e_RAmas_9(i__),e_DEmas_9(i__),e_pmRA_9(i__),e_pmDE_9(i__),
     +  Cmag_9(i__),Rmag_9(i__),Bmag_9(i__),e_Bmag_9(i__),Vmag_9(i__),
     +  e_Vmag_9(i__),Jmag_9(i__),e_Jmag_9(i__),Hmag_9(i__),
     +  e_Hmag_9(i__),Kmag_9(i__),e_Kmag_9(i__),Nobs_9(i__),P_9(i__),
     +  sub_9(i__),r_ID_9(i__),ID_9(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n7500.dat'	! Declination from +75 to +82.5{deg}

C  Format for file interpretation

   11 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n7500.dat')
      write(6,*) '....Loading file: n7500.dat'
      do i__=1,258995
        read(1,'(A166)')ar__10
        read(ar__10,11)
     +  PPMX_10(i__),m_PPMX_10(i__),RAmas_10(i__),DEmas_10(i__),
     +  pmRA_10(i__),pmDE_10(i__),epRA_10(i__),epDE_10(i__),
     +  e_RAmas_10(i__),e_DEmas_10(i__),e_pmRA_10(i__),e_pmDE_10(i__),
     +  Cmag_10(i__),Rmag_10(i__),Bmag_10(i__),e_Bmag_10(i__),
     +  Vmag_10(i__),e_Vmag_10(i__),Jmag_10(i__),e_Jmag_10(i__),
     +  Hmag_10(i__),e_Hmag_10(i__),Kmag_10(i__),e_Kmag_10(i__),
     +  Nobs_10(i__),P_10(i__),sub_10(i__),r_ID_10(i__),ID_10(i__)
        if(ar__10(146:147) .EQ. '') Nobs_10(i__) = iNULL__
c    ..............Just test output...........
        write(6,11)
     +  PPMX_10(i__),m_PPMX_10(i__),RAmas_10(i__),DEmas_10(i__),
     +  pmRA_10(i__),pmDE_10(i__),epRA_10(i__),epDE_10(i__),
     +  e_RAmas_10(i__),e_DEmas_10(i__),e_pmRA_10(i__),e_pmDE_10(i__),
     +  Cmag_10(i__),Rmag_10(i__),Bmag_10(i__),e_Bmag_10(i__),
     +  Vmag_10(i__),e_Vmag_10(i__),Jmag_10(i__),e_Jmag_10(i__),
     +  Hmag_10(i__),e_Hmag_10(i__),Kmag_10(i__),e_Kmag_10(i__),
     +  Nobs_10(i__),P_10(i__),sub_10(i__),r_ID_10(i__),ID_10(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'n8230.dat'	! Declination from +82.5 to +90{deg}

C  Format for file interpretation

   12 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'n8230.dat')
      write(6,*) '....Loading file: n8230.dat'
      do i__=1,73460
        read(1,'(A166)')ar__11
        read(ar__11,12)
     +  PPMX_11(i__),m_PPMX_11(i__),RAmas_11(i__),DEmas_11(i__),
     +  pmRA_11(i__),pmDE_11(i__),epRA_11(i__),epDE_11(i__),
     +  e_RAmas_11(i__),e_DEmas_11(i__),e_pmRA_11(i__),e_pmDE_11(i__),
     +  Cmag_11(i__),Rmag_11(i__),Bmag_11(i__),e_Bmag_11(i__),
     +  Vmag_11(i__),e_Vmag_11(i__),Jmag_11(i__),e_Jmag_11(i__),
     +  Hmag_11(i__),e_Hmag_11(i__),Kmag_11(i__),e_Kmag_11(i__),
     +  Nobs_11(i__),P_11(i__),sub_11(i__),r_ID_11(i__),ID_11(i__)
        if(ar__11(146:147) .EQ. '') Nobs_11(i__) = iNULL__
c    ..............Just test output...........
        write(6,12)
     +  PPMX_11(i__),m_PPMX_11(i__),RAmas_11(i__),DEmas_11(i__),
     +  pmRA_11(i__),pmDE_11(i__),epRA_11(i__),epDE_11(i__),
     +  e_RAmas_11(i__),e_DEmas_11(i__),e_pmRA_11(i__),e_pmDE_11(i__),
     +  Cmag_11(i__),Rmag_11(i__),Bmag_11(i__),e_Bmag_11(i__),
     +  Vmag_11(i__),e_Vmag_11(i__),Jmag_11(i__),e_Jmag_11(i__),
     +  Hmag_11(i__),e_Hmag_11(i__),Kmag_11(i__),e_Kmag_11(i__),
     +  Nobs_11(i__),P_11(i__),sub_11(i__),r_ID_11(i__),ID_11(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's0000.dat'	! Declination from -00 to -07.5{deg}

C  Format for file interpretation

   13 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s0000.dat')
      write(6,*) '....Loading file: s0000.dat'
      do i__=1,970673
        read(1,'(A166)')ar__12
        read(ar__12,13)
     +  PPMX_12(i__),m_PPMX_12(i__),RAmas_12(i__),DEmas_12(i__),
     +  pmRA_12(i__),pmDE_12(i__),epRA_12(i__),epDE_12(i__),
     +  e_RAmas_12(i__),e_DEmas_12(i__),e_pmRA_12(i__),e_pmDE_12(i__),
     +  Cmag_12(i__),Rmag_12(i__),Bmag_12(i__),e_Bmag_12(i__),
     +  Vmag_12(i__),e_Vmag_12(i__),Jmag_12(i__),e_Jmag_12(i__),
     +  Hmag_12(i__),e_Hmag_12(i__),Kmag_12(i__),e_Kmag_12(i__),
     +  Nobs_12(i__),P_12(i__),sub_12(i__),r_ID_12(i__),ID_12(i__)
        if(ar__12(146:147) .EQ. '') Nobs_12(i__) = iNULL__
c    ..............Just test output...........
        write(6,13)
     +  PPMX_12(i__),m_PPMX_12(i__),RAmas_12(i__),DEmas_12(i__),
     +  pmRA_12(i__),pmDE_12(i__),epRA_12(i__),epDE_12(i__),
     +  e_RAmas_12(i__),e_DEmas_12(i__),e_pmRA_12(i__),e_pmDE_12(i__),
     +  Cmag_12(i__),Rmag_12(i__),Bmag_12(i__),e_Bmag_12(i__),
     +  Vmag_12(i__),e_Vmag_12(i__),Jmag_12(i__),e_Jmag_12(i__),
     +  Hmag_12(i__),e_Hmag_12(i__),Kmag_12(i__),e_Kmag_12(i__),
     +  Nobs_12(i__),P_12(i__),sub_12(i__),r_ID_12(i__),ID_12(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's0730.dat'	! Declination from -07.5 to -15{deg}

C  Format for file interpretation

   14 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s0730.dat')
      write(6,*) '....Loading file: s0730.dat'
      do i__=1,970094
        read(1,'(A166)')ar__13
        read(ar__13,14)
     +  PPMX_13(i__),m_PPMX_13(i__),RAmas_13(i__),DEmas_13(i__),
     +  pmRA_13(i__),pmDE_13(i__),epRA_13(i__),epDE_13(i__),
     +  e_RAmas_13(i__),e_DEmas_13(i__),e_pmRA_13(i__),e_pmDE_13(i__),
     +  Cmag_13(i__),Rmag_13(i__),Bmag_13(i__),e_Bmag_13(i__),
     +  Vmag_13(i__),e_Vmag_13(i__),Jmag_13(i__),e_Jmag_13(i__),
     +  Hmag_13(i__),e_Hmag_13(i__),Kmag_13(i__),e_Kmag_13(i__),
     +  Nobs_13(i__),P_13(i__),sub_13(i__),r_ID_13(i__),ID_13(i__)
        if(ar__13(146:147) .EQ. '') Nobs_13(i__) = iNULL__
c    ..............Just test output...........
        write(6,14)
     +  PPMX_13(i__),m_PPMX_13(i__),RAmas_13(i__),DEmas_13(i__),
     +  pmRA_13(i__),pmDE_13(i__),epRA_13(i__),epDE_13(i__),
     +  e_RAmas_13(i__),e_DEmas_13(i__),e_pmRA_13(i__),e_pmDE_13(i__),
     +  Cmag_13(i__),Rmag_13(i__),Bmag_13(i__),e_Bmag_13(i__),
     +  Vmag_13(i__),e_Vmag_13(i__),Jmag_13(i__),e_Jmag_13(i__),
     +  Hmag_13(i__),e_Hmag_13(i__),Kmag_13(i__),e_Kmag_13(i__),
     +  Nobs_13(i__),P_13(i__),sub_13(i__),r_ID_13(i__),ID_13(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's1500.dat'	! Declination from -15 to -22.5{deg}

C  Format for file interpretation

   15 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s1500.dat')
      write(6,*) '....Loading file: s1500.dat'
      do i__=1,1010372
        read(1,'(A166)')ar__14
        read(ar__14,15)
     +  PPMX_14(i__),m_PPMX_14(i__),RAmas_14(i__),DEmas_14(i__),
     +  pmRA_14(i__),pmDE_14(i__),epRA_14(i__),epDE_14(i__),
     +  e_RAmas_14(i__),e_DEmas_14(i__),e_pmRA_14(i__),e_pmDE_14(i__),
     +  Cmag_14(i__),Rmag_14(i__),Bmag_14(i__),e_Bmag_14(i__),
     +  Vmag_14(i__),e_Vmag_14(i__),Jmag_14(i__),e_Jmag_14(i__),
     +  Hmag_14(i__),e_Hmag_14(i__),Kmag_14(i__),e_Kmag_14(i__),
     +  Nobs_14(i__),P_14(i__),sub_14(i__),r_ID_14(i__),ID_14(i__)
        if(ar__14(146:147) .EQ. '') Nobs_14(i__) = iNULL__
c    ..............Just test output...........
        write(6,15)
     +  PPMX_14(i__),m_PPMX_14(i__),RAmas_14(i__),DEmas_14(i__),
     +  pmRA_14(i__),pmDE_14(i__),epRA_14(i__),epDE_14(i__),
     +  e_RAmas_14(i__),e_DEmas_14(i__),e_pmRA_14(i__),e_pmDE_14(i__),
     +  Cmag_14(i__),Rmag_14(i__),Bmag_14(i__),e_Bmag_14(i__),
     +  Vmag_14(i__),e_Vmag_14(i__),Jmag_14(i__),e_Jmag_14(i__),
     +  Hmag_14(i__),e_Hmag_14(i__),Kmag_14(i__),e_Kmag_14(i__),
     +  Nobs_14(i__),P_14(i__),sub_14(i__),r_ID_14(i__),ID_14(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's2230.dat'	! Declination from -22.5 to -30{deg}

C  Format for file interpretation

   16 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s2230.dat')
      write(6,*) '....Loading file: s2230.dat'
      do i__=1,1019597
        read(1,'(A166)')ar__15
        read(ar__15,16)
     +  PPMX_15(i__),m_PPMX_15(i__),RAmas_15(i__),DEmas_15(i__),
     +  pmRA_15(i__),pmDE_15(i__),epRA_15(i__),epDE_15(i__),
     +  e_RAmas_15(i__),e_DEmas_15(i__),e_pmRA_15(i__),e_pmDE_15(i__),
     +  Cmag_15(i__),Rmag_15(i__),Bmag_15(i__),e_Bmag_15(i__),
     +  Vmag_15(i__),e_Vmag_15(i__),Jmag_15(i__),e_Jmag_15(i__),
     +  Hmag_15(i__),e_Hmag_15(i__),Kmag_15(i__),e_Kmag_15(i__),
     +  Nobs_15(i__),P_15(i__),sub_15(i__),r_ID_15(i__),ID_15(i__)
        if(ar__15(146:147) .EQ. '') Nobs_15(i__) = iNULL__
c    ..............Just test output...........
        write(6,16)
     +  PPMX_15(i__),m_PPMX_15(i__),RAmas_15(i__),DEmas_15(i__),
     +  pmRA_15(i__),pmDE_15(i__),epRA_15(i__),epDE_15(i__),
     +  e_RAmas_15(i__),e_DEmas_15(i__),e_pmRA_15(i__),e_pmDE_15(i__),
     +  Cmag_15(i__),Rmag_15(i__),Bmag_15(i__),e_Bmag_15(i__),
     +  Vmag_15(i__),e_Vmag_15(i__),Jmag_15(i__),e_Jmag_15(i__),
     +  Hmag_15(i__),e_Hmag_15(i__),Kmag_15(i__),e_Kmag_15(i__),
     +  Nobs_15(i__),P_15(i__),sub_15(i__),r_ID_15(i__),ID_15(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's3000.dat'	! Declination from -30 to -37.5{deg}

C  Format for file interpretation

   17 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s3000.dat')
      write(6,*) '....Loading file: s3000.dat'
      do i__=1,975900
        read(1,'(A166)')ar__16
        read(ar__16,17)
     +  PPMX_16(i__),m_PPMX_16(i__),RAmas_16(i__),DEmas_16(i__),
     +  pmRA_16(i__),pmDE_16(i__),epRA_16(i__),epDE_16(i__),
     +  e_RAmas_16(i__),e_DEmas_16(i__),e_pmRA_16(i__),e_pmDE_16(i__),
     +  Cmag_16(i__),Rmag_16(i__),Bmag_16(i__),e_Bmag_16(i__),
     +  Vmag_16(i__),e_Vmag_16(i__),Jmag_16(i__),e_Jmag_16(i__),
     +  Hmag_16(i__),e_Hmag_16(i__),Kmag_16(i__),e_Kmag_16(i__),
     +  Nobs_16(i__),P_16(i__),sub_16(i__),r_ID_16(i__),ID_16(i__)
        if(ar__16(146:147) .EQ. '') Nobs_16(i__) = iNULL__
c    ..............Just test output...........
        write(6,17)
     +  PPMX_16(i__),m_PPMX_16(i__),RAmas_16(i__),DEmas_16(i__),
     +  pmRA_16(i__),pmDE_16(i__),epRA_16(i__),epDE_16(i__),
     +  e_RAmas_16(i__),e_DEmas_16(i__),e_pmRA_16(i__),e_pmDE_16(i__),
     +  Cmag_16(i__),Rmag_16(i__),Bmag_16(i__),e_Bmag_16(i__),
     +  Vmag_16(i__),e_Vmag_16(i__),Jmag_16(i__),e_Jmag_16(i__),
     +  Hmag_16(i__),e_Hmag_16(i__),Kmag_16(i__),e_Kmag_16(i__),
     +  Nobs_16(i__),P_16(i__),sub_16(i__),r_ID_16(i__),ID_16(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's3730.dat'	! Declination from -37.5 to -45{deg}

C  Format for file interpretation

   18 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s3730.dat')
      write(6,*) '....Loading file: s3730.dat'
      do i__=1,991495
        read(1,'(A166)')ar__17
        read(ar__17,18)
     +  PPMX_17(i__),m_PPMX_17(i__),RAmas_17(i__),DEmas_17(i__),
     +  pmRA_17(i__),pmDE_17(i__),epRA_17(i__),epDE_17(i__),
     +  e_RAmas_17(i__),e_DEmas_17(i__),e_pmRA_17(i__),e_pmDE_17(i__),
     +  Cmag_17(i__),Rmag_17(i__),Bmag_17(i__),e_Bmag_17(i__),
     +  Vmag_17(i__),e_Vmag_17(i__),Jmag_17(i__),e_Jmag_17(i__),
     +  Hmag_17(i__),e_Hmag_17(i__),Kmag_17(i__),e_Kmag_17(i__),
     +  Nobs_17(i__),P_17(i__),sub_17(i__),r_ID_17(i__),ID_17(i__)
        if(ar__17(146:147) .EQ. '') Nobs_17(i__) = iNULL__
c    ..............Just test output...........
        write(6,18)
     +  PPMX_17(i__),m_PPMX_17(i__),RAmas_17(i__),DEmas_17(i__),
     +  pmRA_17(i__),pmDE_17(i__),epRA_17(i__),epDE_17(i__),
     +  e_RAmas_17(i__),e_DEmas_17(i__),e_pmRA_17(i__),e_pmDE_17(i__),
     +  Cmag_17(i__),Rmag_17(i__),Bmag_17(i__),e_Bmag_17(i__),
     +  Vmag_17(i__),e_Vmag_17(i__),Jmag_17(i__),e_Jmag_17(i__),
     +  Hmag_17(i__),e_Hmag_17(i__),Kmag_17(i__),e_Kmag_17(i__),
     +  Nobs_17(i__),P_17(i__),sub_17(i__),r_ID_17(i__),ID_17(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's4500.dat'	! Declination from -45 to -52.5{deg}

C  Format for file interpretation

   19 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s4500.dat')
      write(6,*) '....Loading file: s4500.dat'
      do i__=1,951351
        read(1,'(A166)')ar__18
        read(ar__18,19)
     +  PPMX_18(i__),m_PPMX_18(i__),RAmas_18(i__),DEmas_18(i__),
     +  pmRA_18(i__),pmDE_18(i__),epRA_18(i__),epDE_18(i__),
     +  e_RAmas_18(i__),e_DEmas_18(i__),e_pmRA_18(i__),e_pmDE_18(i__),
     +  Cmag_18(i__),Rmag_18(i__),Bmag_18(i__),e_Bmag_18(i__),
     +  Vmag_18(i__),e_Vmag_18(i__),Jmag_18(i__),e_Jmag_18(i__),
     +  Hmag_18(i__),e_Hmag_18(i__),Kmag_18(i__),e_Kmag_18(i__),
     +  Nobs_18(i__),P_18(i__),sub_18(i__),r_ID_18(i__),ID_18(i__)
        if(ar__18(146:147) .EQ. '') Nobs_18(i__) = iNULL__
c    ..............Just test output...........
        write(6,19)
     +  PPMX_18(i__),m_PPMX_18(i__),RAmas_18(i__),DEmas_18(i__),
     +  pmRA_18(i__),pmDE_18(i__),epRA_18(i__),epDE_18(i__),
     +  e_RAmas_18(i__),e_DEmas_18(i__),e_pmRA_18(i__),e_pmDE_18(i__),
     +  Cmag_18(i__),Rmag_18(i__),Bmag_18(i__),e_Bmag_18(i__),
     +  Vmag_18(i__),e_Vmag_18(i__),Jmag_18(i__),e_Jmag_18(i__),
     +  Hmag_18(i__),e_Hmag_18(i__),Kmag_18(i__),e_Kmag_18(i__),
     +  Nobs_18(i__),P_18(i__),sub_18(i__),r_ID_18(i__),ID_18(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's5230.dat'	! Declination from -52.5 to -60{deg}

C  Format for file interpretation

   20 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s5230.dat')
      write(6,*) '....Loading file: s5230.dat'
      do i__=1,782446
        read(1,'(A166)')ar__19
        read(ar__19,20)
     +  PPMX_19(i__),m_PPMX_19(i__),RAmas_19(i__),DEmas_19(i__),
     +  pmRA_19(i__),pmDE_19(i__),epRA_19(i__),epDE_19(i__),
     +  e_RAmas_19(i__),e_DEmas_19(i__),e_pmRA_19(i__),e_pmDE_19(i__),
     +  Cmag_19(i__),Rmag_19(i__),Bmag_19(i__),e_Bmag_19(i__),
     +  Vmag_19(i__),e_Vmag_19(i__),Jmag_19(i__),e_Jmag_19(i__),
     +  Hmag_19(i__),e_Hmag_19(i__),Kmag_19(i__),e_Kmag_19(i__),
     +  Nobs_19(i__),P_19(i__),sub_19(i__),r_ID_19(i__),ID_19(i__)
        if(ar__19(146:147) .EQ. '') Nobs_19(i__) = iNULL__
c    ..............Just test output...........
        write(6,20)
     +  PPMX_19(i__),m_PPMX_19(i__),RAmas_19(i__),DEmas_19(i__),
     +  pmRA_19(i__),pmDE_19(i__),epRA_19(i__),epDE_19(i__),
     +  e_RAmas_19(i__),e_DEmas_19(i__),e_pmRA_19(i__),e_pmDE_19(i__),
     +  Cmag_19(i__),Rmag_19(i__),Bmag_19(i__),e_Bmag_19(i__),
     +  Vmag_19(i__),e_Vmag_19(i__),Jmag_19(i__),e_Jmag_19(i__),
     +  Hmag_19(i__),e_Hmag_19(i__),Kmag_19(i__),e_Kmag_19(i__),
     +  Nobs_19(i__),P_19(i__),sub_19(i__),r_ID_19(i__),ID_19(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's6000.dat'	! Declination from -60 to -67.5{deg}

C  Format for file interpretation

   21 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s6000.dat')
      write(6,*) '....Loading file: s6000.dat'
      do i__=1,640425
        read(1,'(A166)')ar__20
        read(ar__20,21)
     +  PPMX_20(i__),m_PPMX_20(i__),RAmas_20(i__),DEmas_20(i__),
     +  pmRA_20(i__),pmDE_20(i__),epRA_20(i__),epDE_20(i__),
     +  e_RAmas_20(i__),e_DEmas_20(i__),e_pmRA_20(i__),e_pmDE_20(i__),
     +  Cmag_20(i__),Rmag_20(i__),Bmag_20(i__),e_Bmag_20(i__),
     +  Vmag_20(i__),e_Vmag_20(i__),Jmag_20(i__),e_Jmag_20(i__),
     +  Hmag_20(i__),e_Hmag_20(i__),Kmag_20(i__),e_Kmag_20(i__),
     +  Nobs_20(i__),P_20(i__),sub_20(i__),r_ID_20(i__),ID_20(i__)
        if(ar__20(146:147) .EQ. '') Nobs_20(i__) = iNULL__
c    ..............Just test output...........
        write(6,21)
     +  PPMX_20(i__),m_PPMX_20(i__),RAmas_20(i__),DEmas_20(i__),
     +  pmRA_20(i__),pmDE_20(i__),epRA_20(i__),epDE_20(i__),
     +  e_RAmas_20(i__),e_DEmas_20(i__),e_pmRA_20(i__),e_pmDE_20(i__),
     +  Cmag_20(i__),Rmag_20(i__),Bmag_20(i__),e_Bmag_20(i__),
     +  Vmag_20(i__),e_Vmag_20(i__),Jmag_20(i__),e_Jmag_20(i__),
     +  Hmag_20(i__),e_Hmag_20(i__),Kmag_20(i__),e_Kmag_20(i__),
     +  Nobs_20(i__),P_20(i__),sub_20(i__),r_ID_20(i__),ID_20(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's6730.dat'	! Declination from -67.5 to -75{deg}

C  Format for file interpretation

   22 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s6730.dat')
      write(6,*) '....Loading file: s6730.dat'
      do i__=1,452628
        read(1,'(A166)')ar__21
        read(ar__21,22)
     +  PPMX_21(i__),m_PPMX_21(i__),RAmas_21(i__),DEmas_21(i__),
     +  pmRA_21(i__),pmDE_21(i__),epRA_21(i__),epDE_21(i__),
     +  e_RAmas_21(i__),e_DEmas_21(i__),e_pmRA_21(i__),e_pmDE_21(i__),
     +  Cmag_21(i__),Rmag_21(i__),Bmag_21(i__),e_Bmag_21(i__),
     +  Vmag_21(i__),e_Vmag_21(i__),Jmag_21(i__),e_Jmag_21(i__),
     +  Hmag_21(i__),e_Hmag_21(i__),Kmag_21(i__),e_Kmag_21(i__),
     +  Nobs_21(i__),P_21(i__),sub_21(i__),r_ID_21(i__),ID_21(i__)
        if(ar__21(146:147) .EQ. '') Nobs_21(i__) = iNULL__
c    ..............Just test output...........
        write(6,22)
     +  PPMX_21(i__),m_PPMX_21(i__),RAmas_21(i__),DEmas_21(i__),
     +  pmRA_21(i__),pmDE_21(i__),epRA_21(i__),epDE_21(i__),
     +  e_RAmas_21(i__),e_DEmas_21(i__),e_pmRA_21(i__),e_pmDE_21(i__),
     +  Cmag_21(i__),Rmag_21(i__),Bmag_21(i__),e_Bmag_21(i__),
     +  Vmag_21(i__),e_Vmag_21(i__),Jmag_21(i__),e_Jmag_21(i__),
     +  Hmag_21(i__),e_Hmag_21(i__),Kmag_21(i__),e_Kmag_21(i__),
     +  Nobs_21(i__),P_21(i__),sub_21(i__),r_ID_21(i__),ID_21(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's7500.dat'	! Declination from -75 to -82.5{deg}

C  Format for file interpretation

   23 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s7500.dat')
      write(6,*) '....Loading file: s7500.dat'
      do i__=1,269692
        read(1,'(A166)')ar__22
        read(ar__22,23)
     +  PPMX_22(i__),m_PPMX_22(i__),RAmas_22(i__),DEmas_22(i__),
     +  pmRA_22(i__),pmDE_22(i__),epRA_22(i__),epDE_22(i__),
     +  e_RAmas_22(i__),e_DEmas_22(i__),e_pmRA_22(i__),e_pmDE_22(i__),
     +  Cmag_22(i__),Rmag_22(i__),Bmag_22(i__),e_Bmag_22(i__),
     +  Vmag_22(i__),e_Vmag_22(i__),Jmag_22(i__),e_Jmag_22(i__),
     +  Hmag_22(i__),e_Hmag_22(i__),Kmag_22(i__),e_Kmag_22(i__),
     +  Nobs_22(i__),P_22(i__),sub_22(i__),r_ID_22(i__),ID_22(i__)
        if(ar__22(146:147) .EQ. '') Nobs_22(i__) = iNULL__
c    ..............Just test output...........
        write(6,23)
     +  PPMX_22(i__),m_PPMX_22(i__),RAmas_22(i__),DEmas_22(i__),
     +  pmRA_22(i__),pmDE_22(i__),epRA_22(i__),epDE_22(i__),
     +  e_RAmas_22(i__),e_DEmas_22(i__),e_pmRA_22(i__),e_pmDE_22(i__),
     +  Cmag_22(i__),Rmag_22(i__),Bmag_22(i__),e_Bmag_22(i__),
     +  Vmag_22(i__),e_Vmag_22(i__),Jmag_22(i__),e_Jmag_22(i__),
     +  Hmag_22(i__),e_Hmag_22(i__),Kmag_22(i__),e_Kmag_22(i__),
     +  Nobs_22(i__),P_22(i__),sub_22(i__),r_ID_22(i__),ID_22(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 's8230.dat'	! Declination from -82.5 to -90{deg}

C  Format for file interpretation

   24 format(
     +  A15,A1,1X,I10,1X,I10,1X,I7,1X,I7,1X,I5,1X,I5,1X,I3,1X,I3,1X,
     +  I3,1X,I3,1X,I5,1X,I5,1X,I5,1X,I3,1X,I5,1X,I3,1X,I5,1X,I3,1X,
     +  I5,1X,I3,1X,I5,1X,I3,1X,I2,1X,A1,1X,A1,1X,A1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'s8230.dat')
      write(6,*) '....Loading file: s8230.dat'
      do i__=1,92360
        read(1,'(A166)')ar__23
        read(ar__23,24)
     +  PPMX_23(i__),m_PPMX_23(i__),RAmas_23(i__),DEmas_23(i__),
     +  pmRA_23(i__),pmDE_23(i__),epRA_23(i__),epDE_23(i__),
     +  e_RAmas_23(i__),e_DEmas_23(i__),e_pmRA_23(i__),e_pmDE_23(i__),
     +  Cmag_23(i__),Rmag_23(i__),Bmag_23(i__),e_Bmag_23(i__),
     +  Vmag_23(i__),e_Vmag_23(i__),Jmag_23(i__),e_Jmag_23(i__),
     +  Hmag_23(i__),e_Hmag_23(i__),Kmag_23(i__),e_Kmag_23(i__),
     +  Nobs_23(i__),P_23(i__),sub_23(i__),r_ID_23(i__),ID_23(i__)
        if(ar__23(146:147) .EQ. '') Nobs_23(i__) = iNULL__
c    ..............Just test output...........
        write(6,24)
     +  PPMX_23(i__),m_PPMX_23(i__),RAmas_23(i__),DEmas_23(i__),
     +  pmRA_23(i__),pmDE_23(i__),epRA_23(i__),epDE_23(i__),
     +  e_RAmas_23(i__),e_DEmas_23(i__),e_pmRA_23(i__),e_pmDE_23(i__),
     +  Cmag_23(i__),Rmag_23(i__),Bmag_23(i__),e_Bmag_23(i__),
     +  Vmag_23(i__),e_Vmag_23(i__),Jmag_23(i__),e_Jmag_23(i__),
     +  Hmag_23(i__),e_Hmag_23(i__),Kmag_23(i__),e_Kmag_23(i__),
     +  Nobs_23(i__),P_23(i__),sub_23(i__),r_ID_23(i__),ID_23(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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