1 files changed, 56 insertions, 0 deletions
diff --git a/src/slalib/kbj.f b/src/slalib/kbj.f
new file mode 100644
index 0000000..6a3ddb4
--- /dev/null
+++ b/src/slalib/kbj.f
@@ -0,0 +1,56 @@
+      SUBROUTINE sla_KBJ (JB, E, K, J)
+*+
+*     - - - -
+*      K B J
+*     - - - -
+*
+*  Select epoch prefix 'B' or 'J'
+*
+*  Given:
+*     JB     int     sla_DBJIN prefix status:  0=none, 1='B', 2='J'
+*     E      dp      epoch - Besselian or Julian
+*
+*  Returned:
+*     K      char    'B' or 'J'
+*     J      int     status:  0=OK
+*
+*  If JB=0, B is assumed for E < 1984D0, otherwise J.
+*
+*  P.T.Wallace   Starlink   31 July 1989
+*
+*  Copyright (C) 1995 Rutherford Appleton Laboratory
+*-
+
+      IMPLICIT NONE
+
+      INTEGER JB
+      DOUBLE PRECISION E
+      CHARACTER K*(*)
+      INTEGER J
+
+*  Preset status
+      J=0
+
+*  If prefix given expressly, use it
+      IF (JB.EQ.1) THEN
+         K='B'
+      ELSE IF (JB.EQ.2) THEN
+         K='J'
+
+*  If no prefix, examine the epoch
+      ELSE IF (JB.EQ.0) THEN
+
+*     If epoch is pre-1984.0, assume Besselian;  otherwise Julian
+         IF (E.LT.1984D0) THEN
+            K='B'
+         ELSE
+            K='J'
+         END IF
+
+*  If illegal prefix, return error status
+      ELSE
+         K=' '
+         J=1
+      END IF
+
+      END