83   PetscReal slab_dip, slab_age, slab_velocity, kappa, z_scale;  member
468   PetscScalar TE = 0.0, TN = 0.0, TS = 0.0, TW = 0.0, dPdx, residual, z_scale;  in XMomentumResidual()  local
472   z_scale = param->z_scale;  in XMomentumResidual()
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()
519   PetscScalar TE = 0.0, TN = 0.0, TS = 0.0, TW = 0.0, dPdz, residual, z_scale;  in ZMomentumResidual()  local
524   z_scale = param->z_scale;  in ZMomentumResidual()
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()
682   PetscScalar epsC = 0.0, etaC, TC, uE, uW, pC, z_scale;  in XNormalStress()  local
686   z_scale = param->z_scale;  in XNormalStress()
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()
722   PetscScalar pC, wN, wS, z_scale;  in ZNormalStress()  local
726   z_scale = param->z_scale;  in ZNormalStress()
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()
885   param->z_scale = param->L * alpha_g_on_cp_units_inverse_km;  in SetParams()
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()