 +      ENDIF
 +      RETURN
 +      END
++C
++      DOUBLE PRECISION FUNCTION BETAEV(IPRN, TTAG, BTTAG, BSMPL, P2, V2)
++C     NAME           BETAEV
++C
++C     PURPOSE        EVALUATE BETA EVOLUTION OVER THE PAST PERIOD
++C
++C     PARAMETERS     DESCRIPTION
++C     IPRN           SV PRN NUMBER
++C     TTAG           OBSERVATION EPOCH TIME TAG (EG SEC OF GPS WEEK)
++C     BTTAG          TIME OF SAMPLED BETA VALUES
++C     BSMPL          SAMPLED BETA VALUES
++C     P2             SQUARE OF THE POSITION VECTOR NORM
++C     V2             SQUARE OF THE VELOCITY VECTOR NORM
++C
++C     RETURN         SIGNED BETA EVOLUTION OVER THE PAST PERIOD,
++C                    OR 0.0D0 IF LESS THAN TWO BETA SAMPLES ARE AVAILABLE
++C
++      IMPLICIT NONE
++      INTEGER   MAXSAT
++      PARAMETER (MAXSAT=136)
++      INTEGER*4 IPRN
++      REAL*8    TTAG, BTTAG(MAXSAT, 2), BSMPL(MAXSAT, 2), P2, V2
++      REAL*8    T01, T20, T21, A, PERIOD, GM, PI, DB, DT
++      REAL*8    SQRT
++      DATA      GM / 3.986004415D14 /
++      DATA      PI / 3.141592653589793D0 /
++C
++      IF (BTTAG(IPRN, 2) .EQ. 0.0D0) THEN
++C         WE HAVE ONLY ONE POINT YET, NO INTERPOLATION
++         BETAEV = 0.0D0
++      ELSE
++C         COMPUTE PERIOD
++C         (WE DON'T USE 360/MURATE BECAUSE MURATE IS THE LOCAL ANGULAR VELOCITY)
++          A      = 1.0D0 / (2.0 / SQRT(P2) - V2 / GM)
++          PERIOD = 2.0D0 * PI * A * SQRT(A / GM)
++C
++          DB     = BSMPL(IPRN, 1) - BSMPL(IPRN, 2)
++          DT     = BTTAG(IPRN, 1) - BTTAG(IPRN, 2)
++          BETAEV = DB * PERIOD / DT
++      ENDIF
++      RETURN
++      END

 C Mar 09, 2016
 C      IF(IPRN.GT.32.AND.ABS(BETADG).LT.ANOON) THEN
        IF(IPRN.GT.32.AND.IPRN.LE.64.AND.ABS(BETADG).LT.ANOON) THEN
@@ -641,7 +652,11 @@
 C           IF(IDIR.GT.0.AND.IPRN.LE.100.AND.ABS(BETADG).LE.0.07D0
          IF(IDIR.GT.0.AND.(IPRN.LE.100.OR.IPRN.EQ.INPRN)

--- eclips_Dec2017.f.orig	2018-10-08 10:47:29.305440637 +0200
+++ eclips_Dec2017.f	2018-10-08 10:49:06.397996483 +0200
@@ -120,6 +120,9 @@
 C                      Jul 12, 2018  UPDATE ESTIMATED ECLIPSE START AND END TIME
 C                                    WHEN THEY ARE ESTIMATED FROM A CLOSER
