[CIG-SEISMO] SPECFEM2D, CUBIT and CPML

[Dimitri Komatitsch]

Hi Andrea,

You can adapt this routine I wrote in the 3D version of the code:

utils/CPML/convert_external_layers_of_a_given_mesh_to_CPML_layers.f90

It assigns CPML flags to the outer layers of an existing CUBIT/TRELIS 
mesh. (note that it is better to have a regular i.e. non distorted mesh 
in these outer layers)

If you adapt it to 2D then please email it to me and I will add it to 
the official Git version of the 2D code. (or commit it to Git using 
https://github.com/geodynamics/specfem3d/wiki/Using-Git-for-SPECFEM  )

For your tests I suggest you use the "devel" version of the code rather 
than the stable one:

git clone --recursive --branch devel 
https://github.com/geodynamics/specfem2d.git

(or git checkout --branch devel  in your existing copy)

To get the 3D program I mentioned:

git clone --recursive --branch devel 
https://github.com/geodynamics/specfem3d.git

Best,
Dimitri.

On 02/09/2015 07:45 PM, Andrea Colombi wrote:
> I would like to use CPML on the 4 boundaries for a p-sv simulation,
> however the implementation remains unclear when one is using CUBIT as
> the mesh generator.
> The only example using CPML seems to be "salt_dome_CUBIT_mesh" but I
> can't figure out how the CPML layers are practically implemented. I have
> tried to implement the PML layer following the procedure in the manual
> (attached python script). Specfem recognises the presence of PML (at
> least the output says so) but the simulations show strong border
> reflections (as surfaces were free). Does anyone have a working example
> of CUBIT and CPML that I could use or have a look?
> I am using the current stable release. Is there any known issue with
> this version and CPML?
> Thanks for the help,
>
> Andrea
>
>
>
> _______________________________________________
> 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