Making it time stable

Author: TatjanaWeiler

src/Setup_geometry.jl

@@ -639,13 +639,9 @@ function add_polygon!(Phase, Temp, Grid::AbstractGeneralGrid;   # required input
     T=nothing, cell=false )                                     # Sets the thermal structure (various functions are available)
 
 
-    xlim = collect(xlim)
-    ylim = collect(ylim)
-    zlim = collect(zlim)
-    
-    xlim = Float64.(xlim)
-    ylim = Float64.(ylim)
-    zlim = Float64.(zlim)
+    xlim_ = Float64.(collect(xlim))
+    ylim_ = Float64.(collect(ylim))
+    zlim_ = Float64.(collect(zlim))
 
 
 # Retrieve 3D data arrays for the grid
@@ -656,8 +652,8 @@ ind_slice = zeros(Bool,size(X[:,1,:]))
 
 # find points within the polygon, only in 2D
 for i = 1:size(Y)[2]
-    if Y[1,i,1] >= ylim[1] && Y[1,i,1]<=ylim[2] 
-        inpolygon!(ind_slice, xlim,zlim, X[:,i,:], Z[:,i,:])
+    if Y[1,i,1] >= ylim_[1] && Y[1,i,1]<=ylim_[2] 
+        inpolygon!(ind_slice, xlim_,zlim_, X[:,i,:], Z[:,i,:])
         ind[:,i,:] = ind_slice
     else
         ind[:,i,:] = zeros(size(X[:,1,:]))

test/test_setup_geometry.jl

@@ -284,6 +284,7 @@ TsHC = HalfspaceCoolingTemp(Tsurface=20.0, Tmantle=1350, Age=120, Adiabat=0.4)
 TsMK = McKenzie_subducting_slab(Tsurface = 20.0, Tmantle = 1350.0, v_cm_yr = 4.0, Adiabat = 0.0)
 T_slab = LinearWeightedTemperature(crit_dist=600, F1=TsHC, F2=TsMK);
 phase = LithosphericPhases(Layers=[5 7 88], Phases = [2 3 4], Tlab=nothing)
+
 t1 = Trench(Start = (400.0,400.0), End = (800.0,800.0),θ_max = 90.0, direction = 1.0, n_seg = 50, Length = 600.0, Thickness = 80.0, Lb = 500.0,d_decoupling = 100.0, type_bending =:Ribe,   WeakzoneThickness=10, WeakzonePhase=9)
 
 add_slab!(Phase,Temp,Cart, t1, phase=phase, T = T_slab)