257   Units     units;  member
301   Units     units = thi->units;  in THIInitialize_HOM_A()  local
304 …p->b     = s - 1000 * units->meter + 500 * units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx)…  in THIInitialize_HOM_A()
311   Units     units = thi->units;  in THIInitialize_HOM_C()  local
314   p->b = s - 1000 * units->meter;  in THIInitialize_HOM_C()
317 …I / thi->Lx) * PetscSinReal(y * 2 * PETSC_PI / thi->Ly)) * units->Pascal * units->year / units->me…  in THIInitialize_HOM_C()
325   Units     units = thi->units;  in THIInitialize_HOM_X()  local
328 …p->b     = s - 1000 * units->meter + 500 * units->meter * PetscSinReal(x + PETSC_PI) * PetscSinRea…  in THIInitialize_HOM_X()
330   p->beta2 = 1000 * (r < 1 ? 2 : 0) * units->Pascal * units->year / units->meter;  in THIInitialize_HOM_X()
336   Units     units = thi->units;  in THIInitialize_HOM_Y()  local
340 …p->b = s - 1000 * units->meter + 500 * units->meter * PetscSinReal(x + PETSC_PI) * PetscSinReal(y …  in THIInitialize_HOM_Y()
341   if (PetscRealPart(p->b) > -700 * units->meter) p->b += 200 * units->meter;  in THIInitialize_HOM_Y()
343 … r) * PetscCosReal(x * 3 / 2) * PetscCosReal(y * 3 / 2)) * units->Pascal * units->year / units->me…  in THIInitialize_HOM_Y()
349   Units     units = thi->units;  in THIInitialize_HOM_Z()  local
353 …p->b     = s - 1000 * units->meter + 500 * units->meter * PetscSinReal(x + PETSC_PI) * PetscSinRea…  in THIInitialize_HOM_Z()
355 … r) * PetscCosReal(x * 3 / 2) * PetscCosReal(y * 3 / 2)) * units->Pascal * units->year / units->me…  in THIInitialize_HOM_Z()
361     Units units           = thi->units;  in THIFriction()  local
362 …  thi->friction.irefgam = 1. / (0.5 * PetscSqr(thi->friction.refvel * units->meter / units->year));  in THIFriction()
363 …thi->friction.eps2    = 0.5 * PetscSqr(thi->friction.epsvel * units->meter / units->year) * thi->f…  in THIFriction()
378     Units           units   = thi->units;  in THIViscosity()  local
381 …A                     = 1.e-16 * PetscPowReal(units->Pascal, -n) / units->year, /* softness parame…  in THIViscosity()
385     thi->viscosity.eps      = 0.5 * PetscSqr(1e-5 / units->year);  in THIViscosity()
424   PetscCall(PetscFree((*thi)->units));  in THIDestroy()
434   Units            units;  in THICreate()  local
444   PetscCall(PetscNew(&thi->units));  in THICreate()
445   units           = thi->units;  in THICreate()
446   units->meter    = 1e-2;  in THICreate()
447   units->second   = 1e-7;  in THICreate()
448   units->kilogram = 1e-12;  in THICreate()
452 …OptionsReal("-units_meter", "1 meter in scaled length units", "", units->meter, &units->meter, NUL…  in THICreate()
453 …OptionsReal("-units_second", "1 second in scaled time units", "", units->second, &units->second, N…  in THICreate()
454 …nsReal("-units_kilogram", "1 kilogram in scaled mass units", "", units->kilogram, &units->kilogram…  in THICreate()
457   units->Pascal = units->kilogram / (units->meter * PetscSqr(units->second));  in THICreate()
458   units->year   = 31556926. * units->second; /* seconds per year */  in THICreate()
542   thi->Lx *= units->meter;  in THICreate()
543   thi->Ly *= units->meter;  in THICreate()
544   thi->Lz *= units->meter;  in THICreate()
551 …units->meter / (3e7 * units->second), gradu = u / (100 * units->meter), eta, deta, rho = 910 * uni…  in THICreate()
552               driving = rho * grav * PetscSinReal(thi->alpha) * 1000 * units->meter;  in THICreate()
556 … %8.2g  kg %8.2g  Pa %8.2g\n", (double)units->meter, (double)units->second, (double)units->kilogra…  in THICreate()
557 …le)thi->Lx, (double)thi->Ly, (double)thi->Lz, (double)(rho * grav * 1e3 * units->meter), (double)d…  in THICreate()
604 …PetscReal Lx = thi->Lx / thi->units->meter, Ly = thi->Ly / thi->units->meter, Lz = thi->Lz / thi->…  in THISetUpDM()
943     nrm2 *= thi->units->year / thi->units->meter;  in THISolveStatistics()
945       min[j] *= thi->units->year / thi->units->meter;  in THISolveStatistics()
946       max[j] *= thi->units->year / thi->units->meter;  in THISolveStatistics()
953       umin *= thi->units->year / thi->units->meter;  in THISolveStatistics()
954       umax *= thi->units->year / thi->units->meter;  in THISolveStatistics()
955       umean *= thi->units->year / thi->units->meter;  in THISolveStatistics()
1361   Units              units = thi->units;  local
1423 … (double)(PetscRealPart(ptr[i]) * units->year / units->meter), (double)(PetscRealPart(ptr[i + 1]) …