241   PetscReal Lx, Ly, Lz; /* Model domain */  member
343 …* units->meter + 500 * units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * …  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()
366 …* units->meter + 100 * units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx); /* * sin(y*2*PETSC…  in THIInitialize_HOM_F()
368 …p->h            = (1 - (atan((x - thi->Lx / 2) / 1.) + PETSC_PI / 2.) / PETSC_PI) * 500 * units->m…  in THIInitialize_HOM_F()
378 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; …  in THIInitialize_HOM_X()
389 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; …  in THIInitialize_HOM_Y()
401 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; …  in THIInitialize_HOM_Z()
524     thi->Lx            = 10.e3;  in THICreate()
543     L = thi->Lx;  in THICreate()
545     if (flg) thi->Lx = thi->Ly = L;  in THICreate()
546     PetscCall(PetscOptionsReal("-thi_Lx", "X Domain size (m)", "", thi->Lx, &thi->Lx, NULL));  in THICreate()
621   thi->Lx *= units->meter;  in THICreate()
636 …(%6.2g,%6.2g,%6.2g), pressure %8.2g, driving stress %8.2g\n", (double)thi->Lx, (double)thi->Ly, (d…  in THICreate()
662 … < info.gys + info.gym; j++) x2[i][j].b += PetscSinReal(thi->alpha) * thi->Lx * (i < 0 ? 1.0 : -1.…  in THIFixGhosts()
678       PetscReal xx = thi->Lx * i / mx, yy = thi->Ly * j / my;  in THIInitializePrm()
706   hx = thi->Lx / mx;  in THIInitial()
761   hx = thi->Lx / info->mz;  in THIFunctionLocal_3D()
1105   hx = thi->Lx / info->mz;  in THIJacobianLocal_Momentum()
1414 …PetscReal xx = thi->Lx * i / mx, yy = thi->Ly * j / my, b = PetscRealPart(y2[i * ym + j][FieldOffs…  in THIDAVecView_VTK_XML()
1560 …PetscReal Lx = thi->Lx / thi->units->meter, Ly = thi->Ly / thi->units->meter, Lz = thi->Lz / thi->…  in main()  local
1564 …nts %3d x %3d x %3d (%8d), size (m) %g x %g x %g\n", i, Lx, Ly, Lz, Mx, My, Mz, Mx * My * Mz, Lx /…  in main()