255   PetscReal Lx, Ly, Lz; /* Model domain */  member
304 …p->b     = s - 1000 * units->meter + 500 * units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx)…  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()
326 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; …  in THIInitialize_HOM_X()
337 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; …  in THIInitialize_HOM_Y()
350 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; …  in THIInitialize_HOM_Z()
460   thi->Lx              = 10.e3;  in THICreate()
473     L = thi->Lx;  in THICreate()
475     if (flg) thi->Lx = thi->Ly = L;  in THICreate()
476     PetscCall(PetscOptionsReal("-thi_Lx", "X Domain size (m)", "", thi->Lx, &thi->Lx, NULL));  in THICreate()
542   thi->Lx *= units->meter;  in THICreate()
557 …(%6.2g,%6.2g,%6.2g), pressure %8.2g, driving stress %8.2g\n", (double)thi->Lx, (double)thi->Ly, (d…  in THICreate()
579       PetscReal xx = thi->Lx * i / mx, yy = thi->Ly * j / my;  in THIInitializePrm()
604 …PetscReal Lx = thi->Lx / thi->units->meter, Ly = thi->Ly / thi->units->meter, Lz = thi->Lz / thi->…  in THISetUpDM()  local
606 … PetscInt_FMT "), size (m) %g x %g\n", level, (double)Lx, (double)Ly, Mx, My, Mx * My, (double)(Lx…  in THISetUpDM()
608 …), size (m) %g x %g x %g\n", level, (double)Lx, (double)Ly, (double)Lz, Mx, My, Mz, Mx * My * Mz, …  in THISetUpDM()
694   hx = thi->Lx / mx;  in THIInitial()
764   hx = thi->Lx / info->mz;  in THIFunctionLocal()
976   hx = thi->Lx / info->my;  in THIJacobianLocal_2D()
1055   hx = thi->Lx / info->mz;
1411             PetscReal xx = thi->Lx * i / mx, yy = thi->Ly * j / my, zz;