diff --git a/src/precondn.f90 b/src/precondn.f90 index 40e416a70..16c2faba2 100644 --- a/src/precondn.f90 +++ b/src/precondn.f90 @@ -73,8 +73,8 @@ SUBROUTINE precondn(lu1, bsq, gsqrt, r12, & t2 = cp25*(xue(l)/shalf(js) + xuo(l))/shalf(js) t3 = cp25*(xue(l-1)/shalf(js) + xuo(l-1))/shalf(js) ax(js,1) = ax(js,1) + ptau(lk)*t1*t1 - ax(js,2) = ax(js,2) + ptau(lk)*(-t1+t3)*(t1+t2) - ax(js,3) = ax(js,3) + ptau(lk)*(t1+t2)*(t1+t2) + ax(js,2) = ax(js,2) + ptau(lk)*(-t1+t3)*( t1+t2) + ax(js,3) = ax(js,3) + ptau(lk)*( t1+t2)*( t1+t2) ax(js,4) = ax(js,4) + ptau(lk)*(-t1+t3)*(-t1+t3) END DO diff --git a/src/scalfor.f90 b/src/scalfor.f90 index 8b008f48c..c7a4cfb3d 100644 --- a/src/scalfor.f90 +++ b/src/scalfor.f90 @@ -98,7 +98,7 @@ SUBROUTINE scalfor(gcx, axm, bxm, axd, bxd, cx, iflag, skip_scalfor_dbg) ! AND zeq IS THE EQUILIBRIUM EDGE VALUE OF Z00 mult_fac = MIN(fac, fac*hs*15.0_dp) - IF (iflag .eq. 1) THEN + IF (iflag .eq. 1) THEN ! this is only active for z ! METHOD 1: SUBTRACT (INSTABILITY) Pedge ~ fac*z/hs FROM PRECONDITIONER AT EDGE dx(ns,0,0) = dx(ns,0,0)*(1.0_dp-mult_fac)/(1.0_dp+edge_pedestal) END IF