[CIG-SEISMO] Homogeneous earth model

[Stefan Wenk]

Dear CIG-Seismo Team,


this is Stefan Wenk from Munich University. Recently we've tried to compare SpecFEM3D with 
our DG Code on global scales. Therefore, the comparison we would like to start 
as simple as possible, using a homogeneous earth model.


For that reason I've tried to change the isotropic PREM routine model_prem.f90 
in the source code and implemented constant material values for rho, vp, and vs 
instead of the polynomials. The compilation worked fine but the code crashed due 
to instability problems. I think this is caused by the constant S-wave velocity 
unequal to zero I've chosen in the outer core.


Does the calculation of the potential you do in the outer core assume vs=0? Is 
there any possibility to circumvent this issue without huge changes in the code?


Best regards,

Stefan

-- 
Department of Earth and Environmental Sciences
Geophysics
Munich University
Theresienstr. 41
80333 Munich / Germany

Room:  441
Phone: +49 (0) 89 2180-4143
Fax:   +49 (0) 89 2180-4205
Mail:  wenk at geophysik.uni-muenchen.de