comments

jonathanschilling · Oct 29, 2023 · d08092e · d08092e
1 parent 40b55e2
commit d08092e
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git 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
@@ -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