You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Creating consistent 3D images of geophysical and geological datasets and turning that into an input model for geodynamic simulations is often challenging. The aim of this package is to help with this, by providing a number of routines to easily import data and create a consistent 3D visualisation from it in the VTK-toolkit format, which can for example be viewed with [Paraview](https://www.paraview.org). In addition, we provide a range of tools that helps to generate input models to perform geodynamic simulations and import the results of such simulations back into julia.
7
+
Creating consistent 3D images of geophysical and geological datasets and turning that into an input model for geodynamic simulations is often challenging. The aim of this package is to help with this, by providing a number of routines to easily import data and create a consistent 3D visualisation from it in the VTK-toolkit format, which can for example be viewed with [Paraview](https://www.paraview.org). In addition, we provide a range of tools that helps to generate input models to perform geodynamic simulations and import the results of such simulations back into julia.
- Grab screenshots of cross-sections or maps in published papers and view them in 3D (together with other data).
27
27
- Create a consistent overview that includes all available data of a certain region.
28
28
- Create initial model setups for the 3D geodynamic code [LaMEM](https://bitbucket.org/bkaus/lamem/src/master/).
29
-
- Import LaMEM timesteps.
29
+
- Import LaMEM timesteps.
30
30
31
31
All data is transformed into either a `GeoData` or a `UTMData` structure which contains info about `longitude/latitude/depth`, `ew/ns/depth` coordinates along with an arbitrary number of scalar/vector datasets, respectively. All data can be exported to Paraview with the `Write_Paraview` routine, which transfers the data to a `ParaviewData` structure (that contains Cartesian Earth-Centered-Earth-Fixed (ECEF) `x/y/z` coordinates, used for plotting)
32
-
33
-
## Usage
32
+
33
+
## Usage
34
34
The best way to learn how to use this is to install the package (see below) and look at the tutorials in the [manual](https://juliageodynamics.github.io/GeophysicalModelGenerator.jl/dev/).
35
35
36
-
## Installation
36
+
## Installation
37
37
First, you need to install julia on your machine. We recommend to use the binaries from [https://julialang.org](https://julialang.org).
38
38
Next, start julia and switch to the julia package manager using `]`, after which you can add the package.
39
39
```julia
@@ -59,7 +59,7 @@ We rely on a number of additional packages, which are all automatically installe
59
59
60
60
61
61
## Visualising Alpine data
62
-
We have used this package to interpret various data sets of the Alps (mostly openly available, sometimes derived from published papers). You can download the resulting paraview files here (using the `*.vts` format), where we also included the julia scripts to do the work (some of which are also described in more detail in the tutorials). Just unzip the files and open the corresponding `*.vts` in Paraview.
62
+
We have used this package to interpret various data sets of the Alps (mostly openly available, sometimes derived from published papers). You can download the resulting paraview files here (using the `*.vts` format), where we also included the julia scripts to do the work (some of which are also described in more detail in the tutorials). Just unzip the files and open the corresponding `*.vts` in Paraview.
@@ -68,12 +68,12 @@ If you want your data be included here as well, give us an email (or even better
68
68
You are very welcome to request new features and point out bugs by opening an issue. You can also help by adding features and creating a pull request.
69
69
70
70
## Development roadmap
71
-
In the pipeline:
71
+
In the pipeline:
72
72
- More tutorials
73
73
- Add more import tools.
74
74
- Compute gravity anomalies for lon/lat datasets, rather than just x/y/z.
75
-
- Provide an interface to [geomIO](https://bitbucket.org/geomio/geomio/wiki/Home) (currently being translated from MATLAB to python) in order to allow creating 3D geometric model setups by drawing in Inkscape.
76
-
- Provide tools to create and export 3D geodynamic model setups.
77
-
75
+
- Provide an interface to [geomIO](https://bitbucket.org/geomio/geomio/wiki/Home) (currently being translated from MATLAB to python) in order to allow creating 3D geometric model setups by drawing in Inkscape.
76
+
- Provide tools to create and export 3D geodynamic model setups.
77
+
78
78
## Funding
79
79
Development of this software package was funded by the German Research Foundation (DFG grants TH2076/7-1 and KA3367/10-1), which are part of the [SPP 2017 4DMB project](http://www.spp-mountainbuilding.de) project as well as by the European Research Council under grant ERC CoG #771143 - [MAGMA](https://magma.uni-mainz.de). The project was initiated at a [TeMaS](https://temas.uni-mainz.de) workshop with researchers from Frankfurt and Mainz where we realized that it is way too timeconsuming to collect available data of a certain region.
Copy file name to clipboardExpand all lines: docs/src/index.md
+1-1Lines changed: 1 addition & 1 deletion
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -9,7 +9,7 @@ Documentation for [GeophysicalModelGenerator](https://github.com/JuliaGeodynamic
9
9
The main purpose of this package is to simplify the process of going from 1D/2D/3D geophysical data to a 3D consistent model of the region. By simplifying the process of plotting the data, it becomes easier to compare different data sets, and generate a 3D models that can be used for other computations such as geodynamic simulations, or forward modelling of gravity anomalies.
10
10
11
11
For this, GeophysicalModelGenerator provides the following functionality:
12
-
- A consistent GeoData structure, that holds the data along with lon/lat/depth information.
12
+
- A consistent GeoData structure, that holds the data along with lon/lat/depth information.
13
13
- Routines to generate VTK files from the GeoData structure in order to visualize results in Paraview.
14
14
- The ability to deal with points, 2D profiles and 3D volumes, for both scalar and vector values.
15
15
- Rapidly import screenshots of published papers and compare them with other data sets in 3D using paraview.
The `LaMEM_grid` structure contains the number of elements in every direction and the number of markers in every cell.
191
191
It also contains `Grid.X`, `Grid.Y`, `Grid.Z`, which are the coordinates of each of the markers in the 3 directions.
192
192
193
-
In a next step we need to give each of these points a `Phase` number (which is an integer, that indicates the type of the rock that point has), as well as the temperature (in Celcius).
193
+
In a next step we need to give each of these points a `Phase` number (which is an integer, that indicates the type of the rock that point has), as well as the temperature (in Celsius).
194
194
195
195
````julia
196
196
Phases =ones(Int32,size(Grid.X))*2;
@@ -298,4 +298,3 @@ If you are interested in doing this, have a look at the LaMEM [wiki](https://bit
298
298
---
299
299
300
300
*This page was generated using [Literate.jl](https://github.com/fredrikekre/Literate.jl).*
Copy file name to clipboardExpand all lines: docs/src/man/installation.md
+10-14Lines changed: 10 additions & 14 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -7,7 +7,7 @@ In order to use then package you obviously need to install julia. We recommend d
7
7
8
8
9
9
### 2. Install Visual Studio Code
10
-
The julia files itself are text files (just like matlab scripts). You may want to edit or modify them at some stage, for which you can use any text editor for that. We prefer to use the freely available [Visual Studio Code](https://code.visualstudio.com) as it has a build-in terminal and is the comes with the (official) julia debugger (install the Julia extension for that).
10
+
The julia files itself are text files (just like matlab scripts). You may want to edit or modify them at some stage, for which you can use any text editor for that. We prefer to use the freely available [Visual Studio Code](https://code.visualstudio.com) as it has a built-in terminal and is the comes with the (official) julia debugger (install the Julia extension for that).
11
11
12
12
### 3. Getting started with julia
13
13
You start julia on the command line with:
@@ -25,19 +25,19 @@ This will start the command-line interface of julia:
25
25
_/ |\__'_|_|_|\__'_| | Official https://julialang.org/ release
26
26
|__/ |
27
27
28
-
julia>
28
+
julia>
29
29
```
30
30
31
31
From the julia prompt, you start the package manager by typing `]`:
32
32
```julia
33
-
(@v1.6) pkg>
33
+
(@v1.6) pkg>
34
34
```
35
35
And you return to the command line with a backspace.
36
36
37
-
Also useful is that julia has a build-in terminal, which you can reach by typing `;` on the command line:
37
+
Also useful is that julia has a built-in terminal, which you can reach by typing `;` on the command line:
38
38
```julia
39
39
julia>;
40
-
shell>
40
+
shell>
41
41
```
42
42
In the shell, you can use the normal commands like listing the content of a directory, or the current path:
43
43
```julia
@@ -48,7 +48,7 @@ shell> pwd
48
48
```
49
49
As before, return to the main command line (called `REPL`) with a backspace.
50
50
51
-
If you want to see help information for any julia function, type `?` followed by the command.
51
+
If you want to see help information for any julia function, type `?` followed by the command.
52
52
An example for`tan` is:
53
53
```julia
54
54
help?> tan
@@ -73,7 +73,7 @@ search: tan tanh tand atan atanh atand instances transpose transcode contains Un
73
73
2×2 Matrix{Float64}:
74
74
-1.09252 -1.09252
75
75
-1.09252 -1.09252
76
-
```
76
+
```
77
77
78
78
If you are in a directory that has a julia file (which have the extension `*.jl`), you can open that file with Visual Studio Code:
79
79
```julia
@@ -115,17 +115,13 @@ julia> using GeophysicalModelGenerator
115
115
```
116
116
117
117
### 5. Other useful packages
118
-
As you will work your way through the tutorials you will see that we often use external packages, for example to load ascii data files into julia. You will find detailed instructions in the respective tutorials.
118
+
As you will work your way through the tutorials you will see that we often use external packages, for example to load ascii data files into julia. You will find detailed instructions in the respective tutorials.
119
119
120
120
If you already want to install some of those, here our favorites. Install them through the package manager:
121
121
122
-
- [CSV](https://github.com/JuliaData/CSV.jl): Read comma-separated data files into julia.
123
-
- [Plots](https://github.com/JuliaPlots/Plots.jl): Create all kinds of plots in julia (quite an extensive package, but very useful to have).
122
+
- [CSV](https://github.com/JuliaData/CSV.jl): Read comma-separated data files into julia.
123
+
- [Plots](https://github.com/JuliaPlots/Plots.jl): Create all kinds of plots in julia (quite an extensive package, but very useful to have).
124
124
- [JLD2](https://github.com/JuliaIO/JLD2.jl): This allows saving julia objects (such as a tomographic model) to a binary file and load it again at a later stage.
125
125
- [Geodesy](https://github.com/JuliaGeo/Geodesy.jl): Convert UTM coordinates to latitude/longitude/altitude.
- [GMT](https://github.com/GenericMappingTools/GMT.jl): A julia interface to the Generic Mapping Tools (GMT), which is a highly popular package to create (geophysical) maps. Note that installing `GMT.jl` is more complicated than installing the other packages listed above, as you first need to have a working version of `GMT` on your machine (it is not yet installed automatically). Installation instructions for Windows/Linux are on their webpage. On a mac, we made the best experiences by downloading the binaries from their webpage and not using a package manager to install GMT.
Copy file name to clipboardExpand all lines: docs/src/man/lamem.md
+2-2Lines changed: 2 additions & 2 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -1,7 +1,7 @@
1
1
# LaMEM
2
2
3
-
In order to generate geodynamic simulations from setups created with `GeophysicalModelGenerator.jl`, we provide a few routines that directly create marker input files for the 3D geodynamic modelling software [LaMEM](https://bitbucket.org/bkaus/lamem), which is an open-source cartesian code that is well-suited to perform crustal and lithospheric-scale simulations.
4
-
If you want to learn how to run LaMEM simulations, please have a look at the [wiki page](https://bitbucket.org/bkaus/lamem/wiki/Home).
3
+
In order to generate geodynamic simulations from setups created with `GeophysicalModelGenerator.jl`, we provide a few routines that directly create marker input files for the 3D geodynamic modelling software [LaMEM](https://bitbucket.org/bkaus/lamem), which is an open-source cartesian code that is well-suited to perform crustal and lithospheric-scale simulations.
4
+
If you want to learn how to run LaMEM simulations, please have a look at the [wiki page](https://bitbucket.org/bkaus/lamem/wiki/Home).
5
5
6
6
The routines provided here have the following functionality:
7
7
- Read LaMEM *.dat files (to get the size of the domain)
0 commit comments