Merge pull request #98 from APiccolo89/bk-geodynamics-tutorial_1

Small additional note to the "Tutorial_NumericalModel_3D.jl" []

Author: boriskaus

src/Setup_geometry.jl

@@ -201,6 +201,13 @@ function AddBox!(Phase, Temp, Grid::AbstractGeneralGrid;                 # requi
         Origin = (xlim[1], ylim[1], zlim[2])  # upper-left corner
     end
 
+    if Origin !== nothing && isa(T, McKenzie_subducting_slab)
+        @warn  "McKenzie temperature does not require the definition of 'Origin' field; if Origin is defined it must be equal to [xmin,ymin,zmax] of the box that has been defined."
+        if Origin[1] != xlim[1] || Origin[2] != ylim[1] || Origin[3] != zlim[2]
+            @error  "Origin is not set up correctly. For fixing the problem Origin can be left blank or Origin = [xmin,ymin,zmax] of the box"
+        end
+    end
+
     # Perform rotation of 3D coordinates:
     Xrot = X .- Origin[1];
     Yrot = Y .- Origin[2];

tutorials/Tutorial_NumericalModel_3D.jl

@@ -33,6 +33,10 @@ Temp = fill(1350.0, nx,ny,nz);
 
 # We use a lithospheric structure, which can be specified with the `LithosphericPhases` structure, where you can indicate `Layers` (the thickness of each lithospheric layer, starting from the top), and `Phases` the phases of the corresponding layers.
 # Note that if the lowermost layer has the same phase as the mantle, you can define `Tlab` as the lithosphere-asthenosphere boundary which will automatically adjust the phase depending on temperature
+
+# Also note that `Origin` must be set at [`xmin`,`ymin`,`zmax`] of the box, or should be left out if you use this along with the `McKenzie_subducting_slab` thermal structure.
+
+
 lith = LithosphericPhases(Layers=[15 45 10], Phases=[0 1 2])
 v_spreading_cm_yr = 3
 AddBox!(Phases, Temp, Grid; xlim=(-800,0.0), ylim=(-400, 400.0), zlim=(-80.0, 0.0), phase = lith,