[Shape opt] water drain response function

applications/ShapeOptimizationApplication/custom_python/add_custom_utilities_to_python.cpp

@@ -37,6 +37,7 @@
 #include "custom_utilities/input_output/universal_file_io.h"
 #include "custom_utilities/search_based_functions.h"
 #include "custom_utilities/response_functions/face_angle_response_function_utility.h"
+#include "custom_utilities/response_functions/water_drain_response_function_utility.h"
 
 // ==============================================================================
 
@@ -294,6 +295,14 @@ void  AddCustomUtilitiesToPython(pybind11::module& m)
         .def("CalculateGradient", &FaceAngleResponseFunctionUtility::CalculateGradient)
         ;
 
+    py::class_<WaterDrainResponseFunctionUtility >(m, "WaterDrainResponseFunctionUtility")
+        .def(py::init<ModelPart&, Parameters>())
+        .def("Initialize", &WaterDrainResponseFunctionUtility::Initialize)
+        .def("InitializeSolutionStep", &WaterDrainResponseFunctionUtility::InitializeSolutionStep)
+        .def("CalculateValue", &WaterDrainResponseFunctionUtility::CalculateValue)
+        .def("CalculateGradient", &WaterDrainResponseFunctionUtility::CalculateGradient)
+        ;
+
     // ========================================================================
     // Additional operations
     // ========================================================================

applications/ShapeOptimizationApplication/custom_python/shape_optimization_python_application.cpp

@@ -120,6 +120,10 @@ PYBIND11_MODULE(KratosShapeOptimizationApplication, m)
     KRATOS_REGISTER_IN_PYTHON_3D_VARIABLE_WITH_COMPONENTS(m, BACKGROUND_NORMAL);
     KRATOS_REGISTER_IN_PYTHON_3D_VARIABLE_WITH_COMPONENTS(m, OUT_OF_PLANE_DELTA);
 
+    // For water drain response
+    KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, WATER_LEVEL);
+    KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, WATER_VOLUMES);
+
     // For sensitivity heatmap
     KRATOS_REGISTER_IN_PYTHON_3D_VARIABLE_WITH_COMPONENTS(m, HEATMAP_DF1DX);
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, HEATMAP_DF1DALPHA);