223   PetscScalar beta2; /* friction */  member
248   PRange    beta2;  member
345 …p->beta2        = -1e-10; /* This value is not used, but it should not be huge because that would …  in THIInitialize_HOM_A()
356 …p->beta2 = 1000 * (1 + PetscSinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * 2 * PETSC_PI / …  in THIInitialize_HOM_C()
370   p->beta2        = -1e-10;  in THIInitialize_HOM_F()
382   p->beta2 = 1000 * (r < 1 ? 2 : 0) * units->Pascal * units->year / units->meter / thi->rhog;  in THIInitialize_HOM_X()
394 …p->beta2 = 1000 * (1. + PetscSinReal(PetscSqrtReal(16 * r)) / PetscSqrtReal(1e-2 + 16 * r) * Petsc…  in THIInitialize_HOM_Y()
405 …p->beta2 = 1000 * (1. + PetscSinReal(PetscSqrtReal(16 * r)) / PetscSqrtReal(1e-2 + 16 * r) * Petsc…  in THIInitialize_HOM_Z()
408 static void THIFriction(THI thi, PetscReal rbeta2, PetscReal gam, PetscReal *beta2, PetscReal *dbet…  in THIFriction()  argument
416     *beta2  = rbeta2;  in THIFriction()
419 …*beta2  = rbeta2 * PetscPowReal(thi->friction.eps2 + gam * thi->friction.irefgam, thi->friction.ex…  in THIFriction()
420 …*dbeta2 = thi->friction.exponent * *beta2 / (thi->friction.eps2 + gam * thi->friction.irefgam) * t…  in THIFriction()
627   PRangeClear(&thi->beta2);  in THICreate()
833               PetscReal       beta2, dbeta2;  in THIFunctionLocal_3D()  local
837                 rbeta2 += phi[l] * pn[l].beta2;  in THIFunctionLocal_3D()
839 …        THIFriction(thi, PetscRealPart(rbeta2), PetscRealPart(u * u + v * v) / 2, &beta2, &dbeta2);  in THIFunctionLocal_3D()
840               RangeUpdate(&beta2min, &beta2max, beta2);  in THIFunctionLocal_3D()
843                 fn[ls + l]->u += pp * jw * beta2 * u;  in THIFunctionLocal_3D()
844                 fn[ls + l]->v += pp * jw * beta2 * v;  in THIFunctionLocal_3D()
854   PetscCall(PRangeMinMax(&thi->beta2, beta2min, beta2max));  in THIFunctionLocal_3D()
888       f[i][j].beta2 = prmdot[i][j].beta2;  in THIFunctionLocal_2D()
981 …                         (double)thi->eta.cmax, (double)thi->beta2.cmin, (double)thi->beta2.cmax));  in THIMatrixStatistics()
1078 …rged range [%g, %g]\n", (double)thi->beta2.min, (double)thi->beta2.max, (double)thi->beta2.cmin, (…  in THISolveStatistics()
1177               PetscReal       beta2, dbeta2;  in THIJacobianLocal_Momentum()  local
1181                 rbeta2 += phi[l] * pn[l].beta2;  in THIJacobianLocal_Momentum()
1183 …        THIFriction(thi, PetscRealPart(rbeta2), PetscRealPart(u * u + v * v) / 2, &beta2, &dbeta2);  in THIJacobianLocal_Momentum()
1188 …               Ke[l * 2 + 0][ll * 2 + 0] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * u * u * ppl;  in THIJacobianLocal_Momentum()
1191 …               Ke[l * 2 + 1][ll * 2 + 1] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * v * v * ppl;  in THIJacobianLocal_Momentum()