FORTRAN Generation
(/./ftp/cats/J/A_A/530/A93)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/530/A93 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-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. J/A+A/530/A93     de Houtman, Kepler and Halley star catalogs   (Verbunt+ 2011)
*================================================================================
*Early star catalogues of the southern sky.
*    Verbunt, F., van Gent, R.H.
*   <Astron. Astrophys. 530, A93 (2011)>
*   =2011A&A...530A..93V
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'houtman.dat'	! Star Catalogue of de Houtman (1603)

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

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

      integer*4     F          (nr__) ! [1/304] Sequence number in de Houtman (1603)
*                                 Catalogue
      integer*4     C          (nr__) ! [1/21] Sequence number of constellation
      character*3   cst        (nr__) ! Abbreviation of constellation name
      integer*4     N          (nr__) ! [1/56] Sequence number of star in constellation
      integer*4     RAd        (nr__) ! (deg) [0/360] Degrees of right ascension (1)
      integer*4     RAm        (nr__) ! (arcmin) Minutes of right ascension (1)
      character*1   DE_        (nr__) ! Sign of declination (1)
      integer*4     DEd        (nr__) ! (deg) Degrees of declination (1)
      integer*4     DEm        (nr__) ! (arcmin) Arcminutes of declination (1)
      integer*4     Vf         (nr__) ! magnitude as given by de Houtman (9=nebulous)
      integer*4     HIP        (nr__) ! Hipparcos number of identification
      integer*4     I          (nr__) ! [1,6] Quality of identification (G3)
      integer*4     v_K        (nr__) ! flag indicating comparison with Knobel (2)
      real*4        Vmag       (nr__) ! Magnitude of Hipparcos identification
      real*4        D_RA       (nr__) ! (arcmin) Difference in RA: Hipparcos-Houtman
      real*4        D_DE       (nr__) ! (arcmin) Difference in DE: Hipparcos-Houtman
      real*4        Del        (nr__) ! (arcmin) Offset between de Houtman star and
*                                   identification
*Note (1): The right ascension is given by RA=RAd+RAm/60;
*     the declination, always negative, is given by DE=-(DEd+DEm/60).
*     The equinox is assumed to be AD 1601.0
*Note (2): Flag for comparison with Knobel (1917MNRAS..77..580K and
*     1917MNRAS..77..414K) as follows:
*      1 = same identification in Knobel
*      2 = identification in Knobel to alternative identification
*      3 = identification in Knobel different from our identification

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

C  Declarations for 'classis.dat'	! Secunda Classis and Tertia Classis Kepler (1627)

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

C  1601 position composed of: Elon.Z Elon.d Elon.m Elat.- Elat.d Elat.m
      real*8        RAdeg_1    (nr__1) ! (deg) Right Ascension 1601
      real*8        DEdeg_1    (nr__1) ! (deg)     Declination 1601
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     K          (nr__1) ! [1005/1438] Sequence number in 
*                                             Classis Catalogue (1)
      integer*4     C_1        (nr__1) ! [1/62] Sequence number of constellation (2)
      character*3   cst_1      (nr__1) ! Abbreviation of constellation name
      integer*4     N_1        (nr__1) ! [1/63] Sequence number of star in 
*                                        constellation (3)
      integer*4     Elon_Z     (nr__1) ! (30deg) [1,12] Zodiakal sign of ecliptic longitude,
*                                        Equinox AD1601 (G1)
      integer*4     Elon_d     (nr__1) ! (deg) [0,30] Degrees of ecliptic longitude (G1)
      real*4        Elon_m     (nr__1) ! (arcmin) Minutes of ecliptic longitude (G1)
      integer*4     Elat_d     (nr__1) ! (deg) Degrees of ecliptic latitude  (G2)
      integer*4     Elat_m     (nr__1) ! (arcmin) Arcminutes of ecliptic latitude (G2)
      character*1   Elat__     (nr__1) ! [AB] Sign of ecliptic latitude, Australis (-)
*                                      or Borealis (+), Equinox AD1601 (G2)
      integer*4     VK         (nr__1) ! [1,9] Magnitude as given by Kepler (4)
      integer*4     HIP_1      (nr__1) ! Hipparcos number of identification
      integer*4     I_1        (nr__1) ! [1,6] Quality of identification (G3)
      real*4        Vmag_1     (nr__1) ! Magnitude of Hipparcos identification
      real*8        Dlon       (nr__1) ! (arcmin) Difference in longitude: Hipparcos-Kepler
      real*4        Dlat       (nr__1) ! (arcmin) Difference in latitude: Hipparcos-Kepler
      real*4        Del_1      (nr__1) ! (arcmin) Offset between Kepler star and identification
      integer*4     P          (nr__1) ! [0/1028] Corresponding sequence number in 
*                                          Ptolemaios (5)
      integer*4     v_P        (nr__1) ! Flag indicating comparison with Ptolemaios (5)
*Note (1):
*  * K < 1005 refers to the Kepler (1627) edition of Brahe's catalogue
*  * 1004 < K < 1303 refers to Secunda Classis, the entries in which
*    we number starting with 1005 (ending with 1302)
*  * K > 1302 refers to Tertia Classis
*Note (2):
*  * C < 47 constellation in Kepler (1627) edition of Brahe catalogue,
*    to which star is added in Secunda Classis
*  * 46 < C < 51 constellation from catalogue of Ptolemaios, not present
*    in Brahe Catalogue, for which Secunda Classis lists stars
*  * C > 50  new constellations from Plancius, in Tertia Classis
*    It sometimes refers to the Secunda Classis. For detail see paper.
*Note (3): Numbering continued for stars in constellation listed in Brahe;
*     starting with 1 for stars in constellations not listed in Brahe
*Note (4): 8 = no magnitude given in catalogue;  9 = nebulous
*Note (5): Comparison with our identification for corresponding star in
*     Toomer's 1998 edition of the star catalogue of Ptolemaios (Princeton
*     Univ. Press " Ptolemy's Almagest", 2nd ed.) as follows:
*      1 = same identification as for star in Classis
*      2 = identification to alternative identification
*      3 = identification different (including identification for new stars)

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

C  Declarations for 'aliter.dat'	! Variants on Secunda Classis (Kepler (1627)

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

C  1601 position composed of: Elon.Z Elon.d Elon.m Elat.- Elat.d Elat.m
      real*8        RAdeg_2    (nr__2) ! (deg) Right Ascension 1601
      real*8        DEdeg_2    (nr__2) ! (deg)     Declination 1601
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     K_1        (nr__2) ! [1005/1438] Sequence number in 
*                                             Classis Catalogue (1)
      integer*4     C_2        (nr__2) ! [1/62] Sequence number of constellation (2)
      character*3   cst_2      (nr__2) ! Abbreviation of constellation name
      integer*4     N_2        (nr__2) ! [1/63] Sequence number of star in 
*                                        constellation (3)
      integer*4     Elon_Z_1   (nr__2) ! (30deg) [1,12] Zodiakal sign of ecliptic longitude,
*                                        Equinox AD1601 (G1)
      integer*4     Elon_d_1   (nr__2) ! (deg) [0,30] Degrees of ecliptic longitude (G1)
      real*4        Elon_m_1   (nr__2) ! (arcmin) Minutes of ecliptic longitude (G1)
      integer*4     Elat_d_1   (nr__2) ! (deg) Degrees of ecliptic latitude  (G2)
      integer*4     Elat_m_1   (nr__2) ! (arcmin) Arcminutes of ecliptic latitude (G2)
      character*1   Elat___1   (nr__2) ! [AB] Sign of ecliptic latitude, Australis (-)
*                                      or Borealis (+), Equinox AD1601 (G2)
      integer*4     VK_1       (nr__2) ! [1,9] Magnitude as given by Kepler (4)
      integer*4     HIP_2      (nr__2) ! Hipparcos number of identification
      integer*4     I_2        (nr__2) ! [1,6] Quality of identification (G3)
      real*4        Vmag_2     (nr__2) ! Magnitude of Hipparcos identification
      real*8        Dlon_1     (nr__2) ! (arcmin) Difference in longitude: Hipparcos-Kepler
      real*4        Dlat_1     (nr__2) ! (arcmin) Difference in latitude: Hipparcos-Kepler
      real*4        Del_2      (nr__2) ! (arcmin) Offset between Kepler star and identification
      integer*4     P_1        (nr__2) ! [0/1028] Corresponding sequence number in 
*                                          Ptolemaios (5)
      integer*4     v_P_1      (nr__2) ! Flag indicating comparison with Ptolemaios (5)
*Note (1):
*  * K < 1005 refers to the Kepler (1627) edition of Brahe's catalogue
*  * 1004 < K < 1303 refers to Secunda Classis, the entries in which
*    we number starting with 1005 (ending with 1302)
*  * K > 1302 refers to Tertia Classis
*Note (2):
*  * C < 47 constellation in Kepler (1627) edition of Brahe catalogue,
*    to which star is added in Secunda Classis
*  * 46 < C < 51 constellation from catalogue of Ptolemaios, not present
*    in Brahe Catalogue, for which Secunda Classis lists stars
*  * C > 50  new constellations from Plancius, in Tertia Classis
*    It sometimes refers to the Secunda Classis. For detail see paper.
*Note (3): Numbering continued for stars in constellation listed in Brahe;
*     starting with 1 for stars in constellations not listed in Brahe
*Note (4): 8 = no magnitude given in catalogue;  9 = nebulous
*Note (5): Comparison with our identification for corresponding star in
*     Toomer's 1998 edition of the star catalogue of Ptolemaios (Princeton
*     Univ. Press " Ptolemy's Almagest", 2nd ed.) as follows:
*      1 = same identification as for star in Classis
*      2 = identification to alternative identification
*      3 = identification different (including identification for new stars)

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

C  Declarations for 'halley.dat'	! Star catalogue of Halley (1679)

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

C  1678 position composed of: Elon.Z Elon.d Elon.m Elat.- Elat.d Elat.m
      real*8        RAdeg_3    (nr__3) ! (deg) Right Ascension 1678
      real*8        DEdeg_3    (nr__3) ! (deg)     Declination 1678
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     E          (nr__3) ! Sequence number in Halley catalogue
      integer*4     C_3        (nr__3) ! Sequence number of constellation
      character*3   cst_3      (nr__3) ! Abbreviation of constellation name
      integer*4     N_3        (nr__3) ! [1/46] Sequence number of star in constellation
      integer*4     Elon_Z_2   (nr__3) ! (30deg) [1,12] Zodiakal sign of ecliptic longitude,
*                                      Equinox AD 1678 (G1)
      integer*4     Elon_d_2   (nr__3) ! (deg) [0,30] Degrees of ecliptic longitude (G1)
      real*4        Elon_m_2   (nr__3) ! (arcmin) Minutes of ecliptic longitude (G1)
      character*1   Elat___2   (nr__3) ! [-] Sign of ecliptic latitude (G2)
      integer*4     Elat_d_2   (nr__3) ! (deg) Degrees of ecliptic latitude, Equinox 1678  (G2)
      real*4        Elat_m_2   (nr__3) ! (arcmin) Arcminutes of ecliptic latitude (G2)
      real*4        VE         (nr__3) ! Magnitude as given by Halley (9=nebulous)
      integer*4     HIP1       (nr__3) ! ? Hipparcos number of first standard star (2)
      integer*4     d1_d       (nr__3) ! (deg) ? Degrees of angle to first standard star (2)
      integer*4     d1_m       (nr__3) ! (arcmin) ? Arcminutes of angle to first standard star (2)
      integer*4     d1_s       (nr__3) ! (arcsec) ? Arcseconds of angle to first standard star (2)
      integer*4     HIP2       (nr__3) ! ? Hipparcos number of second standard star (2)
      integer*4     d2_d       (nr__3) ! (deg) ? Degrees of angle to second standard star (2)
      integer*4     d2_m       (nr__3) ! (arcmin) ? Arcminutes of angle to 2nd standard star (2)
      integer*4     d2_s       (nr__3) ! (arcsec) ? Arcseconds of angle to 2nd standard star (2)
      integer*4     n_HIP2     (nr__3) ! [-3,-1]? only one standard star used (2)
      integer*4     HIP_3      (nr__3) ! ?=0 Hipparcos number of identification
      integer*4     I_3        (nr__3) ! [1,6] Quality of identification (G3)
      real*4        Vmag_3     (nr__3) ! ? Magnitude of Hipparcos identification
      real*4        Dlon_2     (nr__3) ! (arcmin) ? Difference in longitude: Hipparcos-Kepler
      real*4        Dlat_2     (nr__3) ! (arcmin) ? Difference in latitude: Hipparcos-Kepler
      real*4        Del_3      (nr__3) ! (arcmin) ? Offset between Kepler star and identification
      integer*4     HIP3       (nr__3) ! ? Hipparcos number of third standard star (3)
      integer*4     d3_d       (nr__3) ! (deg) ? Degrees of angle to third standard star (3)
      integer*4     d3_m       (nr__3) ! (arcmin) ? Arcminutes of angle to third standard star (3)
      integer*4     d3_s       (nr__3) ! (arcsec) ? Arcseconds of angle to third standard star (3)
      integer*4     HIP4       (nr__3) ! ? Hipparcos number of fourth standard star (3)
      integer*4     d4_d       (nr__3) ! (deg) ? Degrees of angle to fourth standard star (3)
      integer*4     d4_m       (nr__3) ! (arcmin) ? Arcminutes of angle to fourth standard star (3)
      integer*4     d4_s       (nr__3) ! (arcsec) ? Arcseconds of angle to fourth standard star (3)
*Note (2): Halley give measured angles to two standard stars; for some stars
*     more than two standard stars are given, sometimes only one standard
*     star is used, and combined with information from Brahe or Kepler. We
*     indicates these with following n_HIP2 values:
*     -1 = longitude from Brahe
*     -2 = latitude from Brahe
*     -3 = longitude from Kepler
*     These values are announced in the paper as "negative value for HIP"
*Note (3): These parameters give extra distances to standard stars
*    for four sequence numbers, announced in the paper as
*    "at the end of the table, four sequence numbers are repeated".

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

C  Loading file 'houtman.dat'	! Star Catalogue of de Houtman (1603)

C  Format for file interpretation

    1 format(
     +  I3,1X,I2,2X,A3,1X,I2,2X,I3,1X,I2,1X,A1,I2,1X,I2,2X,I1,1X,I6,
     +  2X,I1,1X,I1,1X,F4.1,1X,F6.1,1X,F6.1,1X,F6.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'houtman.dat')
      write(6,*) '....Loading file: houtman.dat'
      do i__=1,304
        read(1,'(A70)')ar__
        read(ar__,1)
     +  F(i__),C(i__),cst(i__),N(i__),RAd(i__),RAm(i__),DE_(i__),
     +  DEd(i__),DEm(i__),Vf(i__),HIP(i__),I(i__),v_K(i__),Vmag(i__),
     +  D_RA(i__),D_DE(i__),Del(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(RAd(i__) .GT. -180) RAdeg(i__)=RAd(i__)
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/60.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  F(i__),C(i__),cst(i__),N(i__),RAd(i__),RAm(i__),DE_(i__),
     +  DEd(i__),DEm(i__),Vf(i__),HIP(i__),I(i__),v_K(i__),Vmag(i__),
     +  D_RA(i__),D_DE(i__),Del(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 'classis.dat'	! Secunda Classis and Tertia Classis Kepler (1627)

C  Format for file interpretation

    2 format(
     +  I4,2X,I2,2X,A3,1X,I2,1X,I2,1X,I2,1X,F4.1,2X,I2,1X,I2,1X,A1,1X,
     +  I1,1X,I6,1X,I1,1X,F4.1,F7.1,1X,F6.1,1X,F6.1,3X,I4,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'classis.dat')
      write(6,*) '....Loading file: classis.dat'
      do i__=1,434
        read(1,'(A82)')ar__1
        read(ar__1,2)
     +  K(i__),C_1(i__),cst_1(i__),N_1(i__),Elon_Z(i__),Elon_d(i__),
     +  Elon_m(i__),Elat_d(i__),Elat_m(i__),Elat__(i__),VK(i__),
     +  HIP_1(i__),I_1(i__),Vmag_1(i__),Dlon(i__),Dlat(i__),
     +  Del_1(i__),P(i__),v_P(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)
c    ..............Just test output...........
        write(6,2)
     +  K(i__),C_1(i__),cst_1(i__),N_1(i__),Elon_Z(i__),Elon_d(i__),
     +  Elon_m(i__),Elat_d(i__),Elat_m(i__),Elat__(i__),VK(i__),
     +  HIP_1(i__),I_1(i__),Vmag_1(i__),Dlon(i__),Dlat(i__),
     +  Del_1(i__),P(i__),v_P(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 'aliter.dat'	! Variants on Secunda Classis (Kepler (1627)

C  Format for file interpretation

    3 format(
     +  I4,2X,I2,2X,A3,1X,I2,1X,I2,1X,I2,1X,F4.1,2X,I2,1X,I2,1X,A1,1X,
     +  I1,1X,I6,1X,I1,1X,F4.1,F7.1,1X,F6.1,1X,F6.1,3X,I4,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'aliter.dat')
      write(6,*) '....Loading file: aliter.dat'
      do i__=1,20
        read(1,'(A82)')ar__2
        read(ar__2,3)
     +  K_1(i__),C_2(i__),cst_2(i__),N_2(i__),Elon_Z_1(i__),
     +  Elon_d_1(i__),Elon_m_1(i__),Elat_d_1(i__),Elat_m_1(i__),
     +  Elat___1(i__),VK_1(i__),HIP_2(i__),I_2(i__),Vmag_2(i__),
     +  Dlon_1(i__),Dlat_1(i__),Del_2(i__),P_1(i__),v_P_1(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)
c    ..............Just test output...........
        write(6,3)
     +  K_1(i__),C_2(i__),cst_2(i__),N_2(i__),Elon_Z_1(i__),
     +  Elon_d_1(i__),Elon_m_1(i__),Elat_d_1(i__),Elat_m_1(i__),
     +  Elat___1(i__),VK_1(i__),HIP_2(i__),I_2(i__),Vmag_2(i__),
     +  Dlon_1(i__),Dlat_1(i__),Del_2(i__),P_1(i__),v_P_1(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 'halley.dat'	! Star catalogue of Halley (1679)

C  Format for file interpretation

    4 format(
     +  I3,1X,I2,2X,A3,1X,I2,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,F4.1,1X,
     +  F3.1,1X,I6,1X,I2,1X,I2,1X,I2,1X,I6,1X,I2,1X,I2,1X,I2,1X,I2,1X,
     +  I6,1X,I1,1X,F4.1,1X,F6.1,1X,F6.1,1X,F6.1,1X,I6,1X,I2,1X,I2,1X,
     +  I2,1X,I5,1X,I2,1X,I2,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'halley.dat')
      write(6,*) '....Loading file: halley.dat'
      do i__=1,341
        read(1,'(A139)')ar__3
        read(ar__3,4)
     +  E(i__),C_3(i__),cst_3(i__),N_3(i__),Elon_Z_2(i__),
     +  Elon_d_2(i__),Elon_m_2(i__),Elat___2(i__),Elat_d_2(i__),
     +  Elat_m_2(i__),VE(i__),HIP1(i__),d1_d(i__),d1_m(i__),d1_s(i__),
     +  HIP2(i__),d2_d(i__),d2_m(i__),d2_s(i__),n_HIP2(i__),
     +  HIP_3(i__),I_3(i__),Vmag_3(i__),Dlon_2(i__),Dlat_2(i__),
     +  Del_3(i__),HIP3(i__),d3_d(i__),d3_m(i__),d3_s(i__),HIP4(i__),
     +  d4_d(i__),d4_m(i__),d4_s(i__)
        if(ar__3(40:45) .EQ. '') HIP1(i__) = iNULL__
        if(ar__3(47:48) .EQ. '') d1_d(i__) = iNULL__
        if(ar__3(50:51) .EQ. '') d1_m(i__) = iNULL__
        if(ar__3(53:54) .EQ. '') d1_s(i__) = iNULL__
        if(ar__3(56:61) .EQ. '') HIP2(i__) = iNULL__
        if(ar__3(63:64) .EQ. '') d2_d(i__) = iNULL__
        if(ar__3(66:67) .EQ. '') d2_m(i__) = iNULL__
        if(ar__3(69:70) .EQ. '') d2_s(i__) = iNULL__
        if(ar__3(72:73) .EQ. '') n_HIP2(i__) = iNULL__
        if(ar__3(84:87) .EQ. '') Vmag_3(i__) = rNULL__
        if(ar__3(89:94) .EQ. '') Dlon_2(i__) = rNULL__
        if(ar__3(96:101) .EQ. '') Dlat_2(i__) = rNULL__
        if(ar__3(103:108) .EQ. '') Del_3(i__) = rNULL__
        if(ar__3(110:115) .EQ. '') HIP3(i__) = iNULL__
        if(ar__3(117:118) .EQ. '') d3_d(i__) = iNULL__
        if(ar__3(120:121) .EQ. '') d3_m(i__) = iNULL__
        if(ar__3(123:124) .EQ. '') d3_s(i__) = iNULL__
        if(ar__3(126:130) .EQ. '') HIP4(i__) = iNULL__
        if(ar__3(132:133) .EQ. '') d4_d(i__) = iNULL__
        if(ar__3(135:136) .EQ. '') d4_m(i__) = iNULL__
        if(ar__3(138:139) .EQ. '') d4_s(i__) = iNULL__
        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)
c    ..............Just test output...........
        write(6,4)
     +  E(i__),C_3(i__),cst_3(i__),N_3(i__),Elon_Z_2(i__),
     +  Elon_d_2(i__),Elon_m_2(i__),Elat___2(i__),Elat_d_2(i__),
     +  Elat_m_2(i__),VE(i__),HIP1(i__),d1_d(i__),d1_m(i__),d1_s(i__),
     +  HIP2(i__),d2_d(i__),d2_m(i__),d2_s(i__),n_HIP2(i__),
     +  HIP_3(i__),I_3(i__),Vmag_3(i__),Dlon_2(i__),Dlat_2(i__),
     +  Del_3(i__),HIP3(i__),d3_d(i__),d3_m(i__),d3_s(i__),HIP4(i__),
     +  d4_d(i__),d4_m(i__),d4_s(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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