[CIG-SEISMO] Can we change the NGLLX?

[Dimitri Komatitsch]

Dear Jungrak,

Yes, you can change it. However, papers (e.g. by Seriani and Oliveira, 
De Basabe and Sen) show that the optimal value is around NGLLX = 5.
2 is way too low (poor accuracy), and 10 will lead to a small time step.

Best regards,
Dimitri.

On 12/02/2014 05:05 AM, Jungrak Son wrote:
> To whom it may concern,
>
> Thanks very much for your time to read this email. Let me briefly ask a
> question about SpecFEM2D.
>
> I know that SpecFEM2D uses 4th order polynomial degree. (NGLLX = 5 in
> constant.h file). *However, can we change the polynomial degree?* Since
> this amazing tool can use irregular mesh, I want to compare the result
> accuracy among different NGLLX settings from 2 to 10.
>
> Thanks for your consideration.
>
> Sincerely,
> Jungrak Son
>
> MS student in Geophysics from Korea
>
>
> _______________________________________________
> 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