[CIG-SEISMO] SPECFEM3D_Cartesian: accuracy in model with water layer

[Dimitri Komatitsch]

Hello,

In such a case, what you need to do is set the "OCEANS" flag to .false. 
and then mesh the ocean i.e. the fluid layer explicitly using CUBIT (or 
any other mesher). This will work fine (if you manage to create a good 
mesh).

I guess we should rename the "OCEANS" flag to "APPROXIMATE_OCEAN_LOAD" 
or something like that to avoid confusion.

Best wishes,
Dimitri.

On 02/26/2013 07:00 PM, Youyi Ruan wrote:
> Hello,
>
> I have been running a high-resolution test case using
> SPECFEM3D_Cartesian (V2.0.2) in a local 3-D model includes a water layer
> and a thin sediment layer. I noticed that in the Par_file of
> waterlayer_examples, the “OCEANS” is set to “.false.”,  what does the
> code really do when it set to “.true.” ? (employ the approximation in
> Komatitsch and Tromp (2002)?)
> My simulations are run in a local 3-D model, the element in water layer
> is about the size of 1km x 1km x 1km.  According to the manual the
> accuracy of the approximation is about 20 s, this is probably not enough
> to my case and I am just wondering how can I do to take into account the
> effects of water phases/effects at period around or shorter than 10 s?
>   Any suggestions?
>
> Thank you,
> Best regards,
>
> Youyi
>
>
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://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