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

[cxh757]

Hi Dimitri

Thanks for your reply.


So If I want to use a force source in specfem3d_globe, I just (1) set USE_FORCE_POINT_SOURCE as .true. in setup/constants.h.in, (2) then copy FORCESOLUTION  in SPECFEM3D (_Cartesian) and paste it in SPECFEM3D_GLOBE.
Is it right?


But I notice that there is "factor force source" flag both in FORCESOLUTION file and setup/constants.h.in file. Which one will be used during calculation?


Yours
Xiaohan Chen

At 2015-06-19 20:09:13, "Dimitri Komatitsch" <komatitsch at lma.cnrs-mrs.fr> wrote:
>
>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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20150619/e9a7024c/attachment.html>