From cb4407dccbfee1b7b4f757b37f365b0ce09f3f1a Mon Sep 17 00:00:00 2001
From: kamischi <karl-michael.schindler@web.de>
Date: Mon, 28 Feb 2022 14:03:49 +0100
Subject: [PATCH] rename undeclared variable T to TN.
---
source/f90/surlos.f90 | 8 ++++----
source/split-f77/surlos.f | 6 +++---
2 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/source/f90/surlos.f90 b/source/f90/surlos.f90
index 7b44126..99ba215 100644
--- a/source/f90/surlos.f90
+++ b/source/f90/surlos.f90
@@ -19,7 +19,7 @@ double precision function surlos(dk, eps, d, layers, nper)
logical :: static, zero, skip
integer :: lstart, n
double precision, allocatable :: arg(:)
- double precision :: argmin, argmax, cn, cnm1, epsmac, sn, snm1, t
+ double precision :: argmin, argmax, cn, cnm1, epsmac, sn, snm1, tn
double complex :: a, b, csi, pnm2, pnm1, pn, pp, qnm2, qnm1, qn, qp, z
common / mulayr / argmin, argmax, epsmac
@@ -148,13 +148,13 @@ double precision function surlos(dk, eps, d, layers, nper)
endif ! n <= 0
if (arg(n) /= 0.0d0) then
- t = tanh(arg(n))
- b = eps(n) + csi * t
+ tn = tanh(arg(n))
+ b = eps(n) + csi * tn
if (zero(b)) then
surlos = 0.0d0
return
endif
- csi = eps(n) * (csi + t * eps(n)) / b
+ csi = eps(n) * (csi + tn * eps(n)) / b
endif
enddo
end function surlos
diff --git a/source/split-f77/surlos.f b/source/split-f77/surlos.f
index 5bad9ff..ba51f93 100644
--- a/source/split-f77/surlos.f
+++ b/source/split-f77/surlos.f
@@ -88,10 +88,10 @@ C
11 N = N-1
IF(N.LE.0) GOTO 15
IF(ARG(N).EQ.0.0D0) GOTO 11
- T = TANH(ARG(N))
- B = EPS(N)+CSI*T
+ TN = TANH(ARG(N))
+ B = EPS(N)+CSI*TN
IF(ZERO(B)) GOTO 13
- CSI = EPS(N)*(CSI+T*EPS(N))/B
+ CSI = EPS(N)*(CSI+TN*EPS(N))/B
GOTO 11
12 N = LSTART
13 IF(N.LE.1) GOTO 14
--
GitLab