239   PetscScalar u, uxx, uyy, vx, vy, avx, avy, vxp, vxm, vyp, vym;  in FormFunctionLocal()  local
335       uyy       = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy;  in FormFunctionLocal()
336       f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx;  in FormFunctionLocal()
341       uyy       = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy;  in FormFunctionLocal()
342       f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy;  in FormFunctionLocal()
347       uyy           = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy;  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()
353       uyy          = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy;  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
524             uyy   = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy;  in NonlinearGS()
525             fu    = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx - bjiu;  in NonlinearGS()
530             uyy   = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy;  in NonlinearGS()
531             fv    = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy - bjiv;  in NonlinearGS()
547             uyy    = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy;  in NonlinearGS()
548 …fomega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy + (vyp *…  in NonlinearGS()
560             uyy   = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy;  in NonlinearGS()
561 …ftemp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)) * hy …  in NonlinearGS()