257   } erosion;  member
446 …etscSqr(vel->v)) / PetscSqr(thi->erosion.refvel), rate = -thi->erosion.rate * PetscPowScalar(magre…  in THIErosion()
449     if (thi->erosion.exponent == 1) {  in THIErosion()
453 …derate->u = 0.5 * thi->erosion.exponent * rate / magref2 * 2. * vel->u / PetscSqr(thi->erosion.ref…  in THIErosion()
454 …derate->v = 0.5 * thi->erosion.exponent * rate / magref2 * 2. * vel->v / PetscSqr(thi->erosion.ref…  in THIErosion()
532   thi->erosion.rate     = 1e-3; /* m/a */  in THICreate()
533   thi->erosion.exponent = 1.;  in THICreate()
534   thi->erosion.refvel   = 1.; /* m/a */  in THICreate()
599 …e to sliding velocity at reference velocity (m/a)", NULL, thi->erosion.rate, &thi->erosion.rate, N…  in THICreate()
600 … sliding velocity appearing in erosion relation", NULL, thi->erosion.exponent, &thi->erosion.expon…  in THICreate()
601 …", "Reference sliding velocity for erosion (m/a)", NULL, thi->erosion.refvel, &thi->erosion.refvel…  in THICreate()
602     thi->erosion.rate *= units->meter / units->year;  in THICreate()
603     thi->erosion.refvel *= units->meter / units->year;  in THICreate()