255   PetscReal Lx, Ly, Lz; /* Model domain */  member
304 …units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * 2 * PETSC_PI / thi->Ly);  in THIInitialize_HOM_A()
317 …SinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * 2 * PETSC_PI / thi->Ly)) * units->Pascal * …  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()
461   thi->Ly              = 10.e3;  in THICreate()
475     if (flg) thi->Lx = thi->Ly = L;  in THICreate()
477     PetscCall(PetscOptionsReal("-thi_Ly", "Y Domain size (m)", "", thi->Ly, &thi->Ly, NULL));  in THICreate()
543   thi->Ly *= units->meter;  in THICreate()
557 …g), pressure %8.2g, driving stress %8.2g\n", (double)thi->Lx, (double)thi->Ly, (double)thi->Lz, (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 …), size (m) %g x %g\n", level, (double)Lx, (double)Ly, Mx, My, Mx * My, (double)(Lx / Mx), (double…  in THISetUpDM()
608 … %g x %g\n", level, (double)Lx, (double)Ly, (double)Lz, Mx, My, Mz, Mx * My * Mz, (double)(Lx / Mx…  in THISetUpDM()
695   hy = thi->Ly / my;  in THIInitial()
765   hy = thi->Ly / info->my;  in THIFunctionLocal()
977   hy = thi->Ly / info->mx;  in THIJacobianLocal_2D()
1056   hy = thi->Ly / info->my;
1411             PetscReal xx = thi->Lx * i / mx, yy = thi->Ly * j / my, zz;