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

adds coupling with wavefield injection from a SPECFEM simulation; adds Vs30 routines; updates setup scripts #1774

Closed
wants to merge 25 commits into from
Closed
Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
25 commits
Select commit Hold shift + click to select a range
4ea2f3a
updates blender scripts
danielpeter Dec 7, 2024
351089e
updates IRIS EMC script (recognizes RHO/VP/VS/.. as parameter names)
danielpeter Dec 7, 2024
6b4ef8c
updates tomography model info output
danielpeter Dec 7, 2024
21efea4
updates scripts; adds script for getting USGS Vs30 data
danielpeter Dec 11, 2024
6811924
updates C-declaration for newer compilers in modified AxiSEM
danielpeter Dec 12, 2024
6e65cfe
adds a "SPECFEM" injection technique (box tomography; simulation coup…
danielpeter Dec 15, 2024
9363912
adds map plot to VS30 script
danielpeter Dec 16, 2024
b9e1ee6
adds an option to output tabulated Modified Mercalli Intensity (MMI) …
danielpeter Dec 16, 2024
6d4e17c
adds routines to impose (USGS) Vs30 model in top 30-m layer
danielpeter Dec 16, 2024
e9ecfac
code cleaning
danielpeter Dec 16, 2024
17aa291
updates source infos (for multiple sources)
danielpeter Dec 17, 2024
8fd4d1a
updates setup scripts to account for UTM zones in Southern hemisphere
danielpeter Dec 17, 2024
d92cfe8
code cleaning (removing white spaces)
danielpeter Dec 18, 2024
8c196c9
adds coupling example for specfem-specfem coupling
danielpeter Dec 18, 2024
1070f98
adds some more comments in the README of the coupling example
danielpeter Dec 18, 2024
352ab97
adds coupling example test
danielpeter Dec 18, 2024
8ccb2db
updates utm2geo python script routine to account for negative zone nu…
danielpeter Dec 18, 2024
78b3bf2
updates test to use less memory for example
danielpeter Dec 18, 2024
c088aba
updates coupling test
danielpeter Dec 19, 2024
fb55620
uses direct access for reading coupling wavefield files
danielpeter Dec 19, 2024
0ebaa4f
updates coupling routine
danielpeter Dec 19, 2024
32c154a
shortens module name (to avoid ifort warnings)
danielpeter Dec 19, 2024
2a5e4d1
updates coupling_points array format
danielpeter Dec 19, 2024
80ed4cc
updates array initialization for scotch partitioning
danielpeter Dec 19, 2024
917e4e0
updates source locations w/ coupling
danielpeter Dec 19, 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
code cleaning (removing white spaces)
  • Loading branch information
danielpeter committed Dec 18, 2024
commit d92cfe8103ad238b60b59d2ea231d196ae5f5aa3
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@
Surface shape functions defined by NGNOD2D = 4 control nodes
Beware! Curvature (i.e. HEX27 elements) is not handled by our internal mesher

velocity model: default
velocity model: default


suppressing UTM projection
Expand All @@ -46,7 +46,7 @@
defined materials : 1
undefined materials : 0
total number of spectral elements: 20736
absorbing boundaries:
absorbing boundaries:
xmin,xmax : 576 576
ymin,ymax : 576 576
bottom,top: 1296 1296
Expand All @@ -60,14 +60,14 @@

create regions:

...allocating arrays
...allocating arrays
NGLLX = 5
NGLLY = 5
NGLLZ = 5
NGNOD = 8
NGNOD2D = 4

main process setup:
main process setup:
nspec = 5136

separating regular/irregular element shapes
Expand All @@ -84,16 +84,16 @@
File DATA/Par_file_faults not found: assuming that there are no faults


...setting up jacobian
...setting up jacobian

...indexing global points
creating ibool indexing : x min/max = 0.00000000 / 134000.000
creating ibool indexing : x min/max = 0.00000000 / 134000.000
creating indirect addressing: nglob = 343525
creating unique point locations

...preparing MPI interfaces
...preparing MPI interfaces
number of interfaces : 2
creating MPI indexing : x min/max = 0.00000000 / 134000.000
creating MPI indexing : x min/max = 0.00000000 / 134000.000
tolerance = 1.3400000000000000E-005

total MPI interface points: 42770
Expand Down Expand Up @@ -124,7 +124,7 @@
90 % time remaining: 2.2891237377624420E-008 s
100 % time remaining: 2.6291509789120710E-010 s

...detecting acoustic-elastic-poroelastic surfaces
...detecting acoustic-elastic-poroelastic surfaces
total acoustic elements : 20736
total elastic elements : 0
total poroelastic elements: 0
Expand All @@ -134,21 +134,21 @@
elastic - poroelastic coupling : total number of faces = 0


...element inner/outer separation
...element inner/outer separation
for overlapping of communications with calculations:
percentage of edge elements 11.2149506 %
percentage of volume elements 88.7850494 %


...element mesh coloring
...element mesh coloring
use coloring = F

...external binary models
no external binary model used
...external binary models
no external binary model used

...creating mass matrix
...creating mass matrix

...setting up mesh adjacency
...setting up mesh adjacency

mesh adjacency:
total number of elements in this slice = 5136
Expand All @@ -174,7 +174,7 @@
(maximum neighbor of neighbors) = 98
total number of neighbors = 512346

Elapsed time for detection of neighbors in seconds = 0.133644789
Elapsed time for detection of neighbors in seconds = 0.133644789


...saving mesh databases
Expand Down Expand Up @@ -209,38 +209,38 @@


********
Model: P velocity min,max = 2800.00000 2800.00000
Model: S velocity min,max = 0.00000000 0.00000000
Model: P velocity min,max = 2800.00000 2800.00000
Model: S velocity min,max = 0.00000000 0.00000000

Model: Poisson's ratio min,max = 0.500000000 0.500000000
Model: Poisson's ratio min,max = 0.500000000 0.500000000
********

*********************************************
*** Verification of simulation parameters ***
*********************************************

*** Xmin and Xmax of the model = 0.00000000 134000.000
*** Ymin and Ymax of the model = 0.00000000 134000.000
*** Zmin and Zmax of the model = -60000.0000 0.00000000
*** Xmin and Xmax of the model = 0.00000000 134000.000
*** Ymin and Ymax of the model = 0.00000000 134000.000
*** Zmin and Zmax of the model = -60000.0000 0.00000000

*** Max GLL point distance = 1227.47656
*** Min GLL point distance = 642.726562
*** Max/min ratio = 1.90979588
*** Max GLL point distance = 1227.47656
*** Min GLL point distance = 642.726562
*** Max/min ratio = 1.90979588

*** Max element size = 3750.00000
*** Min element size = 3722.21875
*** Max/min ratio = 1.00746357
*** Max element size = 3750.00000
*** Min element size = 3722.21875
*** Max/min ratio = 1.00746357

*** Minimum period resolved = 1.67410719
*** Maximum suggested time step = 0.109999999
*** Minimum period resolved = 1.67410719
*** Maximum suggested time step = 0.109999999

Elapsed time for checking mesh resolution in seconds = 1.2908709000001295E-002
saving VTK files for Courant number and minimum period


mesh regions done

min and max of elevation (i.e. height of the upper surface of the mesh) included in mesh in m is 0.0000000000000000 0.0000000000000000
min and max of elevation (i.e. height of the upper surface of the mesh) included in mesh in m is 0.0000000000000000 0.0000000000000000


done mesh setup
Expand Down Expand Up @@ -277,7 +277,7 @@
smallest and largest possible floating-point numbers are: 1.17549435E-38 3.40282347E+38


Elapsed time for mesh generation and buffer creation in seconds = 15.3367910
Elapsed time for mesh generation and buffer creation in seconds = 15.3367910
Elapsed time for mesh generation and buffer creation in hh:mm:ss = 0 h 00 m 15 s

End of mesh generation
Expand Down
Loading