diff --git a/source/split-f77/eels.f b/source/split-f77/eels.f
index ad8c86c868ef42180f288ca807b4c400c203f764..a5b6a908b541a0741976c0bf7a0c2cbac492266e 100644
--- a/source/split-f77/eels.f
+++ b/source/split-f77/eels.f
@@ -193,7 +193,7 @@ C
          WN = WMIN+(IW-1)*DW                                            
          IF(WN.LT.0.0D0) GOTO 45                                        
          IF(WN.EQ.0.0D0) GOTO 40                                        
-         IF(.NOT.USER) CALL SETEPS(NEPS,NOS,OSC,EPSINF,WN,NAME)         
+         IF(.NOT.USER) CALL SETEPS(NEPS,NOS,OSC,JOS,EPSINF,WN,NAME)         
          X = WN/(2.0D0*ENER*PSIA)                                       
          IF(RATION) THEN                                                
             F = QRAT(X)*AIMAG(-2.0/(1.0+EPS(1)))                        
diff --git a/source/split-f77/seteps.f b/source/split-f77/seteps.f
index 36e85638db0f537756b49cf8d9433849bd70d0ec..278c7454b6ab4ce62a57e89c807f5d61f4f73085 100644
--- a/source/split-f77/seteps.f
+++ b/source/split-f77/seteps.f
@@ -1,4 +1,4 @@
-      SUBROUTINE SETEPS(NEPS,NOS,OSC,EPSINF,WN,NAME)                    
+      SUBROUTINE SETEPS(NEPS,NOS,OSC,JOS,EPSINF,WN,NAME)                    
 C    ****************************************************************** 
 C    *                                                                * 
 C    * SET UP LONG-WAVELENGTH DIELECTRIC FUNCTIONS OF THE LAYERS FOR  * 
@@ -9,7 +9,7 @@ C    ******************************************************************
       CHARACTER NAME*10
       DOUBLE PRECISION ARGMIN,ARGMAX,EPSINF,EPSMAC,OSC,THICK,WN,X       
       COMPLEX DENO,EPS,DENO1,DENO2                                      
-      DIMENSION NOS(NEPS),OSC(3,NEPS),EPSINF(NEPS),NAME(NEPS)           
+      DIMENSION NOS(NEPS),OSC(3,JOS),EPSINF(NEPS),NAME(NEPS)           
       COMMON/MULAYR/EPS(LMAX),THICK(LMAX),ARGMIN,ARGMAX,EPSMAC,         
      ,              LAYERS,NPER                                         
       J = 0