[CIG-SEISMO] Floating-point exception - erroneous arithmetic operation

[Dimitri Komatitsch]

Dear Moritz,

It seems you simply forgot to also divide the time step by two to remain 
stable when you increased mesh resolution by two:

  *** Max stability for P wave velocity (must be below about 0.50 or so) 
= 0.972936358376379

Best regards,
Dimitri.

On 05/27/2015 03:36 PM, Moritz Wehler wrote:
> Hey there,
>
> I have some issues with increasing the number of elements along the
> z-axis with specfem2D.
> As soon as i double the number of elements from 111 (in total along the
> z-axis) to 222 (in total), I get the following error message:
>
> "Program received signal SIGFPE: Floating-point exception - erroneous
> arithmetic operation.
>
> Backtrace for this error:
> #0  0x7F62A8D0D7D7
> #1  0x7F62A8D0DDDE
> #2  0x7F62A8449D3F
> #3  0x433DBE in compute_forces_viscoelastic_
> #4  0x4DF021 in MAIN__ at specfem2D.F90:?
> Gleitkomma-Ausnahme (Speicherabzug geschrieben)"
>
> I should add that everything worked perfectly fine BEFORE I increased
> the element number. The corresponding par_file and interface_file are
> attached.
>
> Hope I didn't forget anything important and that someone can help me.
>
> Cheers, Mo
>
>
> _______________________________________________
> 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,
UPR 7051, Marseille, France    http://komatitsch.free.fr