Polygon again

Author: TatjanaWeiler

docs/src/man/tutorial_Polygon_structures.md

@@ -40,12 +40,12 @@ AddBox!(Phase, Temp, Cart; xlim=(0.0,800.0),ylim=(0.0,800.0), zlim=(0.0,50.0), p
  ```
 
 
-#### 2. Add polygon struture
+#### 2. Add polygon structure
 For including the geological structures of a passive margin into the model via the polygon function, depths down to 150 km are focused on. With the polygon function, it is possible to create different shapes. In the example, the sediment basin shows a more trapezial (2D in x-/z-direction) shape, while the thinning of the plate is a more triangular (2D in x-/z-direction). More complex structures are possible to build in the background model due to the non-limitation number of points in the x- and z-direction. In y-direction only the length can be varied and is set by two values. The shape is not changeable. The x- and z-values of the points need to be in the same order for selecting the correct point (P1(1/3), P2(2/2) --> xlim(1,2), ylim(3,2)).
 
 
 ```julia
-# xlim: x-coordiates of the points, same ordering as zlim
+# xlim: x-coordinates of the points, same ordering as zlim
 # zlim: z-coordinates of the points, same ordering as xlim
 # ylim: limits the object within the two ylim values
 # unlimited number of points possible to create the polygon

src/Setup_geometry.jl

@@ -485,9 +485,9 @@ Parameters
 - Phase - Phase array (consistent with Grid)
 - Temp  - Temperature array (consistent with Grid)
 - Grid  - Grid structure (usually obtained with ReadLaMEM_InputFile)
-- xlim  - x-cooridnate of the polygon points, same ordering as zlim, number of points unlimited
-- ylim  - y-cooridnate, limition in length possible (two values (start and stop))
-- zlim  - z-cooridnate of the polygon points, same ordering as xlim, number of points unlimited
+- xlim  - x-coordinate of the polygon points, same ordering as zlim, number of points unlimited
+- ylim  - y-coordinate, limitation in length possible (two values (start and stop))
+- zlim  - z-coordinate of the polygon points, same ordering as xlim, number of points unlimited
 - phase - specifies the phase of the box. See `ConstantPhase()`
 - T - specifies the temperature of the box. See `ConstantTemp()`,`LinearTemp()`,`HalfspaceCoolingTemp()`,`SpreadingRateTemp()`
 
@@ -543,7 +543,7 @@ end
 
 # Compute thermal structure accordingly. See routines below for different options
 if T != nothing
-Temp[ind] = Compute_ThermalStructure(Temp[ind], X[ind], Y[ind], Z[ind], T)
+Temp[ind] = Compute_ThermalStructure(Temp[ind], X[ind], Y[ind], Z[ind], Phase[ind], T)
 end
 
 # Set the phase. Different routines are available for that - see below.

tutorials/Tutorial_polygon_geometry.jl

@@ -21,7 +21,7 @@ Temp    = ones(Float64,size(X))*1350
 AddBox!(Phase, Temp, Cart; xlim=(0.0,800.0),ylim=(0.0,800.0), zlim=(-800.0,0.0), phase = LithosphericPhases(Layers=[15 30 100 800], Phases=[2 3 1 5], Tlab=1300 ), T=LinearTemp(Ttop=20, Tbot=1600) )#T=HalfspaceCoolingTemp(Tsurface=20.0, Tmantle=1350, Age=120, Adiabat=0.4)
 
 
-# xlim: x-coordiates of the points, same ordering as zlim
+# xlim: x-coordinates of the points, same ordering as zlim
 # zlim: z-coordinates of the points, same ordering as xlim
 # ylim: limits the object within the two ylim values
 # unlimited number of points possible to create the polygon