98   PetscReal       dx, dz;  member
313   x  = (i - grid->jlid) * grid->dx;  in HorizVelocity()
330   x  = (i - grid->jlid - 0.5) * grid->dx;  in VertVelocity()
346   x  = (i - grid->jlid - 0.5) * grid->dx;  in Pressure()
412   eps11 = (uE - uW) / grid->dx;  in CalcSecInv()
414   eps12 = 0.5 * ((uN - uS) / grid->dz + (wE - wW) / grid->dx);  in CalcSecInv()
466   PetscScalar dx = grid->dx, dz = grid->dz;  in XMomentumResidual()  local
492   dPdx = (x[j][i + 1].p - x[j][i].p) / dx;  in XMomentumResidual()
493   if (j == jlim) dudzN = etaN * (x[j][i].w - x[j][i + 1].w) / dx;  in XMomentumResidual()
496   dudxE = etaE * (x[j][i + 1].u - x[j][i].u) / dx;  in XMomentumResidual()
497   dudxW = etaW * (x[j][i].u - x[j][i - 1].u) / dx;  in XMomentumResidual()
500            + (dudxE - dudxW) / dx + (dudzN - dudzS) / dz;  in XMomentumResidual()
503     dwdxN = etaN * (x[j][i + 1].w - x[j][i].w) / dx;  in XMomentumResidual()
504     dwdxS = etaS * (x[j - 1][i + 1].w - x[j - 1][i].w) / dx;  in XMomentumResidual()
506     residual += (dudxE - dudxW) / dx + (dwdxN - dwdxS) / dz;  in XMomentumResidual()
517   PetscScalar dx = grid->dx, dz = grid->dz;  in ZMomentumResidual()  local
549   else dwdxE = etaE * (x[j][i + 1].w - x[j][i].w) / dx;  in ZMomentumResidual()
550   dwdxW = 2.0 * etaW * (x[j][i].w - x[j][i - 1].w) / dx;  in ZMomentumResidual()
554            + (dwdzN - dwdzS) / dz + (dwdxE - dwdxW) / dx;  in ZMomentumResidual()
560     residual += (dwdzN - dwdzS) / dz + (dudzE - dudzW) / dx;  in ZMomentumResidual()
574   dudx = (uE - uW) / grid->dx;  in ContinuityResidual()
587   PetscScalar dx = grid->dx, dz = grid->dz;  in EnergyResidual()  local
595   dTdxE = (x[j][i + 1].T - x[j][i].T) / dx;  in EnergyResidual()
596   dTdxW = (x[j][i].T - x[j][i - 1].T) / dx;  in EnergyResidual()
599               (dTdxE - dTdxW) / dx) *  in EnergyResidual()
600              dx * dz / param->peclet;  in EnergyResidual()
623     fN = wN * TN * dx;  in EnergyResidual()
624     fS = wS * TS * dx;  in EnergyResidual()
632 …scAbsScalar(wN) * (-x[j + 2][i].T + 3.0 * (x[j + 1][i].T - x[j][i].T) + x[j - 1][i].T) / 8.0) * dx;  in EnergyResidual()
633 …scAbsScalar(wS) * (-x[j + 1][i].T + 3.0 * (x[j][i].T - x[j - 1][i].T) + x[j - 2][i].T) / 8.0) * dx;  in EnergyResidual()
672   return (uN - uS) / grid->dz + (wE - wW) / grid->dx;  in ShearStress()
680   PetscScalar dx = grid->dx, dz = grid->dz;  in XNormalStress()  local
711   return 2.0 * etaC * (uE - uW) / dx - pC;  in XNormalStress()
784   grid->dx     = param->width / ((PetscReal)(grid->ni - 2)); /* km */  in SetParams()
785   grid->dz     = grid->dx * PetscTanReal(param->slab_dip);   /* km */  in SetParams()
875   grid->dx            = grid->dx / param->L;                       /* dim'less */  in SetParams()
907 …Int_FMT "       [dx,dz] = %g, %g km\n", grid->ni, grid->nj, (double)(grid->dx * param->L), (double…  in ReportParams()
1038       PetscCall(VecSetValue(pars, 2, (PetscScalar)grid->dx, INSERT_VALUES));  in DoOutput()