[CIG-SEISMO] Specfem with an altered crustal model

[Ebru Bozdag]

Hi Elodie,

Are you honouring the crustal thickness in your simulations, i.e., are you using

! to suppress element stretching for 3D moho surface
  logical,parameter :: SUPPRESS_MOHO_STRETCHING = .false.

in your constant.h file? If this is the case perhaps after you add perturbations you may also be distorting the mesh aspect ratio at some locations.

Just try to run the mesher and the solver by setting SUPPRESS_MOHO_STRETCHING = .true. and see if you still have the same problem, if you have not tried it yet.

Best regards,
Ebru


-----------------

Ebru Bozdag
Assistant Professor

COLORADOSCHOOLOFMINES
Department of Geophysics
bozdag at mines.edu<http://www.minesnewsroom.com/email-signature-templates#> | tel: +1-303-273-3578  | fax: +1-303-273-3478







On 06 Jul 2017, at 02:32, Kendall, Elodie <elodie.kendall.14 at ucl.ac.uk<mailto:elodie.kendall.14 at ucl.ac.uk>> wrote:

Hi there,

I was hoping to ask you a question about Specfem3d_globe please? I have been using it a few months now and I am currently trying to introduce crustal thickness perturbations into crust2.0 and analyse the effect on the synthetics. I have set up a new crust2.0 which has 16,200 keys (one for each 2 by 2 degree grid) but the same layers for each grid. I have altered the corresponding subroutines for this new set of keys, no errors with the mesher and solver- same synthetics as before. Now, I add the perturbations to each layer of each grid in my crust2.0 model. I have changed the maximum moho depth from 90 to 115 (the new max.), the mesher runs fine but the solver fails with error "forward simulation became unstable in fluid and blew up". I have triedchanging the time step DT to DT*0.95 and DT*0.8 in the script shared/
get_timestep_and_layers.f90 however the same error appears. I was hoping you could help me please?

Thanks a lot,
Elodie Kendall

_______________________________________________
CIG-SEISMO mailing list
CIG-SEISMO at geodynamics.org<mailto:CIG-SEISMO at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170707/d47a0882/attachment-0001.html>