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

Conversion of standardized ReadMe file for file /./ftp/cats/I/346 into FORTRAN code for reading data files line by line.

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-28
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       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID          ! Serial number of the star
      integer*4     k           ! Number of determinations (plates)
      real*8        RAhour      ! (h) Right Ascension (J2000)
      real*4        e_RAhour    ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_1     ! (deg) Declination (J2000)
      real*4        e_DEdeg     ! (arcsec) Error in Declination
      real*4        Bmag        ! (mag) Johnson Blue magnitude
      real*4        e_Bmag      ! (mag) Error in B magnitude
      real*8        Epoch       ! (yr) Epoch of observations
      real*4        FWHM        ! (1.46pix) Diameter of star image
      real*8        cI          ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh         ! (h) Right ascension (J2000)
      integer*4     RAm         ! (min) Right ascension (J2000)
      real*4        RAs         ! (s) [0/60] Right ascension (J2000)
      character*1   DE_         ! Declination sign (J2000)
      integer*4     DEd         ! (deg) Declination (J2000)
      integer*4     DEm         ! (arcmin) Declination (J2000)
      real*4        DEs         ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag       ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA        ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE        ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag       ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag        ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag        ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag        ! (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     ! (deg) Right Ascension J2000
      real*8        DEdeg_2     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_1        ! Serial number of the star
      integer*4     k_1         ! Number of determinations (plates)
      real*8        RAhour_1    ! (h) Right Ascension (J2000)
      real*4        e_RAhour_1  ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_3     ! (deg) Declination (J2000)
      real*4        e_DEdeg_1   ! (arcsec) Error in Declination
      real*4        Bmag_1      ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_1    ! (mag) Error in B magnitude
      real*8        Epoch_1     ! (yr) Epoch of observations
      real*4        FWHM_1      ! (1.46pix) Diameter of star image
      real*8        cI_1        ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_1       ! (h) Right ascension (J2000)
      integer*4     RAm_1       ! (min) Right ascension (J2000)
      real*4        RAs_1       ! (s) [0/60] Right ascension (J2000)
      character*1   DE__1       ! Declination sign (J2000)
      integer*4     DEd_1       ! (deg) Declination (J2000)
      integer*4     DEm_1       ! (arcmin) Declination (J2000)
      real*4        DEs_1       ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_1     ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_1      ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_1      ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_1     ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_1      ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_1      ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_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     ! (deg) Right Ascension J2000
      real*8        DEdeg_4     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_2        ! Serial number of the star
      integer*4     k_2         ! Number of determinations (plates)
      real*8        RAhour_2    ! (h) Right Ascension (J2000)
      real*4        e_RAhour_2  ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_5     ! (deg) Declination (J2000)
      real*4        e_DEdeg_2   ! (arcsec) Error in Declination
      real*4        Bmag_2      ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_2    ! (mag) Error in B magnitude
      real*8        Epoch_2     ! (yr) Epoch of observations
      real*4        FWHM_2      ! (1.46pix) Diameter of star image
      real*8        cI_2        ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_2       ! (h) Right ascension (J2000)
      integer*4     RAm_2       ! (min) Right ascension (J2000)
      real*4        RAs_2       ! (s) [0/60] Right ascension (J2000)
      character*1   DE__2       ! Declination sign (J2000)
      integer*4     DEd_2       ! (deg) Declination (J2000)
      integer*4     DEm_2       ! (arcmin) Declination (J2000)
      real*4        DEs_2       ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_2     ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_2      ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_2      ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_2     ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_2      ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_2      ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_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     ! (deg) Right Ascension J2000
      real*8        DEdeg_6     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_3        ! Serial number of the star
      integer*4     k_3         ! Number of determinations (plates)
      real*8        RAhour_3    ! (h) Right Ascension (J2000)
      real*4        e_RAhour_3  ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_7     ! (deg) Declination (J2000)
      real*4        e_DEdeg_3   ! (arcsec) Error in Declination
      real*4        Bmag_3      ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_3    ! (mag) Error in B magnitude
      real*8        Epoch_3     ! (yr) Epoch of observations
      real*4        FWHM_3      ! (1.46pix) Diameter of star image
      real*8        cI_3        ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_3       ! (h) Right ascension (J2000)
      integer*4     RAm_3       ! (min) Right ascension (J2000)
      real*4        RAs_3       ! (s) [0/60] Right ascension (J2000)
      character*1   DE__3       ! Declination sign (J2000)
      integer*4     DEd_3       ! (deg) Declination (J2000)
      integer*4     DEm_3       ! (arcmin) Declination (J2000)
      real*4        DEs_3       ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_3     ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_3      ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_3      ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_3     ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_3      ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_3      ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_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     ! (deg) Right Ascension J2000
      real*8        DEdeg_8     ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_4        ! Serial number of the star
      integer*4     k_4         ! Number of determinations (plates)
      real*8        RAhour_4    ! (h) Right Ascension (J2000)
      real*4        e_RAhour_4  ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_9     ! (deg) Declination (J2000)
      real*4        e_DEdeg_4   ! (arcsec) Error in Declination
      real*4        Bmag_4      ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_4    ! (mag) Error in B magnitude
      real*8        Epoch_4     ! (yr) Epoch of observations
      real*4        FWHM_4      ! (1.46pix) Diameter of star image
      real*8        cI_4        ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_4       ! (h) Right ascension (J2000)
      integer*4     RAm_4       ! (min) Right ascension (J2000)
      real*4        RAs_4       ! (s) [0/60] Right ascension (J2000)
      character*1   DE__4       ! Declination sign (J2000)
      integer*4     DEd_4       ! (deg) Declination (J2000)
      integer*4     DEm_4       ! (arcmin) Declination (J2000)
      real*4        DEs_4       ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_4     ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_4      ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_4      ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_4     ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_4      ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_4      ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_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     ! (deg) Right Ascension J2000
      real*8        DEdeg_10    ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ID_5        ! Serial number of the star
      integer*4     k_5         ! Number of determinations (plates)
      real*8        RAhour_5    ! (h) Right Ascension (J2000)
      real*4        e_RAhour_5  ! (arcsec) Error in right ascension, RA*cos(DE)
      real*8        DEdeg_11    ! (deg) Declination (J2000)
      real*4        e_DEdeg_5   ! (arcsec) Error in Declination
      real*4        Bmag_5      ! (mag) Johnson Blue magnitude
      real*4        e_Bmag_5    ! (mag) Error in B magnitude
      real*8        Epoch_5     ! (yr) Epoch of observations
      real*4        FWHM_5      ! (1.46pix) Diameter of star image
      real*8        cI_5        ! Intensity in the center of image
*                                     (in max256 unit)
      integer*4     RAh_5       ! (h) Right ascension (J2000)
      integer*4     RAm_5       ! (min) Right ascension (J2000)
      real*4        RAs_5       ! (s) [0/60] Right ascension (J2000)
      character*1   DE__5       ! Declination sign (J2000)
      integer*4     DEd_5       ! (deg) Declination (J2000)
      integer*4     DEm_5       ! (arcmin) Declination (J2000)
      real*4        DEs_5       ! (arcsec) [0/60] Declination (J2000)
      real*4        f_mag_5     ! (mag) ?=0 UCAC4 f magnitude (1)
      real*8        pmRA_5      ! (mas/yr) ?=0 Proper motion in RA*cos(DE) (1)
      real*8        pmDE_5      ! (mas/yr) ?=0 Proper motion in DE (1)
      real*4        B4mag_5     ! (mag) ?=0 UCAC4 B magnitude (1)
      real*4        Vmag_5      ! (mag) ?=0 UCAC4 V magnitude (1)
      real*4        rmag_5      ! (mag) ?=0 UCAC4 r magnitude (1)
      real*4        imag_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,k,RAhour,e_RAhour,DEdeg_1,e_DEdeg,Bmag,e_Bmag,Epoch,FWHM,
     +  cI,RAh,RAm,RAs,DE_,DEd,DEm,DEs,f_mag,pmRA,pmDE,B4mag,Vmag,
     +  rmag,imag
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAhour .GT. -180) RAdeg=RAhour*15.
        if(DEdeg_1 .GE. 0) DEdeg=DEdeg_1
        if(DEd .GE. 0) DEdeg=DEdeg+DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  ID,k,RAhour,e_RAhour,DEdeg_1,e_DEdeg,Bmag,e_Bmag,Epoch,FWHM,
     +  cI,RAh,RAm,RAs,DE_,DEd,DEm,DEs,f_mag,pmRA,pmDE,B4mag,Vmag,
     +  rmag,imag
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
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,k_1,RAhour_1,e_RAhour_1,DEdeg_3,e_DEdeg_1,Bmag_1,
     +  e_Bmag_1,Epoch_1,FWHM_1,cI_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,
     +  DEm_1,DEs_1,f_mag_1,pmRA_1,pmDE_1,B4mag_1,Vmag_1,rmag_1,imag_1
        RAdeg_1 = rNULL__
        DEdeg_2 = 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 .GT. -180) RAdeg_1=RAhour_1*15.
        if(DEdeg_3 .GE. 0) DEdeg_2=DEdeg_3
        if(DEd_1 .GE. 0) DEdeg_2=DEdeg_2+DEd_1
        if(DEm_1 .GE. 0) DEdeg_2=DEdeg_2+DEm_1/60.
        if(DEs_1 .GE. 0) DEdeg_2=DEdeg_2+DEs_1/3600.
        if(DE__1.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2
c    ..............Just test output...........
        write(6,2)
     +  ID_1,k_1,RAhour_1,e_RAhour_1,DEdeg_3,e_DEdeg_1,Bmag_1,
     +  e_Bmag_1,Epoch_1,FWHM_1,cI_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,
     +  DEm_1,DEs_1,f_mag_1,pmRA_1,pmDE_1,B4mag_1,Vmag_1,rmag_1,imag_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_2
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,k_2,RAhour_2,e_RAhour_2,DEdeg_5,e_DEdeg_2,Bmag_2,
     +  e_Bmag_2,Epoch_2,FWHM_2,cI_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,
     +  DEm_2,DEs_2,f_mag_2,pmRA_2,pmDE_2,B4mag_2,Vmag_2,rmag_2,imag_2
        RAdeg_2 = rNULL__
        DEdeg_4 = 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 .GT. -180) RAdeg_2=RAhour_2*15.
        if(DEdeg_5 .GE. 0) DEdeg_4=DEdeg_5
        if(DEd_2 .GE. 0) DEdeg_4=DEdeg_4+DEd_2
        if(DEm_2 .GE. 0) DEdeg_4=DEdeg_4+DEm_2/60.
        if(DEs_2 .GE. 0) DEdeg_4=DEdeg_4+DEs_2/3600.
        if(DE__2.EQ.'-'.AND.DEdeg_4.GE.0) DEdeg_4=-DEdeg_4
c    ..............Just test output...........
        write(6,3)
     +  ID_2,k_2,RAhour_2,e_RAhour_2,DEdeg_5,e_DEdeg_2,Bmag_2,
     +  e_Bmag_2,Epoch_2,FWHM_2,cI_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,
     +  DEm_2,DEs_2,f_mag_2,pmRA_2,pmDE_2,B4mag_2,Vmag_2,rmag_2,imag_2
        write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_4
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,k_3,RAhour_3,e_RAhour_3,DEdeg_7,e_DEdeg_3,Bmag_3,
     +  e_Bmag_3,Epoch_3,FWHM_3,cI_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,
     +  DEm_3,DEs_3,f_mag_3,pmRA_3,pmDE_3,B4mag_3,Vmag_3,rmag_3,imag_3
        RAdeg_3 = rNULL__
        DEdeg_6 = 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 .GT. -180) RAdeg_3=RAhour_3*15.
        if(DEdeg_7 .GE. 0) DEdeg_6=DEdeg_7
        if(DEd_3 .GE. 0) DEdeg_6=DEdeg_6+DEd_3
        if(DEm_3 .GE. 0) DEdeg_6=DEdeg_6+DEm_3/60.
        if(DEs_3 .GE. 0) DEdeg_6=DEdeg_6+DEs_3/3600.
        if(DE__3.EQ.'-'.AND.DEdeg_6.GE.0) DEdeg_6=-DEdeg_6
c    ..............Just test output...........
        write(6,4)
     +  ID_3,k_3,RAhour_3,e_RAhour_3,DEdeg_7,e_DEdeg_3,Bmag_3,
     +  e_Bmag_3,Epoch_3,FWHM_3,cI_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,
     +  DEm_3,DEs_3,f_mag_3,pmRA_3,pmDE_3,B4mag_3,Vmag_3,rmag_3,imag_3
        write(6,'(6H Pos: 2F8.4)') RAdeg_3,DEdeg_6
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,k_4,RAhour_4,e_RAhour_4,DEdeg_9,e_DEdeg_4,Bmag_4,
     +  e_Bmag_4,Epoch_4,FWHM_4,cI_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,
     +  DEm_4,DEs_4,f_mag_4,pmRA_4,pmDE_4,B4mag_4,Vmag_4,rmag_4,imag_4
        RAdeg_4 = rNULL__
        DEdeg_8 = 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 .GT. -180) RAdeg_4=RAhour_4*15.
        if(DEdeg_9 .GE. 0) DEdeg_8=DEdeg_9
        if(DEd_4 .GE. 0) DEdeg_8=DEdeg_8+DEd_4
        if(DEm_4 .GE. 0) DEdeg_8=DEdeg_8+DEm_4/60.
        if(DEs_4 .GE. 0) DEdeg_8=DEdeg_8+DEs_4/3600.
        if(DE__4.EQ.'-'.AND.DEdeg_8.GE.0) DEdeg_8=-DEdeg_8
c    ..............Just test output...........
        write(6,5)
     +  ID_4,k_4,RAhour_4,e_RAhour_4,DEdeg_9,e_DEdeg_4,Bmag_4,
     +  e_Bmag_4,Epoch_4,FWHM_4,cI_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,
     +  DEm_4,DEs_4,f_mag_4,pmRA_4,pmDE_4,B4mag_4,Vmag_4,rmag_4,imag_4
        write(6,'(6H Pos: 2F8.4)') RAdeg_4,DEdeg_8
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,k_5,RAhour_5,e_RAhour_5,DEdeg_11,e_DEdeg_5,Bmag_5,
     +  e_Bmag_5,Epoch_5,FWHM_5,cI_5,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,
     +  DEm_5,DEs_5,f_mag_5,pmRA_5,pmDE_5,B4mag_5,Vmag_5,rmag_5,imag_5
        RAdeg_5 = rNULL__
        DEdeg_10 = 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 .GT. -180) RAdeg_5=RAhour_5*15.
        if(DEdeg_11 .GE. 0) DEdeg_10=DEdeg_11
        if(DEd_5 .GE. 0) DEdeg_10=DEdeg_10+DEd_5
        if(DEm_5 .GE. 0) DEdeg_10=DEdeg_10+DEm_5/60.
        if(DEs_5 .GE. 0) DEdeg_10=DEdeg_10+DEs_5/3600.
        if(DE__5.EQ.'-'.AND.DEdeg_10.GE.0) DEdeg_10=-DEdeg_10
c    ..............Just test output...........
        write(6,6)
     +  ID_5,k_5,RAhour_5,e_RAhour_5,DEdeg_11,e_DEdeg_5,Bmag_5,
     +  e_Bmag_5,Epoch_5,FWHM_5,cI_5,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,
     +  DEm_5,DEs_5,f_mag_5,pmRA_5,pmDE_5,B4mag_5,Vmag_5,rmag_5,imag_5
        write(6,'(6H Pos: 2F8.4)') RAdeg_5,DEdeg_10
c    .......End.of.Just test output...........
      end do
      close(1)

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