FORTRAN Generation
(/./ftp/cats/III/168)

Conversion of standardized ReadMe file for file /./ftp/cats/III/168 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-Apr-17
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. III/168     General Catalog of S Stars, second edition (Stephenson 1984)
*================================================================================
*General Catalog of S Stars, second edition
*    Stephenson, C.B.
*   <Publ. Warner & Swasey Observatory, 3, no. 1, 1 (1984)>
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'catalog.dat'	! The S Star catalog

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

C  1900 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension 1900
      real*8        DEdeg       ! (deg)     Declination 1900
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     CSS         ! *[1/1347]+ Number
      integer*4     RAh         ! (h) *Right ascension hour (1900)
      integer*4     RAm         ! (min) *Right ascension minute (1900)
      real*4        RAs         ! (s) *Right ascension second (1900)
      character*1   DE_         ! *Sign of declination
      integer*4     DEd         ! (deg) *Declination degree (1900)
      integer*4     DEm         ! (arcmin) *Declination minute (1900)
      integer*4     DEs         ! (arcsec) *Declination second (1900)
      real*4        bmag        ! (mag) *?=0.0 b magnitude
      real*4        vmag        ! (mag) *?=0.0 v magnitude
      real*4        imag        ! (mag) *?=0.0 i magnitude
      real*4        GLON        ! (deg) *Galactic longitude
      real*4        GLAT        ! (deg) *Galactic latitude
      character*5   bmag_pub    ! (mag) *Published b magnitude
      character*5   vmag_pub    ! (mag) *Publsihed v magnitude
      character*5   imag_pub    ! (mag) *Published i magnitude
      character*16  SpType      ! *Spectral type
      character*32  Names       ! *Designations
      character*6   Notes       ! *Notes
      integer*4     PosFlag     ! *Position accuracy flag
*Notes on CSS:
*    The running number in file 1 is newly assigned, but numbers used in
*    the first edition (Catalog <III/60>) are given in the Designations
*    column.
*Note on RAh, RAm, RAs, DE-, DEd, DEm, DEs:
*    Equatorial coordinates, referred to the equinox and equator which is
*    standard in most spectroscopic tabulations to date, given with varying
*    precision. The coordinates are taken from the various finding lists,
*    and averaged by the author in cases of dual or more determinations of
*    nominally comparable precision. In unresolvable cases of disagreement,
*    the quoted number of significant figures, have been downgraded but the
*    original sources also vary in quoted precision. For many of the
*    objective prism surveys -- including all that have been done at the
*    Warner and Swasey Observatory since 1958 -- the coordinates came from
*    measurements of objective prism plates, using methods described in the
*    Warner and Swasey Observatory Publs. Vol. 2, pp. 74-76. Almost all of
*    the stars for which Stephenson is listed among the designations have
*    independent coordinate determinations from his plates, usually to
*    +/-1" - +/-2".
*Note on bmag, vmag, imag:
*    These are, as available, blue (photographic, sometimes B), visual
*    (sometimes red-biased), and infrared magnitudes. The infrared
*    magnitudes refer generally to the 6800- 8800 A region. For known
*    variable stars, the variable star catalog (GCVS) magnitudes at
*    mid-range are cited under the type of magnitude that was given in the
*    GCVS; other quoted magnitude types for known variables are from survey
*    plates taken at unknown phases. No general statements can be made
*    about the survey-magnitudes at large, and the original papers must be
*    consulted for details. For Stephenson's surveys, blue magnitudes are
*    based upon individually-calibrated plates, and visual ones are based
*    upon an average calibration used for the whole survey; the former may
*    carry +/-0.3 mag. probable error, the latter +/-l mag. In all cases, a
*    value of 0.0 implies no value was available. These values are NOT
*    included in the published catalog, but were tabulated along with the
*    published data. Caveat emptor.
*Note on bmag_pub, vmag_pub, imag_pub:
*    These magnitude string values are the values published in the Catalog.
*    Some contain the ":" character that indicates that the value given is
*    more uncertain than most.
*Note on SpType:
*   This column furnishes, as available, spectral types in various systems.
*   Omitting sources that are either rare or also used under Designations,
*   the abbreviations define their sources as follows:
*                 Bid57     Bidelman (1957)
*                 Bid65     Bidelman (1965)
*                 CatchF    Catchpole & Feast (1971)
*                 De        Dean (1972)
*                 F66       Feast (1966)
*                 Houk      HD reclassification
*                 K         Keenan (1954)
*                 KB        Keenan & Boeshaar (1980)
*                 K66       Keenan (1966)
*                 K74       Keenan (1974)
*                 Merr      Merrill (1940)
*                 Sanf41    Sanford (1941)
*                 Sanf49    Sanford (1949)
*                 Sharp     Sharpless (1964)
*                 WPB       Bidelman (1954a)
*    The stars that have been called MS on objective prism plates are in
*    general stronger S stars than those called MS from slit spectrograms,
*    and also -- at least the ones found by Stephenson and colleagues --
*    are stronger than the weaker S stars published by MacConnell and
*    Henize, who used the highest spectral dispersions of any of the modern
*    surveys. Other spectral symbols may be familiar to most, but will be
*    defined here for completeness: e or E, line emission; p, peculiar; r,
*    uncommonly red; wk or vwk, weak or very weak; S! (used only in
*    Stephenson's surveys), outstanding ZrO strength; colon, uncertainty.
*    Paired numerals flanking either a comma, slash, or asterisk are TiO
*    temperature class and a ZrO strength measure (left to right), in
*    various authors' systems.
*
*    C-S or SC refer to stars that spectroscopically are nearly
*    intermediate between S-type and carbon stars, due (among other things)
*    to closely comparable abundances of carbon and oxygen. The C-S type
*    was introduced by Stephenson and means stars which in the red spectral
*    region at about 1000 A/mm are very red with strong CN, like an N-type
*    carbon star, but show no definite C2 although they sometimes give the
*    appearance of showing weak-to-moderate 6474 ZrO, which however is
*    uncertain because of the strong CN and low dispersion. In general,
*    given sufficient spectral resolution they can be classified as weak S
*    or weak carbon, although they share many properties of both groups.
*    The SC nomenclature was introduced slightly later and has been used in
*    two ways: to denote all the C-S stars, or alternately those that just
*    fall on the S side, in which convention the carbon S-like group
*    becomes CS. See the discussion in Keenan and Boeshaar (1980). In the
*    present catalog the intent has been to include only the C-S stars that
*    are likely S types.
*Note on Names:
*    Here are listed the star's occurrence in various spectroscopic (and
*    some other) finding lists. These are mainly surveys of one kind or
*    another; exceptions will be listed shortly. The author has tried to
*    cite first the earliest survey calling the star an S star (or
*    equivalent, in the case of the HD catalog), Not included here are slit
*    spectrograph observations of known or suspected S stars; those go in
*    the "spectrum" column. Independent spectral classifications are cited
*    here (but not repeated, if already under "spectrum") either if they
*    are quantitative or partially so, or if the star has ever been called
*    a non-S star. The HD number, if one exists, is given without the
*    identifying letters "HD", almost always before anything else.
*
*    Non-spectroscopic finding list numbers cited are: the HR number,
*    labelled HR; the B.D. or Co.D. numbers, in an obvious notation; and
*    the GCVS (General Catalog of Variable Stars; catalog <II/139>) or CSV
*    (Catalog of Stars Suspected of Being Variable; catalog <II/140>)
*    numbers. CSV numbers are from the editions prior to the 1982 one; they
*    have not been compared with the newest edition of the CSV, partly
*    because of the different equinoxes of that catalog and this one, and
*    partly because it is never particularly noteworthy that an S star has
*    been suspected of variability.
*
*    Coding of sources: Sources are labelled by the abbreviations given
*    below, with numbers assigned by the corresponding discoverer separated
*    from his abbreviation by a single space. Numbers not separated from
*    the abbreviation are part of the abbreviation, and hyphens if any
*    separate a numbered published table from numbers within the table. The
*    abbreviations mean the following:
*
*   Bid65          Bidelman (1965).
*   BidK           Bidelman and Krumenaker (1972).
*   Bidun          Bidelman, unpublished, usually from a slit spectrogram.
*   BidMac         Bidelman and MacConnell (1973).
*   BlancN         Blanco and Nassau (1957).
*   BM             Blanco and Munch (1955).
*   BSD            Bergedorfer Spektraldurchmusterung (Schwassmann and Van Rhijn
*                     1935).
*   CatchF         Catchpole and Feast (1971).
*   Cra83          Craine et al (1983).
*   D              Lee et al (1943, 1944, 1947).  The present catalog has not
*                     been fully compared with the 44,000 stars of the D lists,
*                     except for the stars called S, Sp, or S?
*   De             Dean (1972).
*   Dol61          Dolidze (1961).
*   Dol62a            "    (1962a).
*   Dol62b            "    (1962b).
*   Dol65             "    (1965).
*   Dol68             "    (1968).
*   Dol70a            "    (1970a).
*   Dol70b            "    (1970b).
*   Dol70c            "    (1970c).
*   Dol71             "    (1971).  In most of the lists of this series the
*                     coordinates seem to be uncertain by 5' - 10'. Charts
*                     for the stars of this series are in Abastumani Bull.
*                     Vol. 47, 10, 1975.
*   DolJim         Dolidze and Jimsheleishvili (1966).
*   Harw           Harwood (1962).
*   Hen60          Henize (1960).
*   Henun          Henize, unpublished.  These stars are from a list dated 1965,
*                     and are in fact published in that they appear in the
*                     first edition of the present catalog.
*   Hetzler        Hetzler (1937).
*   KB             Keenan and Boeshaar (1980).
*   Krum           Krumenaker (1975).
*   MacCon         MacConnell, unpublished.
*   MacCon67           "      (1967).
*   MacCon79           "      (1979).
*   MacCon82           "      (1982).
*   MacRae         MacRae (1952).
*   Merr           Merrill (1940).
*   MSB            Merrill, Sanford, and Burwell (1933,1942).  Without a number,
*                    unpublished other than in the first edition of this catalog.
*   Nass54         Nassau, Blanco, and Morgan (1954).
*   NassCam        Nassau and Cameron (1956).
*   NassSte        Nassau and Stephenson (1961).  Without a number, unpublished
*                     except for the first ed. of the S-star catalog.
*   Nass64         Nassau, Stephenson and Caprioli (1964).
*   Perr59         Perraud (1959).
*   Perr61            "    (1961).
*   Pes            Pesch, unpublished (or only in 1st ed. this catalog).
*   Rust           Rust (1938).
*   S              First edition of this catalog.
*   Sndlk          Sanduleak, unpublished (or only in 1st ed. this catalog).
*   Ste            Stephenson,    "        "   "   "   "   "    "     "
*   Ste65               "    (1965).
*   Ste73               "    (1973).
*   SteTerr        Stephenson and Terrill (1967).
*   Sto            Stock and Wroblewski (1972).
*   S-WS           Two Case lists of S stars: Nassau, Blanco and Morgan (1954);
*                     Blanco and Nassau (1957).
*   The68          The (1968).
*   VB             Blanco, unpublished except for 1st ed. this catalog.
*   Vys            Vyssotsky (1942), Janssen and Vyssotsky (1943), Vyssotsky and
*                    Miller (1946), always with a star number.  Without a number,
*                    Vyssotsky unpublished.
*   VysBalz        Vyssotsky and Balz (1958).
*   West'd         Westerlund (unpublished).
*   Wray           Wray (1966). Numbers are from his combined table of carbon and
*                     S stars.  His type SE means S with H alpha emission.
*
*    The various lists occasionally identify with one another, but, except
*    for the BD numbers, most identifications of the stars of one list with
*    those of another have been made by the author.
*
*    Additional identifications are often given in the Remarks to
*    "notes" file; unlike the first edition, in the present edition the
*    main catalog allots only one line to a given star.
*
*    If the star was in the first edition of this catalog, its number there
*    is always given here, coded as "S".
*Note on GLON, GLAT:
*     Galactic longitude and latitude, computed by Stephenson.
*Note on Notes:
*    Several symbols are used here, as follows. An asterisk means that two
*    or more sources included under "designations" contradict each other as
*    to the star's identity, including cases where a star in an earlier
*    list was published as new in a later list. A single asterisk means
*    that the author's revision of the identification was done without
*    specially examining an objective prism plate, while a double asterisk
*    means the case was specially checked. A plus sign represents a new
*    identification with a GCVS star, done by the author by comparing a
*    plate of the sky with a published identification chart. WPB stands for
*    Bidelman (1954). An R means that there are remarks in the Notes to
*    file 1 in notes.txt.
*
*    References to Case plates in the Notes mean Warner & Swasey
*    Observatory objective prism plates, with the plate examiner almost
*    always being Stephenson
*
*   Stephenson is always an independent source of uncredited remarks, though not
*   always the earliest such source.
*Note on PosFlag:
*   The position accuracy is described by the following codes:
*      1     hh mm ss.s  dd mm ss
*      2     hh mm ss    dd mm.m
*      3     hh mm       dd mm

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

C  Declarations for 'notes.dat'	! Notes to catalog

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

      integer*4     CSS_1       ! [1/1347]+= Number
      character*117 note        ! Notes

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

C  Declarations for 'reject.dat'	! *Suggested S stars rejected

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

      character*10  RAa         ! *Right ascension (1900)
      character*9   DEa         ! *Declination (1900)
      character*23  name        ! Designations
      character*23  note_1      ! *Notes
*Note on RAa, DEa:
*   The position accuracy in this file varies. For stars with the higher
*   accuracy, the minutes of right ascension are in bytes 4-5, the seconds of
*   right ascension are in bytes 7-10 and the seconds of declination are in bytes
*   20-21. For stars whose position accuracy is more crude, the minutes of right
*   ascension are in bytes 7-10 and bytes 20-21 are blank.
*Note on note:
*   Numerals in parentheses have the following meaning:
*      1.  No candidate star could be found by Stephenson on a blue-region
*          objective prism plate.
*      2.  No candidate star could be found by Stephenson on a red-region
*          objective prism plate.
*      3.  No candidate star could be found by Stephenson on a infrared-region
*          objective prism plate.
*      4.  Candidate star classified by Stephenson on a blue-region objective
*          prism plate.
*      5.  Candidate star classified by Stephenson on a red-region objective
*          prism plate.
*      6.  Candidate star classified by Stephenson on a infrared-region objective
*          prism plate.
*      7.  Slit spectrogram by Stephenson.

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

C  Loading file 'catalog.dat'	! The S Star catalog

C  Format for file interpretation

    1 format(
     +  I4,2X,I2,1X,I2,1X,F4.1,2X,A1,I2,1X,I2,1X,I2,2X,F4.1,1X,F4.1,
     +  1X,F4.1,2X,F5.1,1X,F5.1,2X,A5,1X,A5,1X,A5,2X,A16,2X,A32,2X,A6,
     +  2X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,1347
        read(1,'(A138)')ar__
        read(ar__,1)
     +  CSS,RAh,RAm,RAs,DE_,DEd,DEm,DEs,bmag,vmag,imag,GLON,GLAT,
     +  bmag_pub,vmag_pub,imag_pub,SpType,Names,Notes,PosFlag
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) 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)
     +  CSS,RAh,RAm,RAs,DE_,DEd,DEm,DEs,bmag,vmag,imag,GLON,GLAT,
     +  bmag_pub,vmag_pub,imag_pub,SpType,Names,Notes,PosFlag
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'notes.dat'	! Notes to catalog

C  Format for file interpretation

    2 format(I4,3X,A117)

C  Effective file loading

      open(unit=1,status='old',file=
     +'notes.dat')
      write(6,*) '....Loading file: notes.dat'
      do i__=1,440
        read(1,'(A124)')ar__1
        read(ar__1,2)CSS_1,note
c    ..............Just test output...........
        write(6,2)CSS_1,note
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'reject.dat'	! *Suggested S stars rejected

C  Format for file interpretation

    3 format(A10,3X,A9,3X,A23,1X,A23)

C  Effective file loading

      open(unit=1,status='old',file=
     +'reject.dat')
      write(6,*) '....Loading file: reject.dat'
      do i__=1,56
        read(1,'(A72)')ar__2
        read(ar__2,3)RAa,DEa,name,note_1
c    ..............Just test output...........
        write(6,3)RAa,DEa,name,note_1
c    .......End.of.Just test output...........
      end do
      close(1)

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