239   PetscScalar u, uxx, uyy, vx, vy, avx, avy, vxp, vxm, vyp, vym;  in FormFunctionLocal()  local
334       uxx       = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx;  in FormFunctionLocal()
336       f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx;  in FormFunctionLocal()
340       uxx       = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx;  in FormFunctionLocal()
342       f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy;  in FormFunctionLocal()
346       uxx           = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx;  in FormFunctionLocal()
348 …f[j][i].omega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy +…  in FormFunctionLocal()
352       uxx          = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx;  in FormFunctionLocal()
354 …f[j][i].temp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)…  in FormFunctionLocal()
380   PetscScalar   u, uxx, uyy, vx, vy, avx, avy, vxp, vxm, vyp, vym;  in NonlinearGS()  local
523             uxx   = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx;  in NonlinearGS()
525             fu    = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx - bjiu;  in NonlinearGS()
529             uxx   = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx;  in NonlinearGS()
531             fv    = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy - bjiv;  in NonlinearGS()
546             uxx    = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx;  in NonlinearGS()
548 …fomega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy + (vyp *…  in NonlinearGS()
559             uxx   = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx;  in NonlinearGS()
561 …ftemp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)) * hy …  in NonlinearGS()