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