[CIG-SEISMO] Specfem with an altered crustal model

[Kendall, Elodie]

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 tried changing 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170706/76e28caa/attachment.html>