27{
28 TCLAP::CmdLine cmd(
29 "Scales a property of a mesh.\n\n"
30 "OpenGeoSys-6 software, version " +
32 ".\n"
33 "Copyright (c) 2012-2025, OpenGeoSys Community "
34 "(http://www.opengeosys.org)",
36
37 TCLAP::ValueArg<std::string> out_mesh_arg(
38 "o",
39 "out-mesh",
40 "Output (.vtu). The output mesh is stored to a file of this name",
41 true,
42 "",
43 "OUTPUT_FILE");
44 cmd.add(out_mesh_arg);
45
46 TCLAP::ValueArg<std::string> property_arg(
47 "p",
48 "property-name",
49 "the name of the property the values are stored for",
50 true,
51 "",
52 "PROP_NAME");
53 cmd.add(property_arg);
54
55 TCLAP::ValueArg<std::string> mesh_arg(
56 "m", "mesh", "Input (.vtu). The input mesh is read from this file",
57 true, "", "INPUT_FILE");
58 cmd.add(mesh_arg);
59
60 std::vector<std::string> allowed_units{"mm/a", "mm/month", "m/s"};
61 TCLAP::ValuesConstraint<std::string> allowed_units_constraints{
62 allowed_units};
63 TCLAP::ValueArg<std::string> unit_arg("u",
64 "input-unit",
65 "input unit of the data",
66 true,
67 "m/s",
68 &allowed_units_constraints);
69 cmd.add(unit_arg);
70
71 cmd.parse(argc, argv);
72
74
75 std::unique_ptr<MeshLib::Mesh> mesh(
77
79 if (unit_arg.getValue() == "m/s")
80 {
82 }
83 else if (unit_arg.getValue() == "mm/a")
84 {
85 scale = 1e-3 / (365.25 * 86400);
86 }
87 else if (unit_arg.getValue() == "mm/month")
88 {
89 scale = 1e-3 * (12.0 / (365.25 * 86400));
90 }
91
93
95
96 return EXIT_SUCCESS;
97}
GITINFOLIB_EXPORT const std::string ogs_version
void scale(PETScVector &x, PetscScalar const a)
MeshLib::Mesh * readMeshFromFile(const std::string &file_name, bool const compute_element_neighbors)
int writeMeshToFile(const MeshLib::Mesh &mesh, std::filesystem::path const &file_path, std::set< std::string > variable_output_names)
void scaleMeshPropertyVector(MeshLib::Mesh &mesh, std::string const &property_name, double factor)