77   PetscScalar u, v, omega, temp;  member
222       x[j][i].temp  = (grashof > 0) * i * dx;  in FormInitialGuess()
274       f[j][i].temp  = x[j][i].temp - x[j + 1][i].temp;  in FormFunctionLocal()
287       f[j][i].temp  = x[j][i].temp - x[j - 1][i].temp;  in FormFunctionLocal()
300       f[j][i].temp  = x[j][i].temp;  in FormFunctionLocal()
313       f[j][i].temp  = x[j][i].temp - (PetscReal)(grashof > 0);  in FormFunctionLocal()
348 … + vym * (x[j + 1][i].omega - u)) * hx - .5 * grashof * (x[j][i + 1].temp - x[j][i - 1].temp) * hy;  in FormFunctionLocal()
351       u            = x[j][i].temp;  in FormFunctionLocal()
352       uxx          = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx;  in FormFunctionLocal()
353       uyy          = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy;  in FormFunctionLocal()
354 …temp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)) * hy +…  in FormFunctionLocal()
512             bjitemp  = b[j][i].temp;  in NonlinearGS()
548 …ym * (x[j + 1][i].omega - u)) * hx - .5 * grashof * (x[j][i + 1].temp - x[j][i - 1].temp) * hy - b…  in NonlinearGS()
558             u     = x[j][i].temp;  in NonlinearGS()
559             uxx   = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx;  in NonlinearGS()
560             uyy   = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy;  in NonlinearGS()
561 … ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)) * hy + (vyp * (u - x[j - 1][i].tem…  in NonlinearGS()
563             if (PetscRealPart(vx) > 0.0) dftdu = prandtl * (u - x[j][i - 1].temp) * hy;  in NonlinearGS()
564             else dftdu = prandtl * (x[j][i + 1].temp - u) * hy;  in NonlinearGS()
566             if (PetscRealPart(vy) > 0.0) dftdv = prandtl * (u - x[j - 1][i].temp) * hx;  in NonlinearGS()
567             else dftdv = prandtl * (x[j + 1][i].temp - u) * hx;  in NonlinearGS()
583             x[j][i].temp  = x[j][i].temp - yt;  in NonlinearGS()
588             ftemp         = x[j][i].temp - bjitemp;  in NonlinearGS()
592             x[j][i].temp  = x[j][i].temp - ftemp;  in NonlinearGS()
596             ftemp         = x[j][i].temp - (PetscReal)(grashof > 0) - bjitemp;  in NonlinearGS()
600             x[j][i].temp  = x[j][i].temp - ftemp;  in NonlinearGS()
604             ftemp         = x[j][i].temp - x[j + 1][i].temp - bjitemp;  in NonlinearGS()
608             x[j][i].temp  = x[j][i].temp - ftemp;  in NonlinearGS()
612             ftemp         = x[j][i].temp - x[j - 1][i].temp - bjitemp;  in NonlinearGS()
616             x[j][i].temp  = x[j][i].temp - ftemp;  in NonlinearGS()
623 …* x[j][i].u + x[j][i].v * x[j][i].v + x[j][i].omega * x[j][i].omega + x[j][i].temp * x[j][i].temp);  in NonlinearGS()