diff --git a/original/EELS2/eels b/original/EELS2/eels
index ba043095cf1923a8717abe3ffe6c341eee4a33bc..68fd45a863104976f2e9ba4dbb3cd0b8b88f7f91 100755
Binary files a/original/EELS2/eels and b/original/EELS2/eels differ
diff --git a/original/EELS2/eels.for b/original/EELS2/eels.for
index 2ca87a563e81b05697a981b3c3732366a8df655b..0ccb844331a6e7a8c74cb69aa1eef71ddee774ee 100755
--- a/original/EELS2/eels.for
+++ b/original/EELS2/eels.for
@@ -186,7 +186,7 @@ C                                                                       ABTI0209
          WN = WMIN+(IW-1)*DW                                            ABTI0215
          IF(WN.LT.0.0D0) GOTO 45                                        ABTI0216
          IF(WN.EQ.0.0D0) GOTO 40                                        ABTI0217
-         IF(.NOT.USER) CALL SETEPS(NEPS,NOS,OSC,EPSINF,WN)              ABTI0218
+         IF(.NOT.USER) CALL SETEPS(NEPS,NOS,OSC,JOS,EPSINF,WN)          ABTI0218
          X = WN/(2.0D0*ENER*PSIA)                                       ABTI0219
          IF(RATION) THEN                                                ABTI0220
             F = QRAT(X)*AIMAG(-2.0/(1.0+EPS(1)))                         ABTI0221
@@ -483,7 +483,7 @@ C                                                                       ABTI0508
       ERREST = 2.0D0*ERREST                                             ABTI0512
       GO TO 82                                                          ABTI0513
       END                                                               ABTI0514
-      SUBROUTINE SETEPS(NEPS,NOS,OSC,EPSINF,WN)                         ABTI0515
+      SUBROUTINE SETEPS(NEPS,NOS,OSC,JOS,EPSINF,WN)                     ABTI0515
 C    ****************************************************************** ABTI0516
 C    *                                                                * ABTI0517
 C    * SET UP LONG-WAVELENGTH DIELECTRIC FUNCTIONS OF THE LAYERS FOR  * ABTI0518
@@ -493,7 +493,7 @@ C    ****************************************************************** ABTI0521
       PARAMETER(LMAX=100)                                               ABTI0522
       DOUBLE PRECISION ARGMIN,ARGMAX,EPSINF,EPSMAC,OSC,THICK,WN,X       ABTI0523
       COMPLEX DENO,EPS                                                  ABTI0524
-      DIMENSION NOS(NEPS),OSC(3,NEPS),EPSINF(NEPS)                      ABTI0525
+      DIMENSION NOS(NEPS),OSC(3,JOS),EPSINF(NEPS)                       ABTI0525
       COMMON/MULAYR/EPS(LMAX),THICK(LMAX),ARGMIN,ARGMAX,EPSMAC,         ABTI0526
      ,              LAYERS,NPER                                         ABTI0527
       J = 0                                                             ABTI0528
diff --git a/source/f77/eels.for b/source/f77/eels.for
index 92632f004943c4f0d52ae1ff5a0ce978a9d0358e..e5bb6ea25b7f86e58d62dad928ef8764c550f34a 100644
--- a/source/f77/eels.for
+++ b/source/f77/eels.for
@@ -193,7 +193,7 @@ C                                                                       ABTI0209
          WN = WMIN+(IW-1)*DW                                            ABTI0215
          IF(WN.LT.0.0D0) GOTO 45                                        ABTI0216
          IF(WN.EQ.0.0D0) GOTO 40                                        ABTI0217
-         IF(.NOT.USER) CALL SETEPS(NEPS,NOS,OSC,EPSINF,WN,NAME)         ABTI0218
+         IF(.NOT.USER) CALL SETEPS(NEPS,NOS,OSC,JOS,EPSINF,WN,NAME)     ABTI0218
          X = WN/(2.0D0*ENER*PSIA)                                       ABTI0219
          IF(RATION) THEN                                                ABTI0220
             F = QRAT(X)*AIMAG(-2.0/(1.0+EPS(1)))                        ABTI0221
@@ -492,7 +492,7 @@ C                                                                       ABTI0508
       ERREST = 2.0D0*ERREST                                             ABTI0512
       GO TO 82                                                          ABTI0513
       END                                                               ABTI0514
-      SUBROUTINE SETEPS(NEPS,NOS,OSC,EPSINF,WN,NAME)                    ABTI0515
+      SUBROUTINE SETEPS(NEPS,NOS,OSC,JOS,EPSINF,WN,NAME)                ABTI0515
 C    ****************************************************************** ABTI0516
 C    *                                                                * ABTI0517
 C    * SET UP LONG-WAVELENGTH DIELECTRIC FUNCTIONS OF THE LAYERS FOR  * ABTI0518
@@ -503,7 +503,7 @@ C    ****************************************************************** ABTI0521
       CHARACTER NAME*10
       DOUBLE PRECISION ARGMIN,ARGMAX,EPSINF,EPSMAC,OSC,THICK,WN,X       ABTI0523
       COMPLEX DENO,EPS,DENO1,DENO2                                      ABTI0524
-      DIMENSION NOS(NEPS),OSC(3,NEPS),EPSINF(NEPS),NAME(NEPS)           ABTI0525
+      DIMENSION NOS(NEPS),OSC(3,JOS),EPSINF(NEPS),NAME(NEPS)            ABTI0525
       COMMON/MULAYR/EPS(LMAX),THICK(LMAX),ARGMIN,ARGMAX,EPSMAC,         ABTI0526
      ,              LAYERS,NPER                                         ABTI0527
       J = 0                                                             ABTI0528