Index: DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/soil_profile_1d_csv_tutorial.rst
===================================================================
diff -u -r3529 -r3536
--- DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/soil_profile_1d_csv_tutorial.rst	(.../soil_profile_1d_csv_tutorial.rst)	(revision 3529)
+++ DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/soil_profile_1d_csv_tutorial.rst	(.../soil_profile_1d_csv_tutorial.rst)	(revision 3536)
@@ -21,18 +21,24 @@
 
 Then a second loop is created that loops all csv rows with a certain profile id.
 A :py:class:`dampythoninterface.soilprofile1D.Layer1D` class is initialized for each row.
+If the soil type is sand the the IsAquifer property is set to True otherwise is set to False,
+otherwise the Piping calculation cannot be performed.
 This is append to the list of the 1D layers.
 
 .. code-block:: python
 
     for csv_row in profile_csv[profile_csv["soilprofile_id"] == profile_id].index:
+            if "zand" in profile_csv["soil_name"][csv_row]:
+                is_aquifer = True
+            else:
+                is_aquifer = False
             list_of_1D_layers.append(
                 dpi.Layer1D(
                     **{
                         "Name": profile_csv["soil_name"][csv_row],
                         "SoilName": profile_csv["soil_name"][csv_row],
                         "TopLevel": profile_csv["top_level"][csv_row],
-                        "IsAquifer": False,
+                        "IsAquifer": is_aquifer,
                         "WaterpressureInterpolationModel": dpi.WaterpressureInterpolationModelType.Automatic,
                     }
                 )
Index: DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/full_tutorial_snippet.rst
===================================================================
diff -u -r3529 -r3536
--- DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/full_tutorial_snippet.rst	(.../full_tutorial_snippet.rst)	(revision 3529)
+++ DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/full_tutorial_snippet.rst	(.../full_tutorial_snippet.rst)	(revision 3536)
@@ -148,13 +148,17 @@
         for profile_id in profile_csv["soilprofile_id"].unique():
             list_of_1D_layers = []
             for csv_row in profile_csv[profile_csv["soilprofile_id"] == profile_id].index:
+                if "zand" in profile_csv["soil_name"][csv_row]:
+                    is_aquifer = True
+                else:
+                    is_aquifer = False
                 list_of_1D_layers.append(
                     dpi.Layer1D(
                         **{
                             "Name": profile_csv["soil_name"][csv_row],
                             "SoilName": profile_csv["soil_name"][csv_row],
                             "TopLevel": profile_csv["top_level"][csv_row],
-                            "IsAquifer": False,
+                            "IsAquifer": is_aquifer,
                             "WaterpressureInterpolationModel": dpi.WaterpressureInterpolationModelType.Automatic,
                         }
                     )
Index: DamClients/DamPythonInterface/trunk/src/tests/test_dampythonframework.py
===================================================================
diff -u -r3534 -r3536
--- DamClients/DamPythonInterface/trunk/src/tests/test_dampythonframework.py	(.../test_dampythonframework.py)	(revision 3534)
+++ DamClients/DamPythonInterface/trunk/src/tests/test_dampythonframework.py	(.../test_dampythonframework.py)	(revision 3536)
@@ -181,13 +181,17 @@
             for csv_row in profile_csv[
                 profile_csv["soilprofile_id"] == profile_id
             ].index:
+                if "zand" in profile_csv["soil_name"][csv_row]:
+                    is_aquifer = True
+                else:
+                    is_aquifer = False
                 list_of_1D_layers.append(
                     dpi.Layer1D(
                         **{
                             "Name": profile_csv["soil_name"][csv_row],
                             "SoilName": profile_csv["soil_name"][csv_row],
                             "TopLevel": profile_csv["top_level"][csv_row],
-                            "IsAquifer": False,
+                            "IsAquifer": is_aquifer,
                             "WaterpressureInterpolationModel": dpi.WaterpressureInterpolationModelType.Automatic,
                         }
                     )
Index: DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/soil_profile_1d_coders_snippet.rst
===================================================================
diff -u -r3529 -r3536
--- DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/soil_profile_1d_coders_snippet.rst	(.../soil_profile_1d_coders_snippet.rst)	(revision 3529)
+++ DamClients/DamPythonInterface/trunk/doc/source/tutorial_sections/soil_profile_1d_coders_snippet.rst	(.../soil_profile_1d_coders_snippet.rst)	(revision 3536)
@@ -15,13 +15,17 @@
         for profile_id in profile_csv["soilprofile_id"].unique():
             list_of_1D_layers = []
             for csv_row in profile_csv[profile_csv["soilprofile_id"] == profile_id].index:
+                if "zand" in profile_csv["soil_name"][csv_row]:
+                    is_aquifer = True
+                else:
+                    is_aquifer = False
                 list_of_1D_layers.append(
                     dpi.Layer1D(
                         **{
                             "Name": profile_csv["soil_name"][csv_row],
                             "SoilName": profile_csv["soil_name"][csv_row],
                             "TopLevel": profile_csv["top_level"][csv_row],
-                            "IsAquifer": False,
+                            "IsAquifer": is_aquifer,
                             "WaterpressureInterpolationModel": dpi.WaterpressureInterpolationModelType.Automatic,
                         }
                     )