bugfixes

Author: boriskaus

src/Setup_geometry.jl

@@ -1070,7 +1070,7 @@ function Compute_ThermalStructure(Temp, X, Y, Z,Phase, s::McKenzie_subducting_sl
 
     # Thickness of the layer: 
     D0          =   (maximum(Z)-minimum(Z));
-    Zshift      =   Z .- Z[1]       # McKenzie model is defined with Z = 0 at the bottom of the slab
+    Zshift      =   Z .- Z[end]       # McKenzie model is defined with Z = 0 at the bottom of the slab
 
     # Convert subduction velocity from cm/yr -> m/s; 
     convert_velocity = 1/(100.0*365.25*60.0*60.0*24.0);

test/test_setup_geometry.jl

@@ -193,7 +193,7 @@ AddBox!(Phase, Temp, Cart; xlim=(0.0,600.0),ylim=(0.0,600.0), zlim=(-80.0, 0.0),
 # inclined slab
 Temp    = ones(Float64,size(Cart))*1350;
 AddBox!(Phase, Temp, Cart; xlim=(0.0,600.0),ylim=(0.0,600.0), zlim=(-80.0,0),StrikeAngle=0, DipAngle=45, phase = ConstantPhase(5), T=TsMK);
-@test sum(Temp)  ≈ 3.5133669349123573e8
+@test sum(Temp)  ≈ 3.5125017626287365e8