FORTRAN Generation
(/./ftp/cats/II/228A)

Conversion of standardized ReadMe file for file /./ftp/cats/II/228A into FORTRAN code for loading all data files into arrays.

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



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

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

C=============================================================================
Cat. II/228A    DENIS Catalogue toward Magellanic Clouds (DCMC)     (Cioni+ 2000)
*================================================================================
*The DENIS Point Source Catalogue towards the Magellanic Clouds
*    Cioni M.-R., Loup C., Habing H.J., Fouque P., Bertin E., Deul E.,
*    Egret D., Alard C., de  Batz B., Borsenberger J., Dennefeld M.,
*    Epchtein N., Forveille T., Garzon F., Hron J., Kimeswenger S.,
*    Lacombe F., Le  Bertre T., Mamon G.A., Omont A., Paturel G., Persi P.,
*    Robin A., Rouan D., Simon G., Tiphene D., Vauglin I., Wagner S.
*   <Astron. Astrophys. Suppl. Ser., 144, 235 (2000)>
*   =2000A&AS..144..235C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'lmc_ijk.dat'	! LMC stars observed in the 3 bands I, J and K (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC       (nr__) ! Name of the source (based on J2000 position)
      integer*4     RAh        (nr__) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm        (nr__) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs        (nr__) ! (s) Right Ascension J2000 (seconds)
      character*1   DE_        (nr__) ! Declination J2000 (sign)
      integer*4     DEd        (nr__) ! (deg) Declination J2000 (degrees)
      integer*4     DEm        (nr__) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs        (nr__) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr     (nr__) ! (arcsec) Systematic positional error
      real*8        Ixpos      (nr__) ! (pix) Object position along x-axis in I image
      real*8        Iypos      (nr__) ! (pix) Object position along y-axis in I image
      integer*4     Ifield     (nr__) ! ?=999 Field number of I image
      integer*4     Istrip     (nr__) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag       (nr__) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag     (nr__) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag     (nr__) ! [0,7] Extraction flags (2)
      integer*4     n_Imag     (nr__) ! [0,63] Image flags in I band (3)
      real*8        Jxpos      (nr__) ! (pix) Object position along x-axis in J image
      real*8        Jypos      (nr__) ! (pix) Object position along y-axis in J image
      integer*4     Jfield     (nr__) ! ?=999 Field number of J image
      integer*4     Jstrip     (nr__) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag       (nr__) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag     (nr__) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag     (nr__) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag     (nr__) ! [0,127] Image flags in J band (3)
      real*8        Kxpos      (nr__) ! (pix) Object position along x-axis in K image
      real*8        Kypos      (nr__) ! (pix) Object position along y-axis in K image
      integer*4     Kfield     (nr__) ! ?=999 Field number of K image
      integer*4     Kstrip     (nr__) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag       (nr__) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag     (nr__) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag     (nr__) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag     (nr__) ! [0,127] Image flags in J band (3)
      real*4        Bmag       (nr__) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag       (nr__) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'lmc_ij.dat'	! LMC stars observed in I and J bands (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1    (nr__1) ! (deg) Right Ascension J2000
      real*8        DEdeg_1    (nr__1) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC_1     (nr__1) ! Name of the source (based on J2000 position)
      integer*4     RAh_1      (nr__1) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_1      (nr__1) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_1      (nr__1) ! (s) Right Ascension J2000 (seconds)
      character*1   DE__1      (nr__1) ! Declination J2000 (sign)
      integer*4     DEd_1      (nr__1) ! (deg) Declination J2000 (degrees)
      integer*4     DEm_1      (nr__1) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs_1      (nr__1) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr_1   (nr__1) ! (arcsec) Systematic positional error
      real*8        Ixpos_1    (nr__1) ! (pix) Object position along x-axis in I image
      real*8        Iypos_1    (nr__1) ! (pix) Object position along y-axis in I image
      integer*4     Ifield_1   (nr__1) ! ?=999 Field number of I image
      integer*4     Istrip_1   (nr__1) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag_1     (nr__1) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag_1   (nr__1) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag_1   (nr__1) ! [0,7] Extraction flags (2)
      integer*4     n_Imag_1   (nr__1) ! [0,63] Image flags in I band (3)
      real*8        Jxpos_1    (nr__1) ! (pix) Object position along x-axis in J image
      real*8        Jypos_1    (nr__1) ! (pix) Object position along y-axis in J image
      integer*4     Jfield_1   (nr__1) ! ?=999 Field number of J image
      integer*4     Jstrip_1   (nr__1) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag_1     (nr__1) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag_1   (nr__1) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag_1   (nr__1) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag_1   (nr__1) ! [0,127] Image flags in J band (3)
      real*8        Kxpos_1    (nr__1) ! (pix) Object position along x-axis in K image
      real*8        Kypos_1    (nr__1) ! (pix) Object position along y-axis in K image
      integer*4     Kfield_1   (nr__1) ! ?=999 Field number of K image
      integer*4     Kstrip_1   (nr__1) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag_1     (nr__1) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag_1   (nr__1) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag_1   (nr__1) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag_1   (nr__1) ! [0,127] Image flags in J band (3)
      real*4        Bmag_1     (nr__1) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag_1     (nr__1) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'lmc_ik.dat'	! LMC stars observed in I and K bands (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_2    (nr__2) ! (deg) Right Ascension J2000
      real*8        DEdeg_2    (nr__2) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC_2     (nr__2) ! Name of the source (based on J2000 position)
      integer*4     RAh_2      (nr__2) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_2      (nr__2) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_2      (nr__2) ! (s) Right Ascension J2000 (seconds)
      character*1   DE__2      (nr__2) ! Declination J2000 (sign)
      integer*4     DEd_2      (nr__2) ! (deg) Declination J2000 (degrees)
      integer*4     DEm_2      (nr__2) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs_2      (nr__2) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr_2   (nr__2) ! (arcsec) Systematic positional error
      real*8        Ixpos_2    (nr__2) ! (pix) Object position along x-axis in I image
      real*8        Iypos_2    (nr__2) ! (pix) Object position along y-axis in I image
      integer*4     Ifield_2   (nr__2) ! ?=999 Field number of I image
      integer*4     Istrip_2   (nr__2) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag_2     (nr__2) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag_2   (nr__2) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag_2   (nr__2) ! [0,7] Extraction flags (2)
      integer*4     n_Imag_2   (nr__2) ! [0,63] Image flags in I band (3)
      real*8        Jxpos_2    (nr__2) ! (pix) Object position along x-axis in J image
      real*8        Jypos_2    (nr__2) ! (pix) Object position along y-axis in J image
      integer*4     Jfield_2   (nr__2) ! ?=999 Field number of J image
      integer*4     Jstrip_2   (nr__2) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag_2     (nr__2) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag_2   (nr__2) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag_2   (nr__2) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag_2   (nr__2) ! [0,127] Image flags in J band (3)
      real*8        Kxpos_2    (nr__2) ! (pix) Object position along x-axis in K image
      real*8        Kypos_2    (nr__2) ! (pix) Object position along y-axis in K image
      integer*4     Kfield_2   (nr__2) ! ?=999 Field number of K image
      integer*4     Kstrip_2   (nr__2) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag_2     (nr__2) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag_2   (nr__2) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag_2   (nr__2) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag_2   (nr__2) ! [0,127] Image flags in J band (3)
      real*4        Bmag_2     (nr__2) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag_2     (nr__2) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'lmc_jk.dat'	! LMC stars observed in J and K bands (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_3    (nr__3) ! (deg) Right Ascension J2000
      real*8        DEdeg_3    (nr__3) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC_3     (nr__3) ! Name of the source (based on J2000 position)
      integer*4     RAh_3      (nr__3) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_3      (nr__3) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_3      (nr__3) ! (s) Right Ascension J2000 (seconds)
      character*1   DE__3      (nr__3) ! Declination J2000 (sign)
      integer*4     DEd_3      (nr__3) ! (deg) Declination J2000 (degrees)
      integer*4     DEm_3      (nr__3) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs_3      (nr__3) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr_3   (nr__3) ! (arcsec) Systematic positional error
      real*8        Ixpos_3    (nr__3) ! (pix) Object position along x-axis in I image
      real*8        Iypos_3    (nr__3) ! (pix) Object position along y-axis in I image
      integer*4     Ifield_3   (nr__3) ! ?=999 Field number of I image
      integer*4     Istrip_3   (nr__3) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag_3     (nr__3) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag_3   (nr__3) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag_3   (nr__3) ! [0,7] Extraction flags (2)
      integer*4     n_Imag_3   (nr__3) ! [0,63] Image flags in I band (3)
      real*8        Jxpos_3    (nr__3) ! (pix) Object position along x-axis in J image
      real*8        Jypos_3    (nr__3) ! (pix) Object position along y-axis in J image
      integer*4     Jfield_3   (nr__3) ! ?=999 Field number of J image
      integer*4     Jstrip_3   (nr__3) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag_3     (nr__3) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag_3   (nr__3) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag_3   (nr__3) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag_3   (nr__3) ! [0,127] Image flags in J band (3)
      real*8        Kxpos_3    (nr__3) ! (pix) Object position along x-axis in K image
      real*8        Kypos_3    (nr__3) ! (pix) Object position along y-axis in K image
      integer*4     Kfield_3   (nr__3) ! ?=999 Field number of K image
      integer*4     Kstrip_3   (nr__3) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag_3     (nr__3) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag_3   (nr__3) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag_3   (nr__3) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag_3   (nr__3) ! [0,127] Image flags in J band (3)
      real*4        Bmag_3     (nr__3) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag_3     (nr__3) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'smc_ijk.dat'	! SMC stars observed in the 3 bands I, J and K (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_4    (nr__4) ! (deg) Right Ascension J2000
      real*8        DEdeg_4    (nr__4) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC_4     (nr__4) ! Name of the source (based on J2000 position)
      integer*4     RAh_4      (nr__4) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_4      (nr__4) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_4      (nr__4) ! (s) Right Ascension J2000 (seconds)
      character*1   DE__4      (nr__4) ! Declination J2000 (sign)
      integer*4     DEd_4      (nr__4) ! (deg) Declination J2000 (degrees)
      integer*4     DEm_4      (nr__4) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs_4      (nr__4) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr_4   (nr__4) ! (arcsec) Systematic positional error
      real*8        Ixpos_4    (nr__4) ! (pix) Object position along x-axis in I image
      real*8        Iypos_4    (nr__4) ! (pix) Object position along y-axis in I image
      integer*4     Ifield_4   (nr__4) ! ?=999 Field number of I image
      integer*4     Istrip_4   (nr__4) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag_4     (nr__4) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag_4   (nr__4) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag_4   (nr__4) ! [0,7] Extraction flags (2)
      integer*4     n_Imag_4   (nr__4) ! [0,63] Image flags in I band (3)
      real*8        Jxpos_4    (nr__4) ! (pix) Object position along x-axis in J image
      real*8        Jypos_4    (nr__4) ! (pix) Object position along y-axis in J image
      integer*4     Jfield_4   (nr__4) ! ?=999 Field number of J image
      integer*4     Jstrip_4   (nr__4) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag_4     (nr__4) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag_4   (nr__4) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag_4   (nr__4) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag_4   (nr__4) ! [0,127] Image flags in J band (3)
      real*8        Kxpos_4    (nr__4) ! (pix) Object position along x-axis in K image
      real*8        Kypos_4    (nr__4) ! (pix) Object position along y-axis in K image
      integer*4     Kfield_4   (nr__4) ! ?=999 Field number of K image
      integer*4     Kstrip_4   (nr__4) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag_4     (nr__4) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag_4   (nr__4) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag_4   (nr__4) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag_4   (nr__4) ! [0,127] Image flags in J band (3)
      real*4        Bmag_4     (nr__4) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag_4     (nr__4) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'smc_ij.dat'	! SMC stars observed in I and J bands (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_5    (nr__5) ! (deg) Right Ascension J2000
      real*8        DEdeg_5    (nr__5) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC_5     (nr__5) ! Name of the source (based on J2000 position)
      integer*4     RAh_5      (nr__5) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_5      (nr__5) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_5      (nr__5) ! (s) Right Ascension J2000 (seconds)
      character*1   DE__5      (nr__5) ! Declination J2000 (sign)
      integer*4     DEd_5      (nr__5) ! (deg) Declination J2000 (degrees)
      integer*4     DEm_5      (nr__5) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs_5      (nr__5) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr_5   (nr__5) ! (arcsec) Systematic positional error
      real*8        Ixpos_5    (nr__5) ! (pix) Object position along x-axis in I image
      real*8        Iypos_5    (nr__5) ! (pix) Object position along y-axis in I image
      integer*4     Ifield_5   (nr__5) ! ?=999 Field number of I image
      integer*4     Istrip_5   (nr__5) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag_5     (nr__5) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag_5   (nr__5) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag_5   (nr__5) ! [0,7] Extraction flags (2)
      integer*4     n_Imag_5   (nr__5) ! [0,63] Image flags in I band (3)
      real*8        Jxpos_5    (nr__5) ! (pix) Object position along x-axis in J image
      real*8        Jypos_5    (nr__5) ! (pix) Object position along y-axis in J image
      integer*4     Jfield_5   (nr__5) ! ?=999 Field number of J image
      integer*4     Jstrip_5   (nr__5) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag_5     (nr__5) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag_5   (nr__5) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag_5   (nr__5) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag_5   (nr__5) ! [0,127] Image flags in J band (3)
      real*8        Kxpos_5    (nr__5) ! (pix) Object position along x-axis in K image
      real*8        Kypos_5    (nr__5) ! (pix) Object position along y-axis in K image
      integer*4     Kfield_5   (nr__5) ! ?=999 Field number of K image
      integer*4     Kstrip_5   (nr__5) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag_5     (nr__5) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag_5   (nr__5) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag_5   (nr__5) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag_5   (nr__5) ! [0,127] Image flags in J band (3)
      real*4        Bmag_5     (nr__5) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag_5     (nr__5) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'smc_ik.dat'	! SMC stars observed in I and K bands (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_6    (nr__6) ! (deg) Right Ascension J2000
      real*8        DEdeg_6    (nr__6) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC_6     (nr__6) ! Name of the source (based on J2000 position)
      integer*4     RAh_6      (nr__6) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_6      (nr__6) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_6      (nr__6) ! (s) Right Ascension J2000 (seconds)
      character*1   DE__6      (nr__6) ! Declination J2000 (sign)
      integer*4     DEd_6      (nr__6) ! (deg) Declination J2000 (degrees)
      integer*4     DEm_6      (nr__6) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs_6      (nr__6) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr_6   (nr__6) ! (arcsec) Systematic positional error
      real*8        Ixpos_6    (nr__6) ! (pix) Object position along x-axis in I image
      real*8        Iypos_6    (nr__6) ! (pix) Object position along y-axis in I image
      integer*4     Ifield_6   (nr__6) ! ?=999 Field number of I image
      integer*4     Istrip_6   (nr__6) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag_6     (nr__6) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag_6   (nr__6) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag_6   (nr__6) ! [0,7] Extraction flags (2)
      integer*4     n_Imag_6   (nr__6) ! [0,63] Image flags in I band (3)
      real*8        Jxpos_6    (nr__6) ! (pix) Object position along x-axis in J image
      real*8        Jypos_6    (nr__6) ! (pix) Object position along y-axis in J image
      integer*4     Jfield_6   (nr__6) ! ?=999 Field number of J image
      integer*4     Jstrip_6   (nr__6) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag_6     (nr__6) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag_6   (nr__6) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag_6   (nr__6) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag_6   (nr__6) ! [0,127] Image flags in J band (3)
      real*8        Kxpos_6    (nr__6) ! (pix) Object position along x-axis in K image
      real*8        Kypos_6    (nr__6) ! (pix) Object position along y-axis in K image
      integer*4     Kfield_6   (nr__6) ! ?=999 Field number of K image
      integer*4     Kstrip_6   (nr__6) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag_6     (nr__6) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag_6   (nr__6) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag_6   (nr__6) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag_6   (nr__6) ! [0,127] Image flags in J band (3)
      real*4        Bmag_6     (nr__6) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag_6     (nr__6) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'smc_jk.dat'	! SMC stars observed in J and K bands (1)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_7    (nr__7) ! (deg) Right Ascension J2000
      real*8        DEdeg_7    (nr__7) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*19  DCMC_7     (nr__7) ! Name of the source (based on J2000 position)
      integer*4     RAh_7      (nr__7) ! (h) Right Ascension J2000 (hours)
      integer*4     RAm_7      (nr__7) ! (min) Right Ascension J2000 (minutes)
      real*4        RAs_7      (nr__7) ! (s) Right Ascension J2000 (seconds)
      character*1   DE__7      (nr__7) ! Declination J2000 (sign)
      integer*4     DEd_7      (nr__7) ! (deg) Declination J2000 (degrees)
      integer*4     DEm_7      (nr__7) ! (arcmin) Declination J2000 (minutes)
      real*4        DEs_7      (nr__7) ! (arcsec) Declination J2000 (seconds)
      real*4        PosErr_7   (nr__7) ! (arcsec) Systematic positional error
      real*8        Ixpos_7    (nr__7) ! (pix) Object position along x-axis in I image
      real*8        Iypos_7    (nr__7) ! (pix) Object position along y-axis in I image
      integer*4     Ifield_7   (nr__7) ! ?=999 Field number of I image
      integer*4     Istrip_7   (nr__7) ! ?=0 Strip number of I image (in field Ifield)
      real*4        Imag_7     (nr__7) ! (mag) ?=99.0 Magnitude derived within 7" diameter
      real*4        e_Imag_7   (nr__7) ! (mag) ?=99.0 Mean error in Imag (1)
      integer*4     f_Imag_7   (nr__7) ! [0,7] Extraction flags (2)
      integer*4     n_Imag_7   (nr__7) ! [0,63] Image flags in I band (3)
      real*8        Jxpos_7    (nr__7) ! (pix) Object position along x-axis in J image
      real*8        Jypos_7    (nr__7) ! (pix) Object position along y-axis in J image
      integer*4     Jfield_7   (nr__7) ! ?=999 Field number of J image
      integer*4     Jstrip_7   (nr__7) ! ?=0 Strip number of J image (in field Jfield)
      real*4        Jmag_7     (nr__7) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Jmag_7   (nr__7) ! (mag) ?=99.00 Mean error in Jmag (1)
      integer*4     f_Jmag_7   (nr__7) ! [0,7] Extraction flags (2)
      integer*4     n_Jmag_7   (nr__7) ! [0,127] Image flags in J band (3)
      real*8        Kxpos_7    (nr__7) ! (pix) Object position along x-axis in K image
      real*8        Kypos_7    (nr__7) ! (pix) Object position along y-axis in K image
      integer*4     Kfield_7   (nr__7) ! ?=999 Field number of K image
      integer*4     Kstrip_7   (nr__7) ! ?=0 Strip number of K image (in field Kfield)
      real*4        Kmag_7     (nr__7) ! (mag) ?=99.00 Magnitude derived within 7" diameter
      real*4        e_Kmag_7   (nr__7) ! (mag) ?=99.00 Mean error in Kmag (1)
      integer*4     f_Kmag_7   (nr__7) ! [0,7] Extraction flags (2)
      integer*4     n_Kmag_7   (nr__7) ! [0,127] Image flags in J band (3)
      real*4        Bmag_7     (nr__7) ! (mag) ?=99. Blue magnitude from USNO-A2.0 <I/252>
      real*4        Rmag_7     (nr__7) ! (mag) ?=99. Red  magnitude from USNO-A2.0 <I/252>
*Note (1): mean errors are derived from the photon statistics
*    within the aperture
*Note (2): the extractions flags are the summation of the following
*    exclusive numbers:
*    1 = Bright neighbours or bad pixels
*        (more than 10  of the integrated area affected)
*    2 = The object was originally blended with another
*    4 = At least one pixel of the object is saturated
*    Some of the very bright sources with the highest extraction flag values
*    may suffer from large (>10arcsec) astrometric shift.
*Note (3): the image flags are:
*     1 = Regression failed for this pixel
*     2 = Flat derived is less than 0.01 (dead pixels)
*     3 = Abnormally high flat
*     4 = Abnormally high residual after adjustment (5{sigma})
*     5 = Bias more than 5{sigma} of mean
*     6 = Flat more than 5{sigma} of mean
*     9 = Insufficient number of pixels left after iterative clipping
*    16 = Flat < 0.7 and/or flux < -200 ADU
*    32 = Flat > 1.3 and/or flux > 49500 ADU

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

C  Declarations for 'qual.dat'	! Qualities of observation strips (table 8)

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

      integer*4     Strip      (nr__8) ! Strip number
      character*8   Date       (nr__8) ! ("DD/MM/YY") Date of observation
      real*4        Izp        (nr__8) ! (mag) Zero point in I band
      real*4        e_Izp      (nr__8) ! (mag) Standard error on Izp
      real*4        Jzp        (nr__8) ! (mag) Zero point in J band
      real*4        e_Jzp      (nr__8) ! (mag) Standard error on Jzp
      real*4        Kzp        (nr__8) ! (mag) Zero point in K band
      real*4        e_Kzp      (nr__8) ! (mag) Standard error on Kzp
      character*15  Notes      (nr__8) ! Value of the photometric shift, and/or
*                                   default zero point assumed when "ZP_Idef"

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

C  Loading file 'lmc_ijk.dat'	! LMC stars observed in the 3 bands I, J and K (1)

C  Format for file interpretation

    1 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'lmc_ijk.dat')
      write(6,*) '....Loading file: lmc_ijk.dat'
      do i__=1,264347
        read(1,'(A210)')ar__
        read(ar__,1)
     +  DCMC(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),PosErr(i__),Ixpos(i__),Iypos(i__),
     +  Ifield(i__),Istrip(i__),Imag(i__),e_Imag(i__),f_Imag(i__),
     +  n_Imag(i__),Jxpos(i__),Jypos(i__),Jfield(i__),Jstrip(i__),
     +  Jmag(i__),e_Jmag(i__),f_Jmag(i__),n_Jmag(i__),Kxpos(i__),
     +  Kypos(i__),Kfield(i__),Kstrip(i__),Kmag(i__),e_Kmag(i__),
     +  f_Kmag(i__),n_Kmag(i__),Bmag(i__),Rmag(i__)
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  DCMC(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),PosErr(i__),Ixpos(i__),Iypos(i__),
     +  Ifield(i__),Istrip(i__),Imag(i__),e_Imag(i__),f_Imag(i__),
     +  n_Imag(i__),Jxpos(i__),Jypos(i__),Jfield(i__),Jstrip(i__),
     +  Jmag(i__),e_Jmag(i__),f_Jmag(i__),n_Jmag(i__),Kxpos(i__),
     +  Kypos(i__),Kfield(i__),Kstrip(i__),Kmag(i__),e_Kmag(i__),
     +  f_Kmag(i__),n_Kmag(i__),Bmag(i__),Rmag(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'lmc_ij.dat'	! LMC stars observed in I and J bands (1)

C  Format for file interpretation

    2 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'lmc_ij.dat')
      write(6,*) '....Loading file: lmc_ij.dat'
      do i__=1,1045630
        read(1,'(A210)')ar__1
        read(ar__1,2)
     +  DCMC_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
     +  DEd_1(i__),DEm_1(i__),DEs_1(i__),PosErr_1(i__),Ixpos_1(i__),
     +  Iypos_1(i__),Ifield_1(i__),Istrip_1(i__),Imag_1(i__),
     +  e_Imag_1(i__),f_Imag_1(i__),n_Imag_1(i__),Jxpos_1(i__),
     +  Jypos_1(i__),Jfield_1(i__),Jstrip_1(i__),Jmag_1(i__),
     +  e_Jmag_1(i__),f_Jmag_1(i__),n_Jmag_1(i__),Kxpos_1(i__),
     +  Kypos_1(i__),Kfield_1(i__),Kstrip_1(i__),Kmag_1(i__),
     +  e_Kmag_1(i__),f_Kmag_1(i__),n_Kmag_1(i__),Bmag_1(i__),
     +  Rmag_1(i__)
        RAdeg_1(i__) = rNULL__
        DEdeg_1(i__) = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_1 from input data
c  (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
        if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
        if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  DCMC_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
     +  DEd_1(i__),DEm_1(i__),DEs_1(i__),PosErr_1(i__),Ixpos_1(i__),
     +  Iypos_1(i__),Ifield_1(i__),Istrip_1(i__),Imag_1(i__),
     +  e_Imag_1(i__),f_Imag_1(i__),n_Imag_1(i__),Jxpos_1(i__),
     +  Jypos_1(i__),Jfield_1(i__),Jstrip_1(i__),Jmag_1(i__),
     +  e_Jmag_1(i__),f_Jmag_1(i__),n_Jmag_1(i__),Kxpos_1(i__),
     +  Kypos_1(i__),Kfield_1(i__),Kstrip_1(i__),Kmag_1(i__),
     +  e_Kmag_1(i__),f_Kmag_1(i__),n_Kmag_1(i__),Bmag_1(i__),
     +  Rmag_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'lmc_ik.dat'	! LMC stars observed in I and K bands (1)

C  Format for file interpretation

    3 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'lmc_ik.dat')
      write(6,*) '....Loading file: lmc_ik.dat'
      do i__=1,8234
        read(1,'(A210)')ar__2
        read(ar__2,3)
     +  DCMC_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
     +  DEd_2(i__),DEm_2(i__),DEs_2(i__),PosErr_2(i__),Ixpos_2(i__),
     +  Iypos_2(i__),Ifield_2(i__),Istrip_2(i__),Imag_2(i__),
     +  e_Imag_2(i__),f_Imag_2(i__),n_Imag_2(i__),Jxpos_2(i__),
     +  Jypos_2(i__),Jfield_2(i__),Jstrip_2(i__),Jmag_2(i__),
     +  e_Jmag_2(i__),f_Jmag_2(i__),n_Jmag_2(i__),Kxpos_2(i__),
     +  Kypos_2(i__),Kfield_2(i__),Kstrip_2(i__),Kmag_2(i__),
     +  e_Kmag_2(i__),f_Kmag_2(i__),n_Kmag_2(i__),Bmag_2(i__),
     +  Rmag_2(i__)
        RAdeg_2(i__) = rNULL__
        DEdeg_2(i__) = rNULL__
c  Derive coordinates RAdeg_2 and DEdeg_2 from input data
c  (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
        if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
        if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
        if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
        if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
        if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
        if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
        if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  DCMC_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
     +  DEd_2(i__),DEm_2(i__),DEs_2(i__),PosErr_2(i__),Ixpos_2(i__),
     +  Iypos_2(i__),Ifield_2(i__),Istrip_2(i__),Imag_2(i__),
     +  e_Imag_2(i__),f_Imag_2(i__),n_Imag_2(i__),Jxpos_2(i__),
     +  Jypos_2(i__),Jfield_2(i__),Jstrip_2(i__),Jmag_2(i__),
     +  e_Jmag_2(i__),f_Jmag_2(i__),n_Jmag_2(i__),Kxpos_2(i__),
     +  Kypos_2(i__),Kfield_2(i__),Kstrip_2(i__),Kmag_2(i__),
     +  e_Kmag_2(i__),f_Kmag_2(i__),n_Kmag_2(i__),Bmag_2(i__),
     +  Rmag_2(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'lmc_jk.dat'	! LMC stars observed in J and K bands (1)

C  Format for file interpretation

    4 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'lmc_jk.dat')
      write(6,*) '....Loading file: lmc_jk.dat'
      do i__=1,1689
        read(1,'(A210)')ar__3
        read(ar__3,4)
     +  DCMC_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),
     +  DEd_3(i__),DEm_3(i__),DEs_3(i__),PosErr_3(i__),Ixpos_3(i__),
     +  Iypos_3(i__),Ifield_3(i__),Istrip_3(i__),Imag_3(i__),
     +  e_Imag_3(i__),f_Imag_3(i__),n_Imag_3(i__),Jxpos_3(i__),
     +  Jypos_3(i__),Jfield_3(i__),Jstrip_3(i__),Jmag_3(i__),
     +  e_Jmag_3(i__),f_Jmag_3(i__),n_Jmag_3(i__),Kxpos_3(i__),
     +  Kypos_3(i__),Kfield_3(i__),Kstrip_3(i__),Kmag_3(i__),
     +  e_Kmag_3(i__),f_Kmag_3(i__),n_Kmag_3(i__),Bmag_3(i__),
     +  Rmag_3(i__)
        RAdeg_3(i__) = rNULL__
        DEdeg_3(i__) = rNULL__
c  Derive coordinates RAdeg_3 and DEdeg_3 from input data
c  (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
        if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
        if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
        if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
        if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
        if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
        if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
        if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  DCMC_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),
     +  DEd_3(i__),DEm_3(i__),DEs_3(i__),PosErr_3(i__),Ixpos_3(i__),
     +  Iypos_3(i__),Ifield_3(i__),Istrip_3(i__),Imag_3(i__),
     +  e_Imag_3(i__),f_Imag_3(i__),n_Imag_3(i__),Jxpos_3(i__),
     +  Jypos_3(i__),Jfield_3(i__),Jstrip_3(i__),Jmag_3(i__),
     +  e_Jmag_3(i__),f_Jmag_3(i__),n_Jmag_3(i__),Kxpos_3(i__),
     +  Kypos_3(i__),Kfield_3(i__),Kstrip_3(i__),Kmag_3(i__),
     +  e_Kmag_3(i__),f_Kmag_3(i__),n_Kmag_3(i__),Bmag_3(i__),
     +  Rmag_3(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'smc_ijk.dat'	! SMC stars observed in the 3 bands I, J and K (1)

C  Format for file interpretation

    5 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'smc_ijk.dat')
      write(6,*) '....Loading file: smc_ijk.dat'
      do i__=1,69829
        read(1,'(A210)')ar__4
        read(ar__4,5)
     +  DCMC_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),DE__4(i__),
     +  DEd_4(i__),DEm_4(i__),DEs_4(i__),PosErr_4(i__),Ixpos_4(i__),
     +  Iypos_4(i__),Ifield_4(i__),Istrip_4(i__),Imag_4(i__),
     +  e_Imag_4(i__),f_Imag_4(i__),n_Imag_4(i__),Jxpos_4(i__),
     +  Jypos_4(i__),Jfield_4(i__),Jstrip_4(i__),Jmag_4(i__),
     +  e_Jmag_4(i__),f_Jmag_4(i__),n_Jmag_4(i__),Kxpos_4(i__),
     +  Kypos_4(i__),Kfield_4(i__),Kstrip_4(i__),Kmag_4(i__),
     +  e_Kmag_4(i__),f_Kmag_4(i__),n_Kmag_4(i__),Bmag_4(i__),
     +  Rmag_4(i__)
        RAdeg_4(i__) = rNULL__
        DEdeg_4(i__) = rNULL__
c  Derive coordinates RAdeg_4 and DEdeg_4 from input data
c  (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown)
        if(RAh_4(i__) .GT. -180) RAdeg_4(i__)=RAh_4(i__)*15.
        if(RAm_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAm_4(i__)/4.
        if(RAs_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAs_4(i__)/240.
        if(DEd_4(i__) .GE. 0) DEdeg_4(i__)=DEd_4(i__)
        if(DEm_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm_4(i__)/60.
        if(DEs_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs_4(i__)/3600.
        if(DE__4(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_4(i__)
c    ..............Just test output...........
        write(6,5)
     +  DCMC_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),DE__4(i__),
     +  DEd_4(i__),DEm_4(i__),DEs_4(i__),PosErr_4(i__),Ixpos_4(i__),
     +  Iypos_4(i__),Ifield_4(i__),Istrip_4(i__),Imag_4(i__),
     +  e_Imag_4(i__),f_Imag_4(i__),n_Imag_4(i__),Jxpos_4(i__),
     +  Jypos_4(i__),Jfield_4(i__),Jstrip_4(i__),Jmag_4(i__),
     +  e_Jmag_4(i__),f_Jmag_4(i__),n_Jmag_4(i__),Kxpos_4(i__),
     +  Kypos_4(i__),Kfield_4(i__),Kstrip_4(i__),Kmag_4(i__),
     +  e_Kmag_4(i__),f_Kmag_4(i__),n_Kmag_4(i__),Bmag_4(i__),
     +  Rmag_4(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_4(i__),DEdeg_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'smc_ij.dat'	! SMC stars observed in I and J bands (1)

C  Format for file interpretation

    6 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'smc_ij.dat')
      write(6,*) '....Loading file: smc_ij.dat'
      do i__=1,242837
        read(1,'(A210)')ar__5
        read(ar__5,6)
     +  DCMC_5(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),DE__5(i__),
     +  DEd_5(i__),DEm_5(i__),DEs_5(i__),PosErr_5(i__),Ixpos_5(i__),
     +  Iypos_5(i__),Ifield_5(i__),Istrip_5(i__),Imag_5(i__),
     +  e_Imag_5(i__),f_Imag_5(i__),n_Imag_5(i__),Jxpos_5(i__),
     +  Jypos_5(i__),Jfield_5(i__),Jstrip_5(i__),Jmag_5(i__),
     +  e_Jmag_5(i__),f_Jmag_5(i__),n_Jmag_5(i__),Kxpos_5(i__),
     +  Kypos_5(i__),Kfield_5(i__),Kstrip_5(i__),Kmag_5(i__),
     +  e_Kmag_5(i__),f_Kmag_5(i__),n_Kmag_5(i__),Bmag_5(i__),
     +  Rmag_5(i__)
        RAdeg_5(i__) = rNULL__
        DEdeg_5(i__) = rNULL__
c  Derive coordinates RAdeg_5 and DEdeg_5 from input data
c  (RAdeg_5 and DEdeg_5 are set to rNULL__ when unknown)
        if(RAh_5(i__) .GT. -180) RAdeg_5(i__)=RAh_5(i__)*15.
        if(RAm_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAm_5(i__)/4.
        if(RAs_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAs_5(i__)/240.
        if(DEd_5(i__) .GE. 0) DEdeg_5(i__)=DEd_5(i__)
        if(DEm_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEm_5(i__)/60.
        if(DEs_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEs_5(i__)/3600.
        if(DE__5(i__).EQ.'-'.AND.DEdeg_5(i__).GE.0) DEdeg_5(i__)=-DEdeg_5(i__)
c    ..............Just test output...........
        write(6,6)
     +  DCMC_5(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),DE__5(i__),
     +  DEd_5(i__),DEm_5(i__),DEs_5(i__),PosErr_5(i__),Ixpos_5(i__),
     +  Iypos_5(i__),Ifield_5(i__),Istrip_5(i__),Imag_5(i__),
     +  e_Imag_5(i__),f_Imag_5(i__),n_Imag_5(i__),Jxpos_5(i__),
     +  Jypos_5(i__),Jfield_5(i__),Jstrip_5(i__),Jmag_5(i__),
     +  e_Jmag_5(i__),f_Jmag_5(i__),n_Jmag_5(i__),Kxpos_5(i__),
     +  Kypos_5(i__),Kfield_5(i__),Kstrip_5(i__),Kmag_5(i__),
     +  e_Kmag_5(i__),f_Kmag_5(i__),n_Kmag_5(i__),Bmag_5(i__),
     +  Rmag_5(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_5(i__),DEdeg_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'smc_ik.dat'	! SMC stars observed in I and K bands (1)

C  Format for file interpretation

    7 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'smc_ik.dat')
      write(6,*) '....Loading file: smc_ik.dat'
      do i__=1,2302
        read(1,'(A210)')ar__6
        read(ar__6,7)
     +  DCMC_6(i__),RAh_6(i__),RAm_6(i__),RAs_6(i__),DE__6(i__),
     +  DEd_6(i__),DEm_6(i__),DEs_6(i__),PosErr_6(i__),Ixpos_6(i__),
     +  Iypos_6(i__),Ifield_6(i__),Istrip_6(i__),Imag_6(i__),
     +  e_Imag_6(i__),f_Imag_6(i__),n_Imag_6(i__),Jxpos_6(i__),
     +  Jypos_6(i__),Jfield_6(i__),Jstrip_6(i__),Jmag_6(i__),
     +  e_Jmag_6(i__),f_Jmag_6(i__),n_Jmag_6(i__),Kxpos_6(i__),
     +  Kypos_6(i__),Kfield_6(i__),Kstrip_6(i__),Kmag_6(i__),
     +  e_Kmag_6(i__),f_Kmag_6(i__),n_Kmag_6(i__),Bmag_6(i__),
     +  Rmag_6(i__)
        RAdeg_6(i__) = rNULL__
        DEdeg_6(i__) = rNULL__
c  Derive coordinates RAdeg_6 and DEdeg_6 from input data
c  (RAdeg_6 and DEdeg_6 are set to rNULL__ when unknown)
        if(RAh_6(i__) .GT. -180) RAdeg_6(i__)=RAh_6(i__)*15.
        if(RAm_6(i__) .GT. -180) RAdeg_6(i__)=RAdeg_6(i__)+RAm_6(i__)/4.
        if(RAs_6(i__) .GT. -180) RAdeg_6(i__)=RAdeg_6(i__)+RAs_6(i__)/240.
        if(DEd_6(i__) .GE. 0) DEdeg_6(i__)=DEd_6(i__)
        if(DEm_6(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEm_6(i__)/60.
        if(DEs_6(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEs_6(i__)/3600.
        if(DE__6(i__).EQ.'-'.AND.DEdeg_6(i__).GE.0) DEdeg_6(i__)=-DEdeg_6(i__)
c    ..............Just test output...........
        write(6,7)
     +  DCMC_6(i__),RAh_6(i__),RAm_6(i__),RAs_6(i__),DE__6(i__),
     +  DEd_6(i__),DEm_6(i__),DEs_6(i__),PosErr_6(i__),Ixpos_6(i__),
     +  Iypos_6(i__),Ifield_6(i__),Istrip_6(i__),Imag_6(i__),
     +  e_Imag_6(i__),f_Imag_6(i__),n_Imag_6(i__),Jxpos_6(i__),
     +  Jypos_6(i__),Jfield_6(i__),Jstrip_6(i__),Jmag_6(i__),
     +  e_Jmag_6(i__),f_Jmag_6(i__),n_Jmag_6(i__),Kxpos_6(i__),
     +  Kypos_6(i__),Kfield_6(i__),Kstrip_6(i__),Kmag_6(i__),
     +  e_Kmag_6(i__),f_Kmag_6(i__),n_Kmag_6(i__),Bmag_6(i__),
     +  Rmag_6(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_6(i__),DEdeg_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'smc_jk.dat'	! SMC stars observed in J and K bands (1)

C  Format for file interpretation

    8 format(
     +  5X,A19,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.2,1X,
     +  F7.3,1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,
     +  1X,F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,
     +  F7.3,1X,I3,1X,I4,1X,F6.3,1X,F6.3,1X,I3,1X,I3,2X,F5.2,3X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'smc_jk.dat')
      write(6,*) '....Loading file: smc_jk.dat'
      do i__=1,812
        read(1,'(A210)')ar__7
        read(ar__7,8)
     +  DCMC_7(i__),RAh_7(i__),RAm_7(i__),RAs_7(i__),DE__7(i__),
     +  DEd_7(i__),DEm_7(i__),DEs_7(i__),PosErr_7(i__),Ixpos_7(i__),
     +  Iypos_7(i__),Ifield_7(i__),Istrip_7(i__),Imag_7(i__),
     +  e_Imag_7(i__),f_Imag_7(i__),n_Imag_7(i__),Jxpos_7(i__),
     +  Jypos_7(i__),Jfield_7(i__),Jstrip_7(i__),Jmag_7(i__),
     +  e_Jmag_7(i__),f_Jmag_7(i__),n_Jmag_7(i__),Kxpos_7(i__),
     +  Kypos_7(i__),Kfield_7(i__),Kstrip_7(i__),Kmag_7(i__),
     +  e_Kmag_7(i__),f_Kmag_7(i__),n_Kmag_7(i__),Bmag_7(i__),
     +  Rmag_7(i__)
        RAdeg_7(i__) = rNULL__
        DEdeg_7(i__) = rNULL__
c  Derive coordinates RAdeg_7 and DEdeg_7 from input data
c  (RAdeg_7 and DEdeg_7 are set to rNULL__ when unknown)
        if(RAh_7(i__) .GT. -180) RAdeg_7(i__)=RAh_7(i__)*15.
        if(RAm_7(i__) .GT. -180) RAdeg_7(i__)=RAdeg_7(i__)+RAm_7(i__)/4.
        if(RAs_7(i__) .GT. -180) RAdeg_7(i__)=RAdeg_7(i__)+RAs_7(i__)/240.
        if(DEd_7(i__) .GE. 0) DEdeg_7(i__)=DEd_7(i__)
        if(DEm_7(i__) .GE. 0) DEdeg_7(i__)=DEdeg_7(i__)+DEm_7(i__)/60.
        if(DEs_7(i__) .GE. 0) DEdeg_7(i__)=DEdeg_7(i__)+DEs_7(i__)/3600.
        if(DE__7(i__).EQ.'-'.AND.DEdeg_7(i__).GE.0) DEdeg_7(i__)=-DEdeg_7(i__)
c    ..............Just test output...........
        write(6,8)
     +  DCMC_7(i__),RAh_7(i__),RAm_7(i__),RAs_7(i__),DE__7(i__),
     +  DEd_7(i__),DEm_7(i__),DEs_7(i__),PosErr_7(i__),Ixpos_7(i__),
     +  Iypos_7(i__),Ifield_7(i__),Istrip_7(i__),Imag_7(i__),
     +  e_Imag_7(i__),f_Imag_7(i__),n_Imag_7(i__),Jxpos_7(i__),
     +  Jypos_7(i__),Jfield_7(i__),Jstrip_7(i__),Jmag_7(i__),
     +  e_Jmag_7(i__),f_Jmag_7(i__),n_Jmag_7(i__),Kxpos_7(i__),
     +  Kypos_7(i__),Kfield_7(i__),Kstrip_7(i__),Kmag_7(i__),
     +  e_Kmag_7(i__),f_Kmag_7(i__),n_Kmag_7(i__),Bmag_7(i__),
     +  Rmag_7(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_7(i__),DEdeg_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'qual.dat'	! Qualities of observation strips (table 8)

C  Format for file interpretation

    9 format(
     +  I4,1X,A8,4X,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,2X,
     +  A15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'qual.dat')
      write(6,*) '....Loading file: qual.dat'
      do i__=1,187
        read(1,'(A68)')ar__8
        read(ar__8,9)
     +  Strip(i__),Date(i__),Izp(i__),e_Izp(i__),Jzp(i__),e_Jzp(i__),
     +  Kzp(i__),e_Kzp(i__),Notes(i__)
c    ..............Just test output...........
        write(6,9)
     +  Strip(i__),Date(i__),Izp(i__),e_Izp(i__),Jzp(i__),e_Jzp(i__),
     +  Kzp(i__),e_Kzp(i__),Notes(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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