82   PetscReal width, depth, scaled_width, scaled_depth, peclet, potentialT;  member
475     TS = param->potentialT * TInterp(x, i, j - 1) * PetscExpScalar((j - 1.0) * dz * z_scale);  in XMomentumResidual()
477     else TN = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale);  in XMomentumResidual()
478     TW = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale);  in XMomentumResidual()
479     TE = param->potentialT * x[j][i + 1].T * PetscExpScalar((j - 0.5) * dz * z_scale);  in XMomentumResidual()
528     TN = param->potentialT * x[j + 1][i].T * PetscExpScalar((j + 0.5) * dz * z_scale);  in ZMomentumResidual()
529     TS = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale);  in ZMomentumResidual()
530     TW = param->potentialT * TInterp(x, i - 1, j) * PetscExpScalar(j * dz * z_scale);  in ZMomentumResidual()
532     else TE = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale);  in ZMomentumResidual()
690     TC = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale);  in XNormalStress()
701     TC = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale);  in XNormalStress()
730     TC = param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * z_scale);  in ZNormalStress()
740     TC = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale);  in ZNormalStress()
807   param->potentialT    = 1300.0;    /* degrees C */  in SetParams()
813   PetscCall(PetscOptionsGetReal(NULL, NULL, "-potentialT", &param->potentialT, NULL));  in SetParams()
1048       PetscCall(VecSetValue(pars, 12, (PetscScalar)param->potentialT, INSERT_VALUES));  in DoOutput()
1125 …T = PetscRealPart(param->potentialT * x[j][i].T * PetscExpScalar((j - 0.5) * dz * param->z_scale));  in ViscosityField()
1127 …TC = PetscRealPart(param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * param->z_scale));  in ViscosityField()