FORTRAN Generation
(/./ftp/cats/I/334)

Conversion of standardized ReadMe file for file /./ftp/cats/I/334 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 2017-Dec-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. I/334             W1J00 and W2J00 Transit Circle Catalogs      (Rafferty+, 2016)
*================================================================================
*W1J00 results of the Observations made with the Six-inch Transit
*Circle 1977-1982 AND W2J00 results of Pole-to-Pole Observations made with the
*Six-inch and Seven-inch Transit Circles 1985-1996.
*     Rafferty T.J., Holdenried E.R., Urban S.E.
*    <Publ. USNO, 2nd series, vol. XXVII (part 1), Washington (2016)>
*    Holdenried, E.R. and Rafferty, T.J.
*    <Publ. USNO, 2nd series, vol. XXVII (part 2), Washington (2016)>
*   =2016yCat.1334....0R
*   =2016PUSNO..27a...1R
*   +2016PUSNO..27b...1H
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'w1j00sta.dat'	! W1J00 star positions (means) (W1J00_stars.dat)

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

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

      integer*4     W1J00      (nr__) ! W1J00 identifier
      integer*4     RAh        (nr__) ! (h) Right Ascension at epoch EpRA (ICRS) (1)
      integer*4     RAm        (nr__) ! (min) Right Ascension at epoch EpRA (ICRS) (1)
      real*4        RAs        (nr__) ! (s) Right Ascension at epoch EpRA (ICRS) (1)
      integer*4     e_RAs      (nr__) ! (mas) Standard error of the mean, RA
      integer*4     NumRA      (nr__) ! Number of observations used for RA
      real*8        EpRA       (nr__) ! (yr) Mean epoch of RA
      real*4        dRA        (nr__) ! (s) RA difference, Absolute - ICRS (2)
      character*1   DE_        (nr__) ! Declination sign at epoch=EpDE (ICRS) (1)
      integer*4     DEd        (nr__) ! (deg) Declination at epoch EpDE (ICRS) (1)
      integer*4     DEm        (nr__) ! (arcmin) Declination at epoch EpDE (ICRS) (1)
      real*4        DEs        (nr__) ! (arcsec) Declination at epoch EpDE (ICRS) (1)
      integer*4     e_DEs      (nr__) ! (mas) Standard error of the mean, DE
      integer*4     NumDE      (nr__) ! Number of observations used for DE
      real*8        EpDE       (nr__) ! (yr) Mean epoch of DE
      real*4        dDE        (nr__) ! (arcsec) DE difference, Absolute - ICRS (2)
      real*4        Vmag       (nr__) ! (mag) Visual magnitude (3)
      integer*4     HIC        (nr__) ! ? Hipparcos Input Catalog (Cat. I/196)
*                                     number (4)
      character*2   Comp       (nr__) ! Hipparcos Input Catalog component (4)
      integer*4     DM         (nr__) ! ? Durchmusterung (Cat. I/122) number (4)
      integer*4     FK5        (nr__) ! ? FK5 (Cats. I/149 and I/175) number (4)
      character*9   WDS        (nr__) ! Washington Double Star Catalog (Cat. B/wds)
*                                   identifier (4)
*Note (1): Mean positions are on the International Celestial Reference System
*  (ICRS), differentially reduced using the Hipparcos catalog.  The epochs for
*  the right ascension and declination coordinates are found in the
*  EpRA and EpDE fields respectively.
*Note (2): The observing program was designed to be an "absolute" catalog,
*  independent of other catalogs.  Data reductions were carried out according to
*  previous absolute transit circle catalogs. Following the reductions,
*  comparisons with the Hipparcos Catalogue revealed remaining, unaccounted for
*  systematic deviations in the absolute positions. The authors decided to
*  differentially reduce the absolute positions using the Hipparcos data.
*  Those differential positions are listed in the RAh, RAm, RAs, DE-, DEd, DEm
*  and DEs fields. For users who may be interested in the absolute positions,
*  values of the absolute position minus the differential position are provided
*  in the dRA and dDE fields.
*Note (3): Visual magnitude for identification purposes. This was drawn from a
*  variety of sources, but largely from the Hipparcos Input Catalogue.
*Note (4): Cross references to several widely known catalogs are provided for
*   convenience of the user.

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

C  Declarations for 'w1j00sol.dat'	! W1J00 solar system object positions
                                   (W1J00_solsys.dat)

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

      character*7   Obj        (nr__1) ! Solar system object identifier
      character*12  RAa        (nr__1) ! Apparent Right Ascencion at Epoch (1)
      real*4        dRA_1      (nr__1) ! (s) ? RA difference, Absolute - HIP (3)
      character*12  DEa        (nr__1) ! Apparent Declination at Epoch (1)
      real*4        dDE_1      (nr__1) ! (arcsec) ? DE difference, Absolute - HIP  (3)
      real*8        Epoch      (nr__1) ! (d) Julian date (UT1) of observation
      character*2   Obs        (nr__1) ! Observer code (4)
      character*1   Clamp      (nr__1) ! Clamp orientation of instrument (5)
      character*1   RALimb     (nr__1) ! Limbs or center of light measured in RA (6)
      character*1   DELimb     (nr__1) ! Limbs or center of light measured in DE (6)
*Note (1): Positions are apparent places, reduced to be systematically consistent
*  with the Hipparcos Catalogue.  An apparent place is a geocentric direction of
*  an object that takes into account orbital motion, space motion, light-time,
*  light deflection, and annual aberration. Apparent place is given with respect
*  to the true equator and equinox "of date"; in this case, it is the date of
*  observation found in the Epoch field.
* In some cases, an observation is made in only right ascension or declination.
*  If only the right ascension coordinate was observed, then the declination
*  seconds field (DEs) will be blank, as will dDE and DELimb.
*  If only the declination coordinate was observed, then the right ascension
*  seconds field (RAs) will be blank, as will dRA and RALimb.
*Note (3): The observing program was designed to be an "absolute" catalog,
*  independent of other catalogs.  Data reductions were carried out according to
*  previous absolute transit circle catalogs. Following the reductions,
*  comparisons of stars with the Hipparcos Catalogue revealed remaining,
*  unaccounted for systematic deviations in the absolute positions. The authors
*  decided to differentially reduce the absolute positions using the Hipparcos
*  data. Those differential positions are listed in the RAh, RAm, RAs, DE-, DEd,
*  DEm and DEs fields. For users who may be interested in the absolute positions,
*  values of the absolute position minus the differential position are provided
*  in the dRA and dDE fields.
*Note (4): The Six-inch Transit Circle measurements were manual, in the
*    sense that a person measured the position while looking through the
*    instrument. The observer code has the following meaning:
*       RB = R.L. Branham
*       TC = T.E. Corbin
*       SD = S.J. Dick
*       SG = F.S. Gauss
*       EH = E.R. Holdenried
*       EJ = E.S. Jackson
*       BK = B.L. Klock
*       RJ = R.J. Miller
*       DR = M.D. Robinson
*       RR = R.W. Rhynsburger
*       TR = T.J. Rafferty
*       CS = C.A. Smith
*       RS = R.C. Stone
*       ZY = Z.G. Yao
*Note (5): The Six-inch was reversed (rotated 180 degrees in azimuth)
*  interchanging the east and west pivots approximately every 30 days.
*  The orientation of the instrument was referenced to the location of the
*  clamping device, which fixed the altitude of the instrument after it was
*  pointed to a star, and thus referred to either as Clamp West or Clamp East.
*  The flag has the following meaning:
*       W = Clamp West
*       E = Clamp East
*Note (6):  The Sun and planets subtend sizable disks in the
*    instrument.  Primarily depending on the object, either the limbs
*    or the center of light were observed. The flag has the following
*    meaning:
*       L = limb was measured
*       C = center of light was measured
*       R = ring was measured (Saturn only)
*       P = preceding limb was measured (Venus only)
*       F = following limb was measured (Venus only)
*       N = North limb was measured (Venus only)
*       S = South limb was measured (Venus only)
*    Despite how an object was observed, the W1J00 positions refer to
*    the center of the object.

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

C  Declarations for 'w2j00sta.dat'	! W2J00 star positions (means) (W2J00_stars.dat)

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

C  EpRA (ICRS) (1) position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1    (nr__2) ! (deg) Right Ascension EpRA (ICRS) (1)
      real*8        DEdeg_1    (nr__2) ! (deg)     Declination EpRA (ICRS) (1)
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     W2J00      (nr__2) ! W2J00 identifier
      integer*4     RAh_1      (nr__2) ! (h) Right Ascension at epoch EpRA (ICRS) (1)
      integer*4     RAm_1      (nr__2) ! (min) Right Ascension at epoch EpRA (ICRS) (1)
      real*4        RAs_1      (nr__2) ! (s) Right Ascension at epoch EpRA (ICRS) (1)
      integer*4     e_RAs_1    (nr__2) ! (mas) Weighted standard error of the mean, RA
      integer*4     NumRA_1    (nr__2) ! Number of observations used for RA
      real*8        EpRA_1     (nr__2) ! (yr) Mean epoch of RA
      character*1   DE__1      (nr__2) ! Declination sign at epoch EpDE (ICRS) (1)
      integer*4     DEd_1      (nr__2) ! (deg) Declination at epoch EpDE (ICRS) (1)
      integer*4     DEm_1      (nr__2) ! (arcmin) Declination at epoch EpDE (ICRS) (1)
      real*4        DEs_1      (nr__2) ! (arcsec) Declination at epoch EpDE (ICRS) (1)
      integer*4     e_DEs_1    (nr__2) ! (mas) Weighted standard error of the mean, DE
      integer*4     NumDE_1    (nr__2) ! Number of observations used for DE
      real*8        EpDE_1     (nr__2) ! (yr) Mean epoch of DE
      real*4        Vmag_1     (nr__2) ! (mag) Visual magnitude (2)
      integer*4     HIP        (nr__2) ! ? Hipparcos Catalogue (Cat. I/239) number (3)
      integer*4     RA6h       (nr__2) ! (h) ? Right ascension from 6-inch,
*                                   at epoch EpRA6 (ICRS) (4)
      integer*4     RA6m       (nr__2) ! (min) ? Right ascension from 6-inch,
*                                   at epoch EpRA6 (ICRS) (4)
      real*4        RA6s       (nr__2) ! (s) ? Right ascension from 6-inch,
*                                   at epoch EpRA6 (ICRS) (4)
      integer*4     e_RA6s     (nr__2) ! (mas) ? Weighted standard error of the mean, RA6
      integer*4     NumRA6     (nr__2) ! ? Number of observations used for RA6
      real*8        EpRA6      (nr__2) ! (yr) ? Mean epoch of RA6
      character*1   DE6_       (nr__2) ! Declination sign from 6-inch,
*                                 at epoch EpDE6 (ICRS) (4)
      integer*4     DE6d       (nr__2) ! (deg) ? Declination from 6-inch,
*                                 at epoch EpDE6 (ICRS) (4)
      integer*4     DE6m       (nr__2) ! (arcmin) ? Declination from 6-inch,
*                                 at epoch EpDE6 (ICRS) (4)
      real*4        DE6s       (nr__2) ! (arcsec) ? Declination from 6-inch,
*                                 at epoch EpDE6 (ICRS) (4)
      integer*4     e_DE6s     (nr__2) ! (mas) ? Weighted standard error of the mean, DE6
      integer*4     NumDE6     (nr__2) ! ? Number of observations used for DE6
      real*8        EpDE6      (nr__2) ! (yr) ? Mean epoch of DE6
      integer*4     RA7h       (nr__2) ! (h) ? Right ascension from 7-inch,
*                                 at epoch EpRA7 (ICRS) (5)
      integer*4     RA7m       (nr__2) ! (min) ? Right ascension from 7-inch,
*                                 at epoch EpRA7 (ICRS) (5)
      real*4        RA7s       (nr__2) ! (s) ? Right ascension from 7-inch,
*                                 at epoch EpRA7 (ICRS) (5)
      integer*4     e_RA7s     (nr__2) ! (mas) ? Weighted standard error of the mean, RA7
      integer*4     NumRA7     (nr__2) ! ? Number of observations used for RA7
      real*8        EpRA7      (nr__2) ! (yr) ? Mean epoch of RA7
      character*1   DE7_       (nr__2) ! Declination sign from 7-inch,
*                                 at epoch EpDE7 (ICRS) (5)
      integer*4     DE7d       (nr__2) ! (deg) ? Declination from 7-inch,
*                                  at epoch EpDE7 (ICRS) (5)
      integer*4     DE7m       (nr__2) ! (arcmin) ? Declination from 7-inch,
*                                  at epoch EpDE7 (ICRS) (5)
      real*4        DE7s       (nr__2) ! (arcsec) ? Declination from 7-inch,
*                                  at epoch EpDE7 (ICRS) (5)
      integer*4     e_DE7s     (nr__2) ! (mas) ? Weighted standard error of the mean, DE7
      integer*4     NumDE7     (nr__2) ! ? Number of observations used for DE7
      real*8        EpDE7      (nr__2) ! (yr) ? Mean epoch of DE7
*Note (1): Mean positions are on the International Celestial Reference System
*  (ICRS), differentially reduced using the Hipparcos Catalogue. The positions
*  are a weighted mean of the Six-inch Transit Circle and Seven-inch Transit
*  Circle observations, found in columns 79-197. The epochs for the right
*  ascension and declination coordinates are found in the EpRA and EpDE fields
*  respectively.
*Note (2): Visual magnitude for identification purposes.
*Note (3): Cross references to the Hipparcos Catalogue are provided for
*  convenience of the user.
*Note (4): Columns 1-65 give positions that are a weighted mean of the Six-inch
*  Transit Circle and Seven-inch Transit Circle observations.
*  Columns 79-137 give the positions from the Six-inch Transit Circle only.
*Note (5): Columns 1-65 give positions that are a weighted mean of the Six-inch
*  Transit Circle and Seven-inch Transit Circle observations.
*  Columns 139-197 give the positions from the Seven-inch Transit Circle only.

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

C  Declarations for 'w2j00sol.dat'	! W2J00 solar system object positions
                                   (W2J00_solsys.dat)

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

      character*10  Obj_1      (nr__3) ! Solar system object identifier
      character*12  RAa_1      (nr__3) ! Apparent Right Ascencion at Epoch (1)
      character*12  DEa_1      (nr__3) ! Apparent Declination at Epoch (1)
      real*8        Epoch_1    (nr__3) ! (d) Julian date (UT1) of observation
      character*3   Obs_1      (nr__3) ! Observer (3)
      character*1   Clamp_1    (nr__3) ! Clamp orientation of instrument (4)
      character*1   RALimb_1   (nr__3) ! Limbs or center of light measured in RA (5)
      character*1   DELimb_1   (nr__3) ! Limbs or center of light measured in DE (5)
      character*2   Tel        (nr__3) ! Telescope used (6)
*Note (1): Positions are apparent places, reduced to be systematically consistent
*  with the Hipparcos Catalogue.  An apparent place is a geocentric direction of
*  an object that takes into account orbital motion, space motion, light-time,
*  light deflection, and annual aberration. Apparent place is given with respect
*  to the true equator and equinox "of date"; in this case, it is the date of
*  observation found in the Epoch field.
* In some cases, an observation is made in only right ascension or declination.
*  If only the right ascension coordinate was observed, then the declination
*  seconds field (DEs) will be blank, as will the DELimb.
*  If only the declination coordinate was observed, then the right ascension
*  seconds field (RAs) will be blank, as will the RALimb.
*Note (3): The Six-inch Transit Circle measurements were manual, in the sense
*  that a person measured the position while looking through the instrument.
*  The Seven-inch Transit Circle's observations were automated, so the observer
*  is not listed.  The observer code has the following meaning:
*       CSC = C.S. Cole
*       TEC = T.E. Corbin
*       HEC = H.E. Crull
*       SD  = S.J. Dick
*       JCD = J.C. Doty
*       RE  = R.   Etheridge
*       FSG = F.S. Gauss
*       MEG = M.E. Germain
*       DMH = D.M. Hall
*       GSH = G.S. Hennessy
*       JH  = J.L. Hershey
*       RBH = R.B. Hindsely
*       ERH = E.R. Holdenried
*       ESJ = E.S. Jackson
*       IAN = I.   Jordan
*       KAL = V.   Kallarakal
*       JCM = J.C. Martin
*       JMM = J.M. Muse
*       RJM = R.J. Miller
*       MDR = M.D. Robinson
*       TJR = T.J. Rafferty
*       CBS = C.B. Sande
*       DS  = D.K. Scott
*       CAS = C.A. Smith
*       SEU = S.E. Urban
*       GLW = G.L. Wycoff
*       ZGY = Z.G. Yao
*Note (4): Both transit circle telescopes were reversed (rotated 180 degrees in
*  azimuth) interchanging the east and west pivots approximately every 30 days.
*  The orientation of the instrument was referenced to the location of the
*  clamping device, which fixed the altitude of the instrument after it was
*  pointed to a star, and thus referred to either as Clamp West or Clamp East.
*  The flag has the following meaning:
*       W = Clamp West
*       E = Clamp East
*Note (5): Some planets subtend sizable disks in the instrument. Primarily
*  depending on the object, either the limbs or the center-of-light were
*  observed. The flag has the following meaning:
*       L = limb was measured
*       C = center of light was measured
*       R = ring was measured (Saturn only)
*  Despite how an object was observed, the W2J00 positions refer to the center
*   of the object.
*Note (6): The observing program used two telescopes, the Six-inch Transit Circle
*  in Washington DC, and the Seven-inch Transit Circle at Black Birch New
*  Zealand. The "Tel" flag has the following meaning:
*       6" = Six-inch Transit Circle, Washington D.C., USA
*       7" = Seven-inch Transit Circle, Black Birch, New Zealand

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

C  Loading file 'w1j00sta.dat'	! W1J00 star positions (means) (W1J00_stars.dat)

C  Format for file interpretation

    1 format(
     +  I4,1X,I2,1X,I2,1X,F6.3,1X,I3,1X,I3,1X,F8.3,1X,F6.3,1X,A1,I2,
     +  1X,I2,1X,F5.2,1X,I3,1X,I3,1X,F8.3,1X,F5.2,1X,F6.3,1X,I6,A2,1X,
     +  I8,1X,I4,1X,A9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'w1j00sta.dat')
      write(6,*) '....Loading file: w1j00sta.dat'
      do i__=1,7267
        read(1,'(A117)')ar__
        read(ar__,1)
     +  W1J00(i__),RAh(i__),RAm(i__),RAs(i__),e_RAs(i__),NumRA(i__),
     +  EpRA(i__),dRA(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
     +  e_DEs(i__),NumDE(i__),EpDE(i__),dDE(i__),Vmag(i__),HIC(i__),
     +  Comp(i__),DM(i__),FK5(i__),WDS(i__)
        if(ar__(86:91) .EQ. '') HIC(i__) = iNULL__
        if(ar__(95:102) .EQ. '') DM(i__) = iNULL__
        if(ar__(104:107) .EQ. '') FK5(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)
     +  W1J00(i__),RAh(i__),RAm(i__),RAs(i__),e_RAs(i__),NumRA(i__),
     +  EpRA(i__),dRA(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
     +  e_DEs(i__),NumDE(i__),EpDE(i__),dDE(i__),Vmag(i__),HIC(i__),
     +  Comp(i__),DM(i__),FK5(i__),WDS(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 'w1j00son.dat&	) W1J00 solar system object positions
*                                   (W1J00_solsys.dat)

C  Format for file interpretation

    2 format(
     +  A7,1X,A12,1X,F6.3,1X,A12,1X,F5.2,1X,F13.5,1X,A2,1X,A1,1X,A1,
     +  1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'w1j00sol.dat')
      write(6,*) '....Loading file: w1j00sol.dat'
      do i__=1,4383
        read(1,'(A69)')ar__1
        read(ar__1,2)
     +  Obj(i__),RAa(i__),dRA_1(i__),DEa(i__),dDE_1(i__),Epoch(i__),
     +  Obs(i__),Clamp(i__),RALimb(i__),DELimb(i__)
        if(ar__1(22:27) .EQ. '') dRA_1(i__) = rNULL__
        if(ar__1(42:46) .EQ. '') dDE_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Obj(i__),RAa(i__),dRA_1(i__),DEa(i__),dDE_1(i__),Epoch(i__),
     +  Obs(i__),Clamp(i__),RALimb(i__),DELimb(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'w2j00sta.dat'	! W2J00 star positions (means) (W2J00_stars.dat)

C  Format for file interpretation

    3 format(
     +  I5,1X,I2,1X,I2,1X,F6.3,1X,I3,1X,I3,1X,F8.3,1X,A1,I2,1X,I2,1X,
     +  F5.2,1X,I3,1X,I3,1X,F8.3,1X,F4.1,1X,I6,1X,I2,1X,I2,1X,F6.3,1X,
     +  I3,1X,I3,1X,F8.3,1X,A1,I2,1X,I2,1X,F5.2,1X,I3,1X,I3,1X,F8.3,
     +  1X,I2,1X,I2,1X,F6.3,1X,I3,1X,I3,1X,F8.3,1X,A1,I2,1X,I2,1X,
     +  F5.2,1X,I3,1X,I3,1X,F8.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'w2j00sta.dat')
      write(6,*) '....Loading file: w2j00sta.dat'
      do i__=1,44395
        read(1,'(A197)')ar__2
        read(ar__2,3)
     +  W2J00(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),e_RAs_1(i__),
     +  NumRA_1(i__),EpRA_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),e_DEs_1(i__),NumDE_1(i__),EpDE_1(i__),Vmag_1(i__),
     +  HIP(i__),RA6h(i__),RA6m(i__),RA6s(i__),e_RA6s(i__),
     +  NumRA6(i__),EpRA6(i__),DE6_(i__),DE6d(i__),DE6m(i__),
     +  DE6s(i__),e_DE6s(i__),NumDE6(i__),EpDE6(i__),RA7h(i__),
     +  RA7m(i__),RA7s(i__),e_RA7s(i__),NumRA7(i__),EpRA7(i__),
     +  DE7_(i__),DE7d(i__),DE7m(i__),DE7s(i__),e_DE7s(i__),
     +  NumDE7(i__),EpDE7(i__)
        if(ar__2(72:77) .EQ. '') HIP(i__) = iNULL__
        if(ar__2(79:80) .EQ. '') RA6h(i__) = iNULL__
        if(ar__2(82:83) .EQ. '') RA6m(i__) = iNULL__
        if(ar__2(85:90) .EQ. '') RA6s(i__) = rNULL__
        if(ar__2(92:94) .EQ. '') e_RA6s(i__) = iNULL__
        if(ar__2(96:98) .EQ. '') NumRA6(i__) = iNULL__
        if(ar__2(100:107) .EQ. '') EpRA6(i__) = rNULL__
        if(ar__2(110:111) .EQ. '') DE6d(i__) = iNULL__
        if(ar__2(113:114) .EQ. '') DE6m(i__) = iNULL__
        if(ar__2(116:120) .EQ. '') DE6s(i__) = rNULL__
        if(ar__2(122:124) .EQ. '') e_DE6s(i__) = iNULL__
        if(ar__2(126:128) .EQ. '') NumDE6(i__) = iNULL__
        if(ar__2(130:137) .EQ. '') EpDE6(i__) = rNULL__
        if(ar__2(139:140) .EQ. '') RA7h(i__) = iNULL__
        if(ar__2(142:143) .EQ. '') RA7m(i__) = iNULL__
        if(ar__2(145:150) .EQ. '') RA7s(i__) = rNULL__
        if(ar__2(152:154) .EQ. '') e_RA7s(i__) = iNULL__
        if(ar__2(156:158) .EQ. '') NumRA7(i__) = iNULL__
        if(ar__2(160:167) .EQ. '') EpRA7(i__) = rNULL__
        if(ar__2(170:171) .EQ. '') DE7d(i__) = iNULL__
        if(ar__2(173:174) .EQ. '') DE7m(i__) = iNULL__
        if(ar__2(176:180) .EQ. '') DE7s(i__) = rNULL__
        if(ar__2(182:184) .EQ. '') e_DE7s(i__) = iNULL__
        if(ar__2(186:188) .EQ. '') NumDE7(i__) = iNULL__
        if(ar__2(190:197) .EQ. '') EpDE7(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.Edeg1(h__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,3)
     +  W2J00(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),e_RAs_1(i__),
     +  NumRA_1(i__),EpRA_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),e_DEs_1(i__),NumDE_1(i__),EpDE_1(i__),Vmag_1(i__),
     +  HIP(i__),RA6h(i__),RA6m(i__),RA6s(i__),e_RA6s(i__),
     +  NumRA6(i__),EpRA6(i__),DE6_(i__),DE6d(i__),DE6m(i__),
     +  DE6s(i__),e_DE6s(i__),NumDE6(i__),EpDE6(i__),RA7h(i__),
     +  RA7m(i__),RA7s(i__),e_RA7s(i__),NumRA7(i__),EpRA7(i__),
     +  DE7_(i__),DE7d(i__),DE7m(i__),DE7s(i__),e_DE7s(i__),
     +  NumDE7(i__),EpDE7(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 'w2j00sol.dat'	! W2J00 solar system object positions
*                                   (W2J00_solsys.dat)

C  Format for file interpretation

    4 format(A10,1X,A12,1X,A12,1X,F13.5,1X,A3,1X,A1,1X,A1,1X,A1,1X,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'w2j00sol.dat')
      write(6,*) '....Loading file: w2j00sol.dat'
      do i__=1,11566
        read(1,'(A63)')ar__3
        read(ar__3,4)
     +  Obj_1(i__),RAa_1(i__),DEa_1(i__),Epoch_1(i__),Obs_1(i__),
     +  Clamp_1(i__),RALimb_1(i__),DELimb_1(i__),Tel(i__)
c    ..............Just test output...........
        write(6,4)
     +  Obj_1(i__),RAa_1(i__),DEa_1(i__),Epoch_1(i__),Obs_1(i__),
     +  Clamp_1(i__),RALimb_1(i__),DELimb_1(i__),Tel(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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