FORTRAN Generation
(/./ftp/cats/J/ApJ/837/30)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/837/30 into FORTRAN code for reading data files line by line.

Note that special values are assigned to unknown or unspecified numbers (also called NULL numbers); when necessary, the coordinate components making up the right ascension and declination are converted into floating-point numbers representing these angles in degrees.



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-29
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. J/ApJ/837/30   25yrs monitoring of stellar orbits in the GC   (Gillessen+, 2017)
*================================================================================
*An update on monitoring stellar orbits in the Galactic Center.
*    Gillessen S., Plewa P.M., Eisenhauer F., Sari R., Waisberg I., Habibi M.,
*    Pfuhl O., George E., Dexter J., von Fellenberg S., Ott T., Genzel R.
*   <Astrophys. J., 837, 30-30 (2017)>
*   =2017ApJ...837...30G    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Orbital parameters of the 40 stars for which
                                 we were able to determine orbits

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

      character*4   Star        ! Star Name
      character*1   f_Star      ! [a] Flag on Object S111 (1)
      real*8        a           ! (arcsec) [-12.3/4.6] Semi-major axis of the orbit,
*                                 from Sag A*
      real*4        e_a         ! (arcsec) [0.0002/8.4] Uncertainty in a
      real*4        e           ! [0.1/1.1] Eccentricity
      real*4        e_e         ! [0.0003/0.3] Uncertainty in e
      real*4        i           ! (deg) [24.7/171.1] Inclination
      real*4        e_i         ! (deg) [0.06/8.3] Uncertainty in i
      real*4        Omega       ! (deg) [7.9/345] Position angle of the ascending node
      real*4        e_Omega     ! (deg) [0.07/19] Uncertainty in {Omega}
      real*4        w           ! (deg) [17.9/357] Longitude of the pericenter
      real*4        e_w         ! (deg) [0.07/24] Uncertainty in w
      real*8        Tp          ! (yr) [611/2132] Epoch of pericenter passage,
*                                fractional year (G1)
      real*4        e_Tp        ! (yr) [0.01/154] Uncertainty in Tp
      real*8        Per         ! (yr) [12.8/3580]? Orbital period (G1)
      real*8        e_Per       ! (yr) [0.02/2550]? Uncertainty in Per
      character*1   SpT         ! [el] Spectral type: e=early-type stars,
*                                l=late-type stars
      real*4        Kmag        ! (mag) [10/18] K-band magnitude
      real*4        r           ! [0.9/3.4] Global rescaling factor
*Note (1):
*   a = Object S111 formally has a negative semi major axis, indicative for
*       a hyperbolic orbit with e>1 and no solution for an orbital period.

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

C  Declarations for 'table5.dat'	! *Orbital data used in this work

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

      real*8        Date        ! (yr) Observation date, fractional year (G1)
      character*2   Flag        ! [a rv] Observation flag (a=astrometric data;
*                                rv=radial velocity data)
      real*4        oRA_S1      ! (mas) ? Offset in Right Ascension, S1 relative to SgrA*
      real*4        e_oRA_S1    ! (mas) ? Uncertainty in oRA-S1
      real*4        oDE_S1      ! (mas) ? Offset in Declination, S1 relative to SgrA*
      real*4        e_oDE_S1    ! (mas) ? Uncertainty in oDE-S1
      integer*4     RV_S1       ! (km/s) [-1192/-917]? Radial Velocity for S1
      integer*4     e_RV_S1     ! (km/s) [35/175]? Uncertainty in RVel-S1
      real*4        oRA_S2      ! (mas) ? Offset in Right Ascension, S2 relative to SgrA*
      real*4        e_oRA_S2    ! (mas) ? Uncertainty in oRA-S2
      real*4        oDE_S2      ! (mas) ? Offset in Declination, S2 relative to SgrA*
      real*4        e_oDE_S2    ! (mas) ? Uncertainty in oDE-S2
      integer*4     RV_S2       ! (km/s) [-1571/1199]? Radial Velocity for S2
      integer*4     e_RV_S2     ! (km/s) [17/113]? Uncertainty in RVel-S2
      real*4        oRA_S4      ! (mas) ? Offset in Right Ascension, S4 relative to SgrA*
      real*4        e_oRA_S4    ! (mas) ? Uncertainty in oRA-S4
      real*4        oDE_S4      ! (mas) ? Offset in Declination, S4 relative to SgrA*
      real*4        e_oDE_S4    ! (mas) ? Uncertainty in oDE-S4
      integer*4     RV_S4       ! (km/s) [-1015/-525]? Radial Velocity for S4
      integer*4     e_RV_S4     ! (km/s) [25/125]? Uncertainty in RVel-S4
      real*4        oRA_S8      ! (mas) ? Offset in Right Ascension, S8 relative to SgrA*
      real*4        e_oRA_S8    ! (mas) ? Uncertainty in oRA-S8
      real*4        oDE_S8      ! (mas) ? Offset in Declination, S8 relative to SgrA*
      real*4        e_oDE_S8    ! (mas) ? Uncertainty in oDE-S8
      integer*4     RV_S8       ! (km/s) [-306/26]? Radial Velocity for S8
      integer*4     e_RV_S8     ! (km/s) [30/118]? Uncertainty in RVel-S8
      real*4        oRA_S9      ! (mas) ? Offset in Right Ascension, S9 relative to SgrA*
      real*4        e_oRA_S9    ! (mas) ? Uncertainty in oRA-S9
      real*4        oDE_S9      ! (mas) ? Offset in Declination, S9 relative to SgrA*
      real*4        e_oDE_S9    ! (mas) ? Uncertainty in oDE-S9
      integer*4     RV_S9       ! (km/s) [390/805]? Radial Velocity for S9
      integer*4     e_RV_S9     ! (km/s) [41/165]? Uncertainty in RVel-S9
      real*4        oRA_S12     ! (mas) ? Offset in RA, S12 relative to SgrA*
      real*4        e_oRA_S12   ! (mas) ? Uncertainty in oRA-S12
      real*4        oDE_S12     ! (mas) ? Offset in Declination, S12 relative to SgrA*
      real*4        e_oDE_S12   ! (mas) ? Uncertainty in oDE-S12
      integer*4     RV_S12      ! (km/s) [-90/349]? Radial Velocity for S12
      integer*4     e_RV_S12    ! (km/s) [71/272]? Uncertainty in RVel-S12
      real*4        oRA_S13     ! (mas) ? Offset in RA, S13 relative to SgrA*
      real*4        e_oRA_S13   ! (mas) ? Uncertainty in oRA-S13
      real*4        oDE_S13     ! (mas) ? Offset in Declination, S13 relative to SgrA*
      real*4        e_oDE_S13   ! (mas) ? Uncertainty in oDE-S13
      integer*4     RV_S13      ! (km/s) [-409/602]? Radial Velocity for S13
      integer*4     e_RV_S13    ! (km/s) [130/375]? Uncertainty in RVel-S13
      real*4        oRA_S14     ! (mas) ? Offset in RA, S14 relative to SgrA*
      real*4        e_oRA_S14   ! (mas) ? Uncertainty in oRA-S14
      real*4        oDE_S14     ! (mas) ? Offset in Declination, S14 relative to SgrA*
      real*4        e_oDE_S14   ! (mas) ? Uncertainty in oDE-S14
      integer*4     RV_S14      ! (km/s) [131/765]? Radial Velocity for S14
      integer*4     e_RV_S14    ! (km/s) [108/259]? Uncertainty in RVel-S14
      real*4        oRA_S17     ! (mas) ? Offset in RA, S17 relative to SgrA*
      real*4        e_oRA_S17   ! (mas) ? Uncertainty in oRA-S17
      real*4        oDE_S17     ! (mas) ? Offset in Declination, S17 relative to SgrA*
      real*4        e_oDE_S17   ! (mas) ? Uncertainty in oDE-S17
      integer*4     RV_S17      ! (km/s) [-165/864]? Radial Velocity for S17
      integer*4     e_RV_S17    ! (km/s) [68/242]? Uncertainty in RVel-S17
      real*4        oRA_S18     ! (mas) ? Offset in RA, S18 relative to SgrA*
      real*4        e_oRA_S18   ! (mas) ? Uncertainty in oRA-S18
      real*4        oDE_S18     ! (mas) ? Offset in Declination, S18 relative to SgrA*
      real*4        e_oDE_S18   ! (mas) ? Uncertainty in oDE-S18
      integer*4     RV_S18      ! (km/s) [-605/-190]? Radial Velocity for S18
      integer*4     e_RV_S18    ! (km/s) [114/342]? Uncertainty in RVel-S18
      real*4        oRA_S19     ! (mas) ? Offset in RA, S19 relative to SgrA*
      real*4        e_oRA_S19   ! (mas) ? Uncertainty in oRA-S19
      real*4        oDE_S19     ! (mas) ? Offset in Declination, S19 relative to SgrA*
      real*4        e_oDE_S19   ! (mas) ? Uncertainty in oDE-S19
      integer*4     RV_S19      ! (km/s) [-2326/-724]? Radial Velocity for S19
      integer*4     e_RV_S19    ! (km/s) [128/385]? Uncertainty in RVel-S19
      real*4        oRA_S21     ! (mas) ? Offset in RA, S21 relative to SgrA*
      real*4        e_oRA_S21   ! (mas) ? Uncertainty in oRA-S21
      real*4        oDE_S21     ! (mas) ? Offset in Declination, S21 relative to SgrA*
      real*4        e_oDE_S21   ! (mas) ? Uncertainty in oDE-S21
      integer*4     RV_S21      ! (km/s) [373/523]? Radial Velocity for S21
      integer*4     e_RV_S21    ! (km/s) [48/112]? Uncertainty in RVel-S21
      real*4        oRA_S24     ! (mas) ? Offset in RA, S24 relative to SgrA*
      real*4        e_oRA_S24   ! (mas) ? Uncertainty in oRA-S24
      real*4        oDE_S24     ! (mas) ? Offset in Declination, S24 relative to SgrA*
      real*4        e_oDE_S24   ! (mas) ? Uncertainty in oDE-S24
      integer*4     RV_S24      ! (km/s) [-987/-771]? Radial Velocity for S24
      integer*4     e_RV_S24    ! (km/s) [39/154]? Uncertainty in RVel-S24
      real*4        oRA_S31     ! (mas) ? Offset in RA, S31 relative to SgrA*
      real*4        e_oRA_S31   ! (mas) ? Uncertainty in oRA-S31
      real*4        oDE_S31     ! (mas) ? Offset in Declination, S31 relative to SgrA*
      real*4        e_oDE_S31   ! (mas) ? Uncertainty in oDE-S31
      integer*4     RV_S31      ! (km/s) [-480/1020]? Radial Velocity for S31
      integer*4     e_RV_S31    ! (km/s) [79/473]? Uncertainty in RVel-S31
      real*4        oRA_S38     ! (mas) ? Offset in RA, S38 relative to SgrA*
      real*4        e_oRA_S38   ! (mas) ? Uncertainty in oRA-S38
      real*4        oDE_S38     ! (mas) ? Offset in Declination, S38 relative to SgrA*
      real*4        e_oDE_S38   ! (mas) ? Uncertainty in oDE-S38
      integer*4     RV_S38      ! (km/s) [-189/-2]? Radial Velocity for S38
      integer*4     e_RV_S38    ! (km/s) [173/248]? Uncertainty in RVel-S38
      real*4        oRA_S54     ! (mas) ? Offset in RA, S54 relative to SgrA*
      real*4        e_oRA_S54   ! (mas) ? Uncertainty in oRA-S54
      real*4        oDE_S54     ! (mas) ? Offset in Declination, S54 relative to SgrA*
      real*4        e_oDE_S54   ! (mas) ? Uncertainty in oDE-S54
      integer*4     RV_S54      ! (km/s) [-1790/-1046]? Radial Velocity for S54
      integer*4     e_RV_S54    ! (km/s) [182/467]? Uncertainty in RVel-S54
      real*4        oRA_S55     ! (mas) ? Offset in RA, S55 relative to SgrA*
      real*4        e_oRA_S55   ! (mas) ? Uncertainty in oRA-S55
      real*4        oDE_S55     ! (mas) ? Offset in Declination, S55 relative to SgrA*
      real*4        e_oDE_S55   ! (mas) ? Uncertainty in oDE-S55
      integer*4     RV_S55      ! (km/s) [-732/-400]? Radial Velocity for S55
      integer*4     e_RV_S55    ! (km/s) [145/210]? Uncertainty in RVel-S55
      character*6   Inst        ! Observation Source (VLT (201 occurrences),
*                                NTT (11 occurrences), Keck (3 occurrences) or
*                                GEMINI)

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

C  Loading file 'table3.dat'	! Orbital parameters of the 40 stars for which
*                                 we were able to determine orbits

C  Format for file interpretation

    1 format(
     +  A4,1X,A1,1X,F8.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.2,1X,F4.2,1X,
     +  F6.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,F7.2,1X,F6.2,1X,F7.2,1X,F7.2,
     +  1X,A1,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,40
        read(1,'(A118)')ar__
        read(ar__,1)
     +  Star,f_Star,a,e_a,e,e_e,i,e_i,Omega,e_Omega,w,e_w,Tp,e_Tp,Per,
     +  e_Per,SpT,Kmag,r
        if(ar__(91:97) .EQ. '') Per = rNULL__
        if(ar__(99:105) .EQ. '') e_Per = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Star,f_Star,a,e_a,e,e_e,i,e_i,Omega,e_Omega,w,e_w,Tp,e_Tp,Per,
     +  e_Per,SpT,Kmag,r
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! *Orbital data used in this work

C  Format for file interpretation

    2 format(
     +  F8.3,1X,A2,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,I5,1X,I3,1X,
     +  F5.1,1X,F3.1,1X,F5.1,1X,F3.1,1X,I5,1X,I3,1X,F5.1,1X,F3.1,1X,
     +  F5.1,1X,F3.1,1X,I5,1X,I3,1X,F5.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,
     +  I4,1X,I3,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,I3,1X,I3,1X,F5.1,
     +  1X,F4.1,1X,F5.1,1X,F4.1,1X,I3,1X,I3,1X,F6.1,1X,F4.1,1X,F6.1,
     +  1X,F4.1,1X,I4,1X,I3,1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,I3,1X,
     +  I3,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,I4,1X,I3,1X,F6.1,1X,
     +  F4.1,1X,F6.1,1X,F4.1,1X,I4,1X,I3,1X,F6.1,1X,F4.1,1X,F6.1,1X,
     +  F4.1,1X,I5,1X,I3,1X,F6.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,I3,1X,I3,
     +  1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,I4,1X,I3,1X,F6.1,1X,F4.1,
     +  1X,F5.1,1X,F4.1,1X,I4,1X,I3,1X,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,
     +  1X,I4,1X,I3,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,I5,1X,I3,1X,
     +  F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,I4,1X,I3,1X,A6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,216
        read(1,'(A557)')ar__1
        read(ar__1,2)
     +  Date,Flag,oRA_S1,e_oRA_S1,oDE_S1,e_oDE_S1,RV_S1,e_RV_S1,
     +  oRA_S2,e_oRA_S2,oDE_S2,e_oDE_S2,RV_S2,e_RV_S2,oRA_S4,e_oRA_S4,
     +  oDE_S4,e_oDE_S4,RV_S4,e_RV_S4,oRA_S8,e_oRA_S8,oDE_S8,e_oDE_S8,
     +  RV_S8,e_RV_S8,oRA_S9,e_oRA_S9,oDE_S9,e_oDE_S9,RV_S9,e_RV_S9,
     +  oRA_S12,e_oRA_S12,oDE_S12,e_oDE_S12,RV_S12,e_RV_S12,oRA_S13,
     +  e_oRA_S13,oDE_S13,e_oDE_S13,RV_S13,e_RV_S13,oRA_S14,e_oRA_S14,
     +  oDE_S14,e_oDE_S14,RV_S14,e_RV_S14,oRA_S17,e_oRA_S17,oDE_S17,
     +  e_oDE_S17,RV_S17,e_RV_S17,oRA_S18,e_oRA_S18,oDE_S18,e_oDE_S18,
     +  RV_S18,e_RV_S18,oRA_S19,e_oRA_S19,oDE_S19,e_oDE_S19,RV_S19,
     +  e_RV_S19,oRA_S21,e_oRA_S21,oDE_S21,e_oDE_S21,RV_S21,e_RV_S21,
     +  oRA_S24,e_oRA_S24,oDE_S24,e_oDE_S24,RV_S24,e_RV_S24,oRA_S31,
     +  e_oRA_S31,oDE_S31,e_oDE_S31,RV_S31,e_RV_S31,oRA_S38,e_oRA_S38,
     +  oDE_S38,e_oDE_S38,RV_S38,e_RV_S38,oRA_S54,e_oRA_S54,oDE_S54,
     +  e_oDE_S54,RV_S54,e_RV_S54,oRA_S55,e_oRA_S55,oDE_S55,e_oDE_S55,
     +  RV_S55,e_RV_S55,Inst
        if(ar__1(13:18) .EQ. '') oRA_S1 = rNULL__
        if(ar__1(20:23) .EQ. '') e_oRA_S1 = rNULL__
        if(ar__1(25:30) .EQ. '') oDE_S1 = rNULL__
        if(ar__1(32:35) .EQ. '') e_oDE_S1 = rNULL__
        if(ar__1(37:41) .EQ. '') RV_S1 = iNULL__
        if(ar__1(43:45) .EQ. '') e_RV_S1 = iNULL__
        if(ar__1(47:51) .EQ. '') oRA_S2 = rNULL__
        if(ar__1(53:55) .EQ. '') e_oRA_S2 = rNULL__
        if(ar__1(57:61) .EQ. '') oDE_S2 = rNULL__
        if(ar__1(63:65) .EQ. '') e_oDE_S2 = rNULL__
        if(ar__1(67:71) .EQ. '') RV_S2 = iNULL__
        if(ar__1(73:75) .EQ. '') e_RV_S2 = iNULL__
        if(ar__1(77:81) .EQ. '') oRA_S4 = rNULL__
        if(ar__1(83:85) .EQ. '') e_oRA_S4 = rNULL__
        if(ar__1(87:91) .EQ. '') oDE_S4 = rNULL__
        if(ar__1(93:95) .EQ. '') e_oDE_S4 = rNULL__
        if(ar__1(97:101) .EQ. '') RV_S4 = iNULL__
        if(ar__1(103:105) .EQ. '') e_RV_S4 = iNULL__
        if(ar__1(107:111) .EQ. '') oRA_S8 = rNULL__
        if(ar__1(113:115) .EQ. '') e_oRA_S8 = rNULL__
        if(ar__1(117:122) .EQ. '') oDE_S8 = rNULL__
        if(ar__1(124:126) .EQ. '') e_oDE_S8 = rNULL__
        if(ar__1(128:131) .EQ. '') RV_S8 = iNULL__
        if(ar__1(133:135) .EQ. '') e_RV_S8 = iNULL__
        if(ar__1(137:141) .EQ. '') oRA_S9 = rNULL__
        if(ar__1(143:146) .EQ. '') e_oRA_S9 = rNULL__
        if(ar__1(148:153) .EQ. '') oDE_S9 = rNULL__
        if(ar__1(155:158) .EQ. '') e_oDE_S9 = rNULL__
        if(ar__1(160:162) .EQ. '') RV_S9 = iNULL__
        if(ar__1(164:166) .EQ. '') e_RV_S9 = iNULL__
        if(ar__1(168:172) .EQ. '') oRA_S12 = rNULL__
        if(ar__1(174:177) .EQ. '') e_oRA_S12 = rNULL__
        if(ar__1(179:183) .EQ. '') oDE_S12 = rNULL__
        if(ar__1(185:188) .EQ. '') e_oDE_S12 = rNULL__
        if(ar__1(190:192) .EQ. '') RV_S12 = iNULL__
        if(ar__1(194:196) .EQ. '') e_RV_S12 = iNULL__
        if(ar__1(198:203) .EQ. '') oRA_S13 = rNULL__
        if(ar__1(205:208) .EQ. '') e_oRA_S13 = rNULL__
        if(ar__1(210:215) .EQ. '') oDE_S13 = rNULL__
        if(ar__1(217:220) .EQ. '') e_oDE_S13 = rNULL__
        if(ar__1(222:225) .EQ. '') RV_S13 = iNULL__
        if(ar__1(227:229) .EQ. '') e_RV_S13 = iNULL__
        if(ar__1(231:235) .EQ. '') oRA_S14 = rNULL__
        if(ar__1(237:240) .EQ. '') e_oRA_S14 = rNULL__
        if(ar__1(242:246) .EQ. '') oDE_S14 = rNULL__
        if(ar__1(248:251) .EQ. '') e_oDE_S14 = rNULL__
        if(ar__1(253:255) .EQ. '') RV_S14 = iNULL__
        if(ar__1(257:259) .EQ. '') e_RV_S14 = iNULL__
        if(ar__1(261:265) .EQ. '') oRA_S17 = rNULL__
        if(ar__1(267:270) .EQ. '') e_oRA_S17 = rNULL__
        if(ar__1(272:277) .EQ. '') oDE_S17 = rNULL__
        if(ar__1(279:282) .EQ. '') e_oDE_S17 = rNULL__
        if(ar__1(284:287) .EQ. '') RV_S17 = iNULL__
        if(ar__1(289:291) .EQ. '') e_RV_S17 = iNULL__
        if(ar__1(293:298) .EQ. '') oRA_S18 = rNULL__
        if(ar__1(300:303) .EQ. '') e_oRA_S18 = rNULL__
        if(ar__1(305:310) .EQ. '') oDE_S18 = rNULL__
        if(ar__1(312:315) .EQ. '') e_oDE_S18 = rNULL__
        if(ar__1(317:320) .EQ. '') RV_S18 = iNULL__
        if(ar__1(322:324) .EQ. '') e_RV_S18 = iNULL__
        if(ar__1(326:331) .EQ. '') oRA_S19 = rNULL__
        if(ar__1(333:336) .EQ. '') e_oRA_S19 = rNULL__
        if(ar__1(338:343) .EQ. '') oDE_S19 = rNULL__
        if(ar__1(345:348) .EQ. '') e_oDE_S19 = rNULL__
        if(ar__1(350:354) .EQ. '') RV_S19 = iNULL__
        if(ar__1(356:358) .EQ. '') e_RV_S19 = iNULL__
        if(ar__1(360:365) .EQ. '') oRA_S21 = rNULL__
        if(ar__1(367:369) .EQ. '') e_oRA_S21 = rNULL__
        if(ar__1(371:376) .EQ. '') oDE_S21 = rNULL__
        if(ar__1(378:380) .EQ. '') e_oDE_S21 = rNULL__
        if(ar__1(382:384) .EQ. '') RV_S21 = iNULL__
        if(ar__1(386:388) .EQ. '') e_RV_S21 = iNULL__
        if(ar__1(390:395) .EQ. '') oRA_S24 = rNULL__
        if(ar__1(397:400) .EQ. '') e_oRA_S24 = rNULL__
        if(ar__1(402:407) .EQ. '') oDE_S24 = rNULL__
        if(ar__1(409:412) .EQ. '') e_oDE_S24 = rNULL__
        if(ar__1(414:417) .EQ. '') RV_S24 = iNULL__
        if(ar__1(419:421) .EQ. '') e_RV_S24 = iNULL__
        if(ar__1(423:428) .EQ. '') oRA_S31 = rNULL__
        if(ar__1(430:433) .EQ. '') e_oRA_S31 = rNULL__
        if(ar__1(435:439) .EQ. '') oDE_S31 = rNULL__
        if(ar__1(441:444) .EQ. '') e_oDE_S31 = rNULL__
        if(ar__1(446:449) .EQ. '') RV_S31 = iNULL__
        if(ar__1(451:453) .EQ. '') e_RV_S31 = iNULL__
        if(ar__1(455:460) .EQ. '') oRA_S38 = rNULL__
        if(ar__1(462:465) .EQ. '') e_oRA_S38 = rNULL__
        if(ar__1(467:471) .EQ. '') oDE_S38 = rNULL__
        if(ar__1(473:476) .EQ. '') e_oDE_S38 = rNULL__
        if(ar__1(478:481) .EQ. '') RV_S38 = iNULL__
        if(ar__1(483:485) .EQ. '') e_RV_S38 = iNULL__
        if(ar__1(487:491) .EQ. '') oRA_S54 = rNULL__
        if(ar__1(493:496) .EQ. '') e_oRA_S54 = rNULL__
        if(ar__1(498:503) .EQ. '') oDE_S54 = rNULL__
        if(ar__1(505:508) .EQ. '') e_oDE_S54 = rNULL__
        if(ar__1(510:514) .EQ. '') RV_S54 = iNULL__
        if(ar__1(516:518) .EQ. '') e_RV_S54 = iNULL__
        if(ar__1(520:524) .EQ. '') oRA_S55 = rNULL__
        if(ar__1(526:529) .EQ. '') e_oRA_S55 = rNULL__
        if(ar__1(531:536) .EQ. '') oDE_S55 = rNULL__
        if(ar__1(538:541) .EQ. '') e_oDE_S55 = rNULL__
        if(ar__1(543:546) .EQ. '') RV_S55 = iNULL__
        if(ar__1(548:550) .EQ. '') e_RV_S55 = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  Date,Flag,oRA_S1,e_oRA_S1,oDE_S1,e_oDE_S1,RV_S1,e_RV_S1,
     +  oRA_S2,e_oRA_S2,oDE_S2,e_oDE_S2,RV_S2,e_RV_S2,oRA_S4,e_oRA_S4,
     +  oDE_S4,e_oDE_S4,RV_S4,e_RV_S4,oRA_S8,e_oRA_S8,oDE_S8,e_oDE_S8,
     +  RV_S8,e_RV_S8,oRA_S9,e_oRA_S9,oDE_S9,e_oDE_S9,RV_S9,e_RV_S9,
     +  oRA_S12,e_oRA_S12,oDE_S12,e_oDE_S12,RV_S12,e_RV_S12,oRA_S13,
     +  e_oRA_S13,oDE_S13,e_oDE_S13,RV_S13,e_RV_S13,oRA_S14,e_oRA_S14,
     +  oDE_S14,e_oDE_S14,RV_S14,e_RV_S14,oRA_S17,e_oRA_S17,oDE_S17,
     +  e_oDE_S17,RV_S17,e_RV_S17,oRA_S18,e_oRA_S18,oDE_S18,e_oDE_S18,
     +  RV_S18,e_RV_S18,oRA_S19,e_oRA_S19,oDE_S19,e_oDE_S19,RV_S19,
     +  e_RV_S19,oRA_S21,e_oRA_S21,oDE_S21,e_oDE_S21,RV_S21,e_RV_S21,
     +  oRA_S24,e_oRA_S24,oDE_S24,e_oDE_S24,RV_S24,e_RV_S24,oRA_S31,
     +  e_oRA_S31,oDE_S31,e_oDE_S31,RV_S31,e_RV_S31,oRA_S38,e_oRA_S38,
     +  oDE_S38,e_oDE_S38,RV_S38,e_RV_S38,oRA_S54,e_oRA_S54,oDE_S54,
     +  e_oDE_S54,RV_S54,e_RV_S54,oRA_S55,e_oRA_S55,oDE_S55,e_oDE_S55,
     +  RV_S55,e_RV_S55,Inst
c    .......End.of.Just test output...........
      end do
      close(1)

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