Skip to content

Commit 9cab422

Browse files
committed
Fix dsigma for nVarCond=1
1 parent 5631ab0 commit 9cab422

1 file changed

Lines changed: 1 addition & 1 deletion

File tree

src/radial.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -894,7 +894,7 @@ subroutine transportProperties()
894894
a =(-one+con_LambdaOut)/(s1-s2)
895895
c =(s1-s2*con_LambdaOut)/(s1-s2)
896896
sigma(:) = a*tanh(b*(r0-r(:)))+c
897-
dsigma(:) =-a*b/cosh(b*(r0-r(:)))
897+
dsigma(:) =-a*b/(cosh(b*(r0-r(:))))**2
898898
lambda(:) =one/sigma(:)
899899
dLlambda(:)=-dsigma(:)/sigma(:)
900900

0 commit comments

Comments
 (0)