diff --git a/web/content/docs/tutorials/advancing-glacier/glacierclass.py b/web/content/docs/tutorials/advancing-glacier/glacierclass.py
index 82f9cde739c7ce7b0b539450ee53ed0ad5b3f2db..ccd75f959292a7f460a7ccfb561a0c5eb796d692 100644
--- a/web/content/docs/tutorials/advancing-glacier/glacierclass.py
+++ b/web/content/docs/tutorials/advancing-glacier/glacierclass.py
@@ -24,14 +24,14 @@ class glacier:
 
     # analytical function for the glacier's shape
     def local_height(self, x, t):
-        l = self.length(t)
-        if l == 0:
+        length = self.length(t)
+        if length == 0:
             return 0 * x
-        else:
-            xi = (x - self.x_0) / l
-            xi = np.array(xi)
-            xi[xi > 1] = 1.0
-            return self.height(t) * np.sqrt(1 - xi**1)
+
+        xi = (x - self.x_0) / length
+        xi = np.array(xi)
+        xi[xi > 1] = 1.0
+        return self.height(t) * np.sqrt(1 - xi**1)
 
     def height(self, t):
         return self.H_max * (t - self.t_0) / self.t_1
diff --git a/web/content/docs/tutorials/advancing-glacier/mesh_basin.py b/web/content/docs/tutorials/advancing-glacier/mesh_basin.py
index f3f3f50868cd4ad391e6df6e289ffa78fb1bd0d1..5f3b22da946ad594ae480824fc9208c9e1f38887 100644
--- a/web/content/docs/tutorials/advancing-glacier/mesh_basin.py
+++ b/web/content/docs/tutorials/advancing-glacier/mesh_basin.py
@@ -118,8 +118,8 @@ gmsh.model.geo.addCurveLoop([2, 3, 4, 10, 13, 14], 4)
 # loop defines the exterior contour; additional curve loop define holes.
 # (only one here, representing the external contour, since there are no holes
 # --see `t4.py' for an example of a surface with a hole):
-for l in range(1, 5):
-    gmsh.model.geo.addPlaneSurface([l], l)
+for L in range(1, 5):
+    gmsh.model.geo.addPlaneSurface([L], L)
 
 # At this level, Gmsh knows everything to display the surfaces and
 # to mesh it. An optional step is needed if we want to group elementary
diff --git a/web/content/docs/tutorials/advancing-glacier/timeBCs_glacier.py b/web/content/docs/tutorials/advancing-glacier/timeBCs_glacier.py
index 3860321dec9bb8c0dbb5041870ee7b48e19a3de4..a35006ce43a9b2e12a0505aa17c59f8cdd55ba99 100644
--- a/web/content/docs/tutorials/advancing-glacier/timeBCs_glacier.py
+++ b/web/content/docs/tutorials/advancing-glacier/timeBCs_glacier.py
@@ -22,7 +22,7 @@ class BC_Y(OpenGeoSys.BoundaryCondition):
         self.glacier.plotEvolution()
 
     def getFlux(
-        self, t, coords, primary_vars
+        self, t, coords, _primary_vars
     ):  # here Neumann BC: flux of linear momentum
         x, y, z = coords