diff --git a/source/split-f77/seteps.f b/source/split-f77/seteps.f index 278c7454b6ab4ce62a57e89c807f5d61f4f73085..a4f87daf480ea096d7007d4f56bd6757781827a1 100644 --- a/source/split-f77/seteps.f +++ b/source/split-f77/seteps.f @@ -29,12 +29,14 @@ C ****************************************************************** EPS(L) = -OSC(1,J)**2/( WN**2+CMPLX(0.0D0,WN*OSC(3,J)) ) ELSE EPS(L) = EPSINF(L) + do k = 1, nos(l) J = J+1 X = WN/OSC(1,J) DENO = X*CMPLX(X,OSC(3,J)) IF(OSC(2,J).GE.0.0D0) DENO = 1.0-DENO IF(ABS(DENO).EQ.0.0) DENO = EPSMAC EPS(L) = EPS(L)+OSC(2,J)/DENO + enddo ENDIF 2 CONTINUE IF(NEPS.EQ.LAYERS+1) THEN