Adapted error condition in FIC element to react to GEO_DRAINAGE_TYPE

Author: WPK4FEM

applications/GeoMechanicsApplication/custom_elements/U_Pw_small_strain_FIC_element.cpp

@@ -13,6 +13,7 @@
 
 // Application includes
 #include "custom_elements/U_Pw_small_strain_FIC_element.h"
+#include "custom_utilities/constitutive_law_utilities.h"
 #include "custom_utilities/equation_of_motion_utilities.hpp"
 #include "custom_utilities/extrapolation_utilities.h"
 
@@ -121,13 +122,20 @@ int UPwSmallStrainFICElement<TDim, TNumNodes>::Check(const ProcessInfo& rCurrent
     int ierr = UPwSmallStrainElement<TDim, TNumNodes>::Check(rCurrentProcessInfo);
     if (ierr != 0) return ierr;
 
-    const PropertiesType& Prop = this->GetProperties();
+    const PropertiesType& r_properties = this->GetProperties();
 
     // Verify specific properties
-    if (Prop[IGNORE_UNDRAINED])
-        KRATOS_ERROR << "IGNORE_UNDRAINED cannot be used in FIC elements. Use "
-                        "Non FIC elements instead"
-                     << this->Id() << std::endl;
+    auto is_constant_pw_field = false;
+    if (r_properties.Has(GEO_DRAINAGE_TYPE)) {
+        is_constant_pw_field =
+            ConstitutiveLawUtilities::StringToDrainageType(r_properties[GEO_DRAINAGE_TYPE]) ==
+            DrainageType::CONSTANT_WATER_PRESSURE;
+    } else if (r_properties.Has(IGNORE_UNDRAINED)) {
+        is_constant_pw_field = r_properties[IGNORE_UNDRAINED];
+    }
+    KRATOS_ERROR_IF(is_constant_pw_field) << "IGNORE_UNDRAINED cannot be used in FIC elements. Use "
+                                             "Non FIC elements instead"
+                                          << this->Id() << std::endl;
 
     return ierr;