keep column bottom values pinned at basalwater compositions

grahamedwards · Mar 1, 2024 · 9300abf · 9300abf
1 parent 36521c4
commit 9300abf
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/diffuse.jl b/src/diffuse.jl
@@ -125,9 +125,9 @@ function diffuseadvectcolumn!(sc::SedimentColumn, k::Constants)
     end
 
 # set bottom values to penultimate values
-    sc.O.p[k.nz] = sc.O.p[k.penultimate_node]
-    sc.Cl.p[k.nz] = sc.Cl.p[k.penultimate_node]
-    sc.rho.p[k.nz] = sc.rho.p[k.penultimate_node]
+    #sc.O.p[k.nz] = sc.O.p[k.penultimate_node]
+    #sc.Cl.p[k.nz] = sc.Cl.p[k.penultimate_node]
+    #sc.rho.p[k.nz] = sc.rho.p[k.penultimate_node]
 
 # And time steps forward, replacing o with p.
     sc.O.o .= sc.O.p