<!doctype html public "-//W3C//DTD W3 HTML//EN">
<html><head><style type="text/css"><!--
blockquote, dl, ul, ol, li { padding-top: 0 ; padding-bottom: 0 }
--></style><title>SPECFEM3D_GLOBE_V5.0.0 problem</title></head><body>
<div>Hello Jeroen -</div>
<div><br></div>
<div>Greetings from Grenoble. I'm settled in an started trying to run
SPECFEM3D_GLOBE_V5.0.0.</div>
<div><br></div>
<div>I have had several cases crash with the solution blowing up in
the fluid. I started with a 3D case (s29ea) to compare to previous
results (one-chunk, 45-degree width, north korean explosion). I also
tried with a 1D model (1D_iasp91) and it blew up as well.</div>
<div><br></div>
<div>I have incrementally decreased the time step by as much as 25%,
by changing the DT in read_compute_parameters.f90 as below:</div>
<div><br></div>
<div>! takes a 5% safety margin on the maximum stable time step</div>
<div>! which was obtained by trial and error</div>
<div>! DT = DT * (1.d0 - 0.05d0)</div>
<div> DT = DT * (1.d0 - 0.25d0)</div>
<div><font face="Arial"><br></font></div>
<div><font face="Arial">I'm running on atlas, a linux cluster running
CHAOS 4.2 on AMD Opteron cpu's and the code was compiled with
both mpiifort. I had problems when using mpif90 also.</font></div>
<div><font face="Arial"><br></font></div>
<div><font face="Arial">I have not had to make any changes in DT with
the version 4.0.4 of the code.</font></div>
<div><font face="Arial"><br></font></div>
<div>Could there be problems with running a chunk < 90
degrees?</div>
<div><br></div>
<div>Regards,</div>
<div>Artie</div>
</body>
</html>