Merge pull request #831 from danielpeter/devel

updates EMC model routines (1d depth average/interpolation/tapering); initial work to support full gravity
SPECFEM · Feb 25, 2024 · dc089fc · dc089fc
2 parents b6fa4b6 + 3850570
commit dc089fc
Show file tree

Hide file tree

Showing 401 changed files with 9,746 additions and 3,513 deletions.
diff --git a/CITATION.cff b/CITATION.cff
@@ -0,0 +1,201 @@
+cff-version: 1.2.0
+
+message: "If this software turns out to be helpful for your research, please cite both the article from preferred-citation and the software itself. More detailed citation suggestions can be found in the user manual."
+title: 'SPECFEM/specfem3d_globe: SPECFEM3D_GLOBE v8.1.0'
+version: v8.1.0
+date-released: 2023-12-20
+repository-code: https://github.com/SPECFEM/specfem3d_globe/tree/v8.1.0
+url: "https://github.com/SPECFEM/specfem3d_globe"
+type: software
+identifiers:
+  - type: doi
+  - value: 10.5281/zenodo.10411115
+license:
+  - gpl-3.0-or-later
+
+authors:
+- family-names: Komatitsch
+  given-names: Dimitri
+- family-names: Tromp
+  given-names: Jeroen
+- family-names: Peter
+  given-names: Daniel
+- family-names: Sawade
+  given-names: Lucas
+- family-names: Thurin
+  given-names: Julien
+- family-names: Afanasiev
+  given-names: Michael
+- family-names: Alexeev
+  given-names: Dmitry
+- family-names: Ampuero
+  given-names: Jean-Paul
+- family-names: Ando
+  given-names: Kazuto
+- family-names: Bachmann
+  given-names: Etienne
+- family-names: Bai
+  given-names: Kangchen
+- family-names: Bangerth
+  given-names: Wolfgang
+- family-names: Basini
+  given-names: Piero
+- family-names: Beller
+  given-names: Stephen
+- family-names: Blitz
+  given-names: Celine
+- family-names: Bottero
+  given-names: Alexis
+- family-names: Bozdag
+  given-names: Ebru
+- family-names: Casarotti
+  given-names: Emanuele
+- family-names: Charles
+  given-names: Joseph
+- family-names: Chen
+  given-names: Min
+- family-names: Ciardelli
+  given-names: Caio
+- family-names: Cristini
+  given-names: Paul
+- family-names: Cui
+  given-names: Congyue
+- family-names: Durochat
+  given-names: Clement
+- family-names: Espindola Carmona
+  given-names: Armando
+- family-names: Galvez Barron
+  given-names: Percy
+- family-names: Gassmoeller
+  given-names: Rene
+- family-names: Gharti
+  given-names: Hom Nath
+- family-names: Goeddeke
+  given-names: Dominik
+- family-names: Gogar
+  given-names: Sunny
+- family-names: Grinberg
+  given-names: Leopold
+- family-names: Hjoerleifsdottir
+  given-names: Vala
+- family-names: Huang
+  given-names: Quancheng
+- family-names: Kendall
+  given-names: Elodie
+- family-names: Kientz
+  given-names: Sue
+- family-names: Labarta
+  given-names: Jesus
+- family-names: Lanucara
+  given-names: Piero
+- family-names: Le Goff
+  given-names: Nicolas
+- family-names: Le Loher
+  given-names: Pieyre
+- family-names: Lefebvre
+  given-names: Matthieu
+- family-names: Lei
+  given-names: Wenjie
+- family-names: Liu
+  given-names: Qinya
+- family-names: Liu
+  given-names: Youshan
+- family-names: Luet
+  given-names: David
+- family-names: Luo
+  given-names: Yang
+- family-names: Maggi
+  given-names: Alessia
+- family-names: Magnoni
+  given-names: Federica
+- family-names: Martin
+  given-names: Roland
+- family-names: Matzen
+  given-names: Rene
+- family-names: McPherson
+  given-names: Amanda
+- family-names: McRitchie
+  given-names: Dennis
+- family-names: Meenakshi
+  given-names: Yellapragada
+- family-names: Mehaut
+  given-names: Jean-Francois
+- family-names: Meschede
+  given-names: Matthias
+- family-names: Messmer
+  given-names: Peter
+- family-names: Michea
+  given-names: David
+- family-names: Miyoshi
+  given-names: Takayuki
+- family-names: Monteiller
+  given-names: Vadim
+- family-names: Nadh Somala
+  given-names: Surendra
+- family-names: Nissen-Meyer
+  given-names: Tarje
+- family-names: Nowacki
+  given-names: Andy
+- family-names: Orsvuran
+  given-names: Ridvan
+- family-names: Parisi
+  given-names: Laura
+- family-names: Podhorszki
+  given-names: Norbert
+- family-names: Pouget
+  given-names: Kevin
+- family-names: Raghukanth
+  given-names: S.T.G.
+- family-names: Rietmann
+  given-names: Max
+- family-names: Ruggiero
+  given-names: Vittorio
+- family-names: Sales de Andrade
+  given-names: Elliott
+- family-names: Savage
+  given-names: Brian
+- family-names: Schirwon
+  given-names: Malte
+- family-names: Schuberth
+  given-names: Bernhard
+- family-names: Sieminski
+  given-names: Anne
+- family-names: Smith
+  given-names: James
+- family-names: Strand
+  given-names: Leif
+- family-names: Tao
+  given-names: Kai
+- family-names: Tape
+  given-names: Carl
+- family-names: Tsuboi
+  given-names: Seiji
+- family-names: Valero Cano
+  given-names: Eduardo
+- family-names: Videau
+  given-names: Brice
+- family-names: Vilotte
+  given-names: Jean-Pierre
+- family-names: Xie
+  given-names: Zhinan
+- family-names: Zhang
+  given-names: Chang-Hua
+- family-names: Zhu
+  given-names: Hejun
+
+
+preferred-citation:
+  type: article
+  authors:
+  - family-names: Komatitsch
+    given-names: Dimitri
+  - family-names: Tromp
+    given-names: Jeroen
+  title: "Spectral-element simulations of global seismic wave propagation - I. Validation"
+  journal: "Geophys. J. Int."
+  volume: 149
+  issue: 2
+  start: 390  # First page number
+  end: 412    # Last page number
+  year: 2002
+  doi: "10.1046/j.1365-246X.2002.01653.x"
diff --git a/DATA/Par_file b/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .true.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 2.5d0
 

diff --git a/DATA/topo_bathy/display_topography_of_SPECFEM3D_GLOBE_with_GMT.sh b/DATA/topo_bathy/display_topography_of_SPECFEM3D_GLOBE_with_GMT.sh
@@ -22,16 +22,16 @@ gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:m
 # convert the non-regular cubed-sphere grid output of SPECFEM to a regular grid for GMT
  surface observation_grid_long_lat_topo_for_GMT.txt -Gtopography_of_SPECFEM_mesh.grd -Rd -I4m -f0x,1y -V
 
-## DK DK offset of 4.5cm
-## DK DK -Rd is an alias for -R-180/180/-90/90
+## offset of 4.5cm
+## -Rd is an alias for -R-180/180/-90/90
 grdimage topography_of_SPECFEM_mesh.grd -CGMTglobe.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > $ps
 
-## DK DK offset of -1.5 inch
+# offset of -1.5 inch
 pscoast -Rd -JK$central_meridian/9i -B45g30:."Topography of SPECFEM3D_GLOBE mesh": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE calculations by Dimitri Komatitsch" -V -K -O >> $ps
 
 psscale -CGMTglobe.cpt -D12.5/-1.5/16/0.25h -B2000.:"Elevation (m)": -V -O >> $ps
 
-## DK DK convert the final file to PDF
+# convert the final file to PDF
 ps2pdf $ps
 
 # Clean up

diff --git a/EXAMPLES/benchmarks/attenuation_benchmark_GJI_2002_versus_normal_modes/DATA/Par_file b/EXAMPLES/benchmarks/attenuation_benchmark_GJI_2002_versus_normal_modes/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .false.
 ROTATION                        = .false.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 55.0d0
 

diff --git a/EXAMPLES/benchmarks/attenuation_benchmark_GJI_2002_versus_normal_modes/convolve_jeroen.f90 b/EXAMPLES/benchmarks/attenuation_benchmark_GJI_2002_versus_normal_modes/convolve_jeroen.f90
@@ -1,14 +1,14 @@
 
   program convolve_jeroen
 
-!! JT JT January 2001
+!! January 2001, by Jeroen Tromp
 
-!! DK DK added option to convolve with a triangle
-!! DK DK also now reads hdur and number of lines from external file
+!! Dimitri Komatitsch: added option to convolve with a triangle
+!!                     also now reads hdur and number of lines from external file
 
-!! DK DK suppressed last "hdur" seconds in seismograms
-!! DK DK to remove artefact from the convolution at the end of the seismogram
-!! DK DK which was then creating some problems in the "remove_offset" routine
+!! suppressed last "hdur" seconds in seismograms
+!! to remove artefact from the convolution at the end of the seismogram
+!! which was then creating some problems in the "remove_offset" routine
 
   implicit none
 

diff --git a/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file b/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .false.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 50.0d0
 

diff --git a/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file b/EXAMPLES/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .false.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 100.0d0
 

diff --git a/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file b/EXAMPLES/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .false.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 100.0d0
 

diff --git a/EXAMPLES/benchmarks/undo_attenuation_benchmark_GJI_2016/DATA/Par_file b/EXAMPLES/benchmarks/undo_attenuation_benchmark_GJI_2016/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .true.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 90.0d0
 

diff --git a/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file b/EXAMPLES/global_PREM_kernels/amplitude/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .false.
 ROTATION                        = .false.
 ATTENUATION                     = .false.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 15.0d0
 

diff --git a/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file b/EXAMPLES/global_PREM_kernels/traveltime/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .false.
 ROTATION                        = .false.
 ATTENUATION                     = .false.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 15.0d0
 

diff --git a/EXAMPLES/global_s362ani_shakemovie/DATA/Par_file b/EXAMPLES/global_s362ani_shakemovie/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .true.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 100.0d0
 

diff --git a/EXAMPLES/global_small/DATA/Par_file b/EXAMPLES/global_small/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .true.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 0.5d0
 

diff --git a/EXAMPLES/mars_global/DATA/Par_file b/EXAMPLES/mars_global/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .true.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 15.0d0
 

diff --git a/EXAMPLES/mars_regional/DATA/Par_file b/EXAMPLES/mars_regional/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .true.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 2.5d0
 

diff --git a/EXAMPLES/moon_global/DATA/Par_file b/EXAMPLES/moon_global/DATA/Par_file
@@ -73,6 +73,10 @@ GRAVITY                         = .true.
 ROTATION                        = .true.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 1.5d0
 

diff --git a/EXAMPLES/noise_examples/global_long/DATA/Par_file_NOISE_1_attenuation b/EXAMPLES/noise_examples/global_long/DATA/Par_file_NOISE_1_attenuation
@@ -73,6 +73,10 @@ GRAVITY                         = .false.
 ROTATION                        = .false.
 ATTENUATION                     = .true.
 
+# full gravity calculation by solving Poisson's equation for gravity potential instead of using a Cowling approximation
+# (must have also GRAVITY flag set to .true. to become active)
+FULL_GRAVITY                    = .false.
+
 # record length in minutes
 RECORD_LENGTH_IN_MINUTES        = 205.0d0