FORTRAN Generation
(/./ftp/cats/V/61)

Conversion of standardized ReadMe file for file /./ftp/cats/V/61 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. V/61                Almagest  (Ptolemy's Star Catalog)
*================================================================================
*The Almagest: Ptolemy's star catalogue
*     Ptolemy C.: Almagest (years 127-141)
*     Manitius K.: 1913
*    <ed. B.G. Teubner, Leipzig (1913)>
*     Jaschek C.: 1987
*    <Bull. Inform. CDS 33, 145 (1987)>
*    =1987BICDS..33..125J
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'north'	! Northern part (except Zodiacal stars)

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

      character*16  Const      (nr__) ! Constellation name
      integer*4     HR         (nr__) ! [0/9110]? Bright Star Number, 0 if cluster
      character*3   m_HR       (nr__) ! Second HR number for binaries
      integer*4     Seq        (nr__) ! Number in constellation
      integer*4     Lon_d      (nr__) ! (deg) Ptolemean (ecliptic) longitude (degrees) (1)
      integer*4     Lon_m      (nr__) ! (arcmin) Ptolemean (ecliptic) longitude (minutes) (1)
      character*1   Lat_       (nr__) ! Ptolemean (ecliptic) latitude (sign)
      integer*4     Lat_d      (nr__) ! (deg) Ptolemean (ecliptic) latitude (degrees)
      integer*4     Lat_m      (nr__) ! (arcmin) Ptolemean (ecliptic) latitude (minutes)
      real*4        Mag        (nr__) ! (mag) [1/8] Ptolemean magnitude (2)
      character*1   Disag      (nr__) ! [D] 'D' if a disagreement exists
*                                    concerning identification
*Note (1): the Ptolomean longitudes were converted by adding 30deg in Taurus
*    etc (see Description).
*    The epoch of observation is compared to -128 by Grasshoff G., 1990,
*    "The history of Ptolemy's Star Catalogue", ed. G.J. Toomer,
*    Springer-Verlag (ISBN 0-387-97181-5)
*Note (2): the Ptolomean magnitudes are only integer numbers between 1 and 6;
*    the values "7." correspond to "faint", and "8." for "nebula".
*    The decimal value "m.3" (m between 1 and 4) indicate original "(m)-(m+1)"
*    The decimal value "m.7" (m between 1 and 4) indicate original "(m-1)-(m)"
*    Examples:   "2.7" correspond to original  "3-2"
*                "3.3" correspond to original  "3-4"

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

C  Declarations for 'south'	! Southern part (except Zodiacal stars)

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

      character*16  Const_1    (nr__1) ! Constellation name
      integer*4     HR_1       (nr__1) ! [0/9110]? Bright Star Number, 0 if cluster
      character*3   m_HR_1     (nr__1) ! Second HR number for binaries
      integer*4     Seq_1      (nr__1) ! Number in constellation
      integer*4     Lon_d_1    (nr__1) ! (deg) Ptolemean (ecliptic) longitude (degrees) (1)
      integer*4     Lon_m_1    (nr__1) ! (arcmin) Ptolemean (ecliptic) longitude (minutes) (1)
      character*1   Lat__1     (nr__1) ! Ptolemean (ecliptic) latitude (sign)
      integer*4     Lat_d_1    (nr__1) ! (deg) Ptolemean (ecliptic) latitude (degrees)
      integer*4     Lat_m_1    (nr__1) ! (arcmin) Ptolemean (ecliptic) latitude (minutes)
      real*4        Mag_1      (nr__1) ! (mag) [1/8] Ptolemean magnitude (2)
      character*1   Disag_1    (nr__1) ! [D] 'D' if a disagreement exists
*                                    concerning identification
*Note (1): the Ptolomean longitudes were converted by adding 30deg in Taurus
*    etc (see Description).
*    The epoch of observation is compared to -128 by Grasshoff G., 1990,
*    "The history of Ptolemy's Star Catalogue", ed. G.J. Toomer,
*    Springer-Verlag (ISBN 0-387-97181-5)
*Note (2): the Ptolomean magnitudes are only integer numbers between 1 and 6;
*    the values "7." correspond to "faint", and "8." for "nebula".
*    The decimal value "m.3" (m between 1 and 4) indicate original "(m)-(m+1)"
*    The decimal value "m.7" (m between 1 and 4) indicate original "(m-1)-(m)"
*    Examples:   "2.7" correspond to original  "3-2"
*                "3.3" correspond to original  "3-4"

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

C  Declarations for 'zod-n'	! Zodiacal stars, North

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

      character*16  Const_2    (nr__2) ! Constellation name
      integer*4     HR_2       (nr__2) ! [0/9110]? Bright Star Number, 0 if cluster
      character*3   m_HR_2     (nr__2) ! Second HR number for binaries
      integer*4     Seq_2      (nr__2) ! Number in constellation
      integer*4     Lon_d_2    (nr__2) ! (deg) Ptolemean (ecliptic) longitude (degrees) (1)
      integer*4     Lon_m_2    (nr__2) ! (arcmin) Ptolemean (ecliptic) longitude (minutes) (1)
      character*1   Lat__2     (nr__2) ! Ptolemean (ecliptic) latitude (sign)
      integer*4     Lat_d_2    (nr__2) ! (deg) Ptolemean (ecliptic) latitude (degrees)
      integer*4     Lat_m_2    (nr__2) ! (arcmin) Ptolemean (ecliptic) latitude (minutes)
      real*4        Mag_2      (nr__2) ! (mag) [1/8] Ptolemean magnitude (2)
      character*1   Disag_2    (nr__2) ! [D] 'D' if a disagreement exists
*                                    concerning identification
*Note (1): the Ptolomean longitudes were converted by adding 30deg in Taurus
*    etc (see Description).
*    The epoch of observation is compared to -128 by Grasshoff G., 1990,
*    "The history of Ptolemy's Star Catalogue", ed. G.J. Toomer,
*    Springer-Verlag (ISBN 0-387-97181-5)
*Note (2): the Ptolomean magnitudes are only integer numbers between 1 and 6;
*    the values "7." correspond to "faint", and "8." for "nebula".
*    The decimal value "m.3" (m between 1 and 4) indicate original "(m)-(m+1)"
*    The decimal value "m.7" (m between 1 and 4) indicate original "(m-1)-(m)"
*    Examples:   "2.7" correspond to original  "3-2"
*                "3.3" correspond to original  "3-4"

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

C  Declarations for 'zod-s'	! Zodiacal stars, South

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

      character*16  Const_3    (nr__3) ! Constellation name
      integer*4     HR_3       (nr__3) ! [0/9110]? Bright Star Number, 0 if cluster
      character*3   m_HR_3     (nr__3) ! Second HR number for binaries
      integer*4     Seq_3      (nr__3) ! Number in constellation
      integer*4     Lon_d_3    (nr__3) ! (deg) Ptolemean (ecliptic) longitude (degrees) (1)
      integer*4     Lon_m_3    (nr__3) ! (arcmin) Ptolemean (ecliptic) longitude (minutes) (1)
      character*1   Lat__3     (nr__3) ! Ptolemean (ecliptic) latitude (sign)
      integer*4     Lat_d_3    (nr__3) ! (deg) Ptolemean (ecliptic) latitude (degrees)
      integer*4     Lat_m_3    (nr__3) ! (arcmin) Ptolemean (ecliptic) latitude (minutes)
      real*4        Mag_3      (nr__3) ! (mag) [1/8] Ptolemean magnitude (2)
      character*1   Disag_3    (nr__3) ! [D] 'D' if a disagreement exists
*                                    concerning identification
*Note (1): the Ptolomean longitudes were converted by adding 30deg in Taurus
*    etc (see Description).
*    The epoch of observation is compared to -128 by Grasshoff G., 1990,
*    "The history of Ptolemy's Star Catalogue", ed. G.J. Toomer,
*    Springer-Verlag (ISBN 0-387-97181-5)
*Note (2): the Ptolomean magnitudes are only integer numbers between 1 and 6;
*    the values "7." correspond to "faint", and "8." for "nebula".
*    The decimal value "m.3" (m between 1 and 4) indicate original "(m)-(m+1)"
*    The decimal value "m.7" (m between 1 and 4) indicate original "(m-1)-(m)"
*    Examples:   "2.7" correspond to original  "3-2"
*                "3.3" correspond to original  "3-4"

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

C  Loading file 'north'	! Northern part (except Zodiacal stars)

C  Format for file interpretation

    1 format(
     +  A16,1X,I4,A3,1X,I2,1X,I3,1X,I2,1X,A1,I2,1X,I2,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'north')
      write(6,*) '....Loading file: north'
      do i__=1,360
        read(1,'(A47)')ar__
        read(ar__,1)
     +  Const(i__),HR(i__),m_HR(i__),Seq(i__),Lon_d(i__),Lon_m(i__),
     +  Lat_(i__),Lat_d(i__),Lat_m(i__),Mag(i__),Disag(i__)
        if(ar__(18:21) .EQ. '') HR(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Const(i__),HR(i__),m_HR(i__),Seq(i__),Lon_d(i__),Lon_m(i__),
     +  Lat_(i__),Lat_d(i__),Lat_m(i__),Mag(i__),Disag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'south'	! Southern part (except Zodiacal stars)

C  Format for file interpretation

    2 format(
     +  A16,1X,I4,A3,1X,I2,1X,I3,1X,I2,1X,A1,I2,1X,I2,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'south')
      write(6,*) '....Loading file: south'
      do i__=1,317
        read(1,'(A47)')ar__1
        read(ar__1,2)
     +  Const_1(i__),HR_1(i__),m_HR_1(i__),Seq_1(i__),Lon_d_1(i__),
     +  Lon_m_1(i__),Lat__1(i__),Lat_d_1(i__),Lat_m_1(i__),Mag_1(i__),
     +  Disag_1(i__)
        if(ar__1(18:21) .EQ. '') HR_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  Const_1(i__),HR_1(i__),m_HR_1(i__),Seq_1(i__),Lon_d_1(i__),
     +  Lon_m_1(i__),Lat__1(i__),Lat_d_1(i__),Lat_m_1(i__),Mag_1(i__),
     +  Disag_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'zod-n'	! Zodiacal stars, North

C  Format for file interpretation

    3 format(
     +  A16,1X,I4,A3,1X,I2,1X,I3,1X,I2,1X,A1,I2,1X,I2,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'zod-n')
      write(6,*) '....Loading file: zod-n'
      do i__=1,167
        read(1,'(A47)')ar__2
        read(ar__2,3)
     +  Const_2(i__),HR_2(i__),m_HR_2(i__),Seq_2(i__),Lon_d_2(i__),
     +  Lon_m_2(i__),Lat__2(i__),Lat_d_2(i__),Lat_m_2(i__),Mag_2(i__),
     +  Disag_2(i__)
        if(ar__2(18:21) .EQ. '') HR_2(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  Const_2(i__),HR_2(i__),m_HR_2(i__),Seq_2(i__),Lon_d_2(i__),
     +  Lon_m_2(i__),Lat__2(i__),Lat_d_2(i__),Lat_m_2(i__),Mag_2(i__),
     +  Disag_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'zod-s'	! Zodiacal stars, South

C  Format for file interpretation

    4 format(
     +  A16,1X,I4,A3,1X,I2,1X,I3,1X,I2,1X,A1,I2,1X,I2,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'zod-s')
      write(6,*) '....Loading file: zod-s'
      do i__=1,183
        read(1,'(A47)')ar__3
        read(ar__3,4)
     +  Const_3(i__),HR_3(i__),m_HR_3(i__),Seq_3(i__),Lon_d_3(i__),
     +  Lon_m_3(i__),Lat__3(i__),Lat_d_3(i__),Lat_m_3(i__),Mag_3(i__),
     +  Disag_3(i__)
        if(ar__3(18:21) .EQ. '') HR_3(i__) = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  Const_3(i__),HR_3(i__),m_HR_3(i__),Seq_3(i__),Lon_d_3(i__),
     +  Lon_m_3(i__),Lat__3(i__),Lat_d_3(i__),Lat_m_3(i__),Mag_3(i__),
     +  Disag_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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