Skip to content

Commit

Permalink
Update Tutorial_AlpineData.jl
Browse files Browse the repository at this point in the history
  • Loading branch information
mthielma committed Feb 29, 2024
1 parent 6ee0479 commit 3592b9a
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions tutorials/Tutorial_AlpineData.jl
Original file line number Diff line number Diff line change
Expand Up @@ -12,15 +12,15 @@
# To add the GMT package, simply add it with the julia package manager:
# ```julia
# julia> ]
# (@v1.8) pkg> add GMT
# (@v1.10) pkg> add GMT
# ```
# and load both GMG and GMT with:

using GeophysicalModelGenerator, GMT

# When loading both packages, several GMT routines within GMG will be loaded. One of these routines is the function ImportTopo, where one simply has to provide the region for which to download the topographic data and the data source.

Topo = ImportTopo([4,20,37,49], file="@earth_relief_01m.grd")
Topo = ImportTopo([4,20,37,50], file="@earth_relief_01m.grd")
# The data is available in different resolutions; see [here](http://gmt.soest.hawaii.edu/doc/latest/grdimage.html) for an overview. Generally, it is advisable to not use the largest resolution if you have a large area.

# If you have issues with loading the topography with GMT, there is also the alternative to download the data yourself and import it using Rasters.jl. This procedure is explained here. (ADD LINK TO TUTORIAL)
Expand Down Expand Up @@ -149,9 +149,6 @@ for iunit = 1:length(units)

end




# ## 3. Seismicity
# Earthquakes are always interesting, so we will now import the seismicity data from ISC.
# ### 3.1 Download and import
Expand Down Expand Up @@ -279,3 +276,6 @@ Data_GPS_Sanchez = GeoData(Lon,Lat,topo_v,(Velocity_mm_year=(Ve,Vn,Vz),V_north=V
# And as always, we'll save everything in VTK format and in jld2 format
Write_Paraview(Data_GPS_Sanchez, "GPS_Sanchez")
save_GMG("GPS_Sanchez",Data_GPS_Sanchez)


# Paffrath

0 comments on commit 3592b9a

Please sign in to comment.