Merge branch 'ap-mckenzie-temp' of https://github.com/APiccolo89/GeophysicalModelGenerator.jl into ap-mckenzie-temp

Author: boriskaus

docs/src/man/geodynamic_setups.md

@@ -20,4 +20,5 @@ GeophysicalModelGenerator.LithosphericTemp
 GeophysicalModelGenerator.ConstantPhase
 GeophysicalModelGenerator.Compute_Phase
 GeophysicalModelGenerator.LithosphericPhases
+GeophysicalModelGenerator.McKenzie_subducting_slab
 ```

src/Setup_geometry.jl

@@ -1063,7 +1063,18 @@ end
 """ 
     Compute_ThermalStructure(Temp, X, Y, Z, s::McKenzie_subducting_slab)
 
-Compute the temperature field of a `McKenzie_subducting_slab` scenario 
+Compute the temperature field of a `McKenzie_subducting_slab`. Uses the analytical solution
+of McKenzie (1969) ["Speculations on the consequences and causes of plate motions"]. The functions assumes
+that the bottom of the slab is the coordinate Z=0. Internally the function shifts the coordinate. 
+Parameters
+=============================
+Temp Temperature array
+X    X Array 
+Y    Y Array 
+Z    Z Array 
+Phase Phase array 
+s    McKenzie_subducting_slab
+
 """
 function Compute_ThermalStructure(Temp, X, Y, Z,Phase, s::McKenzie_subducting_slab)
     @unpack Tsurface, Tmantle, Adiabat, v_cm_yr, κ, it = s