Fix indexing in free_surface_bcs! (#135)

* fix indexing in free_surface_bcs! * fix another incorrect index
PTsolvers · Apr 27, 2024 · 2117f77 · 2117f77
1 parent f4518bc
commit 2117f77
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/boundaryconditions/BoundaryConditions.jl b/src/boundaryconditions/BoundaryConditions.jl
@@ -212,7 +212,7 @@ end
 function free_surface_bcs!(
     stokes, bcs::FlowBoundaryConditions, η, rheology, phase_ratios, dt, di
 )
-    indices_range(::Any, Vy) = @idx (size(Vy, 2) - 1)
+    indices_range(::Any, Vy) = @idx (size(Vy, 1) - 2)
     indices_range(::Any, ::Any, Vz) = @idx (size(Vz, 1) - 2, size(Vz, 2) - 2)
 
     V = @velocity(stokes)