[cig-commits] [commit] devel, master: Added a routine moho_stretching.f90 to slightly squish the mesh underneath the oceans and stretch it underneath the continents to get better sampling of the crust. This is accomplished based upon a degree-12 Moho map. Once paraview reads the new file structure we should make some mesh crss sections with the wave speeds superimposed to make sure we want to use this stretching. Tests suggests the stretched mesh gives slightly better dispersion. (7202119)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d_globe

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

>---------------------------------------------------------------

commit 7202119f17b475c4d0ce7968c6be79b555026003
Author: Jeroen Tromp <jtromp at princeton.edu>
Date:   Mon Jul 30 10:59:18 2007 +0000

    Added a routine moho_stretching.f90 to slightly squish the mesh underneath the oceans and
    stretch it underneath the continents to get better sampling of the crust. This is accomplished
    based upon a degree-12 Moho map. Once paraview reads the new file structure we should make
    some mesh crss sections with the wave speeds superimposed to make sure we want to use this
    stretching. Tests suggests the stretched mesh gives slightly better dispersion.


>---------------------------------------------------------------

7202119f17b475c4d0ce7968c6be79b555026003
 CMTSOLUTION |  4 ++--
 Par_file    | 12 ++++++------
 2 files changed, 8 insertions(+), 8 deletions(-)

diff --git a/CMTSOLUTION b/CMTSOLUTION
index 65ab6cb..cc874e6 100644
--- a/CMTSOLUTION
+++ b/CMTSOLUTION
@@ -1,10 +1,10 @@
 PDE 1999 11 26 13 21 15.60 -16.4200  168.2100  33.0 6.4 7.3 VANUATU ISLANDS  
 event name:     112699G
 time shift:      0.0000
-half duration:  15.0000
+half duration:  18.0000
 latitude:      -16.0800
 longitude:     168.3100
-depth:         515.0000
+depth:          15.0000
 Mrr:       1.310000e+27
 Mtt:       1.200000e+26
 Mpp:      -1.420000e+27
diff --git a/Par_file b/Par_file
index c18a965..b790490 100644
--- a/Par_file
+++ b/Par_file
@@ -15,12 +15,12 @@ GAMMA_ROTATION_AZIMUTH        = 20.d0
 
 # number of elements at the surface along the two sides of the first chunk
 # (must be 16 * multiple of NPROC below)
-NEX_XI                          = 160  
-NEX_ETA                         = 160 
+NEX_XI                          = 256  
+NEX_ETA                         = 256 
 
 # number of MPI processors along the two sides of the first chunk
-NPROC_XI                        = 5
-NPROC_ETA                       = 5
+NPROC_XI                        = 8
+NPROC_ETA                       = 8
 
 # 1D models with real structure:
 # 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref
@@ -31,7 +31,7 @@ NPROC_ETA                       = 5
 # fully 3D models:
 # transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
 # s20rts, s362ani, s362wmani, s362ani_prem, s29ea
-MODEL                           = s20rts
+MODEL                           = s362ani
 
 # parameters describing the Earth model
 OCEANS                          = .true.
@@ -45,7 +45,7 @@ ATTENUATION                     = .true.
 ABSORBING_CONDITIONS            = .false.
 
 # record length in minutes
-RECORD_LENGTH_IN_MINUTES        =  0.5d0
+RECORD_LENGTH_IN_MINUTES        = 1.0d0
 
 # save AVS or OpenDX movies
 MOVIE_SURFACE                   = .false.