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

Conversion of standardized ReadMe file for file /./ftp/cats/V/105/ 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/105               SKY2000 Catalog, Version 3         (Myers+ 2000)
*================================================================================
*SKY2000 Master Catalog, Version 3 (Rev. 1)
*     Myers J.R., Sande C.B., Miller A.C., Warren Jr. W.H., Tracewell D.A.
*    <Goddard Space Flight Center, Flight Dynamics Division (2000)>
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'sky2v3r1.dat'	! The SkyMap 2000 Version 3 Catalog Rev. 1

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

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

      character*19  SKY2000    (nr__) ! *Identifier based on J2000 position
      integer*4     ID         (nr__) ! Skymap number
      integer*4     HD         (nr__) ! ?Henry Draper <III/135> number
      character*1   m_HD       (nr__) ! *[1239]? HD duplicity indication
      character*1   u_HD       (nr__) ! HD identification uncertain
      integer*4     SAO        (nr__) ! ? SAO <I/131> number
      character*1   m_SAO      (nr__) ! SAO component
      character*11  DM         (nr__) ! Durchmusterung (BD <I/122>; SD <I/119>;
*                                        CD <I/114>;  CP <I/108>)
      character*1   m_DM       (nr__) ! Durchmusterung supplement letter
      character*1   u_DM       (nr__) ! [: ] DM identification uncertain
      integer*4     HR         (nr__) ! ?Harvard Revised <V/50> num. (=BS)
      character*10  WDS        (nr__) ! Washington Double Stars <I/237> number
      character*5   m_WDS      (nr__) ! WDS components
      character*1   u_WDS      (nr__) ! [: ] WDS identification uncertain
      integer*4     PPM        (nr__) ! ?Position and Proper Motion number
*                                         (<I/146>, <I/193>, <I/208>)
      character*1   u_PPM      (nr__) ! [: ] PPM identification uncertain
      integer*4     ID_merg    (nr__) ! ?Skymap num. of last skymap entry merged
*                                       with this star
      character*10  Name       (nr__) ! Star name (or AGK3 number)
      character*10  Vname      (nr__) ! Variable star name (or
*                                        doubtful variability)
      integer*4     RAh        (nr__) ! (h) Right ascension (J2000) hours
      integer*4     RAm        (nr__) ! (min) Right ascension (J2000) minutes
      real*8        RAs        (nr__) ! (s) Right ascension (J2000) seconds
      character*1   DE_        (nr__) ! Declination sign
      integer*4     DEd        (nr__) ! (deg) Declination degrees (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Declination minutes (J2000)
      real*4        DEs        (nr__) ! (arcsec) Declination seconds (J2000)
      real*4        e_pos      (nr__) ! (arcsec) Position uncertainty
      character*1   f_pos      (nr__) ! [b] Blended position flag
      integer*4     r_pos      (nr__) ! Source of position
      real*8        pmRA       (nr__) ! (s/a) Proper motion in RA (J2000)
      real*8        pmDE       (nr__) ! (arcsec/a) Proper motion in Dec (J2000)
      integer*4     r_pm       (nr__) ! ?Source of proper motion data
      real*4        RV         (nr__) ! (km/s) ?Radial velocity
      integer*4     r_RV       (nr__) ! ?Source of radial velocity data
      real*8        Plx        (nr__) ! (arcsec) ?Trigonometric parallax
      real*8        e_Plx      (nr__) ! (arcsec) ?Trigonometric parallax uncertainty
      integer*4     r_Plx      (nr__) ! ?Source of trigonometric parallax data
      real*8        GCI_X      (nr__) ! *GCI unit vector in X (J2000)
      real*8        GCI_Y      (nr__) ! *GCI unit vector in Y (J2000)
      real*8        GCI_Z      (nr__) ! *GCI unit vector in Z (J2000)
      real*4        GLON       (nr__) ! (deg) Galactic longitude (B1950)
      real*4        GLAT       (nr__) ! (deg) Galactic latitude (B1950)
      real*4        Vmag       (nr__) ! (mag) ?Observed visual magnitude (V or v)
      real*4        Vder       (nr__) ! (mag) ?Derived visual magnitude
      real*4        e_Vmag     (nr__) ! (mag) ?Derived v or observed visual magnitude
*                                        uncertainty
      character*1   f_Vmag     (nr__) ! [b] Blended visual magnitude flag
      integer*4     r_Vmag     (nr__) ! ?Source of visual magnitude
      integer*4     n_Vmag     (nr__) ! ?V magnitude derivation flag
      real*4        Bmag       (nr__) ! (mag) ?B-magnitude (observed)
      real*4        B_V        (nr__) ! (mag) ?B-V color (observed)
      real*4        e_Bmag     (nr__) ! (mag) ?B or (B-V) magnitude uncertainty
      character*1   f_Bmag     (nr__) ! [b] Blended b-magnitude flag
      integer*4     r_Bmag     (nr__) ! ?Source of b-magnitude
      real*4        Umag       (nr__) ! (mag) ?U-magnitude (observed)
      real*4        U_B        (nr__) ! (mag) ?U-B color (observed)
      real*4        e_Umag     (nr__) ! (mag) ?U or (U-B) magnitude uncertainty
      character*1   n_Umag     (nr__) ! Blended u-magnitude flag
      integer*4     r_Umag     (nr__) ! ?Source of u-magnitude
      real*4        Ptv        (nr__) ! (mag) ?Photovisual magnitude (observed)
      integer*4     r_Ptv      (nr__) ! ?Source of ptv magnitudes
      real*4        Ptg        (nr__) ! (mag) ?Photographic magnitude (observed)
      integer*4     r_Ptg      (nr__) ! ?Source of ptg magnitudes
      character*30  SpMK       (nr__) ! Morgan-Keenan (MK) spectral type
      integer*4     r_SpMK     (nr__) ! ?Source of MK spectral type data
      character*3   Sp         (nr__) ! *One-dimensional spectral class
      integer*4     r_Sp       (nr__) ! ?Source of one-dimen. spectral class
      real*8        sep        (nr__) ! (arcsec) ?Separation of brightest and second
*                                        brightest components
      real*4        Dmag       (nr__) ! (mag) ?Magnitude difference of the brightest
*                                        and second brightest components
      real*8        orbPer     (nr__) ! (yr) *?Orbital period
      integer*4     PA         (nr__) ! (deg) ?Position angle
      real*8        date       (nr__) ! (yr) ?Year of observation (AD)
      integer*4     r_dup      (nr__) ! ?Source of multiplicity data
      character*1   n_Dmag     (nr__) ! Passband of multiple star mag. dif.
      real*8        dist1      (nr__) ! (deg) ?Distance to nearest neighboring star in
*                                        the master catalog
      real*8        dist2      (nr__) ! (deg) ?Dist. to nearest neighboring master
*                                        cat. star no more than 2 mag. fainter
      integer*4     ID_A       (nr__) ! ?Skymap number of primary component
      integer*4     ID_B       (nr__) ! ?Skymap number of second component
      integer*4     ID_C       (nr__) ! ?Skymap number of third component
      real*4        magMax     (nr__) ! (mag) *?Maximum variable magnitude
      real*4        magMin     (nr__) ! (mag) *?Minimum variable magnitude
      real*4        varAmp     (nr__) ! (mag) ?Variability amplitude
      character*1   n_varAmp   (nr__) ! Passband of variability amplitude
      real*8        varPer     (nr__) ! (d) ?Period of variability
      real*8        varEp      (nr__) ! (d) ?Epoch of variability (JD-2400000)
      integer*4     varTyp     (nr__) ! ?Type of variable star
      integer*4     r_var      (nr__) ! ?Source of variability data
      real*4        mag1       (nr__) ! (mag) ?Passband #1-magnitude (observed)
      real*4        v_mag1     (nr__) ! (mag) ?v - passband #1 color
      real*4        e_mag1     (nr__) ! (mag) ?Passband #1 uncertainty in mag. or col.
      character*1   n_mag1     (nr__) ! *[RJC] Passband #1 photometric system
      character*1   p_mag1     (nr__) ! *[R] Passband #1
      integer*4     r_mag1     (nr__) ! ?Source of passband #1: mag. or color
      real*4        mag2       (nr__) ! (mag) ?Passband #2-magnitude (observed)
      real*4        v_mag2     (nr__) ! (mag) ?v - passband #2 color
      real*4        e_mag2     (nr__) ! (mag) ?Passband #2 uncertainty in mag. or col.
      character*1   n_mag2     (nr__) ! *[JEC] Passband #2 photometric system
      character*1   p_mag2     (nr__) ! *[I] Passband #2
      integer*4     r_mag2     (nr__) ! ?Source of passband #2: mag. or color
      real*4        ci1_2      (nr__) ! (mag) ?Passband #1 - passband #2 color
      character*1   f_mag1     (nr__) ! [b] Blended passband #1 mag/color flag
      character*1   f_mag2     (nr__) ! [b] Blended passband #2 mag/color flag
      real*4        mag3       (nr__) ! (mag) *?Passband #3-magnitude (observed)
      real*4        v_mag3     (nr__) ! (mag) *?v - passband #3 color
      real*4        e_mag3     (nr__) ! (mag) *?Passband #3 uncertainty in mag. or col.
      character*1   n_mag3     (nr__) ! *Passband #3 photometric system
      character*1   p_mag3     (nr__) ! *[X] Passband #3
      integer*4     r_mag3     (nr__) ! *?Source of passband #3: mag. or color
      character*1   f_mag3     (nr__) ! *[b] Blended passband #3 mag/color flag
*Note on SKY2000:
*   Field added to SKY2000 Version 1 for SKY2000 Version 2.
*   Note that there are 91 stars for which the SKY2000 name (based on the
*   J2000 position) differs from their actual position, due to corrections
*   performed on some positions; the SKY2000 names were intentionally kept
*   identical to the original SKY2000 version. The list of these 91 stars
*   can be found in the "List of stars with differences" section below.
*Note on m_HD:
*   1. indicates this the brighter component of a system with a companion
*      >= 0.3 mag.
*   2. indicates this is the fainter component of a system with a companion
*      <= 0.3 mag.
*   3. for a third component
*   9 indicates that this is a blend of two HD stars.  The HD number is the
*      lower of the two.
*Note on GCI_X:
*   cos(Dec)*cos(RA)
*Note on GCI_Y:
*   cos(Dec)*sin(RA)
*Note on GCI_Z:
*   sin(Dec)
*Note on Sp:
*   i.e. HD, AGK3, or SAO
*Note on orbPer, magMax and magMin:
*   Data for these fields have not been added to the SKY2000 Master Star
*      Catalog.
*Note on n_mag1:
*   Photometry is on the RI system. This byte indicates whether it is
*      Johnson or Russian)
*Note on p_mag1:
*   Filter used (currently, only R)
*Note on n_mag2:
*   Photometry is on the Johnson, Eggen, or Cousins system. This byte
*      indicates which.
*Note on p_mag2:
*   Filter used (currently, only I)
*Note on mag3, v-mag3, e_mag3, n_mag3, p_mag3, r_mag3 and f_mag3:
*   Fields added to SKY2000 Version 1 for SKY2000 Version 2.
*   The column p_mag3 contains X, which refers to the star trackers
*   aboard the RXTE satellite. The response curve of these trackers
*   covers the 0.4 to 1.1um range, with its peak sensitivity covering
*   the R and I bands (.65 to .9um)

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

C  Declarations for 'refs.dat'	! Reference table

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

      integer*4     RefNo      (nr__1) ! [1/99]+= Reference number
      character*4   n_RefNo    (nr__1) ! *Note on the addition of the reference
      character*73  Text       (nr__1) ! *Text of reference
*Note on n_RefNo:
*    (98) = via SKYMAP Master Star Catalog, Version 3.7 (Cat. <V/77>)
*     .0  = new sources for Version 4.0
*     .0a = new sources for Version 4.0a
*     +2  = new sources for SKY2000 Version 2
*    All other sources added for Versions 4.1 and 4.2
*Note on Text:
*   The first line contains generally the title, followed by an asterisk (*)
*   when item referenced in the paper, and the number(s) assigned to the
*   catalog in the CDS / ADC Archives within < >
*   The following lines provide the full reference.
*--------------------------------------------------------------------------------
*
*Acknowledgements;
*    It is a pleasure to thank Chris Sande and Wayne H. Warren Jr
*    for giving a copy of their files.

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

C  Loading file 'sky2v3r1.dat'	! The SkyMap 2000 Version 3 Catalog Rev. 1

C  Format for file interpretation

    1 format(
     +  8X,A19,I8,I6,A1,A1,I6,A1,A11,A1,A1,I4,A10,A5,A1,I6,A1,I8,A10,
     +  A10,I2,I2,F7.4,A1,I2,I2,F6.3,F6.4,A1,I2,F8.5,F8.4,I2,F6.1,I2,
     +  F8.5,F8.6,I2,F9.6,F9.6,F9.6,F6.2,F6.2,F6.3,F5.2,F5.3,A1,I2,I1,
     +  F6.3,F6.3,F5.3,A1,I2,F6.3,F6.3,F5.3,A1,I2,F4.1,I2,F4.1,I2,A30,
     +  I2,A3,I2,F7.3,F5.2,F7.2,I3,F7.2,I2,A1,F7.4,F7.4,I8,I8,I8,F5.2,
     +  F5.2,F5.2,A1,F8.2,F8.2,I3,I2,F6.3,F6.3,F5.3,A1,A1,I2,F6.3,
     +  F6.3,F5.3,A1,A1,I2,F6.3,A1,A1,F6.3,F6.3,F5.3,A1,A1,I2,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'sky2v3r1.dat')
      write(6,*) '....Loading file: sky2v3r1.dat'
      do i__=1,299161
        read(1,'(A520)')ar__
        read(ar__,1)
     +  SKY2000(i__),ID(i__),HD(i__),m_HD(i__),u_HD(i__),SAO(i__),
     +  m_SAO(i__),DM(i__),m_DM(i__),u_DM(i__),HR(i__),WDS(i__),
     +  m_WDS(i__),u_WDS(i__),PPM(i__),u_PPM(i__),ID_merg(i__),
     +  Name(i__),Vname(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),e_pos(i__),f_pos(i__),r_pos(i__),
     +  pmRA(i__),pmDE(i__),r_pm(i__),RV(i__),r_RV(i__),Plx(i__),
     +  e_Plx(i__),r_Plx(i__),GCI_X(i__),GCI_Y(i__),GCI_Z(i__),
     +  GLON(i__),GLAT(i__),Vmag(i__),Vder(i__),e_Vmag(i__),
     +  f_Vmag(i__),r_Vmag(i__),n_Vmag(i__),Bmag(i__),B_V(i__),
     +  e_Bmag(i__),f_Bmag(i__),r_Bmag(i__),Umag(i__),U_B(i__),
     +  e_Umag(i__),n_Umag(i__),r_Umag(i__),Ptv(i__),r_Ptv(i__),
     +  Ptg(i__),r_Ptg(i__),SpMK(i__),r_SpMK(i__),Sp(i__),r_Sp(i__),
     +  sep(i__),Dmag(i__),orbPer(i__),PA(i__),date(i__),r_dup(i__),
     +  n_Dmag(i__),dist1(i__),dist2(i__),ID_A(i__),ID_B(i__),
     +  ID_C(i__),magMax(i__),magMin(i__),varAmp(i__),n_varAmp(i__),
     +  varPer(i__),varEp(i__),varTyp(i__),r_var(i__),mag1(i__),
     +  v_mag1(i__),e_mag1(i__),n_mag1(i__),p_mag1(i__),r_mag1(i__),
     +  mag2(i__),v_mag2(i__),e_mag2(i__),n_mag2(i__),p_mag2(i__),
     +  r_mag2(i__),ci1_2(i__),f_mag1(i__),f_mag2(i__),mag3(i__),
     +  v_mag3(i__),e_mag3(i__),n_mag3(i__),p_mag3(i__),r_mag3(i__),
     +  f_mag3(i__)
        if(ar__(36:41) .EQ. '') HD(i__) = iNULL__
        if(ar__(44:49) .EQ. '') SAO(i__) = iNULL__
        if(ar__(64:67) .EQ. '') HR(i__) = iNULL__
        if(ar__(84:89) .EQ. '') PPM(i__) = iNULL__
        if(ar__(91:98) .EQ. '') ID_merg(i__) = iNULL__
        if(ar__(166:167) .EQ. '') r_pm(i__) = iNULL__
        if(ar__(168:173) .EQ. '') RV(i__) = rNULL__
        if(ar__(174:175) .EQ. '') r_RV(i__) = iNULL__
        if(ar__(176:183) .EQ. '') Plx(i__) = rNULL__
        if(ar__(184:191) .EQ. '') e_Plx(i__) = rNULL__
        if(ar__(192:193) .EQ. '') r_Plx(i__) = iNULL__
        if(ar__(233:238) .EQ. '') Vmag(i__) = rNULL__
        if(ar__(239:243) .EQ. '') Vder(i__) = rNULL__
        if(ar__(244:248) .EQ. '') e_Vmag(i__) = rNULL__
        if(ar__(250:251) .EQ. '') r_Vmag(i__) = iNULL__
        if(ar__(252:252) .EQ. '') n_Vmag(i__) = iNULL__
        if(ar__(253:258) .EQ. '') Bmag(i__) = rNULL__
        if(ar__(259:264) .EQ. '') B_V(i__) = rNULL__
        if(ar__(265:269) .EQ. '') e_Bmag(i__) = rNULL__
        if(ar__(271:272) .EQ. '') r_Bmag(i__) = iNULL__
        if(ar__(273:278) .EQ. '') Umag(i__) = rNULL__
        if(ar__(279:284) .EQ. '') U_B(i__) = rNULL__
        if(ar__(285:289) .EQ. '') e_Umag(i__) = rNULL__
        if(ar__(291:292) .EQ. '') r_Umag(i__) = iNULL__
        if(ar__(293:296) .EQ. '') Ptv(i__) = rNULL__
        if(ar__(297:298) .EQ. '') r_Ptv(i__) = iNULL__
        if(ar__(299:302) .EQ. '') Ptg(i__) = rNULL__
        if(ar__(303:304) .EQ. '') r_Ptg(i__) = iNULL__
        if(ar__(335:336) .EQ. '') r_SpMK(i__) = iNULL__
        if(ar__(340:341) .EQ. '') r_Sp(i__) = iNULL__
        if(ar__(342:348) .EQ. '') sep(i__) = rNULL__
        if(ar__(349:353) .EQ. '') Dmag(i__) = rNULL__
        if(ar__(354:360) .EQ. '') orbPer(i__) = rNULL__
        if(ar__(361:363) .EQ. '') PA(i__) = iNULL__
        if(ar__(364:370) .EQ. '') date(i__) = rNULL__
        if(ar__(371:372) .EQ. '') r_dup(i__) = iNULL__
        if(ar__(374:380) .EQ. '') dist1(i__) = rNULL__
        if(ar__(381:387) .EQ. '') dist2(i__) = rNULL__
        if(ar__(388:395) .EQ. '') ID_A(i__) = iNULL__
        if(ar__(396:403) .EQ. '') ID_B(i__) = iNULL__
        if(ar__(404:411) .EQ. '') ID_C(i__) = iNULL__
        if(ar__(412:416) .EQ. '') magMax(i__) = rNULL__
        if(ar__(417:421) .EQ. '') magMin(i__) = rNULL__
        if(ar__(422:426) .EQ. '') varAmp(i__) = rNULL__
        if(ar__(428:435) .EQ. '') varPer(i__) = rNULL__
        if(ar__(436:443) .EQ. '') varEp(i__) = rNULL__
        if(ar__(444:446) .EQ. '') varTyp(i__) = iNULL__
        if(ar__(447:448) .EQ. '') r_var(i__) = iNULL__
        if(ar__(449:454) .EQ. '') mag1(i__) = rNULL__
        if(ar__(455:460) .EQ. '') v_mag1(i__) = rNULL__
        if(ar__(461:465) .EQ. '') e_mag1(i__) = rNULL__
        if(ar__(468:469) .EQ. '') r_mag1(i__) = iNULL__
        if(ar__(470:475) .EQ. '') mag2(i__) = rNULL__
        if(ar__(476:481) .EQ. '') v_mag2(i__) = rNULL__
        if(ar__(482:486) .EQ. '') e_mag2(i__) = rNULL__
        if(ar__(489:490) .EQ. '') r_mag2(i__) = iNULL__
        if(ar__(491:496) .EQ. '') ci1_2(i__) = rNULL__
        if(ar__(499:504) .EQ. '') mag3(i__) = rNULL__
        if(ar__(505:510) .EQ. '') v_mag3(i__) = rNULL__
        if(ar__(511:515) .EQ. '') e_mag3(i__) = rNULL__
        if(ar__(518:519) .EQ. '') r_mag3(i__) = iNULL__
        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)
     +  SKY2000(i__),ID(i__),HD(i__),m_HD(i__),u_HD(i__),SAO(i__),
     +  m_SAO(i__),DM(i__),m_DM(i__),u_DM(i__),HR(i__),WDS(i__),
     +  m_WDS(i__),u_WDS(i__),PPM(i__),u_PPM(i__),ID_merg(i__),
     +  Name(i__),Vname(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),e_pos(i__),f_pos(i__),r_pos(i__),
     +  pmRA(i__),pmDE(i__),r_pm(i__),RV(i__),r_RV(i__),Plx(i__),
     +  e_Plx(i__),r_Plx(i__),GCI_X(i__),GCI_Y(i__),GCI_Z(i__),
     +  GLON(i__),GLAT(i__),Vmag(i__),Vder(i__),e_Vmag(i__),
     +  f_Vmag(i__),r_Vmag(i__),n_Vmag(i__),Bmag(i__),B_V(i__),
     +  e_Bmag(i__),f_Bmag(i__),r_Bmag(i__),Umag(i__),U_B(i__),
     +  e_Umag(i__),n_Umag(i__),r_Umag(i__),Ptv(i__),r_Ptv(i__),
     +  Ptg(i__),r_Ptg(i__),SpMK(i__),r_SpMK(i__),Sp(i__),r_Sp(i__),
     +  sep(i__),Dmag(i__),orbPer(i__),PA(i__),date(i__),r_dup(i__),
     +  n_Dmag(i__),dist1(i__),dist2(i__),ID_A(i__),ID_B(i__),
     +  ID_C(i__),magMax(i__),magMin(i__),varAmp(i__),n_varAmp(i__),
     +  varPer(i__),varEp(i__),varTyp(i__),r_var(i__),mag1(i__),
     +  v_mag1(i__),e_mag1(i__),n_mag1(i__),p_mag1(i__),r_mag1(i__),
     +  mag2(i__),v_mag2(i__),e_mag2(i__),n_mag2(i__),p_mag2(i__),
     +  r_mag2(i__),ci1_2(i__),f_mag1(i__),f_mag2(i__),mag3(i__),
     +  v_mag3(i__),e_mag3(i__),n_mag3(i__),p_mag3(i__),r_mag3(i__),
     +  f_mag3(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 'refs.dat'	! Reference table

C  Format for file interpretation

    2 format(1X,I2,A4,A73)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,190
        read(1,'(A80)')ar__1
        read(ar__1,2)RefNo(i__),n_RefNo(i__),Text(i__)
c    ..............Just test output...........
        write(6,2)RefNo(i__),n_RefNo(i__),Text(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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