update tutorials

Author: boriskaus

docs/src/man/Tutorial_Basic.md

@@ -204,10 +204,10 @@ In creating this image, I used the `Clip` tool of Paraview to only show topograp
 ### 4. Cartesian data
 As you can see, the curvature or the Earth is taken into account here. Yet, for many applications it is more convenient to work in Cartesian coordinates (kilometers) rather then in geographic coordinates.
 `GeophysicalModelGenerator` has a number of tools for this.
-First we need do define a `ProjectionPoint`  around which we project the data
+First we need do define a `projectionPoint`  around which we project the data
 
 ```julia
-proj = ProjectionPoint(Lon=12.0,Lat =43)
+proj = projectionPoint(Lon=12.0,Lat =43)
 
 Topo_cart = convert2CartData(Topo_Alps, proj)
 ```

docs/src/man/Tutorial_Jura.md

@@ -100,10 +100,10 @@ Moreover, the resolution of the grids is different. Whereas the `TopoGeology` ha
 It is often useful to have them on exactly the same size grid
 
 We can do this in two steps:
-First, we define a `ProjectionPoint` along which we perform the projection
+First, we define a `projectionPoint` along which we perform the projection
 
 ```julia
-proj = ProjectionPoint(Lon=6, Lat=46.5)
+proj = projectionPoint(Lon=6, Lat=46.5)
 ```
 
 We can simply transfer the TopoGeology map to Cartesian values with:

docs/src/man/Tutorial_LaPalma.md

@@ -65,7 +65,7 @@ In order to create model setups, it is helpful to first transfer the data to Car
 This requires us to first determine a *projection point*, that is fixed. Often, it is helpful to use the center of the topography for this. In the present example, we will center the model around La Palma itself:
 
 ```julia
-proj = ProjectionPoint(Lon=-17.84, Lat=28.56)
+proj = projectionPoint(Lon=-17.84, Lat=28.56)
 ```
 
 Once this is done you can convert the topographic data to the cartesian reference frame

docs/src/man/Tutorial_NumericalModel_2D.md

@@ -17,7 +17,7 @@ using GeophysicalModelGenerator
 nx,nz = 512,128
 x = range(-1000,1000, nx);
 z = range(-660,0,    nz);
-Grid2D = CartData(XYZGrid(x,0,z))
+Grid2D = CartData(xyzGrid(x,0,z))
 ```
 
 ````
@@ -47,13 +47,13 @@ Temp = fill(1350.0, nx, 1, nz);
 We will start with a simple subduction setup, which consists of a horizontal part:
 
 ```julia
-AddBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = ConstantPhase(1));
+addBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = ConstantPhase(1));
 ```
 
 And with the inclined part:
 
 ```julia
-AddBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = ConstantPhase(1), DipAngle=30);
+addBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = ConstantPhase(1), DipAngle=30);
 ```
 
 Add them to the `CartData` dataset:
@@ -75,7 +75,7 @@ CartData
 Which looks like
 
 ```julia
-Write_Paraview(Grid2D,"Grid2D_SubductionMechanical");
+write_Paraview(Grid2D,"Grid2D_SubductionMechanical");
 ```
 
 ````
@@ -100,15 +100,15 @@ LithosphericPhases([15 55], [1 2], nothing)
 and set the slab again:
 
 ```julia
-AddBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith);
-AddBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = lith, DipAngle=30);
+addBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith);
+addBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = lith, DipAngle=30);
 ```
 
 Which looks like:
 
 ```julia
 Grid2D = addField(Grid2D,(;Phases, Temp))
-Write_Paraview(Grid2D,"Grid2D_SubductionMechanicalLayered");
+write_Paraview(Grid2D,"Grid2D_SubductionMechanicalLayered");
 ```
 
 ````
@@ -124,15 +124,15 @@ We can do that by specifying a thermal structure. For example, we can use the ha
 
 ```julia
 therm = HalfspaceCoolingTemp(Age=40)
-AddBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith, T=therm);
-AddBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = lith, T = therm, DipAngle=30);
+addBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith, T=therm);
+addBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = lith, T = therm, DipAngle=30);
 ```
 
 Which looks like:
 
 ```julia
 Grid2D = addField(Grid2D,(;Phases, Temp))
-Write_Paraview(Grid2D,"Grid2D_SubductionHalfspaceCooling");
+write_Paraview(Grid2D,"Grid2D_SubductionHalfspaceCooling");
 ```
 
 ````
@@ -159,13 +159,13 @@ a spreading velocity (note that this simply relates to the thermal structure and
 
 ```julia
 lith = LithosphericPhases(Layers=[15 55], Phases=[1 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
+addBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
 ```
 
 For the subduction we use a thermal structure of a slab heated by hot asthenosphere
 
 ```julia
-AddBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = lith, T = McKenzie_subducting_slab(Tsurface=0,v_cm_yr=3), DipAngle=30);
+addBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 0.0), phase = lith, T = McKenzie_subducting_slab(Tsurface=0,v_cm_yr=3), DipAngle=30);
 ```
 
 We can set the mantle lithosphere that is hotter > 1250 C to mantle:
@@ -175,7 +175,7 @@ ind = findall(Temp .> 1250 .&& Phases .==2);
 Phases[ind] .= 0;
 
 Grid2D = addField(Grid2D,(;Phases, Temp))
-Write_Paraview(Grid2D,"Grid2D_SubductionRidge");
+write_Paraview(Grid2D,"Grid2D_SubductionRidge");
 ```
 
 ````
@@ -190,21 +190,21 @@ Ok, lets add an overriding slab as well. For this, we use the `AddLayer!` functi
 
 ```julia
 lith = LithosphericPhases(Layers=[15 20 55], Phases=[3 4 5], Tlab=1250)
-AddBox!(Phases, Temp, Grid2D; xlim=(0,1000), zlim=(-80.0, 0.0), phase = lith, T=HalfspaceCoolingTemp(Age=80));
+addBox!(Phases, Temp, Grid2D; xlim=(0,1000), zlim=(-80.0, 0.0), phase = lith, T=HalfspaceCoolingTemp(Age=80));
 ```
 
 The oceanic plate is as before
 
 ```julia
 lith = LithosphericPhases(Layers=[15 55], Phases=[1 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
+addBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
 ```
 
 For the inclined part, we set a layer above the slab (the "weak" layer to facilitate subduction initiation )
 
 ```julia
 lith = LithosphericPhases(Layers=[10 15 55], Phases=[6 1 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 10.0), phase = lith, T = McKenzie_subducting_slab(Tsurface=0,v_cm_yr=3), DipAngle=30);
+addBox!(Phases, Temp, Grid2D; xlim=(0,300), zlim=(-80.0, 10.0), phase = lith, T = McKenzie_subducting_slab(Tsurface=0,v_cm_yr=3), DipAngle=30);
 ```
 
 Lithosphere-asthenosphere boundary:
@@ -214,7 +214,7 @@ ind = findall(Temp .> 1250 .&& Phases .==2);
 Phases[ind] .= 0;
 
 Grid2D = addField(Grid2D,(;Phases, Temp))
-Write_Paraview(Grid2D,"Grid2D_SubductionOverriding");
+write_Paraview(Grid2D,"Grid2D_SubductionOverriding");
 ```
 
 ````
@@ -225,18 +225,18 @@ Saved file: Grid2D_SubductionOverriding.vts
 ![Mechanical2D_Tutorial_5](../assets/img/Mechanical2D_Tutorial_5.png)
 
 #### Curved slab - mechanics
-So far, the subducting part of the slab was always straight. We can also create a curved slab by using the `AddSlab!` function. This uses a parametric representation of the slab and is a bit more involved than the `AddBox!` function.
+So far, the subducting part of the slab was always straight. We can also create a curved slab by using the `addSlab!` function. This uses a parametric representation of the slab and is a bit more involved than the `addBox!` function.
 
 We start with the horizontal part:
 
 ```julia
 nx,nz = 512,128
 x = range(-1000,1000, nx);
 z = range(-660,0,    nz);
-Grid2D = CartData(XYZGrid(x,0,z))
+Grid2D = CartData(xyzGrid(x,0,z))
 Phases = zeros(Int64, nx, 1, nz);
 Temp = fill(1350.0, nx, 1, nz);
-AddBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = ConstantPhase(1));
+addBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-80.0, 0.0), phase = ConstantPhase(1));
 ```
 
 Next, we should define a `Trench` structure, which contains info about the trench which goes in 3D from `Start` - `End` coordinates (`x`,`y`)-coordinates respectively. As we are dealing with a 2D model, we set the `y`-coordinates to -100.0 and 100.0 respectively.
@@ -251,7 +251,7 @@ Add them to the `CartData` dataset:
 
 ```julia
 Grid2D = addField(Grid2D,(;Phases, Temp))
-Write_Paraview(Grid2D,"Grid2D_SubductionCurvedMechanical");
+write_Paraview(Grid2D,"Grid2D_SubductionCurvedMechanical");
 ```
 
 ````
@@ -271,7 +271,7 @@ Our starting basis is the example above with ridge and overriding slab
 nx,nz = 512,128
 x = range(-1000,1000, nx);
 z = range(-660,0,    nz);
-Grid2D = CartData(XYZGrid(x,0,z))
+Grid2D = CartData(xyzGrid(x,0,z))
 Phases = zeros(Int64, nx, 1, nz);
 Temp = fill(1350.0, nx, 1, nz);
 lith = LithosphericPhases(Layers=[15 20 55], Phases=[3 4 5], Tlab=1250)
@@ -284,16 +284,16 @@ LithosphericPhases([15 20 55], [3 4 5], 1250)
 Lets start with defining the horizontal part of the overriding plate. Note that we define this twice with different thickness to deal with the bending subduction area:
 
 ```julia
-AddBox!(Phases, Temp, Grid2D; xlim=(200,1000), zlim=(-150.0, 0.0), phase = lith, T=HalfspaceCoolingTemp(Age=80));
-AddBox!(Phases, Temp, Grid2D; xlim=(0,200), zlim=(-50.0, 0.0), phase = lith, T=HalfspaceCoolingTemp(Age=80));
+addBox!(Phases, Temp, Grid2D; xlim=(200,1000), zlim=(-150.0, 0.0), phase = lith, T=HalfspaceCoolingTemp(Age=80));
+addBox!(Phases, Temp, Grid2D; xlim=(0,200), zlim=(-50.0, 0.0), phase = lith, T=HalfspaceCoolingTemp(Age=80));
 ```
 
 The horizontal part of the oceanic plate is as before:
 
 ```julia
 v_spread_cm_yr = 3      #spreading velocity
 lith = LithosphericPhases(Layers=[15 55], Phases=[1 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-150.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=v_spread_cm_yr));
+addBox!(Phases, Temp, Grid2D; xlim=(-800,0.0), zlim=(-150.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=v_spread_cm_yr));
 ```
 
 Yet, now we add a trench as well. The starting thermal age at the trench is that of the horizontal part of the oceanic plate:
@@ -331,7 +331,7 @@ ind = findall(Temp .> 1250 .&& (Phases.==2 .|| Phases.==5));
 Phases[ind] .= 0;
 
 Grid2D = addField(Grid2D,(;Phases, Temp))
-Write_Paraview(Grid2D,"Grid2D_SubductionCurvedOverriding");
+write_Paraview(Grid2D,"Grid2D_SubductionCurvedOverriding");
 ```
 
 ````

docs/src/man/Tutorial_NumericalModel_3D.md

@@ -18,7 +18,7 @@ nx,ny,nz = 512,512,128
 x = range(-1000,1000, nx);
 y = range(-1000,1000, ny);
 z = range(-660,0,    nz);
-Grid = CartData(XYZGrid(x,y,z));
+Grid = CartData(xyzGrid(x,y,z));
 ```
 
 Now we create an integer array that will hold the `Phases` information (which usually refers to the material or rock type in the simulation)
@@ -43,15 +43,15 @@ Note that if the lowermost layer has the same phase as the mantle, you can defin
 
 ```julia
 lith = LithosphericPhases(Layers=[15 45 10], Phases=[0 1 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid; xlim=(-800,0.0), ylim=(-400, 400.0), zlim=(-80.0, 0.0), phase = lith,
+addBox!(Phases, Temp, Grid; xlim=(-800,0.0), ylim=(-400, 400.0), zlim=(-80.0, 0.0), phase = lith,
         Origin=(-0,0,0),
         T=SpreadingRateTemp(SpreadingVel=3, MORside="right"), StrikeAngle=30);
 ```
 
 And an an inclined part:
 
 ```julia
-AddBox!(Phases, Temp, Grid; xlim=(0,300), ylim=(-400, 400.0), zlim=(-80.0, 0.0), phase = lith,
+addBox!(Phases, Temp, Grid; xlim=(0,300), ylim=(-400, 400.0), zlim=(-80.0, 0.0), phase = lith,
         Origin=(-0,0,0),
         T=McKenzie_subducting_slab(Tsurface=0,v_cm_yr=3), DipAngle=30, StrikeAngle=30);
 ```
@@ -75,7 +75,7 @@ CartData
 Which looks like
 
 ```julia
-Write_Paraview(Grid,"Grid3D_FreeSubduction");
+write_Paraview(Grid,"Grid3D_FreeSubduction");
 ```
 
 ````
@@ -94,7 +94,7 @@ nx,ny,nz = 512,512,128
 x = range(-1000,1000, nx);
 y = range(-1000,1000, ny);
 z = range(-660,0,    nz);
-Grid = CartData(XYZGrid(x,y,z));
+Grid = CartData(xyzGrid(x,y,z));
 
 Phases = fill(2,nx,ny,nz);
 Temp = fill(1350.0, nx,ny,nz);
@@ -104,27 +104,27 @@ Overriding plate with a 30 km crust and mantle lithosphere that where T<1250 cel
 
 ```julia
 lith_cont = LithosphericPhases(Layers=[30 200 50], Phases=[3 4 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid; xlim=(400,1000), ylim=(-1000, 0.0), zlim=(-240.0, 0.0), phase = lith_cont, T=HalfspaceCoolingTemp(Age=150));
-AddBox!(Phases, Temp, Grid; xlim=(200,1000), ylim=(-1000, 0.0), zlim=(-80.0, 0.0), phase = lith_cont,  T=HalfspaceCoolingTemp(Age=150));
+addBox!(Phases, Temp, Grid; xlim=(400,1000), ylim=(-1000, 0.0), zlim=(-240.0, 0.0), phase = lith_cont, T=HalfspaceCoolingTemp(Age=150));
+addBox!(Phases, Temp, Grid; xlim=(200,1000), ylim=(-1000, 0.0), zlim=(-80.0, 0.0), phase = lith_cont,  T=HalfspaceCoolingTemp(Age=150));
 
 lith_cont = LithosphericPhases(Layers=[30 200 10], Phases=[5 6 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid; xlim=(400,1000), ylim=(0, 1000),    zlim=(-240.0, 0.0), phase = lith_cont, T=HalfspaceCoolingTemp(Age=200));
-AddBox!(Phases, Temp, Grid; xlim=(200,1000), ylim=(0, 1000),    zlim=( -80.0, 0.0), phase = lith_cont, T=HalfspaceCoolingTemp(Age=200));
+addBox!(Phases, Temp, Grid; xlim=(400,1000), ylim=(0, 1000),    zlim=(-240.0, 0.0), phase = lith_cont, T=HalfspaceCoolingTemp(Age=200));
+addBox!(Phases, Temp, Grid; xlim=(200,1000), ylim=(0, 1000),    zlim=( -80.0, 0.0), phase = lith_cont, T=HalfspaceCoolingTemp(Age=200));
 ```
 
 Define an oceanic plate with ridge
 
 ```julia
 v_spread_cm_yr = 3      #spreading velocity
 lith = LithosphericPhases(Layers=[15 45 10], Phases=[0 1 2], Tlab=1250)
-AddBox!(Phases, Temp, Grid; xlim=(-800 , 200), ylim=(-1000, -400.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
-AddBox!(Phases, Temp, Grid; xlim=(-1000,-800), ylim=(-1000, -400.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3,MORside="right"));
+addBox!(Phases, Temp, Grid; xlim=(-800 , 200), ylim=(-1000, -400.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
+addBox!(Phases, Temp, Grid; xlim=(-1000,-800), ylim=(-1000, -400.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3,MORside="right"));
 
-AddBox!(Phases, Temp, Grid; xlim=(-700,  200), ylim=(-400, 200.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
-AddBox!(Phases, Temp, Grid; xlim=(-1000,-700), ylim=(-400, 200.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3,MORside="right"));
+addBox!(Phases, Temp, Grid; xlim=(-700,  200), ylim=(-400, 200.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
+addBox!(Phases, Temp, Grid; xlim=(-1000,-700), ylim=(-400, 200.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3,MORside="right"));
 
-AddBox!(Phases, Temp, Grid; xlim=(-650,  200), ylim=(200, 1000.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
-AddBox!(Phases, Temp, Grid; xlim=(-1000,-650), ylim=(200, 1000.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3,MORside="right"));
+addBox!(Phases, Temp, Grid; xlim=(-650,  200), ylim=(200, 1000.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3));
+addBox!(Phases, Temp, Grid; xlim=(-1000,-650), ylim=(200, 1000.0), zlim=(-80.0, 0.0), phase = lith, T=SpreadingRateTemp(SpreadingVel=3,MORside="right"));
 ```
 
 Subducting parts of the oceanic plate
@@ -163,7 +163,7 @@ Finally, we can add all this to the `CartData` dataset:
 
 ```julia
 Grid = addField(Grid,(;Phases, Temp))
-Write_Paraview(Grid,"Grid3D_Ridges");
+write_Paraview(Grid,"Grid3D_Ridges");
 ```
 
 ````

docs/src/man/datastructures.md

@@ -1,7 +1,7 @@
 # Data structures
 
 The main data structure used in GeophysicalModelGenerator.jl is `GeoData`, which contains info about the `longitude`,`latitude`, and `depth` of a data set, as well as several data sets itself.
-We also provide a `UTMData`, which is essentially the same but with UTM coordinates, and a `CartData` structure, which has Cartesian coordinates in kilometers (as used in many geodynamic codes). If one wishes to transfer `GeoData` to `CartData`, one needs to provide a `ProjectionPoint`.
+We also provide a `UTMData`, which is essentially the same but with UTM coordinates, and a `CartData` structure, which has Cartesian coordinates in kilometers (as used in many geodynamic codes). If one wishes to transfer `GeoData` to `CartData`, one needs to provide a `projectionPoint`.
 For plotting, we transfer this into the `ParaviewData` structure, which has cartesian coordinates centered around the center of the Earth. We employ the `wgs84` reference ellipsoid as provided by the [Geodesy.jl](https://github.com/JuliaGeo/Geodesy.jl) package to perform this transformation. 
 
 ```@docs
@@ -11,5 +11,5 @@ UTMData
 ParaviewData
 lonlatdepthGrid
 xyzGrid
-ProjectionPoint
+projectionPoint
 ```

docs/src/man/projection.md

@@ -4,7 +4,7 @@ Typically, you load a dataset by reading it into julia and either generating a `
 
 If you write the data to `Paraview`, it is internally converted to a Paraview structure (which involves `x,y,z` Cartesian Earth-Centered-Earth-Fixed (ECEF) coordinates using the `wgs84` ellipsoid). 
 
-Yet, if you do geodynamic calculations the chances are that the geodynamic code does not operate in spherical coordinates, but rather use cartesian ones. In that case you should transfer your data to the `CartData` structure, which requires you to specify a `ProjectionPoint` that is a point on the map that will later have the coordinates `(0,0)` in the `CartData` structure.
+Yet, if you do geodynamic calculations the chances are that the geodynamic code does not operate in spherical coordinates, but rather use cartesian ones. In that case you should transfer your data to the `CartData` structure, which requires you to specify a `projectionPoint` that is a point on the map that will later have the coordinates `(0,0)` in the `CartData` structure.
 
 
 #### 1. Converting
@@ -41,8 +41,8 @@ As the area is large, it covers a range of `UTM` zones (and every point has a UT
 Yet, what we could do instead is show all data with respect to a single UTM zone. For this, we have to select a point around which we project (in this case more or less in the center):
 
 ```julia
-julia> p=ProjectionPoint(Lon=17.3, Lat=37.5)
-ProjectionPoint(37.5, 17.3, 703311.4380385976, 4.152826288024972e6, 33, true)
+julia> p=projectionPoint(Lon=17.3, Lat=37.5)
+projectionPoint(37.5, 17.3, 703311.4380385976, 4.152826288024972e6, 33, true)
 ```
 
 Projecting the `GeoData` set using this projection point is done with:

docs/src/man/visualise.md

@@ -36,8 +36,8 @@ Which will look like:
 
 This is an example where we used `GeoData` to visualize results. Alternatively, we can also visualize results in km (often more useful for numerical modelling setups). For `Visualize` to work with this, we however need orthogonal cartesian data, which can be obtained by projecting both the data.
 ```julia
-julia> p=ProjectionPoint(Lon=10, Lat=45)
-ProjectionPoint(45.0, 10.0, 578815.302916711, 4.983436768349297e6, 32, true)
+julia> p=projectionPoint(Lon=10, Lat=45)
+projectionPoint(45.0, 10.0, 578815.302916711, 4.983436768349297e6, 32, true)
 julia> Data_Cart  = CartData(xyzGrid(-600:10:600,-600:10:600,-1000:10:-1));
 julia> Topo_Cart  = CartData(xyzGrid(-600:10:600,-600:10:600,0));
 julia> Topo_Cart  = projectCartData(Topo_Cart, Topo, p)

src/Setup_geometry.jl

@@ -174,7 +174,7 @@ julia> write_Paraview(Model3D,"LaMEM_ModelSetup")           # Save model to para
 
 Example 2) Box with halfspace cooling profile
 ```julia
-julia> Grid = CartData(XYZGrid(-1000:10:1000,0,-660:10:0))
+julia> Grid = CartData(xyzGrid(-1000:10:1000,0,-660:10:0))
 julia> Phases = zeros(Int32,   size(Grid));
 julia> Temp   = zeros(Float64, size(Grid));
 julia> addBox!(Phases,Temp,Grid, xlim=(0,500), zlim=(-50,0), phase=ConstantPhase(3), DipAngle=10, T=HalfspaceCoolingTemp(Age=30))