239   PetscScalar beta2; /* friction */  member
261   PRange    beta2;  member
306   p->beta2 = 1e30;  in THIInitialize_HOM_A()
317 …p->beta2 = 1000 * (1 + PetscSinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * 2 * PETSC_PI / …  in THIInitialize_HOM_C()
330   p->beta2 = 1000 * (r < 1 ? 2 : 0) * units->Pascal * units->year / units->meter;  in THIInitialize_HOM_X()
343 …p->beta2 = 1000 * (1. + PetscSinReal(PetscSqrtReal(16 * r)) / PetscSqrtReal(1e-2 + 16 * r) * Petsc…  in THIInitialize_HOM_Y()
355 …p->beta2 = 1000 * (1. + PetscSinReal(PetscSqrtReal(16 * r)) / PetscSqrtReal(1e-2 + 16 * r) * Petsc…  in THIInitialize_HOM_Z()
358 static void THIFriction(THI thi, PetscReal rbeta2, PetscReal gam, PetscReal *beta2, PetscReal *dbet…  in THIFriction()  argument
366     *beta2  = rbeta2;  in THIFriction()
369 …*beta2  = rbeta2 * PetscPowReal(thi->friction.eps2 + gam * thi->friction.irefgam, thi->friction.ex…  in THIFriction()
370 …*dbeta2 = thi->friction.exponent * *beta2 / (thi->friction.eps2 + gam * thi->friction.irefgam) * t…  in THIFriction()
548   PRangeClear(&thi->beta2);  in THICreate()
829               PetscReal       beta2, dbeta2;  in THIFunctionLocal()  local
833                 rbeta2 += phi[l] * pn[l].beta2;  in THIFunctionLocal()
835 …        THIFriction(thi, PetscRealPart(rbeta2), PetscRealPart(u * u + v * v) / 2, &beta2, &dbeta2);  in THIFunctionLocal()
836               RangeUpdate(&beta2min, &beta2max, beta2);  in THIFunctionLocal()
839                 fn[ls + l]->u += pp * jw * beta2 * u;  in THIFunctionLocal()
840                 fn[ls + l]->v += pp * jw * beta2 * v;  in THIFunctionLocal()
852   PetscCall(PRangeMinMax(&thi->beta2, beta2min, beta2max));  in THIFunctionLocal()
871 …  (double)thi->eta.cmin, (double)thi->eta.cmax, (double)thi->beta2.cmin, (double)thi->beta2.cmax));  in THIMatrixStatistics()
960 …rged range %g to %g\n", (double)thi->beta2.min, (double)thi->beta2.max, (double)thi->beta2.cmin, (…  in THISolveStatistics()
991         PetscReal   phi[4], dphi[4][2], jw, eta, deta, beta2, dbeta2;  in THIJacobianLocal_2D()  local
998           rbeta2 += phi[l] * pn[l].beta2;  in THIJacobianLocal_2D()
1002         THIFriction(thi, PetscRealPart(rbeta2), PetscRealPart(u * u + v * v) / 2, &beta2, &dbeta2);  in THIJacobianLocal_2D()
1011 …= dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + pp * jw * (beta2 / h) * ppl * thi->…  in THIJacobianLocal_2D()
1014 …= dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + pp * jw * (beta2 / h) * ppl * thi->…  in THIJacobianLocal_2D()
1151               PetscReal       beta2, dbeta2;  local
1155                 rbeta2 += phi[l] * pn[l].beta2;
1157 …        THIFriction(thi, PetscRealPart(rbeta2), PetscRealPart(u * u + v * v) / 2, &beta2, &dbeta2);
1162 …               Ke[l * 2 + 0][ll * 2 + 0] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * u * u * ppl;
1165 …               Ke[l * 2 + 1][ll * 2 + 1] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * v * v * ppl;