Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

use seissol v1.3/Rusanov fluxes/ vtkhdf/ extend tpv13 #52

Merged
merged 24 commits into from
Nov 11, 2024
Merged
Changes from 1 commit
Commits
Show all changes
24 commits
Select commit Hold shift + click to select a range
cc84e30
use rusanov
Thomas-Ulrich Nov 4, 2024
e1cddd3
fix Dockerfile_jupyterlab
Thomas-Ulrich Nov 4, 2024
4051c03
small improvements on tpv13_training.geo
Thomas-Ulrich Nov 7, 2024
010b1f0
update to msh4
Thomas-Ulrich Nov 7, 2024
4e9b92e
switch to main
Thomas-Ulrich Nov 7, 2024
2e007f3
update both dockers
Thomas-Ulrich Nov 7, 2024
5abb3c7
use latest releases of SeisSol and PUMGen
Thomas-Ulrich Nov 7, 2024
9e2929c
use vtk output for fault
Thomas-Ulrich Nov 7, 2024
65cfc19
add rupture speed exercise
Thomas-Ulrich Nov 7, 2024
fbbf8b9
combine subsidence plots and fix typo in geo
Thomas-Ulrich Nov 9, 2024
637d10e
close plotter to avoid memory leak
Thomas-Ulrich Nov 9, 2024
05c2e87
precise exercise and input file
Thomas-Ulrich Nov 9, 2024
e727737
disable vtk fault output for now
Thomas-Ulrich Nov 9, 2024
8d7c653
add small exercice related with the initial stress
Thomas-Ulrich Nov 10, 2024
b366ca6
add sympy
Thomas-Ulrich Nov 10, 2024
e73d4ed
fix dictionnary
Thomas-Ulrich Nov 10, 2024
933bcaf
more precise instructions
Thomas-Ulrich Nov 10, 2024
358069b
fix missing import module and format files naemes and variables with
Thomas-Ulrich Nov 10, 2024
a5ef412
add subtitle for ground deformation exercise
Thomas-Ulrich Nov 10, 2024
03de600
improve exercise with varying dip
Thomas-Ulrich Nov 10, 2024
d6cad50
add pandas (for moment rate release)
Thomas-Ulrich Nov 10, 2024
cfc93eb
fix bug intoduced previously
Thomas-Ulrich Nov 11, 2024
7345b35
change task title
Thomas-Ulrich Nov 11, 2024
459abea
fix clim for Vr
Thomas-Ulrich Nov 11, 2024
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
combine subsidence plots and fix typo in geo
Thomas-Ulrich committed Nov 9, 2024
commit fbbf8b93040668c050a0d40dd603d59aaf52ce23
26 changes: 16 additions & 10 deletions tpv13/tpv13.ipynb
Original file line number Diff line number Diff line change
@@ -266,16 +266,22 @@
" p.camera.zoom(3)\n",
"p.show(jupyter_backend='panel')\n",
" \n",
"# plot subsidence/uplift profile\n",
"for i, pref in enumerate(prefixList):\n",
" load_mesh(pref).plot_over_line(\n",
" a,\n",
" b,\n",
" resolution=10000,\n",
" title=f\"subsidence/uplift profile for {pref}\",\n",
" ylabel=\"vertical displacement (m)\",\n",
" scalars='u3'\n",
" )"
"def plot_combined_profiles(prefixList, a, b, resolution=10000):\n",
" plt.figure(figsize=(6, 4))\n",
" for prefix in prefixList:\n",
" mesh = load_mesh(prefix)\n",
" sampled_data = mesh.sample_over_line(a, b, resolution=resolution)\n",
" distance = sampled_data['Distance']\n",
" displacement = sampled_data['u3']\n",
" plt.plot(distance/1000, displacement, label=prefix, linewidth=2)\n",
" plt.title(\"Comparison of Subsidence/Uplift Profiles\")\n",
" plt.xlabel(\"Distance (km)\")\n",
" plt.ylabel(\"Vertical Displacement (m)\")\n",
" plt.legend()\n",
" plt.tight_layout()\n",
" plt.show()\n",
"\n",
"plot_combined_profiles(prefixList, a, b)"
]
},
{
2 changes: 1 addition & 1 deletion tpv13/tpv13_training.geo
Original file line number Diff line number Diff line change
@@ -32,7 +32,7 @@ Z0 = -42e3;
// Create the domain as a box
domain = newv; Box(domain) = {X0, Y0, Z0, X1-X0, Y1-Y0, -Z0};

// Create the fault as a in the x-y plane, centered in x at the hypocenter
// Create the fault as a rectangle in the x-y plane, centered in x at the hypocenter
fault = news; Rectangle(fault) = {-l_f/2, -w_f, 0, l_f, w_f};

// Create the nucleation patch as a smaller rectangle