[CIG-SEISMO] Can SPECEM3D do gravity inversion?

[Dimitri Komatitsch]

Dear Ming,

Perfect! Please let us know if you have any feedback (and if you improve 
it please commit your changes to Git or email them to me and I will add 
them to the official version of the code and will list your name in the 
list of authors of the code).

Thanks,
Cheers,

Dimitri.

On 03/18/2017 06:00 PM, mzhao at cea-igp.ac.cn wrote:
> Dear all,
>
> Coooooool!!I can help to have a try.Many thanks!
>
> Cheers
> Ming
> ------------------------------------------------------------------------
> mzhao at cea-igp.ac.cn
>
> *From:* Dimitri Komatitsch <mailto:komatitsch at lma.cnrs-mrs.fr>
> *Date:* 2017-03-19 00:46
> *To:* mzhao at cea-igp.ac.cn <mailto:mzhao at cea-igp.ac.cn>;
> cig-seismo at geodynamics.org <mailto:cig-seismo at geodynamics.org>; Roland
> Martin <mailto:roland.martin at get.omp.eu>; Sébastien Chevrot
> <mailto:sebastien.chevrot at get.omp.eu>; Jeroen Tromp
> <mailto:jtromp at princeton.edu>
> *Subject:* Re: [CIG-SEISMO] Can SPECEM3D do gravity inversion?
>
> Hi all,
>
> Nice!
>
> For gravity modeling, in the current version of SPECFEM3D one needs to
> change the following flag to .true. in setup/constants.h.in :
>
> !!-----------------------------------------------------------
> !!
> !! gravity integral calculations
> !!
> !!-----------------------------------------------------------
>
> !! DK DK for gravity integrals
>    logical, parameter :: GRAVITY_INTEGRALS = .false. !!! .true.
>
>
> Best regards,
> Dimitri.
>
>
> On 03/18/2017 05:32 PM, Jeroen Tromp wrote:
>> Dear all,
>>
>> This is something Hom Nath Gharti has worked on as well, combining a
>> spectral-element method with an infinite-element method to solve
>> Laplace's equation in outer space. See his recent CIG webinar:
>> https://www.youtube.com/watch?v=uIuvv7tAcfk
>>
>> Best regards,
>>
>> Jeroen Tromp
>>
>>
>> On 3/18/17 12:29 PM, Dimitri Komatitsch wrote:
>>>
>>> Hi Ming,
>>>
>>> Thanks for your message! This is indeed an interesting topic. My
>>> friend Roland Martin has a nice paper that has just appeared on this:
>>> http://komatitsch.free.fr/preprints/GJI_Martin_gravimetry_2017.pdf
>>>
>>> thus let me cc Roland to see what the current version of SPECFEM3D can
>>> do (it can do 3D modeling for sure, but I think the gravity inversion
>>> part was performed outside of SPECFEM using scripts for now (?);
>>> Roland can tell us more about this, he is the expert on that).
>>>
>>> Thanks!
>>> Cheers,
>>> Dimitri.
>>>
>>> On 03/18/2017 04:08 PM, mzhao at cea-igp.ac.cn wrote:
>>>> Hi Dimitri,
>>>>
>>>> I read a very interesting paper recently which use spectral element
>>>> method to do gravity inversion,I wonder if the SPECFEM3D can do gravity
>>>> inversion in the future? We have a lot of seismic data and gravity
>>>> data in
>>>> our group,and I'm using spectral element method all the time,if it can
>>>> do both adjoint inversion and gravity inversion, that would be exciting!
>>>>
>>>> Cheers
>>>> Ming Zhao
>>>> ------------------------------------------------------------------------
>>>> mzhao at cea-igp.ac.cn
>>>
>>
>> _______________________________________________
>> 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, Marseille, France
> http://komatitsch.free.fr
>

-- 
Dimitri Komatitsch, CNRS Research Director (DR CNRS)
Laboratory of Mechanics and Acoustics, Marseille, France
http://komatitsch.free.fr