[CIG-SEISMO] Could I use Force Point Source in specfem3d_globe?

[Dimitri Komatitsch]

Hi Xiaohan Chen, Hi all,

I do not think that is fixed/closed yet, since the Git issue is still 
open, however it is probably easy to test and to fix if needed because 
adding a force source is very easy in the SEM.

I think the current code is probably OK, it is just untested and thus we 
would need somebody to test it carefully once and then close the Git issue.

Thanks,
Dimitri.

On 06/19/2015 12:30 PM, cxh757 wrote:
> Hi all
>
> There is USE_FORCE_POINT_SOURCE flag in Par_file of specfem3d, and we
> can turn this flag to use force point source instead of a moment-tensor
> source.
>
> But I did not find this flag in Par_file of specfem3d_globe. It seems
> there are some related parameters in setup/constants.h.in
> of specfem3d_globe:
>
> ! use a force source located exactly at a grid point instead of a
> CMTSOLUTION source
> ! this can be useful e.g. for asteroid impact simulations
> ! in which the source is a vertical force, normal force, impact etc.
>    logical, parameter :: USE_FORCE_POINT_SOURCE = .false.
>    double precision, parameter :: FACTOR_FORCE_SOURCE = 1.d15
>    integer, parameter :: COMPONENT_FORCE_SOURCE = 3  ! takes direction
> in comp E/N/Z = 1/2/3
>
> I am not sure whether I can use these parameters. According to this
> github issue:https://github.com/geodynamics/specfem3d_globe/issues/270,
> using these parameters is unsafe in specfem3d_globe.
>
> Does anyone could tell me how to use force point source safely
> in specfem3d_globe?
>
> And I really suggest USE_FORCE_POINT_SOURCE flag should be added in
> Par_file of specfem3d_globe, similar to specfem3d
> (https://github.com/geodynamics/specfem3d/issues/54).
>
> Thanks
>
> Yours
> Xiaohan Chen
>
>
> _______________________________________________
> 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