290 static void PrmHexGetZ(const PrmNode pn[], PetscInt k, PetscInt zm, PetscReal zn[])  in PrmHexGetZ()  argument
292 …const PetscScalar zm1 = zm - 1, znl[8] = {pn[0].b + pn[0].h * (PetscScalar)k / zm1,       pn[1].b …  in PrmHexGetZ()
681   PetscInt  i, j, k, xs, xm, ys, ym, zs, zm, mx, my;  in THIInitial()  local
691   PetscCall(DMDAGetCorners(da, &zs, &ys, &xs, &zm, &ym, &xm));  in THIInitial()
698       for (k = zs; k < zs + zm; k++) {  in THIInitial()
699 …const PetscScalar zm1 = zm - 1, drivingx = thi->rhog * (prm[i + 1][j].b + prm[i + 1][j].h - prm[i …  in THIInitial()
754   PetscInt  xs, ys, xm, ym, zm, i, j, k, q, l;  in THIFunctionLocal()  local
761   xm = info->zm;  in THIFunctionLocal()
763   zm = info->xm;  in THIFunctionLocal()
778       for (k = 0; k < zm - 1; k++) {  in THIFunctionLocal()
782         PrmHexGetZ(pn, k, zm, zn);  in THIFunctionLocal()
821             const PetscReal   hz    = PetscRealPart(pn[0].h) / (zm - 1.);  in THIFunctionLocal()
879   PetscInt    i, j, xs, ys, zs, xm, ym, zm, mx, my, mz;  in THISurfaceStatistics()  local
886   PetscCall(DMDAGetCorners(da, &zs, &ys, &xs, &zm, &ym, &xm));  in THISurfaceStatistics()
887   PetscCheck(zs == 0 && zm == mz, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Unexpected decomposition");  in THISurfaceStatistics()
891       PetscReal u = PetscRealPart(x[i][j][zm - 1].u);  in THISurfaceStatistics()
1045   PetscInt  xs, ys, xm, ym, zm, i, j, k, q, l, ll;  local
1052   xm = info->zm;
1054   zm = info->xm;
1066       for (k = 0; k < zm - 1; k++) {
1072         PrmHexGetZ(pn, k, zm, zn);
1143             const PetscReal   hz    = PetscRealPart(pn[0].h) / (zm - 1);
1299     PetscInt i, j, k, xs, ys, zs, xm, ym, zm, mx, my, mz, rstart, cstart;  local
1303     PetscCall(DMDAGetCorners(daf, &zs, &ys, &xs, &zm, &ym, &xm));
1306     PetscCall(MatSetSizes(B, xm * ym * zm, xm * ym, mx * my * mz, mx * my));
1315         for (k = zs; k < zs + zm; k++) {
1316           PetscInt    i2 = i * ym + j, i3 = i2 * zm + k;
1333   PetscInt               xm, ym, zm, dim, dof = 2, starts[3], dims[3];  local
1339   PetscCall(DMDAGetCorners(da, 0, 0, 0, &zm, &ym, &xm));
1342 …PetscCall(MatSetSizes(A, dof * xm * ym * zm, dof * xm * ym * zm, PETSC_DETERMINE, PETSC_DETERMINE)…
1386       PetscInt           i, j, k, xs, xm, ys, ym, zs, zm;  local
1393       zm = range[3];
1396       PetscCheck(xm * ym * zm * dof <= nmax, PETSC_COMM_SELF, PETSC_ERR_PLIB, "should not happen");
1400         PetscCheck(nn == xm * ym * zm * dof, PETSC_COMM_SELF, PETSC_ERR_PLIB, "should not happen");
1403 …MT " %" PetscInt_FMT " %" PetscInt_FMT " %" PetscInt_FMT "\">\n", zs, zs + zm - 1, ys, ys + ym - 1…
1409           for (k = zs; k < zs + zm; k++) {