INTEGRALS: Fix editting error.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2020-12-25 11:29:51 -03:00
parent 4ece753176
commit 62d9a6f97b

View File

@ -743,7 +743,7 @@ module integrals
tmp(:,1,:) = 0.5d+00 * (pdata%q(iby,nbu:neu,nb ,nb:ne) & tmp(:,1,:) = 0.5d+00 * (pdata%q(iby,nbu:neu,nb ,nb:ne) &
- pdata%q(iby,nbl:nel,nb ,nb:ne)) / dh(1) & - pdata%q(iby,nbl:nel,nb ,nb:ne)) / dh(1) &
- (pdata%q(ibx,nb :ne ,nb ,nb:ne) & - (pdata%q(ibx,nb :ne ,nb ,nb:ne) &
- pdata%quibx,nb :ne ,nbl,nb:ne)) / dh(2) - pdata%q(ibx,nb :ne ,nbl,nb:ne)) / dh(2)
inarr(15) = inarr(15) - resistivity * sum(sign(tmp(:,1,:), & inarr(15) = inarr(15) - resistivity * sum(sign(tmp(:,1,:), &
pdata%q(ibx,nb:ne,nb ,nb:ne))) * dxz pdata%q(ibx,nb:ne,nb ,nb:ne))) * dxz