1 files changed, 30 insertions, 0 deletions
diff --git a/src/slalib/epb.f b/src/slalib/epb.f
new file mode 100644
index 0000000..4cc4d49
--- /dev/null
+++ b/src/slalib/epb.f
@@ -0,0 +1,30 @@
+      DOUBLE PRECISION FUNCTION sla_EPB (DATE)
+*+
+*     - - - -
+*      E P B
+*     - - - -
+*
+*  Conversion of Modified Julian Date to Besselian Epoch
+*  (double precision)
+*
+*  Given:
+*     DATE     dp       Modified Julian Date (JD - 2400000.5)
+*
+*  The result is the Besselian Epoch.
+*
+*  Reference:
+*     Lieske,J.H., 1979. Astron.Astrophys.,73,282.
+*
+*  P.T.Wallace   Starlink   February 1984
+*
+*  Copyright (C) 1995 Rutherford Appleton Laboratory
+*-
+
+      IMPLICIT NONE
+
+      DOUBLE PRECISION DATE
+
+
+      sla_EPB = 1900D0 + (DATE-15019.81352D0)/365.242198781D0
+
+      END