FORTRAN Generation
(/./ftp/cats/II/250/)

Conversion of standardized ReadMe file for file /./ftp/cats/II/250/ 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-Apr-16
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. II/250              Combined General Catalogue of Variable Stars   (Samus+ 2004)
*================================================================================
*Combined General Catalog of Variable Stars (GCVS4.2, 2004 Ed.)
*     Samus N.N., Durlevich O.V., et al.
*    <Institute of Astronomy of Russian Academy of Sciences and Sternberg
*     State Astronomical Institute of the Moscow State University (2004)>
*     =2004yCat.2250....0S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'gcvs_cat.dat'	! The GCVS4.2 Catalog (Vol. I-III)

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

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

      integer*4     VarNum     (nr__) ! Numeric designation, made from constellation
*                                 and star names (codes in file "constel.txt")
      character*1   m_VarNum   (nr__) ! Component identification (1)
      character*10  GCVS       (nr__) ! Variable star designation (G1)
      character*1   n_GCVS     (nr__) ! [*] indicates a remark in "gcvs_rem.dat" file
      integer*4     RA1950h    (nr__) ! (h) ? Hours RA, equinox 1950.0 (2)
      integer*4     RA1950m    (nr__) ! (min) ? Minutes RA, equinox 1950.0 (2)
      real*4        RA1950s    (nr__) ! (s) ? Seconds RA, equinox 1950.0 (2)
      character*1   DE1950_    (nr__) ! ? Sign Dec, equinox 1950.0 (2)
      integer*4     DE1950d    (nr__) ! (deg) ? Degrees Dec, equinox 1950.0 (2)
      integer*4     DE1950m    (nr__) ! (arcmin) ? Minutes Dec, equinox 1950.0 (2)
      integer*4     DE1950s    (nr__) ! (arcsec) ? Seconds Dec, equinox 1950.0 (2)
      character*1   u_DE1950s  (nr__) ! [:*] position accuracy flags (G2)
      integer*4     RAh        (nr__) ! (h) ? Hours RA, equinox J2000.0 (2)
      integer*4     RAm        (nr__) ! (min) ? Minutes RA, equinox J2000.0 (2)
      real*4        RAs        (nr__) ! (s) ? Seconds RA, equinox J2000.0 (2)
      character*1   DE_        (nr__) ! Sign Dec, equinox J2000.0 (2)
      integer*4     DEd        (nr__) ! (deg) ? Degrees Dec, equinox J2000.0 (2)
      integer*4     DEm        (nr__) ! (arcmin) ? Minutes Dec, equinox J2000.0 (2)
      integer*4     DEs        (nr__) ! (arcsec) ? Seconds Dec, equinox J2000.0 (2)
      character*1   u_DEs      (nr__) ! [:*] position accuracy flags (G2)
      character*1   n_Pos      (nr__) ! [*] Position computed from B1950 (3)
      character*10  VarType    (nr__) ! Type of variability (see file "vartype.txt")
      character*1   l_magMax   (nr__) ! [<>(] Limit or amplitude symbol on magMax (G3)
      real*4        magMax     (nr__) ! (mag) ? Magnitude at maximum brightness
      character*1   u_magMax   (nr__) ! Uncertainty flag (:) on magMax
      character*1   l_magMin   (nr__) ! [<(] Limit or amplitude symbol on magMin (G3)
      real*4        magMin     (nr__) ! (mag) ? Minimum magnitude or amplitude
      character*1   u_magMin   (nr__) ! Uncertainty flag (:) on magMin
      character*2   n_magMin   (nr__) ! Alternative photometric system for magMin (G4)
      character*2   magCode    (nr__) ! The photometric system for magnitudes (G4)
      real*8        Epoch      (nr__) ! (d) ? Epoch for maximum light, Julian days (G5)
      character*1   u_Epoch    (nr__) ! [:+-] Uncertainty flag on Epoch (4)
      character*4   Year       (nr__) ! Year of outburst for nova or supernova
      character*1   u_Year     (nr__) ! [:] Uncertainty flag on Year of outburst
      character*1   l_Period   (nr__) ! [<>(] Code for upper or lower limits (5)
      real*8        Period     (nr__) ! (d) ? Period of the variable star
      character*3   u_Period   (nr__) ! [*/N)2: ] Uncertainties on Period (6)
      character*3   M_m_D      (nr__) ! (%) ? Rising time or duration of eclipse (G6)
      character*1   u_M_m_D    (nr__) ! Uncertainty flag (:) on M-m/D
      character*1   n_M_m_D    (nr__) ! [*] Note for eclipsing variable (G6)
      character*17  SpType     (nr__) ! MK Spectral type
      character*5   Ref1       (nr__) ! Reference to a study of the star (G7)
      character*5   Ref2       (nr__) ! Reference to a chart or photograph (G7)
      character*1   f_GCVS     (nr__) ! [=N+] "N" if the star does not exist (7)
      character*11  VarName    (nr__) ! Alternative name of the variable (G8)
*Note (1): designations of components of double/multiple stars :
*   A,B; 1,2,3  ... (tau 1). Usually only one component of the binary
*   is variable, but there exist two cases (CE Cas A,B; EQ Peg A,B) when
*   both components are variable and enter the catalogue as separate lines.
*Note (2):
*    the position is not listed when the variable is equivalent to another
*    catalogued variable (see column VarName)
*Note (3):
*    the asterisk (*) means that the coordinates that could not be improved
*    were recalculated from old rough coordinates to the equinox J2000.0
*    with allowance for precession.
*Note (4): the following codes are used:
*    ":" are used if the epoch is uncertain;
*    "+" epoch may be later than that reported;
*    "-" epoch may be earlier than that reported.
*Note (5):
*    ">" or "<" if period is a lower or upper limit, respectively.
*    "(" indicates that the period is the mean cycle time of a U Gem
*        or recurrent nova (a closing bracket exists in u_Period)
*Note (6): the codes are:
*  "*N" if the period may be a multiple of the quoted value;
*  "/N" if the period may be a fraction of the quoted value;
*   ")" closing bracket of l_Period for U Gem-type or recurrent novae.
*Note (7): Besides "N", the codes are:
*   "=" if the star is also known under another variable name
*   "+" if the star is in other catalogues (e.g. Hipparcos or Name Lists)

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

C  Declarations for 'gcvs_rem.dat'	! Remarks on GCVS stars

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

      character*10  GCVS_1     (nr__1) ! Variable star designation
      character*81  Text       (nr__1) ! Text of note.

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

C  Declarations for 'crossid.dat'	! Cross-identifications (Vol. IV)

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

      character*4   Cat        (nr__2) ! Abbreviated name of catalog or list
*                                  (details in the file "abbrev.txt")
      character*25  Seq        (nr__2) ! Star number or name in the respective catalog
      character*12  VarName_1  (nr__2) ! Designation in GCVS or NSV (G8)
      character*1   f_VarName  (nr__2) ! [:?] Flag for doubtful identifications

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

C  Declarations for 'nsv_cat.dat'	! The Suspected Variable stars and Supplement

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

C  J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1    (nr__3) ! (deg) Right Ascension J2000.0
      real*8        DEdeg_1    (nr__3) ! (deg)     Declination J2000.0
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     NSV        (nr__3) ! NSV number (above 15000 for the Supplement)
      character*1   m_NSV      (nr__3) ! NSV letter suffix (1)
      character*1   u_NSV      (nr__3) ! [-] Dubitancy flag if a star's variability
*                                 seems doubtful or erroneous to the compilers
      character*1   n_NSV      (nr__3) ! [*] indicate a remark in file "nsv.rem.dat"
      integer*4     RA1950h_1  (nr__3) ! (h) ? Hours RA, equinox 1950.0 (2)
      integer*4     RA1950m_1  (nr__3) ! (min) ? Minutes RA, equinox 1950.0 (2)
      real*4        RA1950s_1  (nr__3) ! (s) ? Seconds RA, equinox 1950.0 (2)
      character*1   DE1950__1  (nr__3) ! ? Sign Dec, equinox 1950.0 (2)
      integer*4     DE1950d_1  (nr__3) ! (deg) ? Degrees Dec, equinox 1950.0 (2)
      integer*4     DE1950m_1  (nr__3) ! (arcmin) ? Minutes Dec, equinox 1950.0 (2)
      integer*4     DE1950s_1  (nr__3) ! (arcsec) ? Seconds Dec, equinox 1950.0 (2)
      character*1   u_DE1950s_1(nr__3) ! [:*] accuracy flags (G2)
      integer*4     RAh_1      (nr__3) ! (h) ? Hours RA, equinox J2000.0 (2)
      integer*4     RAm_1      (nr__3) ! (min) ? Minutes RA, equinox J2000.0 (2)
      real*4        RAs_1      (nr__3) ! (s) ? Seconds RA, equinox J2000.0 (2)
      character*1   DE__1      (nr__3) ! ? Sign Dec, equinox J2000.0 (2)
      integer*4     DEd_1      (nr__3) ! (deg) ? Degrees Dec, equinox J2000.0 (2)
      integer*4     DEm_1      (nr__3) ! (arcmin) ? Minutes Dec, equinox J2000.0 (2)
      integer*4     DEs_1      (nr__3) ! (arcsec) ? Seconds Dec, equinox J2000.0 (2)
      character*1   u_DEs_1    (nr__3) ! [:*] accuracy flags (G2)
      character*6   VarType_1  (nr__3) ! Type of variability (see file "vartype.txt")
      character*1   l_magMax_1 (nr__3) ! [><(] Limit or amplitude symbol on magMax (G3)
      real*4        magMax_1   (nr__3) ! (mag) ? Magnitude at maximum brightness (3)
      character*1   u_magMax_1 (nr__3) ! [:] Uncertainty flag on magMax
      character*1   f_magMin   (nr__3) ! [(] indicates that magMin is an amplitude
      character*1   l_magMin_1 (nr__3) ! [><] Limit symbol on magMin (G3)
      real*4        magMin_1   (nr__3) ! (mag) ? Minimum magnitude or amplitude (3)
      character*1   u_magMin_1 (nr__3) ! [st:] Uncertainty or steps on magMin (4)
      character*1   n_magMin_1 (nr__3) ! Photometric system for magMin (G4)
      character*2   magCode_1  (nr__3) ! Photometric system for magnitudes (G4)
      character*5   Ref1_1     (nr__3) ! Reference to a study of the star (G7)
      character*5   Ref2_1     (nr__3) ! Chart reference (G7)
      character*15  Desig      (nr__3) ! Designation in paper referenced by Ref1 (5)
      character*15  SpType_1   (nr__3) ! Spectral type (6)
      character*1   f_NSV      (nr__3) ! [=] for named variables.
      character*11  VarName_2  (nr__3) ! Designation in GCVS (G8)
*Note (1):
*    The letter "A" is present for one star, NSV 10360A, which is a
*    completely different object from NSV 10360.
*Note (2):
*   The position is not listed when the variable is equivalent to a GCVS
*   star (see column VarName).
*Note (3):
*    Magnitudes are reported to hundredths if the observations are
*    photoelectric, to tenths or whole magnitudes if they are not.
*    If only an amplitude has been measured photoelectrically, then
*    the maximum magnitude is generally given to tenths only and the
*    minimum is reported to hundredths.
*Note (4):
*    "st" indicates that a light amplitude is reported for the minimum
*    and is given in steps (instrumental measure).
*Note (5): (empty in the NSV Supplement part):
*    identification of the suspected variable in the paper referenced in
*    Ref1. DM numbers are given without prefix, the standard naming
*    convention of "The Henry Draper Catalogue" being used (BD for zones
*    +89 to -22; CD for -23 to -51; CP for -52 to -89). Identifications
*    are sometimes made of coordinates, a 6-digit number consisting of
*    hours, minutes, and seconds (or tenths of a minute) of time and
*    degrees of declination with sign included. An asterisk (*) signifies
*    that a designation is given in the remarks to the published catalog.
*Note (6):
*  HD spectral types are generally in brackets, e.g. (A0).
*  Some special symbols are used:
*       AF  for  A-F stars
*       FG  for  F-G stars
*       KM  for  K-M stars

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

C  Declarations for 'nsv_rem.dat'	! Remarks on NSV (and Supplement) stars

      integer*4 nr__4
      parameter (nr__4=10776)	! Number of records
      character*80 ar__4  	! Full-size record

      integer*4     NSV_1      (nr__4) ! NSV number (above 15000 for the Supplement)
      character*75  Text_1     (nr__4) ! Text of remark

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

C  Declarations for 'evs_cat.dat'	! Extragalactic Variable Stars. Catalogue (Vol. V)

      integer*4 nr__5
      parameter (nr__5=10979)	! Number of records
      character*160 ar__5  	! Full-size record

C  1950.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_2    (nr__5) ! (deg) Right Ascension 1950.0
      real*8        DEdeg_2    (nr__5) ! (deg)     Declination 1950.0
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     VarNum_1   (nr__5) ! [890001/1230042] Variable number (1)
      character*12  Name       (nr__5) ! Variable star designation (3)
      character*1   n_Name     (nr__5) ! [*] indicates a remark in file evs_rem.dat
      integer*4     RAh_2      (nr__5) ! (h) Right ascension, equinox 1950.0
      integer*4     RAm_2      (nr__5) ! (min) (minutes)
      real*4        RAs_2      (nr__5) ! (s) (seconds)
      character*1   DE__2      (nr__5) ! sign of declination, 1950.0
      integer*4     DEd_2      (nr__5) ! (deg) Declination, equinox 1950.0
      integer*4     DEm_2      (nr__5) ! (arcmin) (arc minutes)
      real*4        DEs_2      (nr__5) ! (arcsec) (arc seconds)
      character*1   u_DEs_2    (nr__5) ! [":*!] accuracy flags on position (G2)
      character*8   VarType_2  (nr__5) ! Type of variability (see file vartype.txt)
      real*4        magMax_2   (nr__5) ! (mag) ? Magnitude at maximum brightness
      character*1   u_magMax_2 (nr__5) ! Uncertainty flag (:) on magMax
      character*1   l_magMin_2 (nr__5) ! [<(] Limit or amplitude symbol on magMax (G3)
      real*4        magMin_2   (nr__5) ! (mag) ? Magnitude at minimum brightness or amplitude
      character*1   u_magMin_2 (nr__5) ! Uncertainty flag (:) on magMin
      character*1   magCode_2  (nr__5) ! The photometric system for magnitudes (G4)
      real*8        Epoch_1    (nr__5) ! (d) ? Epoch for maximum light, Julian days (G5)
      character*1   u_Epoch_1  (nr__5) ! Uncertainty flag on Epoch
      character*1   l_P        (nr__5) ! [(] if mean cycle time of U Gem-type variable
      real*8        P          (nr__5) ! (d) ? Period (mean for U Gem-type)
      character*1   u_P        (nr__5) ! Uncertainty flag (:) on P
      integer*4     M_m_D_1    (nr__5) ! (%) ? Rising time or duration of eclipse (G6)
      character*1   u_M_m_D_1  (nr__5) ! Uncertainty flag (:) on M-m/D
      character*1   n_M_m_D_1  (nr__5) ! [*] Note for eclipsing variable (G6)
      character*6   SpType_2   (nr__5) ! Spectral type (4)
      character*5   Ref1_2     (nr__5) ! Reference for study of the star (G7)
      character*5   Ref2_2     (nr__5) ! Reference for chart or photograph (G7)
      character*14  Desig1     (nr__5) ! designation in paper Ref1
      character*12  Desig2     (nr__5) ! designation in paper Ref2
      character*9   VarName_3  (nr__5) ! Designation in GCVS (G8)
      character*2   nonMem     (nr__5) ! [n: ] Non-member note (5)
      integer*4     Year_1     (nr__5) ! (yr) ? Year of outburst for novae or supernovae
      character*1   u_Year_1   (nr__5) ! [:] Uncertainty flag on YearNova
*Note (1): Number of the variable star.
*    This numbering system is made of a 3-digit number identifying the
*    galaxy, followed by a 4-digit number identifying the variable in the
*    galaxy. The numbering system for galaxies begins with 089, the first
*    88 codes being reserved for constellations (cf. column "VarNum" in
*    file "gcvs_cat.dat")
*
*    If in the literature for the considered galaxy there exist several
*    different variable star numbering systems, we have introduced new
*    numbers in the order of increasing right ascension for the equinox
*    1950.0. But if the variable star numbering system for the galaxy is
*    practically of general use, we have retained it - if necessary, with
*    additional numbers ordered by increasing right ascension.
*Note (3): Variable star Designation.
*    The name of the variable in the adopted system of variable star
*    designations: name of the galaxy + "V" + variable star No. within the
*    galaxy (see column VarNo)
*Note (4):
*    If there exist several spectral type determinations for a star, we
*    preferred more recent ones, taking into account reliability of each
*    determination. We converted symbols "d" (dwarfs) and "g" (giants) into
*    luminosity classes V and III. If a nova showed spectral features
*    typical of novae during or after the outburst, a symbol (NOVA) is
*    given in this column. If the spectrum shows features characteristic of
*    U Gem variables, this column contains the symbol (UG). Continuous
*    spectrum is designated "cont", the symbol "e" means emission spectrum,
*    lower case letters "ea", "eb" in this column stand for "e alpha", "e
*    beta". For RR Lyrae stars, spectral types from hydrogen lines are
*    given. A plus sign (+) between two spectral type values means that
*    spectra of two components of a spectroscopic binary are observed. Two
*    spectral type values separated with a minus sign (-) mean the range of
*    spectral type variations in the process of brightness changes.
*Note (5):
*    Symbols "n" or "n:" for "nonmembers" or "possible nonmembers" of the
*    galaxy in question. For the stars that enter the GCVS vols. I-III
*    (file gcvs_cat.dat) or the NSV catalogue (file "nsv_cat.dat") and are
*    not members of the corresponding galaxy, the present volume gives only
*    the running number, coordinates, the symbol "n", and the GCVS or NSV
*    name.

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

C  Declarations for 'evs_rem.dat'	! Extragalactic Variable Stars. Remarks

      integer*4 nr__6
      parameter (nr__6=1936)	! Number of records
      character*80 ar__6  	! Full-size record

      character*12  Name_1     (nr__6) ! Variable star designation
      character*67  Text_2     (nr__6) ! Text of remark

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

C  Declarations for 'sn_cat.dat'	! Extragalactic supernovae. Catalogue

      integer*4 nr__7
      parameter (nr__7=984)	! Number of records
      character*132 ar__7  	! Full-size record

C  1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_3    (nr__7) ! (deg) Right Ascension 1950
      real*8        DEdeg_3    (nr__7) ! (deg)     Declination 1950
C  ---------------------------------- ! (position vector(s) in degrees)

      character*6   SN         (nr__7) ! Supernova designation
      character*1   f_SN       (nr__7) ! [-?] Doubtful (?), or rejected (-) SN
      character*1   n_SN       (nr__7) ! [*] indicates a remark in file "sn_rem.dat"
      character*10  Gal        (nr__7) ! Parent galaxy designation (1)
      integer*4     RAh_3      (nr__7) ! (h) ? RA 1950 (hours) of Parent galaxy
      integer*4     RAm_3      (nr__7) ! (min) ? Right Ascension 1950 (minutes)
      real*4        RAs_3      (nr__7) ! (s) ? Right Ascension 1950 (seconds)
      character*1   DE__3      (nr__7) ! ? Declination 1950 (sign)
      integer*4     DEd_3      (nr__7) ! (deg) ? Declination Deg. 1950 of Parent galaxy
      integer*4     DEm_3      (nr__7) ! (arcmin) ? Declination 1950 (minutes)
      integer*4     DEs_3      (nr__7) ! (arcsec) ? Declination 1950 (seconds)
      character*1   u_DEs_3    (nr__7) ! [":'] Accuracy flags (G2)
      character*6   MType      (nr__7) ! Morphological Type of the parent galaxy
      real*4        BTmag      (nr__7) ! (mag) ? Integrated magnitude of parent galaxy
      character*1   u_BTmag    (nr__7) ! Uncertainty flag (:) on BTmag
      character*1   n_BTmag    (nr__7) ! [BJV] Band indicator, blank for photographic
      character*8   maxDate    (nr__7) ! Date of the SN's maximum or of discovery (2)
      real*4        magMax_3   (nr__7) ! (mag) ? Magnitude at maximum (light) of SN
      character*1   u_magMax_3 (nr__7) ! Uncertainty flag (:) on magMax
      character*1   n_magMax   (nr__7) ! [BJVRIKr] Band indicator (blank=photographic)
      character*1   f_magMax   (nr__7) ! [*] if maximum not observed (2)
      real*4        dRA        (nr__7) ! (arcsec) ? Distance of SN from galaxy nucleus (3)
      character*2   oRA        (nr__7) ! [NSEW:? ] Orientation of dRA
      real*4        dDE        (nr__7) ! (arcsec) ? Distance of SN from galaxy nucleus (4)
      character*2   oDE        (nr__7) ! [NSEW:? ] Orientation of dDE
      integer*4     RASN_h     (nr__7) ! (h) ? Right Ascension 1950 (hours) of SN
      integer*4     RASN_m     (nr__7) ! (min) ? Right Ascension 1950 (minutes)
      real*4        RASN_s     (nr__7) ! (s) ? Right Ascension 1950 (seconds)
      character*1   DESN__     (nr__7) ! ? Declination 1950 (sign)
      integer*4     DESN_d     (nr__7) ! (deg) ? Declination 1950 (degrees) of SN
      integer*4     DESN_m     (nr__7) ! (arcmin) ? Declination 1950 (minutes)
      integer*4     DESN_s     (nr__7) ! (arcsec) ? Declination 1950 (seconds)
      character*1   u_DESN_s   (nr__7) ! [":'] accuracy flags on SN position (G2)
      character*6   SNType     (nr__7) ! SN classification (5)
      character*16  Discov     (nr__7) ! Name of the SN discoverer (6)
      character*5   Ref1_3     (nr__7) ! Reference for study of the star (G7)
      character*5   Ref2_3     (nr__7) ! Reference for chart or photograph (G7)
*Note (1):
*   "N"  means NGC catalogue numbers,
*   "I"  refers to the IC catalogue,
*   "U"  to the UGC catalogue,
*   "M"  to the MCG catalogue,
*   "E"  to ESO survey lists.
*   "An" means an anonymous galaxy
*   "Intergal" is for cases when the supernova cannot be reliably
*   attributed to a particular galaxy and is considered intergalactic.
*Note (2):
*   The date of the SN's maximum or of discovery, or of an isolated
*   observation at high brightness. If maxDate or f_magMax contains an
*   asterisk (*), it means that the date and/or the magnitude refer not
*   to the maximum (not covered with observations), but to the EPOCH OF
*   DISCOVERY or to an isolated observation at high brightness. As a
*   rule, information about maxima was not included if it was based on
*   considerable extrapolation.
*Note (3):
*   In line 393, 1970K, the offset is 0.3E
*Note (4):
*   In line 732, 1988T, the offset is 0.4N
*Note (5):
*   Limited to types I (I-A, I-B) and II (II-l, II-pl), with possible
*   indications of peculiarities (pec).
*Note (6):
*   In case of several discoverers, the name of the first author of the
*   discovery is given followed by "+" symbol. As a rule, information on
*   independent rediscoveries is not presented in the remarks, for the
*   sake of brevity. The abbreviation "BASST" in this column means that
*   the supernova was discovered by the Berkeley Automated Supernova
*   Search Team.

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

C  Declarations for 'sn_rem.dat'	! Extragalactic supernovae. Remarks

      integer*4 nr__8
      parameter (nr__8=277)	! Number of records
      character*80 ar__8  	! Full-size record

      character*6   SN_1       (nr__8) ! Supernova designation, repeated if more than 1 line
      character*1   n_SN_1     (nr__8) ! [-?] Doubtful (?), or rejected (-) SN
      character*72  Text_3     (nr__8) ! Text of remarks

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

C  Declarations for 'refs.dat'	! References

      integer*4 nr__9
      parameter (nr__9=21004)	! Number of records
      character*80 ar__9  	! Full-size record

      character*5   Ref        (nr__9) ! Reference number (G7)
      character*74  Text_4     (nr__9) ! Text of reference

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

C  Loading file 'gcvs_cat.dat'	! The GCVS4.2 Catalog (Vol. I-III)

C  Format for file interpretation

    1 format(
     +  I6,A1,1X,A10,A1,1X,I2,I2,F4.1,A1,I2,I2,I2,A1,1X,I2,I2,F4.1,A1,
     +  I2,I2,I2,A1,A1,A10,1X,A1,F6.3,A1,2X,A1,1X,F6.3,A1,A2,2X,A2,1X,
     +  F13.5,A1,1X,A4,A1,1X,A1,F16.10,A3,1X,A3,A1,A1,1X,A17,1X,A5,1X,
     +  A5,1X,A1,A11)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gcvs_cat.dat')
      write(6,*) '....Loading file: gcvs_cat.dat'
      do i__=1,38624
        read(1,'(A181)')ar__
        read(ar__,1)
     +  VarNum(i__),m_VarNum(i__),GCVS(i__),n_GCVS(i__),RA1950h(i__),
     +  RA1950m(i__),RA1950s(i__),DE1950_(i__),DE1950d(i__),
     +  DE1950m(i__),DE1950s(i__),u_DE1950s(i__),RAh(i__),RAm(i__),
     +  RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),u_DEs(i__),
     +  n_Pos(i__),VarType(i__),l_magMax(i__),magMax(i__),
     +  u_magMax(i__),l_magMin(i__),magMin(i__),u_magMin(i__),
     +  n_magMin(i__),magCode(i__),Epoch(i__),u_Epoch(i__),Year(i__),
     +  u_Year(i__),l_Period(i__),Period(i__),u_Period(i__),
     +  M_m_D(i__),u_M_m_D(i__),n_M_m_D(i__),SpType(i__),Ref1(i__),
     +  Ref2(i__),f_GCVS(i__),VarName(i__)
        if(ar__(21:22) .EQ. '') RA1950h(i__) = iNULL__
        if(ar__(23:24) .EQ. '') RA1950m(i__) = iNULL__
        if(ar__(25:28) .EQ. '') RA1950s(i__) = rNULL__
        if(ar__(30:31) .EQ. '') DE1950d(i__) = iNULL__
        if(ar__(32:33) .EQ. '') DE1950m(i__) = iNULL__
        if(ar__(34:35) .EQ. '') DE1950s(i__) = iNULL__
        if(ar__(38:39) .EQ. '') RAh(i__) = iNULL__
        if(ar__(40:41) .EQ. '') RAm(i__) = iNULL__
        if(ar__(42:45) .EQ. '') RAs(i__) = rNULL__
        if(ar__(47:48) .EQ. '') DEd(i__) = iNULL__
        if(ar__(49:50) .EQ. '') DEm(i__) = iNULL__
        if(ar__(51:52) .EQ. '') DEs(i__) = iNULL__
        if(ar__(67:72) .EQ. '') magMax(i__) = rNULL__
        if(ar__(78:83) .EQ. '') magMin(i__) = rNULL__
        if(ar__(92:104) .EQ. '') Epoch(i__) = rNULL__
        if(ar__(114:129) .EQ. '') Period(i__) = rNULL__
        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(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  VarNum(i__),m_VarNum(i__),GCVS(i__),n_GCVS(i__),RA1950h(i__),
     +  RA1950m(i__),RA1950s(i__),DE1950_(i__),DE1950d(i__),
     +  DE1950m(i__),DE1950s(i__),u_DE1950s(i__),RAh(i__),RAm(i__),
     +  RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),u_DEs(i__),
     +  n_Pos(i__),VarType(i__),l_magMax(i__),magMax(i__),
     +  u_magMax(i__),l_magMin(i__),magMin(i__),u_magMin(i__),
     +  n_magMin(i__),magCode(i__),Epoch(i__),u_Epoch(i__),Year(i__),
     +  u_Year(i__),l_Period(i__),Period(i__),u_Period(i__),
     +  M_m_D(i__),u_M_m_D(i__),n_M_m_D(i__),SpType(i__),Ref1(i__),
     +  Ref2(i__),f_GCVS(i__),VarName(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 'gcvs_rem.dat'	! Remarks on GCVS stars

C  Format for file interpretation

    2 format(A10,3X,A81)

C  Effective file loading

      open(unit=1,status='old',file=
     +'gcvs_rem.dat')
      write(6,*) '....Loading file: gcvs_rem.dat'
      do i__=1,14441
        read(1,'(A94)')ar__1
        read(ar__1,2)GCVS_1(i__),Text(i__)
c    ..............Just test output...........
        write(6,2)GCVS_1(i__),Text(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'crossid.dat'	! Cross-identifications (Vol. IV)

C  Format for file interpretation

    3 format(A4,1X,A25,3X,A12,2X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'crossid.dat')
      write(6,*) '....Loading file: crossid.dat'
      do i__=1,119894
        read(1,'(A48)')ar__2
        read(ar__2,3)Cat(i__),Seq(i__),VarName_1(i__),f_VarName(i__)
c    ..............Just test output...........
        write(6,3)Cat(i__),Seq(i__),VarName_1(i__),f_VarName(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'nsv_cat.dat'	! The Suspected Variable stars and Supplement

C  Format for file interpretation

    4 format(
     +  I5,A1,A1,A1,1X,I2,I2,F4.1,A1,I2,I2,I2,A1,1X,I2,I2,F4.1,A1,I2,
     +  I2,I2,A1,1X,A6,1X,A1,F5.2,A1,1X,A1,A1,F6.3,A1,A1,2X,A2,1X,A5,
     +  1X,A5,1X,A15,1X,A15,A1,A11)

C  Effective file loading

      open(unit=1,status='old',file=
     +'nsv_cat.dat')
      write(6,*) '....Loading file: nsv_cat.dat'
      do i__=1,26017
        read(1,'(A128)')ar__3
        read(ar__3,4)
     +  NSV(i__),m_NSV(i__),u_NSV(i__),n_NSV(i__),RA1950h_1(i__),
     +  RA1950m_1(i__),RA1950s_1(i__),DE1950__1(i__),DE1950d_1(i__),
     +  DE1950m_1(i__),DE1950s_1(i__),u_DE1950s_1(i__),RAh_1(i__),
     +  RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),u_DEs_1(i__),VarType_1(i__),l_magMax_1(i__),
     +  magMax_1(i__),u_magMax_1(i__),f_magMin(i__),l_magMin_1(i__),
     +  magMin_1(i__),u_magMin_1(i__),n_magMin_1(i__),magCode_1(i__),
     +  Ref1_1(i__),Ref2_1(i__),Desig(i__),SpType_1(i__),f_NSV(i__),
     +  VarName_2(i__)
        if(ar__3(10:11) .EQ. '') RA1950h_1(i__) = iNULL__
        if(ar__3(12:13) .EQ. '') RA1950m_1(i__) = iNULL__
        if(ar__3(14:17) .EQ. '') RA1950s_1(i__) = rNULL__
        if(ar__3(19:20) .EQ. '') DE1950d_1(i__) = iNULL__
        if(ar__3(21:22) .EQ. '') DE1950m_1(i__) = iNULL__
        if(ar__3(23:24) .EQ. '') DE1950s_1(i__) = iNULL__
        if(ar__3(27:28) .EQ. '') RAh_1(i__) = iNULL__
        if(ar__3(29:30) .EQ. '') RAm_1(i__) = iNULL__
        if(ar__3(31:34) .EQ. '') RAs_1(i__) = rNULL__
        if(ar__3(36:37) .EQ. '') DEd_1(i__) = iNULL__
        if(ar__3(38:39) .EQ. '') DEm_1(i__) = iNULL__
        if(ar__3(40:41) .EQ. '') DEs_1(i__) = iNULL__
        if(ar__3(52:56) .EQ. '') magMax_1(i__) = rNULL__
        if(ar__3(61:66) .EQ. '') magMin_1(i__) = rNULL__
        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)
        if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
        if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,4)
     +  NSV(i__),m_NSV(i__),u_NSV(i__),n_NSV(i__),RA1950h_1(i__),
     +  RA1950m_1(i__),RA1950s_1(i__),DE1950__1(i__),DE1950d_1(i__),
     +  DE1950m_1(i__),DE1950s_1(i__),u_DE1950s_1(i__),RAh_1(i__),
     +  RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),u_DEs_1(i__),VarType_1(i__),l_magMax_1(i__),
     +  magMax_1(i__),u_magMax_1(i__),f_magMin(i__),l_magMin_1(i__),
     +  magMin_1(i__),u_magMin_1(i__),n_magMin_1(i__),magCode_1(i__),
     +  Ref1_1(i__),Ref2_1(i__),Desig(i__),SpType_1(i__),f_NSV(i__),
     +  VarName_2(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 'nsv_rem.dat'	! Remarks on NSV (and Supplement) stars

C  Format for file interpretation

    5 format(I5,A75)

C  Effective file loading

      open(unit=1,status='old',file=
     +'nsv_rem.dat')
      write(6,*) '....Loading file: nsv_rem.dat'
      do i__=1,10776
        read(1,'(A80)')ar__4
        read(ar__4,5)NSV_1(i__),Text_1(i__)
c    ..............Just test output...........
        write(6,5)NSV_1(i__),Text_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'evs_cat.dat'	! Extragalactic Variable Stars. Catalogue (Vol. V)

C  Format for file interpretation

    6 format(
     +  I7,1X,A12,A1,I2,I2,F5.2,A1,I2,I2,F4.1,A1,A8,F5.2,A1,A1,F6.3,
     +  A1,1X,A1,F14.6,A1,A1,F13.8,A1,1X,I2,A1,A1,A6,1X,A5,A5,1X,A14,
     +  A12,A9,1X,A2,I4,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'evs_cat.dat')
      write(6,*) '....Loading file: evs_cat.dat'
      do i__=1,10979
        read(1,'(A160)')ar__5
        read(ar__5,6)
     +  VarNum_1(i__),Name(i__),n_Name(i__),RAh_2(i__),RAm_2(i__),
     +  RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
     +  u_DEs_2(i__),VarType_2(i__),magMax_2(i__),u_magMax_2(i__),
     +  l_magMin_2(i__),magMin_2(i__),u_magMin_2(i__),magCode_2(i__),
     +  Epoch_1(i__),u_Epoch_1(i__),l_P(i__),P(i__),u_P(i__),
     +  M_m_D_1(i__),u_M_m_D_1(i__),n_M_m_D_1(i__),SpType_2(i__),
     +  Ref1_2(i__),Ref2_2(i__),Desig1(i__),Desig2(i__),
     +  VarName_3(i__),nonMem(i__),Year_1(i__),u_Year_1(i__)
        if(ar__5(49:53) .EQ. '') magMax_2(i__) = rNULL__
        if(ar__5(56:61) .EQ. '') magMin_2(i__) = rNULL__
        if(ar__5(65:78) .EQ. '') Epoch_1(i__) = rNULL__
        if(ar__5(81:93) .EQ. '') P(i__) = rNULL__
        if(ar__5(96:97) .EQ. '') M_m_D_1(i__) = iNULL__
        if(ar__5(156:159) .EQ. '') Year_1(i__) = iNULL__
        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)
        if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
        if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
        if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
        if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
        if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
        if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
        if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c    ..............Just test output...........
        write(6,6)
     +  VarNum_1(i__),Name(i__),n_Name(i__),RAh_2(i__),RAm_2(i__),
     +  RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
     +  u_DEs_2(i__),VarType_2(i__),magMax_2(i__),u_magMax_2(i__),
     +  l_magMin_2(i__),magMin_2(i__),u_magMin_2(i__),magCode_2(i__),
     +  Epoch_1(i__),u_Epoch_1(i__),l_P(i__),P(i__),u_P(i__),
     +  M_m_D_1(i__),u_M_m_D_1(i__),n_M_m_D_1(i__),SpType_2(i__),
     +  Ref1_2(i__),Ref2_2(i__),Desig1(i__),Desig2(i__),
     +  VarName_3(i__),nonMem(i__),Year_1(i__),u_Year_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 'evs_rem.dat'	! Extragalactic Variable Stars. Remarks

C  Format for file interpretation

    7 format(A12,1X,A67)

C  Effective file loading

      open(unit=1,status='old',file=
     +'evs_rem.dat')
      write(6,*) '....Loading file: evs_rem.dat'
      do i__=1,1936
        read(1,'(A80)')ar__6
        read(ar__6,7)Name_1(i__),Text_2(i__)
c    ..............Just test output...........
        write(6,7)Name_1(i__),Text_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'sn_cat.dat'	! Extragalactic supernovae. Catalogue

C  Format for file interpretation

    8 format(
     +  A6,A1,A1,1X,A10,1X,I2,I2,F4.1,A1,I2,I2,I2,A1,A6,F5.2,A1,A1,1X,
     +  A8,1X,F5.2,A1,A1,A1,1X,F5.1,A2,1X,F5.1,A2,I2,I2,F4.1,A1,I2,I2,
     +  I2,A1,A6,1X,A16,A5,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'sn_cat.dat')
      write(6,*) '....Loading file: sn_cat.dat'
      do i__=1,984
        read(1,'(A132)')ar__7
        read(ar__7,8)
     +  SN(i__),f_SN(i__),n_SN(i__),Gal(i__),RAh_3(i__),RAm_3(i__),
     +  RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),
     +  u_DEs_3(i__),MType(i__),BTmag(i__),u_BTmag(i__),n_BTmag(i__),
     +  maxDate(i__),magMax_3(i__),u_magMax_3(i__),n_magMax(i__),
     +  f_magMax(i__),dRA(i__),oRA(i__),dDE(i__),oDE(i__),RASN_h(i__),
     +  RASN_m(i__),RASN_s(i__),DESN__(i__),DESN_d(i__),DESN_m(i__),
     +  DESN_s(i__),u_DESN_s(i__),SNType(i__),Discov(i__),Ref1_3(i__),
     +  Ref2_3(i__)
        if(ar__7(21:22) .EQ. '') RAh_3(i__) = iNULL__
        if(ar__7(23:24) .EQ. '') RAm_3(i__) = iNULL__
        if(ar__7(25:28) .EQ. '') RAs_3(i__) = rNULL__
        if(ar__7(30:31) .EQ. '') DEd_3(i__) = iNULL__
        if(ar__7(32:33) .EQ. '') DEm_3(i__) = iNULL__
        if(ar__7(34:35) .EQ. '') DEs_3(i__) = iNULL__
        if(ar__7(43:47) .EQ. '') BTmag(i__) = rNULL__
        if(ar__7(60:64) .EQ. '') magMax_3(i__) = rNULL__
        if(ar__7(69:73) .EQ. '') dRA(i__) = rNULL__
        if(ar__7(77:81) .EQ. '') dDE(i__) = rNULL__
        if(ar__7(84:85) .EQ. '') RASN_h(i__) = iNULL__
        if(ar__7(86:87) .EQ. '') RASN_m(i__) = iNULL__
        if(ar__7(88:91) .EQ. '') RASN_s(i__) = rNULL__
        if(ar__7(93:94) .EQ. '') DESN_d(i__) = iNULL__
        if(ar__7(95:96) .EQ. '') DESN_m(i__) = iNULL__
        if(ar__7(97:98) .EQ. '') DESN_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)
        if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
        if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
        if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
        if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
        if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
        if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
        if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c    ..............Just test output...........
        write(6,8)
     +  SN(i__),f_SN(i__),n_SN(i__),Gal(i__),RAh_3(i__),RAm_3(i__),
     +  RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),
     +  u_DEs_3(i__),MType(i__),BTmag(i__),u_BTmag(i__),n_BTmag(i__),
     +  maxDate(i__),magMax_3(i__),u_magMax_3(i__),n_magMax(i__),
     +  f_magMax(i__),dRA(i__),oRA(i__),dDE(i__),oDE(i__),RASN_h(i__),
     +  RASN_m(i__),RASN_s(i__),DESN__(i__),DESN_d(i__),DESN_m(i__),
     +  DESN_s(i__),u_DESN_s(i__),SNType(i__),Discov(i__),Ref1_3(i__),
     +  Ref2_3(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'sn_rem.dat'	! Extragalactic supernovae. Remarks

C  Format for file interpretation

    9 format(A6,A1,1X,A72)

C  Effective file loading

      open(unit=1,status='old',file=
     +'sn_rem.dat')
      write(6,*) '....Loading file: sn_rem.dat'
      do i__=1,277
        read(1,'(A80)')ar__8
        read(ar__8,9)SN_1(i__),n_SN_1(i__),Text_3(i__)
c    ..............Just test output...........
        write(6,9)SN_1(i__),n_SN_1(i__),Text_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

   10 format(A5,1X,A74)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,21004
        read(1,'(A80)')ar__9
        read(ar__9,10)Ref(i__),Text_4(i__)
c    ..............Just test output...........
        write(6,10)Ref(i__),Text_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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