[CIG-SEISMO] Is SH plane-wave implemented in specfem2d?

[Dimitri Komatitsch]

Hi Seokho,

I do not think that is implemented yet, but it is easy to add.
Just type this to see where current P-SV plane waves are implemented:

grep -i plane src/*/*90

and when you see lines implementing them (at least in files 
src/specfem2D/prepare_initial_field.F90 and 
src/specfem2D/compute_Bielak_conditions.f90, maybe more) change them to:

if (p_sv) then ! P-SV waves
   leave the current statements here
else ! SH waves
   add your SH plane wave here; apply it to component number 2 of the 
displacement / velocity / acceleration arrays, since if I remember 
correctly in the 2D code components 1 and 3 are for P-SV and component 2 
is for SH (please doublecheck)
endif

If you implement that please submit it to the official source code 
(following 
https://github.com/geodynamics/specfem3d/wiki/Using-Git-for-SPECFEM ), 
it will be a useful contribution (and add your name to the AUTHORS file 
and to authors.tex in the users manual).

Thanks,
Dimitri.

On 08/04/2015 05:55 AM, Seokho Jeong wrote:
> Hi all,
>
> I was trying to run a simulation with SH plane wave as source. I was
> under the impression that this is possible, but when I ran the
> simulation it only produced zero displacement everywhere.
>
> SV plane wave and SH point source appear to be working fine.
>
> Is this not doable yet? or did I make a mistake in my input files?
>
> Best regards,
> Seokho
>
>
>
> _______________________________________________
> 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