Add missing terms for y-grad calculation

LLNL · Sep 22, 2023 · 4b6df21 · 4b6df21
1 parent 596213c
commit 4b6df21
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/bbb/oderhs.m b/bbb/oderhs.m
@@ -2480,6 +2480,8 @@ ccc      if (isupgon .eq. 1 .and. zi(ifld) .eq. 0.0) call neudif
      .                  (ex(ix,iy)+ex(ix1,iy))/gx(ix,iy)
                 seec(ix,iy) = seec(ix,iy) + t1*vol(ix,iy) - t2
             else
+                iyp1 = min(iy+1,ny+1)
+                iym1 = max(iy-1,0)
                 t1 = .5*cvgp*( vex(ix,iy)*
      .                  ave(gx(ix,iy),gx(ix2,iy))*gpex(ix,iy)/gxf(ix,iy) +
      .                  vex(ix1,iy)*