From 8363e232d48c6a68bf1db5490a903384f5c72736 Mon Sep 17 00:00:00 2001
From: kamischi <karl-michael.schindler@web.de>
Date: Sun, 9 Jun 2024 08:19:33 +0200
Subject: [PATCH] reestablish jos for the size of osc
the usage of neps for the size of osc is wrong and gives a range check error.
---
source/split-f77/eels.f | 2 +-
source/split-f77/seteps.f | 4 ++--
2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/source/split-f77/eels.f b/source/split-f77/eels.f
index ad8c86c..a5b6a90 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 36e8563..278c745 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
--
GitLab