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

Conversion of standardized ReadMe file for file /./ftp/cats/I/346 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-Mar-29
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/346     FON Astrographic Catalogue Southern Part (FONAC-S)  (Yuldoshev+, 2017)
*================================================================================
*Catalog of Equatorial Coordinates and B- Magnitudes of Stars of the Kitab Part
*of the FON Project.
*    Yuldoshev Q.X., Ehgamberdiev Sh.A., Muminov M.M., Protsyuk Yu.I., Relke H.,
*    Andruk V.M.
*    <Kinem. and Physics of Celest. Bodies. 33, 250 (2017)>
*    =2017KPCB...33..250Y
*    =2017KFNT...33..250Y
*    =2018yCat.1346....0Y
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'fs00.dat'	! Star data file (zone from +02.0 to -01.9 degrees)

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

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

      integer*4     ID         (nr__) ! Serial number of the star
      integer*4     k          (nr__) ! Number of determinations (plates)
      real*8        RAhour     (nr__) ! (h) Right Ascension (J2000)
      real*4        e_RAhour   (nr__) ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_1    (nr__) ! (deg) Declination (J2000)
      real*4        e_DEdeg    (nr__) ! (arcsec) Error in Declination
      real*4        Bmag       (nr__) ! (mag) Johnson Blue magnitude
      real*4        e_Bmag     (nr__) ! (mag) Error in B magnitude
      real*8        Epoch      (nr__) ! (yr) Epoch of observations
      real*4        FWHM       (nr__) ! (1.46pix) Diameter of star image
      real*8        cI         (nr__) ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh        (nr__) ! (h) Right ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Right ascension (J2000)
      real*4        RAs        (nr__) ! (s) [0/60] Right ascension (J2000)
      character*1   DE_        (nr__) ! Declination sign (J2000)
      integer*4     DEd        (nr__) ! (deg) Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Declination (J2000)
      real*4        DEs        (nr__) ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag      (nr__) ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA       (nr__) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE       (nr__) ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag      (nr__) ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag       (nr__) ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag       (nr__) ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag       (nr__) ! (mag) ?=0 UCAC4 i magnitude (1)
*Note (1):  Proper motions and magnitudes were taken from UCAC4 catalogue.
*           Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z,
*           Cat. I/322).

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

C  Declarations for 'fs04.dat'	! Star data file (zone from -02.0 to -05.9 degrees)

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

C  J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs
      real*8        RAdeg_1    (nr__1) ! (deg) Right Ascension J2000
      real*8        DEdeg_2    (nr__1) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_1       (nr__1) ! Serial number of the star
      integer*4     k_1        (nr__1) ! Number of determinations (plates)
      real*8        RAhour_1   (nr__1) ! (h) Right Ascension (J2000)
      real*4        e_RAhour_1 (nr__1) ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_3    (nr__1) ! (deg) Declination (J2000)
      real*4        e_DEdeg_1  (nr__1) ! (arcsec) Error in Declination
      real*4        Bmag_1     (nr__1) ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_1   (nr__1) ! (mag) Error in B magnitude
      real*8        Epoch_1    (nr__1) ! (yr) Epoch of observations
      real*4        FWHM_1     (nr__1) ! (1.46pix) Diameter of star image
      real*8        cI_1       (nr__1) ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_1      (nr__1) ! (h) Right ascension (J2000)
      integer*4     RAm_1      (nr__1) ! (min) Right ascension (J2000)
      real*4        RAs_1      (nr__1) ! (s) [0/60] Right ascension (J2000)
      character*1   DE__1      (nr__1) ! Declination sign (J2000)
      integer*4     DEd_1      (nr__1) ! (deg) Declination (J2000)
      integer*4     DEm_1      (nr__1) ! (arcmin) Declination (J2000)
      real*4        DEs_1      (nr__1) ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_1    (nr__1) ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_1     (nr__1) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_1     (nr__1) ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_1    (nr__1) ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_1     (nr__1) ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_1     (nr__1) ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_1     (nr__1) ! (mag) ?=0 UCAC4 i magnitude (1)
*Note (1):  Proper motions and magnitudes were taken from UCAC4 catalogue.
*           Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z,
*           Cat. I/322).

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

C  Declarations for 'fs08.dat'	! Star data file (zone from -06.0 to -09.9 degrees)

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

C  J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs
      real*8        RAdeg_2    (nr__2) ! (deg) Right Ascension J2000
      real*8        DEdeg_4    (nr__2) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_2       (nr__2) ! Serial number of the star
      integer*4     k_2        (nr__2) ! Number of determinations (plates)
      real*8        RAhour_2   (nr__2) ! (h) Right Ascension (J2000)
      real*4        e_RAhour_2 (nr__2) ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_5    (nr__2) ! (deg) Declination (J2000)
      real*4        e_DEdeg_2  (nr__2) ! (arcsec) Error in Declination
      real*4        Bmag_2     (nr__2) ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_2   (nr__2) ! (mag) Error in B magnitude
      real*8        Epoch_2    (nr__2) ! (yr) Epoch of observations
      real*4        FWHM_2     (nr__2) ! (1.46pix) Diameter of star image
      real*8        cI_2       (nr__2) ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_2      (nr__2) ! (h) Right ascension (J2000)
      integer*4     RAm_2      (nr__2) ! (min) Right ascension (J2000)
      real*4        RAs_2      (nr__2) ! (s) [0/60] Right ascension (J2000)
      character*1   DE__2      (nr__2) ! Declination sign (J2000)
      integer*4     DEd_2      (nr__2) ! (deg) Declination (J2000)
      integer*4     DEm_2      (nr__2) ! (arcmin) Declination (J2000)
      real*4        DEs_2      (nr__2) ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_2    (nr__2) ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_2     (nr__2) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_2     (nr__2) ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_2    (nr__2) ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_2     (nr__2) ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_2     (nr__2) ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_2     (nr__2) ! (mag) ?=0 UCAC4 i magnitude (1)
*Note (1):  Proper motions and magnitudes were taken from UCAC4 catalogue.
*           Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z,
*           Cat. I/322).

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

C  Declarations for 'fs12.dat'	! Star data file (zone from -10.0 to -13.9 degrees)

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

C  J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs
      real*8        RAdeg_3    (nr__3) ! (deg) Right Ascension J2000
      real*8        DEdeg_6    (nr__3) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_3       (nr__3) ! Serial number of the star
      integer*4     k_3        (nr__3) ! Number of determinations (plates)
      real*8        RAhour_3   (nr__3) ! (h) Right Ascension (J2000)
      real*4        e_RAhour_3 (nr__3) ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_7    (nr__3) ! (deg) Declination (J2000)
      real*4        e_DEdeg_3  (nr__3) ! (arcsec) Error in Declination
      real*4        Bmag_3     (nr__3) ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_3   (nr__3) ! (mag) Error in B magnitude
      real*8        Epoch_3    (nr__3) ! (yr) Epoch of observations
      real*4        FWHM_3     (nr__3) ! (1.46pix) Diameter of star image
      real*8        cI_3       (nr__3) ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_3      (nr__3) ! (h) Right ascension (J2000)
      integer*4     RAm_3      (nr__3) ! (min) Right ascension (J2000)
      real*4        RAs_3      (nr__3) ! (s) [0/60] Right ascension (J2000)
      character*1   DE__3      (nr__3) ! Declination sign (J2000)
      integer*4     DEd_3      (nr__3) ! (deg) Declination (J2000)
      integer*4     DEm_3      (nr__3) ! (arcmin) Declination (J2000)
      real*4        DEs_3      (nr__3) ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_3    (nr__3) ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_3     (nr__3) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_3     (nr__3) ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_3    (nr__3) ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_3     (nr__3) ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_3     (nr__3) ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_3     (nr__3) ! (mag) ?=0 UCAC4 i magnitude (1)
*Note (1):  Proper motions and magnitudes were taken from UCAC4 catalogue.
*           Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z,
*           Cat. I/322).

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

C  Declarations for 'fs16.dat'	! Star data file (zone from -14.0 to -17.9 degrees)

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

C  J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs
      real*8        RAdeg_4    (nr__4) ! (deg) Right Ascension J2000
      real*8        DEdeg_8    (nr__4) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_4       (nr__4) ! Serial number of the star
      integer*4     k_4        (nr__4) ! Number of determinations (plates)
      real*8        RAhour_4   (nr__4) ! (h) Right Ascension (J2000)
      real*4        e_RAhour_4 (nr__4) ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_9    (nr__4) ! (deg) Declination (J2000)
      real*4        e_DEdeg_4  (nr__4) ! (arcsec) Error in Declination
      real*4        Bmag_4     (nr__4) ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_4   (nr__4) ! (mag) Error in B magnitude
      real*8        Epoch_4    (nr__4) ! (yr) Epoch of observations
      real*4        FWHM_4     (nr__4) ! (1.46pix) Diameter of star image
      real*8        cI_4       (nr__4) ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_4      (nr__4) ! (h) Right ascension (J2000)
      integer*4     RAm_4      (nr__4) ! (min) Right ascension (J2000)
      real*4        RAs_4      (nr__4) ! (s) [0/60] Right ascension (J2000)
      character*1   DE__4      (nr__4) ! Declination sign (J2000)
      integer*4     DEd_4      (nr__4) ! (deg) Declination (J2000)
      integer*4     DEm_4      (nr__4) ! (arcmin) Declination (J2000)
      real*4        DEs_4      (nr__4) ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_4    (nr__4) ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_4     (nr__4) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_4     (nr__4) ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_4    (nr__4) ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_4     (nr__4) ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_4     (nr__4) ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_4     (nr__4) ! (mag) ?=0 UCAC4 i magnitude (1)
*Note (1):  Proper motions and magnitudes were taken from UCAC4 catalogue.
*           Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z,
*           Cat. I/322).

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

C  Declarations for 'fs20.dat'	! Star data file (zone from -18.0 to -20.9 degrees)

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

C  J2000 position composed of: RAhour DE- DEdeg DEd DEm DEs
      real*8        RAdeg_5    (nr__5) ! (deg) Right Ascension J2000
      real*8        DEdeg_10   (nr__5) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_5       (nr__5) ! Serial number of the star
      integer*4     k_5        (nr__5) ! Number of determinations (plates)
      real*8        RAhour_5   (nr__5) ! (h) Right Ascension (J2000)
      real*4        e_RAhour_5 (nr__5) ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_11   (nr__5) ! (deg) Declination (J2000)
      real*4        e_DEdeg_5  (nr__5) ! (arcsec) Error in Declination
      real*4        Bmag_5     (nr__5) ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_5   (nr__5) ! (mag) Error in B magnitude
      real*8        Epoch_5    (nr__5) ! (yr) Epoch of observations
      real*4        FWHM_5     (nr__5) ! (1.46pix) Diameter of star image
      real*8        cI_5       (nr__5) ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_5      (nr__5) ! (h) Right ascension (J2000)
      integer*4     RAm_5      (nr__5) ! (min) Right ascension (J2000)
      real*4        RAs_5      (nr__5) ! (s) [0/60] Right ascension (J2000)
      character*1   DE__5      (nr__5) ! Declination sign (J2000)
      integer*4     DEd_5      (nr__5) ! (deg) Declination (J2000)
      integer*4     DEm_5      (nr__5) ! (arcmin) Declination (J2000)
      real*4        DEs_5      (nr__5) ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_5    (nr__5) ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_5     (nr__5) ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_5     (nr__5) ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_5    (nr__5) ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_5     (nr__5) ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_5     (nr__5) ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_5     (nr__5) ! (mag) ?=0 UCAC4 i magnitude (1)
*Note (1):  Proper motions and magnitudes were taken from UCAC4 catalogue.
*           Zacharias N., Finch C.T., Girard T.M., et al. 2013AJ....145...44Z,
*           Cat. I/322).

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

C  Loading file 'fs00.dat'	! Star data file (zone from +02.0 to -01.9 degrees)

C  Format for file interpretation

    1 format(
     +  I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X,
     +  F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3,
     +  2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fs00.dat')
      write(6,*) '....Loading file: fs00.dat'
      do i__=1,2550586
        read(1,'(A164)')ar__
        read(ar__,1)
     +  ID(i__),k(i__),RAhour(i__),e_RAhour(i__),DEdeg_1(i__),
     +  e_DEdeg(i__),Bmag(i__),e_Bmag(i__),Epoch(i__),FWHM(i__),
     +  cI(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
     +  DEs(i__),f_mag(i__),pmRA(i__),pmDE(i__),B4mag(i__),Vmag(i__),
     +  rmag(i__),imag(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(RAhour(i__) .GT. -180) RAdeg(i__)=RAhour(i__)*15.
        if(DEdeg_1(i__) .GE. 0) DEdeg(i__)=DEdeg_1(i__)
        if(DEd(i__) .GE. 0) DEdeg(i__)=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)
     +  ID(i__),k(i__),RAhour(i__),e_RAhour(i__),DEdeg_1(i__),
     +  e_DEdeg(i__),Bmag(i__),e_Bmag(i__),Epoch(i__),FWHM(i__),
     +  cI(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
     +  DEs(i__),f_mag(i__),pmRA(i__),pmDE(i__),B4mag(i__),Vmag(i__),
     +  rmag(i__),imag(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 'fs04.dat'	! Star data file (zone from -02.0 to -05.9 degrees)

C  Format for file interpretation

    2 format(
     +  I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X,
     +  F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3,
     +  2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fs04.dat')
      write(6,*) '....Loading file: fs04.dat'
      do i__=1,2233270
        read(1,'(A164)')ar__1
        read(ar__1,2)
     +  ID_1(i__),k_1(i__),RAhour_1(i__),e_RAhour_1(i__),DEdeg_3(i__),
     +  e_DEdeg_1(i__),Bmag_1(i__),e_Bmag_1(i__),Epoch_1(i__),
     +  FWHM_1(i__),cI_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),f_mag_1(i__),
     +  pmRA_1(i__),pmDE_1(i__),B4mag_1(i__),Vmag_1(i__),rmag_1(i__),
     +  imag_1(i__)
        RAdeg_1(i__) = rNULL__
        DEdeg_2(i__) = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_2 from input data
c  (RAdeg_1 and DEdeg_2 are set to rNULL__ when unknown)
        if(RAhour_1(i__) .GT. -180) RAdeg_1(i__)=RAhour_1(i__)*15.
        if(DEdeg_3(i__) .GE. 0) DEdeg_2(i__)=DEdeg_3(i__)
        if(DEd_1(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c    ..............Just test output...........
        write(6,2)
     +  ID_1(i__),k_1(i__),RAhour_1(i__),e_RAhour_1(i__),DEdeg_3(i__),
     +  e_DEdeg_1(i__),Bmag_1(i__),e_Bmag_1(i__),Epoch_1(i__),
     +  FWHM_1(i__),cI_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),f_mag_1(i__),
     +  pmRA_1(i__),pmDE_1(i__),B4mag_1(i__),Vmag_1(i__),rmag_1(i__),
     +  imag_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fs08.dat'	! Star data file (zone from -06.0 to -09.9 degrees)

C  Format for file interpretation

    3 format(
     +  I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X,
     +  F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3,
     +  2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fs08.dat')
      write(6,*) '....Loading file: fs08.dat'
      do i__=1,2486306
        read(1,'(A164)')ar__2
        read(ar__2,3)
     +  ID_2(i__),k_2(i__),RAhour_2(i__),e_RAhour_2(i__),DEdeg_5(i__),
     +  e_DEdeg_2(i__),Bmag_2(i__),e_Bmag_2(i__),Epoch_2(i__),
     +  FWHM_2(i__),cI_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),f_mag_2(i__),
     +  pmRA_2(i__),pmDE_2(i__),B4mag_2(i__),Vmag_2(i__),rmag_2(i__),
     +  imag_2(i__)
        RAdeg_2(i__) = rNULL__
        DEdeg_4(i__) = rNULL__
c  Derive coordinates RAdeg_2 and DEdeg_4 from input data
c  (RAdeg_2 and DEdeg_4 are set to rNULL__ when unknown)
        if(RAhour_2(i__) .GT. -180) RAdeg_2(i__)=RAhour_2(i__)*15.
        if(DEdeg_5(i__) .GE. 0) DEdeg_4(i__)=DEdeg_5(i__)
        if(DEd_2(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEd_2(i__)
        if(DEm_2(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm_2(i__)/60.
        if(DEs_2(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs_2(i__)/3600.
        if(DE__2(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_4(i__)
c    ..............Just test output...........
        write(6,3)
     +  ID_2(i__),k_2(i__),RAhour_2(i__),e_RAhour_2(i__),DEdeg_5(i__),
     +  e_DEdeg_2(i__),Bmag_2(i__),e_Bmag_2(i__),Epoch_2(i__),
     +  FWHM_2(i__),cI_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),f_mag_2(i__),
     +  pmRA_2(i__),pmDE_2(i__),B4mag_2(i__),Vmag_2(i__),rmag_2(i__),
     +  imag_2(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fs12.dat'	! Star data file (zone from -10.0 to -13.9 degrees)

C  Format for file interpretation

    4 format(
     +  I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X,
     +  F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3,
     +  2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fs12.dat')
      write(6,*) '....Loading file: fs12.dat'
      do i__=1,2416612
        read(1,'(A164)')ar__3
        read(ar__3,4)
     +  ID_3(i__),k_3(i__),RAhour_3(i__),e_RAhour_3(i__),DEdeg_7(i__),
     +  e_DEdeg_3(i__),Bmag_3(i__),e_Bmag_3(i__),Epoch_3(i__),
     +  FWHM_3(i__),cI_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),
     +  DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),f_mag_3(i__),
     +  pmRA_3(i__),pmDE_3(i__),B4mag_3(i__),Vmag_3(i__),rmag_3(i__),
     +  imag_3(i__)
        RAdeg_3(i__) = rNULL__
        DEdeg_6(i__) = rNULL__
c  Derive coordinates RAdeg_3 and DEdeg_6 from input data
c  (RAdeg_3 and DEdeg_6 are set to rNULL__ when unknown)
        if(RAhour_3(i__) .GT. -180) RAdeg_3(i__)=RAhour_3(i__)*15.
        if(DEdeg_7(i__) .GE. 0) DEdeg_6(i__)=DEdeg_7(i__)
        if(DEd_3(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEd_3(i__)
        if(DEm_3(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEm_3(i__)/60.
        if(DEs_3(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEs_3(i__)/3600.
        if(DE__3(i__).EQ.'-'.AND.DEdeg_6(i__).GE.0) DEdeg_6(i__)=-DEdeg_6(i__)
c    ..............Just test output...........
        write(6,4)
     +  ID_3(i__),k_3(i__),RAhour_3(i__),e_RAhour_3(i__),DEdeg_7(i__),
     +  e_DEdeg_3(i__),Bmag_3(i__),e_Bmag_3(i__),Epoch_3(i__),
     +  FWHM_3(i__),cI_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),
     +  DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),f_mag_3(i__),
     +  pmRA_3(i__),pmDE_3(i__),B4mag_3(i__),Vmag_3(i__),rmag_3(i__),
     +  imag_3(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fs16.dat'	! Star data file (zone from -14.0 to -17.9 degrees)

C  Format for file interpretation

    5 format(
     +  I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X,
     +  F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3,
     +  2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fs16.dat')
      write(6,*) '....Loading file: fs16.dat'
      do i__=1,2575447
        read(1,'(A164)')ar__4
        read(ar__4,5)
     +  ID_4(i__),k_4(i__),RAhour_4(i__),e_RAhour_4(i__),DEdeg_9(i__),
     +  e_DEdeg_4(i__),Bmag_4(i__),e_Bmag_4(i__),Epoch_4(i__),
     +  FWHM_4(i__),cI_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),
     +  DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),f_mag_4(i__),
     +  pmRA_4(i__),pmDE_4(i__),B4mag_4(i__),Vmag_4(i__),rmag_4(i__),
     +  imag_4(i__)
        RAdeg_4(i__) = rNULL__
        DEdeg_8(i__) = rNULL__
c  Derive coordinates RAdeg_4 and DEdeg_8 from input data
c  (RAdeg_4 and DEdeg_8 are set to rNULL__ when unknown)
        if(RAhour_4(i__) .GT. -180) RAdeg_4(i__)=RAhour_4(i__)*15.
        if(DEdeg_9(i__) .GE. 0) DEdeg_8(i__)=DEdeg_9(i__)
        if(DEd_4(i__) .GE. 0) DEdeg_8(i__)=DEdeg_8(i__)+DEd_4(i__)
        if(DEm_4(i__) .GE. 0) DEdeg_8(i__)=DEdeg_8(i__)+DEm_4(i__)/60.
        if(DEs_4(i__) .GE. 0) DEdeg_8(i__)=DEdeg_8(i__)+DEs_4(i__)/3600.
        if(DE__4(i__).EQ.'-'.AND.DEdeg_8(i__).GE.0) DEdeg_8(i__)=-DEdeg_8(i__)
c    ..............Just test output...........
        write(6,5)
     +  ID_4(i__),k_4(i__),RAhour_4(i__),e_RAhour_4(i__),DEdeg_9(i__),
     +  e_DEdeg_4(i__),Bmag_4(i__),e_Bmag_4(i__),Epoch_4(i__),
     +  FWHM_4(i__),cI_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),
     +  DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),f_mag_4(i__),
     +  pmRA_4(i__),pmDE_4(i__),B4mag_4(i__),Vmag_4(i__),rmag_4(i__),
     +  imag_4(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_4(i__),DEdeg_8(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fs20.dat'	! Star data file (zone from -18.0 to -20.9 degrees)

C  Format for file interpretation

    6 format(
     +  I7,1X,I2,1X,F12.8,1X,F5.3,2X,F11.7,1X,F5.3,2X,F5.2,1X,F4.2,2X,
     +  F8.3,2X,F6.2,1X,F7.2,2X,I2,I2,F6.3,1X,A1,I2,I2,F5.2,2X,F6.3,
     +  2X,F7.1,2X,F7.1,2X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fs20.dat')
      write(6,*) '....Loading file: fs20.dat'
      do i__=1,1151047
        read(1,'(A164)')ar__5
        read(ar__5,6)
     +  ID_5(i__),k_5(i__),RAhour_5(i__),e_RAhour_5(i__),
     +  DEdeg_11(i__),e_DEdeg_5(i__),Bmag_5(i__),e_Bmag_5(i__),
     +  Epoch_5(i__),FWHM_5(i__),cI_5(i__),RAh_5(i__),RAm_5(i__),
     +  RAs_5(i__),DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),
     +  f_mag_5(i__),pmRA_5(i__),pmDE_5(i__),B4mag_5(i__),Vmag_5(i__),
     +  rmag_5(i__),imag_5(i__)
        RAdeg_5(i__) = rNULL__
        DEdeg_10(i__) = rNULL__
c  Derive coordinates RAdeg_5 and DEdeg_10 from input data
c  (RAdeg_5 and DEdeg_10 are set to rNULL__ when unknown)
        if(RAhour_5(i__) .GT. -180) RAdeg_5(i__)=RAhour_5(i__)*15.
        if(DEdeg_11(i__) .GE. 0) DEdeg_10(i__)=DEdeg_11(i__)
        if(DEd_5(i__) .GE. 0) DEdeg_10(i__)=DEdeg_10(i__)+DEd_5(i__)
        if(DEm_5(i__) .GE. 0) DEdeg_10(i__)=DEdeg_10(i__)+DEm_5(i__)/60.
        if(DEs_5(i__) .GE. 0) DEdeg_10(i__)=DEdeg_10(i__)+DEs_5(i__)/3600.
        if(DE__5(i__).EQ.'-'.AND.DEdeg_10(i__).GE.0) DEdeg_10(i__)=-DEdeg_10(i__)
c    ..............Just test output...........
        write(6,6)
     +  ID_5(i__),k_5(i__),RAhour_5(i__),e_RAhour_5(i__),
     +  DEdeg_11(i__),e_DEdeg_5(i__),Bmag_5(i__),e_Bmag_5(i__),
     +  Epoch_5(i__),FWHM_5(i__),cI_5(i__),RAh_5(i__),RAm_5(i__),
     +  RAs_5(i__),DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),
     +  f_mag_5(i__),pmRA_5(i__),pmDE_5(i__),B4mag_5(i__),Vmag_5(i__),
     +  rmag_5(i__),imag_5(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_5(i__),DEdeg_10(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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