[CIG-SEISMO] specfem3d_globe-7.0.0 overflow

[Dimitri Komatitsch]

Dear Tsuboi-san,

It is probably the time step computed automatically by the code that is 
slightly above the CFL stability limit in that case. Thus, you probably 
just need to reduce the value of variable DT a bit (say 10% or so).

Best regards,
Dimitri.

On 07/02/2018 09:41 AM, Seiji Tsuboi wrote:
> Dear Sir,
> 
> I have used specfem3d_globe-7.0.0 on intel cluster and have overflow 
> problems.
> If I set NPROC=16 and NEX=768 (1536 cores), I could calculate 
> seismograms correctly. However if I increase NEX, for example NPROC=17, 
> NEX=816 (1734 cores), xspecfem3D stops with overflow in the very 
> beginning of computation. ATTENUATION is set as true and it seems that 
> the program overflows in "multiply_accel_elastic" routine. I tried both 
> ifort and gfortran but the situation is the same. I would be grateful if 
> someone would let me know how I should solve this issue.
> 
> Best regards,
> Seiji Tsuboi
> 
> 
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
> 

-- 
Dimitri Komatitsch, CNRS Research Director (DR CNRS)
Laboratory of Mechanics and Acoustics, Marseille, France
http://komatitsch.free.fr