Creating Cross Sections

In scientific publications that deal with seismic tomographies, the results are often shown in 2D horizontal and vertical cross sections (slices and profiles). The Geophysical Model Generator therefore also provides tools do exactly this.

We will use the data that we have previously imported to show how to create cross sections. If you haven’t closed julia, these datasets should still be in memory, otherwise you have to load the corresponding jld2 file:

julia> load("Paffrath_Europe.jld2")

1. Creating horizontal cross sections (slices)

Once you have loaded the required data, creating a horizontal slice is quite easy using the GMG tool CrossSection. To extract a horizontal slice at a given depth, simply call the function as in the example below:

julia> Data_Paffrath_d100  =   CrossSection(Data_Paffrath, Depth_level=-100km)  

where the first argument to the function is the GMG data structure and the second argument denotes the kind of cross section to be extracted. The key word Depth_level denotes that a horizontal slice is to be extracted, with its depth value being -100 km (note that depth is negative).

2. Creating vertical cross sections (profiles)

Creating a vertical cross section is similarly straightforward, the only difference being that in this case, the starting and ending point of the cross section has to be provided:

julia> Data_Paffrath_A  =  CrossSection(Data_Paffrath, Start=(4.65,45.73), End=(17.23,43.80))

Where the start and end points have to be given as (longitude,latitude) pairs.

• Previous: Data Import

• Next: Exporting data