From xiangtian2589 at yahoo.com  Fri Jun  2 07:05:10 2017
From: xiangtian2589 at yahoo.com (xiang tian)
Date: Fri, 2 Jun 2017 14:05:10 +0000 (UTC)
Subject: [CIG-SEISMO] problem with compiling xspecfem3D
References: <319688862.392312.1496412310296.ref@mail.yahoo.com>
Message-ID: <319688862.392312.1496412310296@mail.yahoo.com>

 I also meet a problem while I compile the SPECFEM3D Cartesian.
I run the configure shell script as follow:./configure  FC=gfortran CC=gcc MPIFC=mpif90 --with-mpi
I can make xcheck_mesh_quality  xdecompose_mesh  xgenerate_databases successfully, but I meet a problem while I make xspecfem3D.
The error is as follow:
gfortran -g -O2  -std=f2003 -fimplicit-none -frange-check -fmax-errors=10 -pedantic -pedantic-errors -Waliasing -Wampersand -Wcharacter-truncation -Wline-truncation -Wsurprising -Wno-tabs -Wunderflow -ffpe-trap=invalid,zero,overflow -Wunused -Werror -O3  -J./obj -I./obj -I.  -I./setup  -c -o obj/compute_forces_viscoelastic.spec.o src/specfem3D/compute_forces_viscoelastic.f90
src/specfem3D/compute_forces_viscoelastic.f90:35.57:                        ATTENUATION,deltat,PML_CONDITIONS, &
                                                         1
Warning: Unused dummy argument 'pml_conditions' at (1)
make: *** [obj/compute_forces_viscoelastic.spec.o] Error 1
Does anyone also meet this problem? Does anyone would like to tell us how to deal with this problem?
xiang tian
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170602/0e283bda/attachment.html>

From hgharti at princeton.edu  Fri Jun  2 08:15:20 2017
From: hgharti at princeton.edu (Hom Nath Gharti)
Date: Fri, 2 Jun 2017 11:15:20 -0400
Subject: [CIG-SEISMO] Some problems
In-Reply-To: <f40947bd0b73434283be01f10f8e6e3e@CSGHUB208W.pu.win.princeton.edu>
References: <f40947bd0b73434283be01f10f8e6e3e@CSGHUB208W.pu.win.princeton.edu>
Message-ID: <CAC+G=sNbDE9t5RUQffY3=-SKCtFHWu7XHhsXJ-ZTTMPiZN1dbA@mail.gmail.com>

Dear Tu Xiang,

It seems that there is an unused variable in the
compute_forces_viscoelastic.f90
of the master branch. That is fixed in the development branch.

You can either try development branch

OR

temporarily fix your problem by deleting the word "Wunused" from
flags.guess file located in the parent folder. You need to reconfigure.

Thanks,
Hom Nath

On Tue, May 30, 2017 at 3:15 AM, tu xiang <tuxiang2016 at outlook.com> wrote:

>
> Good afternoon
>
>
> I alse cannot download SPECFEM3D Cartesian from the CIG website.
>
> I have tried this command :
>
> git clone --recursive https://github.com/geodynamics/specfem3d.git
>
>
> The output of this command is :
>
> $ git clone --recursive https://github.com/geodynamics/specfem3d.git
> Initialized empty Git repository in /home/specfem3d/.git/
> remote: Counting objects: 38298, done.
> remote: Compressing objects: 100% (25/25), done.
> remote: Total 38298 (delta 12), reused 23 (delta 11), pack-reused 38261
> Receiving objects: 100% (38298/38298), 372.06 MiB | 2.38 MiB/s, done.
> Resolving deltas: 100% (29156/29156), done.
> Submodule 'm4' (https://github.com/geodynamics/autoconf_cig.git)
> registered for path 'm4'
> Submodule 'utils/ADJOINT_TOMOGRAPHY_TOOLS/flexwin' (https://github.com/
> geodynamics/flexwin.git) registered for path 'utils/ADJOINT_TOMOGRAPHY_
> TOOLS/flexwin'
> Initialized empty Git repository in /home/specfem3d/m4/.git/
> remote: Counting objects: 403, done.
> remote: Total 403 (delta 0), reused 0 (delta 0), pack-reused 403
> Receiving objects: 100% (403/403), 127.11 KiB | 201 KiB/s, done.
> Resolving deltas: 100% (224/224), done.
> Submodule path 'm4': checked out '9912ff1c7f7e66577d849160c4169c
> 0c3be4806f'
> Initialized empty Git repository in /home/specfem3d/utils/ADJOINT_
> TOMOGRAPHY_TOOLS/flexwin/.git/
> remote: Counting objects: 830, done.
> remote: Total 830 (delta 0), reused 0 (delta 0), pack-reused 830
> Receiving objects: 100% (830/830), 2.40 MiB | 577 KiB/s, done.
> Resolving deltas: 100% (517/517), done.
> Submodule path 'utils/ADJOINT_TOMOGRAPHY_TOOLS/flexwin': checked out '
> fde3e0993c1853d98bd564e0cf13d796afe16cbe'
>
>
> I don't know why, and I tried many times on different computers.
>
>
> Then I download one version of SPECFEM3D Cartesian from github
>
> I do not have ifort compiler on my computer, so I want to compile the
> software with gfortran.
>
> I installed gfortran and openmpi.
>
> Then, I run next command:
>
> ./configure FC=gfortran MPIFC=mpif90
>
>
> Next, I run:
>
> $ make xcheck_mesh_quality_CUBIT_Abaqus
> output:
> make: *** No rule to make target `xcheck_mesh_quality_CUBIT_Abaqus'.
> Stop.
>
> But the command:
>
> make xdecompose_mesh and make xgenerate_databases
>
> can run successfully.
>
>
> However, the command:
>
> make xspecfem3D
>
> cannot, the errors are:
>
>
> ATTENUATION,deltat,PML_CONDITIONS, &
>                                                          1
> Warning: Unused dummy argument 'pml_conditions' at (1)
> make: *** [obj/compute_forces_viscoelastic.spec.o] Error 1
>
>
> I do not know how to deal with these problems.
>
> Hope someone give me a hand.
>
>
> Best regards
> Yours sincerely,
> Tu Xiang
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170602/53a1d99f/attachment.html>

From komatitsch at lma.cnrs-mrs.fr  Fri Jun  2 13:44:05 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Fri, 2 Jun 2017 22:44:05 +0200
Subject: [CIG-SEISMO] Some problems
In-Reply-To: <CAC+G=sNbDE9t5RUQffY3=-SKCtFHWu7XHhsXJ-ZTTMPiZN1dbA@mail.gmail.com>
References: <f40947bd0b73434283be01f10f8e6e3e@CSGHUB208W.pu.win.princeton.edu>
 <CAC+G=sNbDE9t5RUQffY3=-SKCtFHWu7XHhsXJ-ZTTMPiZN1dbA@mail.gmail.com>
Message-ID: <ffc6009f-9925-8e30-a324-242150db6f2f@lma.cnrs-mrs.fr>


Dear all,

Thanks!

I have synchronized the main branch with the devel branch, thus both are 
now OK.

Best regards,
Dimitri.

On 06/02/2017 05:15 PM, Hom Nath Gharti wrote:
> Dear Tu Xiang,
> 
> It seems that there is an unused variable in the 
> compute_forces_viscoelastic.f90 of the master branch. That is fixed in 
> the development branch.
> 
> You can either try development branch
> 
> OR
> 
> temporarily fix your problem by deleting the word "Wunused" from 
> flags.guess file located in the parent folder. You need to reconfigure.
> 
> Thanks,
> Hom Nath
> 
> On Tue, May 30, 2017 at 3:15 AM, tu xiang <tuxiang2016 at outlook.com 
> <mailto:tuxiang2016 at outlook.com>> wrote:
> 
> 
>     Good afternoon
> 
> 
>     I alse cannot download SPECFEM3D Cartesian from the CIG website.
> 
>     I have tried this command :
> 
>     git clone --recursive https://github.com/geodynamics/specfem3d.git
>     <https://github.com/geodynamics/specfem3d.git>
> 
> 
>     The output of this command is :
> 
>     $ git clone --recursive https://github.com/geodynamics/specfem3d.git
>     <https://github.com/geodynamics/specfem3d.git>
>     Initialized empty Git repository in /home/specfem3d/.git/
>     remote: Counting objects: 38298, done.
>     remote: Compressing objects: 100% (25/25), done.
>     remote: Total 38298 (delta 12), reused 23 (delta 11), pack-reused 38261
>     Receiving objects: 100% (38298/38298), 372.06 MiB | 2.38 MiB/s, done.
>     Resolving deltas: 100% (29156/29156), done.
>     Submodule 'm4' (https://github.com/geodynamics/autoconf_cig.git
>     <https://github.com/geodynamics/autoconf_cig.git>) registered for
>     path 'm4'
>     Submodule 'utils/ADJOINT_TOMOGRAPHY_TOOLS/flexwin'
>     (https://github.com/geodynamics/flexwin.git
>     <https://github.com/geodynamics/flexwin.git>) registered for path
>     'utils/ADJOINT_TOMOGRAPHY_TOOLS/flexwin'
>     Initialized empty Git repository in /home/specfem3d/m4/.git/
>     remote: Counting objects: 403, done.
>     remote: Total 403 (delta 0), reused 0 (delta 0), pack-reused 403
>     Receiving objects: 100% (403/403), 127.11 KiB | 201 KiB/s, done.
>     Resolving deltas: 100% (224/224), done.
>     Submodule path 'm4': checked out
>     '9912ff1c7f7e66577d849160c4169c0c3be4806f'
>     Initialized empty Git repository in
>     /home/specfem3d/utils/ADJOINT_TOMOGRAPHY_TOOLS/flexwin/.git/
>     remote: Counting objects: 830, done.
>     remote: Total 830 (delta 0), reused 0 (delta 0), pack-reused 830
>     Receiving objects: 100% (830/830), 2.40 MiB | 577 KiB/s, done.
>     Resolving deltas: 100% (517/517), done.
>     Submodule path 'utils/ADJOINT_TOMOGRAPHY_TOOLS/flexwin': checked out
>     'fde3e0993c1853d98bd564e0cf13d796afe16cbe'
> 
> 
>     I don't know why, and I tried many times on different computers.
> 
> 
>     Then I download one version of SPECFEM3D Cartesian from github
> 
>     I do not have ifort compiler on my computer, so I want to compile
>     the software with gfortran.
> 
>     I installed gfortran and openmpi.
> 
>     Then, I run next command:
> 
>     ./configure FC=gfortran MPIFC=mpif90
> 
> 
>     Next, I run:
> 
>     $ make xcheck_mesh_quality_CUBIT_Abaqus
>     output:
>     make: *** No rule to make target
>     `xcheck_mesh_quality_CUBIT_Abaqus'.  Stop.
> 
>     But the command:
> 
>     make xdecompose_mesh and make xgenerate_databases
> 
>     can run successfully.
> 
> 
>     However, the command:
> 
>     make xspecfem3D
> 
>     cannot, the errors are:
> 
> 
>     ATTENUATION,deltat,PML_CONDITIONS, &
>                                                               1
>     Warning: Unused dummy argument 'pml_conditions' at (1)
>     make: *** [obj/compute_forces_viscoelastic.spec.o] Error 1
> 
> 
>     I do not know how to deal with these problems.
> 
>     Hope someone give me a hand.
> 
> 
>     Best regards
>     Yours sincerely,
>     Tu Xiang
> 
> 
> 
> 
> 
> _______________________________________________
> 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

From mzhao at cea-igp.ac.cn  Sat Jun  3 05:05:32 2017
From: mzhao at cea-igp.ac.cn (=?gbk?B?1dTD9w==?=)
Date: Sat, 3 Jun 2017 20:05:32 +0800
Subject: [CIG-SEISMO] suspicious bug of the gravity calculations
Message-ID: <6e9d3efdc7d4454dcf71f8c728de87c1@cea-igp.ac.cn>

Hi Everyone,
I tried the gravity field calculation of SPECFEM3D using the given example: homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides.However,I found the results seem to be abnormal at some observations.
Since it is a homogenous model,the gravity value should be the same to the observation on the surface.However,look at the toto.ps,there are some place (very likely to be the boundary of some elements) the gravity is abnormal. 
Best regardsMing Zhao

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170603/4fb0910f/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: toto.ps
Type: application/postscript
Size: 58900 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170603/4fb0910f/attachment-0001.ps>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: results_norm_of_g_for_GMT.txt
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170603/4fb0910f/attachment-0002.txt>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: observation_grid_to_use_for_gravity.txt
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170603/4fb0910f/attachment-0003.txt>

From komatitsch at lma.cnrs-mrs.fr  Sat Jun  3 07:09:32 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Sat, 3 Jun 2017 16:09:32 +0200
Subject: [CIG-SEISMO] suspicious bug of the gravity calculations
In-Reply-To: <6e9d3efdc7d4454dcf71f8c728de87c1@cea-igp.ac.cn>
References: <6e9d3efdc7d4454dcf71f8c728de87c1@cea-igp.ac.cn>
Message-ID: <e329afdb-985a-f3a0-e2b0-29dd3f0609ce@lma.cnrs-mrs.fr>


Dear Ming, dear all,

Thanks! Let me cc Roland Martin, who is the expert on that.
However, it seems you used an observation surface that is located in z = 
0 i.e. at elevation zero, but the gravity formulas are undefined there 
(they involve a term with 1/r , in which r is the distance to the 
surface). Thus, could you try again using a non-zero elevation? 
everything should work fine then.

Also, I think the results should be symmetric, but not necessarily 
identical at all points (since the mesh has a finite size, you should 
see different values near the edges than in the center, in a symmetric 
fashion (I think, but I am not an expert; Roland will confirm).

Thanks,
Cheers,

Dimitri.

On 06/03/2017 02:05 PM, 赵明 wrote:
> Hi Everyone,
> 
> I tried the gravity field calculation of SPECFEM3D using the given 
> example: 
> homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides.However,I 
> found the results seem to be abnormal at some observations.
> 
> Since it is a homogenous model,the gravity value should be the same to 
> the observation on the surface.However,look at the toto.ps,there are 
> some place (very likely to be the boundary of some elements) the gravity 
> is abnormal.
> 
> Best regards
> Ming Zhao
> 
> 
> 
> 
> _______________________________________________
> 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

From mzhao at cea-igp.ac.cn  Sat Jun  3 23:04:08 2017
From: mzhao at cea-igp.ac.cn (=?utf-8?B?6LW15piO?=)
Date: Sun, 4 Jun 2017 14:04:08 +0800
Subject: [CIG-SEISMO] =?utf-8?b?5Zue5aSNOiBGd2Q6IFJFOiBzdXNwaWNpb3VzIGJ1?=
	=?utf-8?q?g_of_the_gravity_calculations?=
In-Reply-To: <fe18142d-1f62-2749-bb47-bce83ddcfdce@lma.cnrs-mrs.fr>
References: <A90E99DBCF207C4EA878BE64195F509711B2D810@CNMB03WVP.core-res.rootcore.local>, 
 <fe18142d-1f62-2749-bb47-bce83ddcfdce@lma.cnrs-mrs.fr>
Message-ID: <36f0fad090439ce166e286bbc57b1aa8@cea-igp.ac.cn>

Hi Dimitri,Roland,
Thank you very much for your reply and suggestion. the observation can't be at the surface of the model,I tried once again for two depths of the observations plane:z=-10m and z=-122 m. The -122m seems to be quite normal(toto_minus_122.ps), while z=-10m have abnormal values at the down left corner.
Cheers
Ming--------- 原始邮件信息 ---------
发件人: "Dimitri Komatitsch" <komatitsch at lma.cnrs-mrs.fr>
收件人: <mzhao at cea-igp.ac.cn>, "Roland Martin" <roland.martin at get.omp.eu>
主题: Fwd: RE:[CIG-SEISMO] suspicious bug of the gravity calculations
日期: 2017/06/03 23:46:11 (Sat)

Hi Ming,

I am forwarding you this email from Roland, because when he sends it to 
you from his account he gets an error message.

Roland, attached is the file toto.ps

Cheers,
Dimitri.

-------- Forwarded Message --------
Subject: RE:[CIG-SEISMO] suspicious bug of the gravity calculations
Date: Sat, 3 Jun 2017 14:21:15 +0000
From: MARTIN Roland <Roland.MARTIN at Get.omp.eu>
To: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>, 赵明 
<mzhao at cea-igp.ac.cn>

Hi Ming,

It is not possible to put obsevation points on the topography solution 
will become usntable there.
THis version is in revision and we need to work a little bit more on it 
before the release. I can not see the toto.ps because I did nit receive it.

Cheers
Roland
________________________________________
De : Dimitri Komatitsch [komatitsch at lma.cnrs-mrs.fr]
Envoyé : samedi 3 juin 2017 16:09
À : 赵明; cig-seismo at geodynamics.org; 赵明; MARTIN Roland; CHEVROT 
Sebastien; Vadim Monteiller
Objet : Re: [CIG-SEISMO] suspicious bug of the gravity calculations

Dear Ming, dear all,

Thanks! Let me cc Roland Martin, who is the expert on that.
However, it seems you used an observation surface that is located in z =
0 i.e. at elevation zero, but the gravity formulas are undefined there
(they involve a term with 1/r , in which r is the distance to the
surface). Thus, could you try again using a non-zero elevation?
everything should work fine then.

Also, I think the results should be symmetric, but not necessarily
identical at all points (since the mesh has a finite size, you should
see different values near the edges than in the center, in a symmetric
fashion (I think, but I am not an expert; Roland will confirm).

Thanks,
Cheers,

Dimitri.

On 06/03/2017 02:05 PM, 赵明 wrote:
> Hi Everyone,
>
> I tried the gravity field calculation of SPECFEM3D using the given
> example:
> homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides.However,I
> found the results seem to be abnormal at some observations.
>
> Since it is a homogenous model,the gravity value should be the same to
> the observation on the surface.However,look at the toto.ps,there are
> some place (very likely to be the boundary of some elements) the gravity
> is abnormal.
>
> Best regards
> Ming Zhao
>
>
>
>
> _______________________________________________
> 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170604/dae5ac51/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: toto_minus_122.ps
Type: application/postscript
Size: 95055 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170604/dae5ac51/attachment-0002.ps>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: toto_minus_10.ps
Type: application/postscript
Size: 80897 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170604/dae5ac51/attachment-0003.ps>

From komatitsch at lma.cnrs-mrs.fr  Mon Jun  5 04:39:35 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Mon, 5 Jun 2017 13:39:35 +0200
Subject: [CIG-SEISMO]
 =?utf-8?b?5Zue5aSNOiBGd2Q6IFJFOiBzdXNwaWNpb3VzIGJ1?=
 =?utf-8?q?g_of_the_gravity_calculations?=
In-Reply-To: <36f0fad090439ce166e286bbc57b1aa8@cea-igp.ac.cn>
References: <A90E99DBCF207C4EA878BE64195F509711B2D810@CNMB03WVP.core-res.rootcore.local>
 <fe18142d-1f62-2749-bb47-bce83ddcfdce@lma.cnrs-mrs.fr>
 <36f0fad090439ce166e286bbc57b1aa8@cea-igp.ac.cn>
Message-ID: <1b7d197f-c644-e2cf-96b3-c715f98bc22a@lma.cnrs-mrs.fr>


Hi Ming,

Thanks; regarding Z=-10m it depends how large your model is, the minimum 
Z value to use so that 1/z is not singular should probably be a value 
that is not too small compared to the typical size of your model, 
otherwise 1/z is OK mathematically but may be problematic from a 
numerical point of view (?).

Also, in most cases people compute the field outside of the model rather 
than inside, i.e. for positive values of Z. Please check e.g. 
http://komatitsch.free.fr/preprints/GJI_Roland_Martin_2013.pdf or 
references therein to see if the formulas remain OK for points inside 
the model rather than outside (I do not remember; Roland knows more 
about this than I do).

thanks
Cheers,
Dimitri.

On 06/04/2017 08:04 AM, 赵明 wrote:
> Hi Dimitri,Roland,
> 
> Thank you very much for your reply and suggestion. the observation can't 
> be at the surface of the model,I tried once again for two depths of the 
> observations plane:z=-10m and z=-122 m. The -122m seems to be quite 
> normal(toto_minus_122.ps), while z=-10m have abnormal values at the down 
> left corner.
> 
> Cheers
> 
> Ming
> --------- 原始邮件信息 ---------
> *发件人:* "Dimitri Komatitsch" <komatitsch at lma.cnrs-mrs.fr>
> *收件人:* <mzhao at cea-igp.ac.cn>, "Roland Martin" <roland.martin at get.omp.eu>
> *主题:* Fwd: RE:[CIG-SEISMO] suspicious bug of the gravity calculations
> *日期:* 2017/06/03 23:46:11 (Sat)
> 
> Hi Ming,
> 
> I am forwarding you this email from Roland, because when he sends it to
> you from his account he gets an error message.
> 
> Roland, attached is the file toto.ps
> 
> Cheers,
> Dimitri.
> 
> -------- Forwarded Message --------
> Subject: RE:[CIG-SEISMO] suspicious bug of the gravity calculations
> Date: Sat, 3 Jun 2017 14:21:15 +0000
> From: MARTIN Roland <Roland.MARTIN at Get.omp.eu>
> To: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>, 赵明
> <mzhao at cea-igp.ac.cn>
> 
> Hi Ming,
> 
> It is not possible to put obsevation points on the topography solution
> will become usntable there.
> THis version is in revision and we need to work a little bit more on it
> before the release. I can not see the toto.ps because I did nit receive it.
> 
> Cheers
> Roland
> ________________________________________
> De : Dimitri Komatitsch [komatitsch at lma.cnrs-mrs.fr]
> Envoyé : samedi 3 juin 2017 16:09
> À : 赵明; cig-seismo at geodynamics.org; 赵明; MARTIN Roland; CHEVROT
> Sebastien; Vadim Monteiller
> Objet : Re: [CIG-SEISMO] suspicious bug of the gravity calculations
> 
> Dear Ming, dear all,
> 
> Thanks! Let me cc Roland Martin, who is the expert on that.
> However, it seems you used an observation surface that is located in z =
> 0 i.e. at elevation zero, but the gravity formulas are undefined there
> (they involve a term with 1/r , in which r is the distance to the
> surface). Thus, could you try again using a non-zero elevation?
> everything should work fine then.
> 
> Also, I think the results should be symmetric, but not necessarily
> identical at all points (since the mesh has a finite size, you should
> see different values near the edges than in the center, in a symmetric
> fashion (I think, but I am not an expert; Roland will confirm).
> 
> Thanks,
> Cheers,
> 
> Dimitri.
> 
> On 06/03/2017 02:05 PM, 赵明 wrote:
>> Hi Everyone,
>>
>> I tried the gravity field calculation of SPECFEM3D using the given
>> example:
>> homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides.However,I
>> found the results seem to be abnormal at some observations.
>>
>> Since it is a homogenous model,the gravity value should be the same to
>> the observation on the surface.However,look at the toto.ps,there are
>> some place (very likely to be the boundary of some elements) the gravity
>> is abnormal.
>>
>> Best regards
>> Ming Zhao
>>
>>
>>
>>
>> _______________________________________________
>> 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

From komatitsch at lma.cnrs-mrs.fr  Tue Jun  6 09:26:33 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Tue, 6 Jun 2017 18:26:33 +0200
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
Message-ID: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>


Hi all,

Vadim Monteiller here in our group has created a nice integrated FWI 
tool for SPECFEM3D (based on L-BFGS and following 
http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf Sections 3 and 
4). The advantage is that it uses no I/Os to disk, all the inversion 
routines (computing sensitivity kernels, smoothing them, summing, 
iterating) are called from a single (parallel) calling program and thus 
only memory is used.
We have designed a set of examples, which we will put in BuildBot.

We will release it before the summer, Vadim is finishing a few final tests.

Best regards,
Dimitri.

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

From surendra at iith.ac.in  Tue Jun  6 09:31:09 2017
From: surendra at iith.ac.in (Surendra Nadh Somala)
Date: Tue, 6 Jun 2017 22:01:09 +0530
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
In-Reply-To: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
Message-ID: <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>

Dear Dimitri,

Are these tomographic inversions or can SPECFEM3D also be used for source
inversions?

Thanks,
Surendra


___________
*Surendra Nadh Somala*
Assistant Professor
Department of Civil Engineering (Block-E, 208)
Indian Institute of Technology (IIT) Hyderabad
Kandi, Sangareddy
Telangana, India - 502285
Phone : +91 (0)40 2301-8457
http://www.iith.ac.in/~surendra/

On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch <
komatitsch at lma.cnrs-mrs.fr> wrote:

>
> Hi all,
>
> Vadim Monteiller here in our group has created a nice integrated FWI tool
> for SPECFEM3D (based on L-BFGS and following
> http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf Sections 3 and
> 4). The advantage is that it uses no I/Os to disk, all the inversion
> routines (computing sensitivity kernels, smoothing them, summing,
> iterating) are called from a single (parallel) calling program and thus
> only memory is used.
> We have designed a set of examples, which we will put in BuildBot.
>
> We will release it before the summer, Vadim is finishing a few final tests.
>
> Best regards,
> Dimitri.
>
> --
> Dimitri Komatitsch, CNRS Research Director (DR CNRS)
> Laboratory of Mechanics and Acoustics, Marseille, France
> http://komatitsch.free.fr
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170606/98e01f15/attachment.html>

From tuxiang2016 at outlook.com  Tue Jun 13 19:59:38 2017
From: tuxiang2016 at outlook.com (tu xiang)
Date: Wed, 14 Jun 2017 02:59:38 +0000
Subject: [CIG-SEISMO] question about GEOCUBIT
Message-ID: <SG2PR0401MB1501A408511B840840D5BD13A1C30@SG2PR0401MB1501.apcprd04.prod.outlook.com>

Hi, everyone

I have successfully compiled the SPECFEM3D, and I have constructed my own mesh model with Trelis/CUBIT software in python script. However, while I export the mesh model to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT, I met some problems.


After I successfully compiled the xcheck_mesh_quality, xconvolve_source_timefunction, xdecompose_mesh,  xgenerate_databases,  xmeshfem3D, xspecfem3D code in SPECFEM3D and installed the Trelis software, then I added the path to my .bashrc file:

export PYTHONPATH=$PYTHONPATH:/home/soft/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib
export PATH=$PATH:/home/soft/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib


But when I export the mesh model with:

cubit2specfem3d.export2SPECFEM3D('MESH')

command,

it output errors:

Trelis>
Traceback (most recent call last):
File "<string>", line 1, in <module>
File "/home/bei/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py", line 1464, in export2SPECFEM3D
sem_mesh = mesh(hex27, cpml, cpml_size, top_absorbing)
File "/home/bei/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py", line 467, in __init__
self.block_definition()
File "/home/bei/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py", line 554, in block_definition
if qk < 0 or qmu < 0:
UnboundLocalError: local variable 'qk' referenced before assignment
Trelis>
Is anybody can help me?


I look forward to your response.
Best regards,

tuxiang




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170614/d21ab3de/attachment.html>

From komatitsch at lma.cnrs-mrs.fr  Sat Jun 17 03:32:32 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Sat, 17 Jun 2017 12:32:32 +0200
Subject: [CIG-SEISMO] question about GEOCUBIT
In-Reply-To: <SG2PR0401MB1501A408511B840840D5BD13A1C30@SG2PR0401MB1501.apcprd04.prod.outlook.com>
References: <SG2PR0401MB1501A408511B840840D5BD13A1C30@SG2PR0401MB1501.apcprd04.prod.outlook.com>
Message-ID: <14758bb1-b857-df34-652a-d33ad5619b2c@lma.cnrs-mrs.fr>


Hi Tuxiang,

Could you try again? the problem should now be fixed (otherwise please 
let us know; Emanuele is the main contact for GEOCUBIT). To get the new 
version of the code, just type "git pull".

Best regards,
Dimitri.

On 06/14/2017 04:59 AM, tu xiang wrote:
> Hi, everyone
> 
> I have successfully compiled the SPECFEM3D, and I have constructed my 
> own mesh model with Trelis/CUBIT software in python script. However, 
> while I export the mesh model to SPECFEM3D format using 
> cubit2specfem3d.py of GEOCUBIT, I met some problems.
> 
> 
> After I successfully compiled the xcheck_mesh_quality, 
> xconvolve_source_timefunction, xdecompose_mesh,  xgenerate_databases, 
>   xmeshfem3D, xspecfem3D code in SPECFEM3D and installed the Trelis 
> software, then I added the path to my .bashrc file:
> 
> export 
> PYTHONPATH=$PYTHONPATH:/home/soft/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib
> export PATH=$PATH:/home/soft/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib
> 
> But when I export the mesh model with:
> 
> cubit2specfem3d.export2SPECFEM3D('MESH')
> 
> command,
> 
> it output errors:
> 
> Trelis>
> Traceback (most recent call last):
> File "<string>", line 1, in <module>
> File 
> "/home/bei/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py", 
> line 1464, in export2SPECFEM3D
> sem_mesh = mesh(hex27, cpml, cpml_size, top_absorbing)
> File 
> "/home/bei/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py", 
> line 467, in __init__
> self.block_definition()
> File 
> "/home/bei/specfem3d-master/CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py", 
> line 554, in block_definition
> if qk < 0 or qmu < 0:
> UnboundLocalError: local variable 'qk' referenced before assignment
> Trelis>
> Is anybody can help me?
> 
> 
> I look forward to your response.
> Best regards,
> 
> tuxiang
> 
> 
> 
> 
> 
> 
> 
> _______________________________________________
> 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

From komatitsch at lma.cnrs-mrs.fr  Mon Jun 19 10:48:41 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Mon, 19 Jun 2017 19:48:41 +0200
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
In-Reply-To: <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
 <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
Message-ID: <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>


Hi Surendra, Hi all,

Vadim's tool currently inverts for the model only, mostly because we 
developed it for oil industry problems, in which case the location of 
the active source is known. Here is what SPECFEM3D users classically do 
to relocate earthquake sources (paragraph taken from a text written by 
Federica Magnoni and Emanuele Casarotti for Italy):

"As initial source parameters of the chosen events we use Time Domain 
Moment Tensor solutions (TDMT, Dreger and Helmberger, 1993; Scognamiglio 
et al. 2009) that are calculated by inverting full three-component 
traces of regional broad-band stations using 1D wave speed models.
Before starting the tomographic inversion, we then exploit the initial 
3D model of Italy of Di Stefano and Ciaccio (2014) to recalculate the 
source parameters of the considered events. The TDMT solutions, based on 
1D wave speed models, are considered as initial solutions and we invert 
for the six moment tensor components and the event location by using the 
technique presented in Liu et al. (2004). We use SPECFEM3D to 
numerically calculate the Fréchet derivatives with respect to the source 
parameters based on 3D models and we minimize a waveform misfit function 
to invert for source solutions. If the variance reduction of the 
solution with the 3D model is significant, this is considered as the new 
moment tensor solution of the event."

which comes from this paper by Qinya Liu: 
http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf and in 
which the derivatives are computed using finite-differences.

Emanuele and Federica (and Qinya), do you have scripts to do that?
if so, are they in the package somewhere (or should we consider 
committing them?); that would be very useful.

Thanks,
Best regards,
Dimitri.

On 06/06/2017 06:31 PM, Surendra Nadh Somala wrote:
> Dear Dimitri,
> 
> Are these tomographic inversions or can SPECFEM3D also be used for 
> source inversions?
> 
> Thanks,
> Surendra
> 
> 
> ___________
> *Surendra Nadh Somala*
> Assistant Professor
> Department of Civil Engineering (Block-E, 208)
> Indian Institute of Technology (IIT) Hyderabad
> Kandi, Sangareddy
> Telangana, India - 502285
> Phone : +91 (0)40 2301-8457
> http://www.iith.ac.in/~surendra/
> 
> On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch 
> <komatitsch at lma.cnrs-mrs.fr <mailto:komatitsch at lma.cnrs-mrs.fr>> wrote:
> 
> 
>     Hi all,
> 
>     Vadim Monteiller here in our group has created a nice integrated FWI
>     tool for SPECFEM3D (based on L-BFGS and following
>     http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>     <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf> Sections 3
>     and 4). The advantage is that it uses no I/Os to disk, all the
>     inversion routines (computing sensitivity kernels, smoothing them,
>     summing, iterating) are called from a single (parallel) calling
>     program and thus only memory is used.
>     We have designed a set of examples, which we will put in BuildBot.
> 
>     We will release it before the summer, Vadim is finishing a few final
>     tests.
> 
>     Best regards,
>     Dimitri.
> 
>     -- 
>     Dimitri Komatitsch, CNRS Research Director (DR CNRS)
>     Laboratory of Mechanics and Acoustics, Marseille, France
>     http://komatitsch.free.fr
>     _______________________________________________
>     CIG-SEISMO mailing list
>     CIG-SEISMO at geodynamics.org <mailto:CIG-SEISMO at geodynamics.org>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>
> 
> 
> 
> 
> _______________________________________________
> 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

From emanuele.casarotti at gmail.com  Mon Jun 19 10:53:42 2017
From: emanuele.casarotti at gmail.com (emanuele casarotti)
Date: Mon, 19 Jun 2017 17:53:42 +0000
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
In-Reply-To: <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
 <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
 <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>
Message-ID: <CA+6MF__H4_bwT+_pr8v1PLmiNGpX2pJ9gj8Bsa-KKLpmChYARw@mail.gmail.com>

You should consider pycmt3d https://github.com/wjlei1990/pycmt3d

Il giorno lun 19 giu 2017 alle 19:48 Dimitri Komatitsch <
komatitsch at lma.cnrs-mrs.fr> ha scritto:

>
> Hi Surendra, Hi all,
>
> Vadim's tool currently inverts for the model only, mostly because we
> developed it for oil industry problems, in which case the location of
> the active source is known. Here is what SPECFEM3D users classically do
> to relocate earthquake sources (paragraph taken from a text written by
> Federica Magnoni and Emanuele Casarotti for Italy):
>
> "As initial source parameters of the chosen events we use Time Domain
> Moment Tensor solutions (TDMT, Dreger and Helmberger, 1993; Scognamiglio
> et al. 2009) that are calculated by inverting full three-component
> traces of regional broad-band stations using 1D wave speed models.
> Before starting the tomographic inversion, we then exploit the initial
> 3D model of Italy of Di Stefano and Ciaccio (2014) to recalculate the
> source parameters of the considered events. The TDMT solutions, based on
> 1D wave speed models, are considered as initial solutions and we invert
> for the six moment tensor components and the event location by using the
> technique presented in Liu et al. (2004). We use SPECFEM3D to
> numerically calculate the Fréchet derivatives with respect to the source
> parameters based on 3D models and we minimize a waveform misfit function
> to invert for source solutions. If the variance reduction of the
> solution with the 3D model is significant, this is considered as the new
> moment tensor solution of the event."
>
> which comes from this paper by Qinya Liu:
> http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf and in
> which the derivatives are computed using finite-differences.
>
> Emanuele and Federica (and Qinya), do you have scripts to do that?
> if so, are they in the package somewhere (or should we consider
> committing them?); that would be very useful.
>
> Thanks,
> Best regards,
> Dimitri.
>
> On 06/06/2017 06:31 PM, Surendra Nadh Somala wrote:
> > Dear Dimitri,
> >
> > Are these tomographic inversions or can SPECFEM3D also be used for
> > source inversions?
> >
> > Thanks,
> > Surendra
> >
> >
> > ___________
> > *Surendra Nadh Somala*
> > Assistant Professor
> > Department of Civil Engineering (Block-E, 208)
> > Indian Institute of Technology (IIT) Hyderabad
> > Kandi, Sangareddy
> > Telangana, India - 502285
> > Phone : +91 (0)40 2301-8457
> > http://www.iith.ac.in/~surendra/
> >
> > On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch
> > <komatitsch at lma.cnrs-mrs.fr <mailto:komatitsch at lma.cnrs-mrs.fr>> wrote:
> >
> >
> >     Hi all,
> >
> >     Vadim Monteiller here in our group has created a nice integrated FWI
> >     tool for SPECFEM3D (based on L-BFGS and following
> >     http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
> >     <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf> Sections 3
> >     and 4). The advantage is that it uses no I/Os to disk, all the
> >     inversion routines (computing sensitivity kernels, smoothing them,
> >     summing, iterating) are called from a single (parallel) calling
> >     program and thus only memory is used.
> >     We have designed a set of examples, which we will put in BuildBot.
> >
> >     We will release it before the summer, Vadim is finishing a few final
> >     tests.
> >
> >     Best regards,
> >     Dimitri.
> >
> >     --
> >     Dimitri Komatitsch, CNRS Research Director (DR CNRS)
> >     Laboratory of Mechanics and Acoustics, Marseille, France
> >     http://komatitsch.free.fr
> >     _______________________________________________
> >     CIG-SEISMO mailing list
> >     CIG-SEISMO at geodynamics.org <mailto:CIG-SEISMO at geodynamics.org>
> >     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
> >     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>
> >
> >
> >
> >
> > _______________________________________________
> > 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170619/5ca516f0/attachment.html>

From komatitsch at lma.cnrs-mrs.fr  Tue Jun 20 06:40:51 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Tue, 20 Jun 2017 15:40:51 +0200
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
In-Reply-To: <CAEX8VZ=3CjE4FkHFOQYMbgeXHpdMSfvFB=uLXyUrZdMX+Z9B5Q@mail.gmail.com>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
 <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
 <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>
 <CA+6MF__H4_bwT+_pr8v1PLmiNGpX2pJ9gj8Bsa-KKLpmChYARw@mail.gmail.com>
 <CAEX8VZ=3CjE4FkHFOQYMbgeXHpdMSfvFB=uLXyUrZdMX+Z9B5Q@mail.gmail.com>
Message-ID: <39b6f258-8166-4de2-1903-0f30106bf63f@lma.cnrs-mrs.fr>


Hi Federica and Emanuele, Hi all,

Thanks for your answers. Very useful. I am going to cut and paste them 
in the manual of SPECFEM and put a copy of cmt3d and of pycmt3d in the 
Git version of SPECFEM (they are compatible in terms of license, pycmt3d 
is LGPLv3 and SPECFEM is GPLv2+, which makes the combination GPLv3 
according to Wikipedia).

best wishes,
Dimitri.

On 06/20/2017 10:54 AM, Federica Magnoni wrote:
> Hi all,
> 
> Initially we used the fortran code cmt3d 
> (https://github.com/UTCompSeismo/GRD_CMT3D/tree/master/cmt3d) and 
> recently, as Emanuele said, we switched to the python code pycmt3d.
> 
> Best regards,
> Federica
> 
> 2017-06-19 19:53 GMT+02:00 emanuele casarotti 
> <emanuele.casarotti at gmail.com <mailto:emanuele.casarotti at gmail.com>>:
> 
>     You should consider pycmt3d https://github.com/wjlei1990/pycmt3d
>     <https://github.com/wjlei1990/pycmt3d>
> 
> 
>     Il giorno lun 19 giu 2017 alle 19:48 Dimitri Komatitsch
>     <komatitsch at lma.cnrs-mrs.fr <mailto:komatitsch at lma.cnrs-mrs.fr>> ha
>     scritto:
> 
> 
>         Hi Surendra, Hi all,
> 
>         Vadim's tool currently inverts for the model only, mostly because we
>         developed it for oil industry problems, in which case the
>         location of
>         the active source is known. Here is what SPECFEM3D users
>         classically do
>         to relocate earthquake sources (paragraph taken from a text
>         written by
>         Federica Magnoni and Emanuele Casarotti for Italy):
> 
>         "As initial source parameters of the chosen events we use Time
>         Domain
>         Moment Tensor solutions (TDMT, Dreger and Helmberger, 1993;
>         Scognamiglio
>         et al. 2009) that are calculated by inverting full three-component
>         traces of regional broad-band stations using 1D wave speed models.
>         Before starting the tomographic inversion, we then exploit the
>         initial
>         3D model of Italy of Di Stefano and Ciaccio (2014) to
>         recalculate the
>         source parameters of the considered events. The TDMT solutions,
>         based on
>         1D wave speed models, are considered as initial solutions and we
>         invert
>         for the six moment tensor components and the event location by
>         using the
>         technique presented in Liu et al. (2004). We use SPECFEM3D to
>         numerically calculate the Fréchet derivatives with respect to
>         the source
>         parameters based on 3D models and we minimize a waveform misfit
>         function
>         to invert for source solutions. If the variance reduction of the
>         solution with the 3D model is significant, this is considered as
>         the new
>         moment tensor solution of the event."
> 
>         which comes from this paper by Qinya Liu:
>         http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf <http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf>
>         and in
>         which the derivatives are computed using finite-differences.
> 
>         Emanuele and Federica (and Qinya), do you have scripts to do that?
>         if so, are they in the package somewhere (or should we consider
>         committing them?); that would be very useful.
> 
>         Thanks,
>         Best regards,
>         Dimitri.
> 
>         On 06/06/2017 06:31 PM, Surendra Nadh Somala wrote:
>          > Dear Dimitri,
>          >
>          > Are these tomographic inversions or can SPECFEM3D also be
>         used for
>          > source inversions?
>          >
>          > Thanks,
>          > Surendra
>          >
>          >
>          > ___________
>          > *Surendra Nadh Somala*
>          > Assistant Professor
>          > Department of Civil Engineering (Block-E, 208)
>          > Indian Institute of Technology (IIT) Hyderabad
>          > Kandi, Sangareddy
>          > Telangana, India - 502285
>          > Phone : +91 (0)40 2301-8457 <tel:+91%2040%202301%208457>
>          > http://www.iith.ac.in/~surendra/
>         <http://www.iith.ac.in/~surendra/>
>          >
>          > On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch
>          > <komatitsch at lma.cnrs-mrs.fr
>         <mailto:komatitsch at lma.cnrs-mrs.fr>
>         <mailto:komatitsch at lma.cnrs-mrs.fr
>         <mailto:komatitsch at lma.cnrs-mrs.fr>>> wrote:
>          >
>          >
>          >     Hi all,
>          >
>          >     Vadim Monteiller here in our group has created a nice
>         integrated FWI
>          >     tool for SPECFEM3D (based on L-BFGS and following
>          > http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>
>          >     <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>>
>         Sections 3
>          >     and 4). The advantage is that it uses no I/Os to disk,
>         all the
>          >     inversion routines (computing sensitivity kernels,
>         smoothing them,
>          >     summing, iterating) are called from a single (parallel)
>         calling
>          >     program and thus only memory is used.
>          >     We have designed a set of examples, which we will put in
>         BuildBot.
>          >
>          >     We will release it before the summer, Vadim is finishing
>         a few final
>          >     tests.
>          >
>          >     Best regards,
>          >     Dimitri.
>          >
>          >     --
>          >     Dimitri Komatitsch, CNRS Research Director (DR CNRS)
>          >     Laboratory of Mechanics and Acoustics, Marseille, France
>          > http://komatitsch.free.fr
>          >     _______________________________________________
>          >     CIG-SEISMO mailing list
>          > CIG-SEISMO at geodynamics.org
>         <mailto:CIG-SEISMO at geodynamics.org>
>         <mailto:CIG-SEISMO at geodynamics.org
>         <mailto:CIG-SEISMO at geodynamics.org>>
>          >
>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>         <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>
>          >   
>           <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>>
>          >
>          >
>          >
>          >
>          > _______________________________________________
>          > CIG-SEISMO mailing list
>          > CIG-SEISMO at geodynamics.org <mailto:CIG-SEISMO at geodynamics.org>
>          >
>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>         <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
> 
> 
> 
> 
> -- 
> Federica Magnoni
> Istituto Nazionale di Geofisica e Vulcanologia (INGV)
> Via di Vigna Murata, 605
> 00143 ROMA, Italy
> e-mail: federica.magnoni at ingv.it <mailto:federica.magnoni at ingv.it>
> 

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

From komatitsch at lma.cnrs-mrs.fr  Tue Jun 20 06:43:38 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Tue, 20 Jun 2017 15:43:38 +0200
Subject: [CIG-SEISMO] changing the license of SPECFEM from GPLv2+ to GPLv3
In-Reply-To: <39b6f258-8166-4de2-1903-0f30106bf63f@lma.cnrs-mrs.fr>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
 <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
 <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>
 <CA+6MF__H4_bwT+_pr8v1PLmiNGpX2pJ9gj8Bsa-KKLpmChYARw@mail.gmail.com>
 <CAEX8VZ=3CjE4FkHFOQYMbgeXHpdMSfvFB=uLXyUrZdMX+Z9B5Q@mail.gmail.com>
 <39b6f258-8166-4de2-1903-0f30106bf63f@lma.cnrs-mrs.fr>
Message-ID: <3db11cd6-f4ba-5d95-ab42-08520d01472c@lma.cnrs-mrs.fr>


Hi all,

By the way it would make sense to switch all three SPECFEM packages from 
GPL2+ to GPL3. If nobody disagrees I will do that in 7 days.

(most new packages and Python tools are released under version 3, and we 
use or will use several of them: AxiSEM, ObsPy, pycmt3d etc).

Best regards,
Dimitri.

On 06/20/2017 03:40 PM, Dimitri Komatitsch wrote:
> 
> Hi Federica and Emanuele, Hi all,
> 
> Thanks for your answers. Very useful. I am going to cut and paste them 
> in the manual of SPECFEM and put a copy of cmt3d and of pycmt3d in the 
> Git version of SPECFEM (they are compatible in terms of license, pycmt3d 
> is LGPLv3 and SPECFEM is GPLv2+, which makes the combination GPLv3 
> according to Wikipedia).
> 
> best wishes,
> Dimitri.
> 
> On 06/20/2017 10:54 AM, Federica Magnoni wrote:
>> Hi all,
>>
>> Initially we used the fortran code cmt3d 
>> (https://github.com/UTCompSeismo/GRD_CMT3D/tree/master/cmt3d) and 
>> recently, as Emanuele said, we switched to the python code pycmt3d.
>>
>> Best regards,
>> Federica
>>
>> 2017-06-19 19:53 GMT+02:00 emanuele casarotti 
>> <emanuele.casarotti at gmail.com <mailto:emanuele.casarotti at gmail.com>>:
>>
>>     You should consider pycmt3d https://github.com/wjlei1990/pycmt3d
>>     <https://github.com/wjlei1990/pycmt3d>
>>
>>
>>     Il giorno lun 19 giu 2017 alle 19:48 Dimitri Komatitsch
>>     <komatitsch at lma.cnrs-mrs.fr <mailto:komatitsch at lma.cnrs-mrs.fr>> ha
>>     scritto:
>>
>>
>>         Hi Surendra, Hi all,
>>
>>         Vadim's tool currently inverts for the model only, mostly 
>> because we
>>         developed it for oil industry problems, in which case the
>>         location of
>>         the active source is known. Here is what SPECFEM3D users
>>         classically do
>>         to relocate earthquake sources (paragraph taken from a text
>>         written by
>>         Federica Magnoni and Emanuele Casarotti for Italy):
>>
>>         "As initial source parameters of the chosen events we use Time
>>         Domain
>>         Moment Tensor solutions (TDMT, Dreger and Helmberger, 1993;
>>         Scognamiglio
>>         et al. 2009) that are calculated by inverting full 
>> three-component
>>         traces of regional broad-band stations using 1D wave speed 
>> models.
>>         Before starting the tomographic inversion, we then exploit the
>>         initial
>>         3D model of Italy of Di Stefano and Ciaccio (2014) to
>>         recalculate the
>>         source parameters of the considered events. The TDMT solutions,
>>         based on
>>         1D wave speed models, are considered as initial solutions and we
>>         invert
>>         for the six moment tensor components and the event location by
>>         using the
>>         technique presented in Liu et al. (2004). We use SPECFEM3D to
>>         numerically calculate the Fréchet derivatives with respect to
>>         the source
>>         parameters based on 3D models and we minimize a waveform misfit
>>         function
>>         to invert for source solutions. If the variance reduction of the
>>         solution with the 3D model is significant, this is considered as
>>         the new
>>         moment tensor solution of the event."
>>
>>         which comes from this paper by Qinya Liu:
>>         
>> http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf 
>> <http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf>
>>         and in
>>         which the derivatives are computed using finite-differences.
>>
>>         Emanuele and Federica (and Qinya), do you have scripts to do 
>> that?
>>         if so, are they in the package somewhere (or should we consider
>>         committing them?); that would be very useful.
>>
>>         Thanks,
>>         Best regards,
>>         Dimitri.
>>
>>         On 06/06/2017 06:31 PM, Surendra Nadh Somala wrote:
>>          > Dear Dimitri,
>>          >
>>          > Are these tomographic inversions or can SPECFEM3D also be
>>         used for
>>          > source inversions?
>>          >
>>          > Thanks,
>>          > Surendra
>>          >
>>          >
>>          > ___________
>>          > *Surendra Nadh Somala*
>>          > Assistant Professor
>>          > Department of Civil Engineering (Block-E, 208)
>>          > Indian Institute of Technology (IIT) Hyderabad
>>          > Kandi, Sangareddy
>>          > Telangana, India - 502285
>>          > Phone : +91 (0)40 2301-8457 <tel:+91%2040%202301%208457>
>>          > http://www.iith.ac.in/~surendra/
>>         <http://www.iith.ac.in/~surendra/>
>>          >
>>          > On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch
>>          > <komatitsch at lma.cnrs-mrs.fr
>>         <mailto:komatitsch at lma.cnrs-mrs.fr>
>>         <mailto:komatitsch at lma.cnrs-mrs.fr
>>         <mailto:komatitsch at lma.cnrs-mrs.fr>>> wrote:
>>          >
>>          >
>>          >     Hi all,
>>          >
>>          >     Vadim Monteiller here in our group has created a nice
>>         integrated FWI
>>          >     tool for SPECFEM3D (based on L-BFGS and following
>>          > http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>
>>          >     <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>>
>>         Sections 3
>>          >     and 4). The advantage is that it uses no I/Os to disk,
>>         all the
>>          >     inversion routines (computing sensitivity kernels,
>>         smoothing them,
>>          >     summing, iterating) are called from a single (parallel)
>>         calling
>>          >     program and thus only memory is used.
>>          >     We have designed a set of examples, which we will put in
>>         BuildBot.
>>          >
>>          >     We will release it before the summer, Vadim is finishing
>>         a few final
>>          >     tests.
>>          >
>>          >     Best regards,
>>          >     Dimitri.
>>          >
>>          >     --
>>          >     Dimitri Komatitsch, CNRS Research Director (DR CNRS)
>>          >     Laboratory of Mechanics and Acoustics, Marseille, France
>>          > http://komatitsch.free.fr
>>          >     _______________________________________________
>>          >     CIG-SEISMO mailing list
>>          > CIG-SEISMO at geodynamics.org
>>         <mailto:CIG-SEISMO at geodynamics.org>
>>         <mailto:CIG-SEISMO at geodynamics.org
>>         <mailto:CIG-SEISMO at geodynamics.org>>
>>          >
>>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>>         
>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>
>>          >           
>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo 
>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>>
>>          >
>>          >
>>          >
>>          >
>>          > _______________________________________________
>>          > CIG-SEISMO mailing list
>>          > CIG-SEISMO at geodynamics.org <mailto:CIG-SEISMO at geodynamics.org>
>>          >
>>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>>         
>> <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
>>
>>
>>
>>
>> -- 
>> Federica Magnoni
>> Istituto Nazionale di Geofisica e Vulcanologia (INGV)
>> Via di Vigna Murata, 605
>> 00143 ROMA, Italy
>> e-mail: federica.magnoni at ingv.it <mailto:federica.magnoni at ingv.it>
>>
> 

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

From federica.magnoni at gmail.com  Tue Jun 20 01:54:27 2017
From: federica.magnoni at gmail.com (Federica Magnoni)
Date: Tue, 20 Jun 2017 10:54:27 +0200
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
In-Reply-To: <CA+6MF__H4_bwT+_pr8v1PLmiNGpX2pJ9gj8Bsa-KKLpmChYARw@mail.gmail.com>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
 <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
 <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>
 <CA+6MF__H4_bwT+_pr8v1PLmiNGpX2pJ9gj8Bsa-KKLpmChYARw@mail.gmail.com>
Message-ID: <CAEX8VZ=3CjE4FkHFOQYMbgeXHpdMSfvFB=uLXyUrZdMX+Z9B5Q@mail.gmail.com>

Hi all,

Initially we used the fortran code cmt3d (
https://github.com/UTCompSeismo/GRD_CMT3D/tree/master/cmt3d) and recently,
as Emanuele said, we switched to the python code pycmt3d.

Best regards,
Federica

2017-06-19 19:53 GMT+02:00 emanuele casarotti <emanuele.casarotti at gmail.com>
:

> You should consider pycmt3d https://github.com/wjlei1990/pycmt3d
>
>
> Il giorno lun 19 giu 2017 alle 19:48 Dimitri Komatitsch <
> komatitsch at lma.cnrs-mrs.fr> ha scritto:
>
>>
>> Hi Surendra, Hi all,
>>
>> Vadim's tool currently inverts for the model only, mostly because we
>> developed it for oil industry problems, in which case the location of
>> the active source is known. Here is what SPECFEM3D users classically do
>> to relocate earthquake sources (paragraph taken from a text written by
>> Federica Magnoni and Emanuele Casarotti for Italy):
>>
>> "As initial source parameters of the chosen events we use Time Domain
>> Moment Tensor solutions (TDMT, Dreger and Helmberger, 1993; Scognamiglio
>> et al. 2009) that are calculated by inverting full three-component
>> traces of regional broad-band stations using 1D wave speed models.
>> Before starting the tomographic inversion, we then exploit the initial
>> 3D model of Italy of Di Stefano and Ciaccio (2014) to recalculate the
>> source parameters of the considered events. The TDMT solutions, based on
>> 1D wave speed models, are considered as initial solutions and we invert
>> for the six moment tensor components and the event location by using the
>> technique presented in Liu et al. (2004). We use SPECFEM3D to
>> numerically calculate the Fréchet derivatives with respect to the source
>> parameters based on 3D models and we minimize a waveform misfit function
>> to invert for source solutions. If the variance reduction of the
>> solution with the 3D model is significant, this is considered as the new
>> moment tensor solution of the event."
>>
>> which comes from this paper by Qinya Liu:
>> http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf and in
>> which the derivatives are computed using finite-differences.
>>
>> Emanuele and Federica (and Qinya), do you have scripts to do that?
>> if so, are they in the package somewhere (or should we consider
>> committing them?); that would be very useful.
>>
>> Thanks,
>> Best regards,
>> Dimitri.
>>
>> On 06/06/2017 06:31 PM, Surendra Nadh Somala wrote:
>> > Dear Dimitri,
>> >
>> > Are these tomographic inversions or can SPECFEM3D also be used for
>> > source inversions?
>> >
>> > Thanks,
>> > Surendra
>> >
>> >
>> > ___________
>> > *Surendra Nadh Somala*
>> > Assistant Professor
>> > Department of Civil Engineering (Block-E, 208)
>> > Indian Institute of Technology (IIT) Hyderabad
>> > Kandi, Sangareddy
>> > Telangana, India - 502285
>> > Phone : +91 (0)40 2301-8457 <+91%2040%202301%208457>
>> > http://www.iith.ac.in/~surendra/
>> >
>> > On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch
>> > <komatitsch at lma.cnrs-mrs.fr <mailto:komatitsch at lma.cnrs-mrs.fr>> wrote:
>> >
>> >
>> >     Hi all,
>> >
>> >     Vadim Monteiller here in our group has created a nice integrated FWI
>> >     tool for SPECFEM3D (based on L-BFGS and following
>> >     http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>> >     <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf> Sections
>> 3
>> >     and 4). The advantage is that it uses no I/Os to disk, all the
>> >     inversion routines (computing sensitivity kernels, smoothing them,
>> >     summing, iterating) are called from a single (parallel) calling
>> >     program and thus only memory is used.
>> >     We have designed a set of examples, which we will put in BuildBot.
>> >
>> >     We will release it before the summer, Vadim is finishing a few final
>> >     tests.
>> >
>> >     Best regards,
>> >     Dimitri.
>> >
>> >     --
>> >     Dimitri Komatitsch, CNRS Research Director (DR CNRS)
>> >     Laboratory of Mechanics and Acoustics, Marseille, France
>> >     http://komatitsch.free.fr
>> >     _______________________________________________
>> >     CIG-SEISMO mailing list
>> >     CIG-SEISMO at geodynamics.org <mailto:CIG-SEISMO at geodynamics.org>
>> >     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>> >     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>
>> >
>> >
>> >
>> >
>> > _______________________________________________
>> > 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
>>
>


-- 
Federica Magnoni
Istituto Nazionale di Geofisica e Vulcanologia (INGV)
Via di Vigna Murata, 605
00143 ROMA, Italy
e-mail: federica.magnoni at ingv.it
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170620/4d0d5166/attachment.html>

From komatitsch at lma.cnrs-mrs.fr  Thu Jun 22 09:18:55 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Thu, 22 Jun 2017 18:18:55 +0200
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
In-Reply-To: <39b6f258-8166-4de2-1903-0f30106bf63f@lma.cnrs-mrs.fr>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
 <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
 <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>
 <CA+6MF__H4_bwT+_pr8v1PLmiNGpX2pJ9gj8Bsa-KKLpmChYARw@mail.gmail.com>
 <CAEX8VZ=3CjE4FkHFOQYMbgeXHpdMSfvFB=uLXyUrZdMX+Z9B5Q@mail.gmail.com>
 <39b6f258-8166-4de2-1903-0f30106bf63f@lma.cnrs-mrs.fr>
Message-ID: <30170089-9577-9a73-85db-71641521b672@lma.cnrs-mrs.fr>


Hi all,

Done!
(added to specfem3d/src/inverse_problem_for_source;
Vadim's unified FWI framework for SPECFEM3D will go in 
specfem3d/src/inverse_problem_for_model, we will commit it on Monday).

Best regards,
Thanks again to Brad for the information on how to add git submodules!

Dimitri.

On 06/20/2017 03:40 PM, Dimitri Komatitsch wrote:
> 
> Hi Federica and Emanuele, Hi all,
> 
> Thanks for your answers. Very useful. I am going to cut and paste them 
> in the manual of SPECFEM and put a copy of cmt3d and of pycmt3d in the 
> Git version of SPECFEM (they are compatible in terms of license, pycmt3d 
> is LGPLv3 and SPECFEM is GPLv2+, which makes the combination GPLv3 
> according to Wikipedia).
> 
> best wishes,
> Dimitri.
> 
> On 06/20/2017 10:54 AM, Federica Magnoni wrote:
>> Hi all,
>>
>> Initially we used the fortran code cmt3d 
>> (https://github.com/UTCompSeismo/GRD_CMT3D/tree/master/cmt3d) and 
>> recently, as Emanuele said, we switched to the python code pycmt3d.
>>
>> Best regards,
>> Federica
>>
>> 2017-06-19 19:53 GMT+02:00 emanuele casarotti 
>> <emanuele.casarotti at gmail.com <mailto:emanuele.casarotti at gmail.com>>:
>>
>>     You should consider pycmt3d https://github.com/wjlei1990/pycmt3d
>>     <https://github.com/wjlei1990/pycmt3d>
>>
>>
>>     Il giorno lun 19 giu 2017 alle 19:48 Dimitri Komatitsch
>>     <komatitsch at lma.cnrs-mrs.fr <mailto:komatitsch at lma.cnrs-mrs.fr>> ha
>>     scritto:
>>
>>
>>         Hi Surendra, Hi all,
>>
>>         Vadim's tool currently inverts for the model only, mostly 
>> because we
>>         developed it for oil industry problems, in which case the
>>         location of
>>         the active source is known. Here is what SPECFEM3D users
>>         classically do
>>         to relocate earthquake sources (paragraph taken from a text
>>         written by
>>         Federica Magnoni and Emanuele Casarotti for Italy):
>>
>>         "As initial source parameters of the chosen events we use Time
>>         Domain
>>         Moment Tensor solutions (TDMT, Dreger and Helmberger, 1993;
>>         Scognamiglio
>>         et al. 2009) that are calculated by inverting full 
>> three-component
>>         traces of regional broad-band stations using 1D wave speed 
>> models.
>>         Before starting the tomographic inversion, we then exploit the
>>         initial
>>         3D model of Italy of Di Stefano and Ciaccio (2014) to
>>         recalculate the
>>         source parameters of the considered events. The TDMT solutions,
>>         based on
>>         1D wave speed models, are considered as initial solutions and we
>>         invert
>>         for the six moment tensor components and the event location by
>>         using the
>>         technique presented in Liu et al. (2004). We use SPECFEM3D to
>>         numerically calculate the Fréchet derivatives with respect to
>>         the source
>>         parameters based on 3D models and we minimize a waveform misfit
>>         function
>>         to invert for source solutions. If the variance reduction of the
>>         solution with the 3D model is significant, this is considered as
>>         the new
>>         moment tensor solution of the event."
>>
>>         which comes from this paper by Qinya Liu:
>>         
>> http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf 
>> <http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf>
>>         and in
>>         which the derivatives are computed using finite-differences.
>>
>>         Emanuele and Federica (and Qinya), do you have scripts to do 
>> that?
>>         if so, are they in the package somewhere (or should we consider
>>         committing them?); that would be very useful.
>>
>>         Thanks,
>>         Best regards,
>>         Dimitri.
>>
>>         On 06/06/2017 06:31 PM, Surendra Nadh Somala wrote:
>>          > Dear Dimitri,
>>          >
>>          > Are these tomographic inversions or can SPECFEM3D also be
>>         used for
>>          > source inversions?
>>          >
>>          > Thanks,
>>          > Surendra
>>          >
>>          >
>>          > ___________
>>          > *Surendra Nadh Somala*
>>          > Assistant Professor
>>          > Department of Civil Engineering (Block-E, 208)
>>          > Indian Institute of Technology (IIT) Hyderabad
>>          > Kandi, Sangareddy
>>          > Telangana, India - 502285
>>          > Phone : +91 (0)40 2301-8457 <tel:+91%2040%202301%208457>
>>          > http://www.iith.ac.in/~surendra/
>>         <http://www.iith.ac.in/~surendra/>
>>          >
>>          > On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch
>>          > <komatitsch at lma.cnrs-mrs.fr
>>         <mailto:komatitsch at lma.cnrs-mrs.fr>
>>         <mailto:komatitsch at lma.cnrs-mrs.fr
>>         <mailto:komatitsch at lma.cnrs-mrs.fr>>> wrote:
>>          >
>>          >
>>          >     Hi all,
>>          >
>>          >     Vadim Monteiller here in our group has created a nice
>>         integrated FWI
>>          >     tool for SPECFEM3D (based on L-BFGS and following
>>          > http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>
>>          >     <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>>
>>         Sections 3
>>          >     and 4). The advantage is that it uses no I/Os to disk,
>>         all the
>>          >     inversion routines (computing sensitivity kernels,
>>         smoothing them,
>>          >     summing, iterating) are called from a single (parallel)
>>         calling
>>          >     program and thus only memory is used.
>>          >     We have designed a set of examples, which we will put in
>>         BuildBot.
>>          >
>>          >     We will release it before the summer, Vadim is finishing
>>         a few final
>>          >     tests.
>>          >
>>          >     Best regards,
>>          >     Dimitri.
>>          >
>>          >     --
>>          >     Dimitri Komatitsch, CNRS Research Director (DR CNRS)
>>          >     Laboratory of Mechanics and Acoustics, Marseille, France
>>          > http://komatitsch.free.fr
>>          >     _______________________________________________
>>          >     CIG-SEISMO mailing list
>>          > CIG-SEISMO at geodynamics.org
>>         <mailto:CIG-SEISMO at geodynamics.org>
>>         <mailto:CIG-SEISMO at geodynamics.org
>>         <mailto:CIG-SEISMO at geodynamics.org>>
>>          >
>>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>>         
>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>
>>          >           
>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo 
>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>>
>>          >
>>          >
>>          >
>>          >
>>          > _______________________________________________
>>          > CIG-SEISMO mailing list
>>          > CIG-SEISMO at geodynamics.org <mailto:CIG-SEISMO at geodynamics.org>
>>          >
>>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>>         
>> <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
>>
>>
>>
>>
>> -- 
>> Federica Magnoni
>> Istituto Nazionale di Geofisica e Vulcanologia (INGV)
>> Via di Vigna Murata, 605
>> 00143 ROMA, Italy
>> e-mail: federica.magnoni at ingv.it <mailto:federica.magnoni at ingv.it>
>>
> 

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

From komatitsch at lma.cnrs-mrs.fr  Thu Jun 22 10:47:29 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Thu, 22 Jun 2017 19:47:29 +0200
Subject: [CIG-SEISMO] changing the license of SPECFEM from GPLv2+ to
	GPLv3
In-Reply-To: <3db11cd6-f4ba-5d95-ab42-08520d01472c@lma.cnrs-mrs.fr>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
 <CAK6kgGhc0Z6qMesd0Xs8x59EYxAq1UtH1-QBK==pofUN-6jusw@mail.gmail.com>
 <a631d2f1-e2c1-ada3-f3b8-da41d08c004a@lma.cnrs-mrs.fr>
 <CA+6MF__H4_bwT+_pr8v1PLmiNGpX2pJ9gj8Bsa-KKLpmChYARw@mail.gmail.com>
 <CAEX8VZ=3CjE4FkHFOQYMbgeXHpdMSfvFB=uLXyUrZdMX+Z9B5Q@mail.gmail.com>
 <39b6f258-8166-4de2-1903-0f30106bf63f@lma.cnrs-mrs.fr>
 <3db11cd6-f4ba-5d95-ab42-08520d01472c@lma.cnrs-mrs.fr>
Message-ID: <eb9b5214-d31d-0346-863f-fcc81a12ffc5@lma.cnrs-mrs.fr>


Hi all,

Done!

Dimitri.

On 06/20/2017 03:43 PM, Dimitri Komatitsch wrote:
> 
> Hi all,
> 
> By the way it would make sense to switch all three SPECFEM packages from 
> GPL2+ to GPL3. If nobody disagrees I will do that in 7 days.
> 
> (most new packages and Python tools are released under version 3, and we 
> use or will use several of them: AxiSEM, ObsPy, pycmt3d etc).
> 
> Best regards,
> Dimitri.
> 
> On 06/20/2017 03:40 PM, Dimitri Komatitsch wrote:
>>
>> Hi Federica and Emanuele, Hi all,
>>
>> Thanks for your answers. Very useful. I am going to cut and paste them 
>> in the manual of SPECFEM and put a copy of cmt3d and of pycmt3d in the 
>> Git version of SPECFEM (they are compatible in terms of license, 
>> pycmt3d is LGPLv3 and SPECFEM is GPLv2+, which makes the combination 
>> GPLv3 according to Wikipedia).
>>
>> best wishes,
>> Dimitri.
>>
>> On 06/20/2017 10:54 AM, Federica Magnoni wrote:
>>> Hi all,
>>>
>>> Initially we used the fortran code cmt3d 
>>> (https://github.com/UTCompSeismo/GRD_CMT3D/tree/master/cmt3d) and 
>>> recently, as Emanuele said, we switched to the python code pycmt3d.
>>>
>>> Best regards,
>>> Federica
>>>
>>> 2017-06-19 19:53 GMT+02:00 emanuele casarotti 
>>> <emanuele.casarotti at gmail.com <mailto:emanuele.casarotti at gmail.com>>:
>>>
>>>     You should consider pycmt3d https://github.com/wjlei1990/pycmt3d
>>>     <https://github.com/wjlei1990/pycmt3d>
>>>
>>>
>>>     Il giorno lun 19 giu 2017 alle 19:48 Dimitri Komatitsch
>>>     <komatitsch at lma.cnrs-mrs.fr <mailto:komatitsch at lma.cnrs-mrs.fr>> ha
>>>     scritto:
>>>
>>>
>>>         Hi Surendra, Hi all,
>>>
>>>         Vadim's tool currently inverts for the model only, mostly 
>>> because we
>>>         developed it for oil industry problems, in which case the
>>>         location of
>>>         the active source is known. Here is what SPECFEM3D users
>>>         classically do
>>>         to relocate earthquake sources (paragraph taken from a text
>>>         written by
>>>         Federica Magnoni and Emanuele Casarotti for Italy):
>>>
>>>         "As initial source parameters of the chosen events we use Time
>>>         Domain
>>>         Moment Tensor solutions (TDMT, Dreger and Helmberger, 1993;
>>>         Scognamiglio
>>>         et al. 2009) that are calculated by inverting full 
>>> three-component
>>>         traces of regional broad-band stations using 1D wave speed 
>>> models.
>>>         Before starting the tomographic inversion, we then exploit the
>>>         initial
>>>         3D model of Italy of Di Stefano and Ciaccio (2014) to
>>>         recalculate the
>>>         source parameters of the considered events. The TDMT solutions,
>>>         based on
>>>         1D wave speed models, are considered as initial solutions and we
>>>         invert
>>>         for the six moment tensor components and the event location by
>>>         using the
>>>         technique presented in Liu et al. (2004). We use SPECFEM3D to
>>>         numerically calculate the Fréchet derivatives with respect to
>>>         the source
>>>         parameters based on 3D models and we minimize a waveform misfit
>>>         function
>>>         to invert for source solutions. If the variance reduction of the
>>>         solution with the 3D model is significant, this is considered as
>>>         the new
>>>         moment tensor solution of the event."
>>>
>>>         which comes from this paper by Qinya Liu:
>>> http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf 
>>> <http://komatitsch.free.fr/preprints/bssa_Qinya_inversion_2004.pdf>
>>>         and in
>>>         which the derivatives are computed using finite-differences.
>>>
>>>         Emanuele and Federica (and Qinya), do you have scripts to do 
>>> that?
>>>         if so, are they in the package somewhere (or should we consider
>>>         committing them?); that would be very useful.
>>>
>>>         Thanks,
>>>         Best regards,
>>>         Dimitri.
>>>
>>>         On 06/06/2017 06:31 PM, Surendra Nadh Somala wrote:
>>>          > Dear Dimitri,
>>>          >
>>>          > Are these tomographic inversions or can SPECFEM3D also be
>>>         used for
>>>          > source inversions?
>>>          >
>>>          > Thanks,
>>>          > Surendra
>>>          >
>>>          >
>>>          > ___________
>>>          > *Surendra Nadh Somala*
>>>          > Assistant Professor
>>>          > Department of Civil Engineering (Block-E, 208)
>>>          > Indian Institute of Technology (IIT) Hyderabad
>>>          > Kandi, Sangareddy
>>>          > Telangana, India - 502285
>>>          > Phone : +91 (0)40 2301-8457 <tel:+91%2040%202301%208457>
>>>          > http://www.iith.ac.in/~surendra/
>>>         <http://www.iith.ac.in/~surendra/>
>>>          >
>>>          > On Tue, Jun 6, 2017 at 9:56 PM, Dimitri Komatitsch
>>>          > <komatitsch at lma.cnrs-mrs.fr
>>>         <mailto:komatitsch at lma.cnrs-mrs.fr>
>>>         <mailto:komatitsch at lma.cnrs-mrs.fr
>>>         <mailto:komatitsch at lma.cnrs-mrs.fr>>> wrote:
>>>          >
>>>          >
>>>          >     Hi all,
>>>          >
>>>          >     Vadim Monteiller here in our group has created a nice
>>>         integrated FWI
>>>          >     tool for SPECFEM3D (based on L-BFGS and following
>>>          > http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>>>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>
>>>          >     <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf
>>>         <http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf>>
>>>         Sections 3
>>>          >     and 4). The advantage is that it uses no I/Os to disk,
>>>         all the
>>>          >     inversion routines (computing sensitivity kernels,
>>>         smoothing them,
>>>          >     summing, iterating) are called from a single (parallel)
>>>         calling
>>>          >     program and thus only memory is used.
>>>          >     We have designed a set of examples, which we will put in
>>>         BuildBot.
>>>          >
>>>          >     We will release it before the summer, Vadim is finishing
>>>         a few final
>>>          >     tests.
>>>          >
>>>          >     Best regards,
>>>          >     Dimitri.
>>>          >
>>>          >     --
>>>          >     Dimitri Komatitsch, CNRS Research Director (DR CNRS)
>>>          >     Laboratory of Mechanics and Acoustics, Marseille, France
>>>          > http://komatitsch.free.fr
>>>          >     _______________________________________________
>>>          >     CIG-SEISMO mailing list
>>>          > CIG-SEISMO at geodynamics.org
>>>         <mailto:CIG-SEISMO at geodynamics.org>
>>>         <mailto:CIG-SEISMO at geodynamics.org
>>>         <mailto:CIG-SEISMO at geodynamics.org>>
>>>          >
>>>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>
>>>          > 
>>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo 
>>> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo>>
>>>          >
>>>          >
>>>          >
>>>          >
>>>          > _______________________________________________
>>>          > CIG-SEISMO mailing list
>>>          > CIG-SEISMO at geodynamics.org 
>>> <mailto:CIG-SEISMO at geodynamics.org>
>>>          >
>>>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>>> <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
>>>
>>>
>>>
>>>
>>> -- 
>>> Federica Magnoni
>>> Istituto Nazionale di Geofisica e Vulcanologia (INGV)
>>> Via di Vigna Murata, 605
>>> 00143 ROMA, Italy
>>> e-mail: federica.magnoni at ingv.it <mailto:federica.magnoni at ingv.it>
>>>
>>
> 

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

From komatitsch at lma.cnrs-mrs.fr  Sun Jun 25 09:12:59 2017
From: komatitsch at lma.cnrs-mrs.fr (Dimitri Komatitsch)
Date: Sun, 25 Jun 2017 18:12:59 +0200
Subject: [CIG-SEISMO] integrated inversion tool for SPECFEM3D (will be
 released before the summer)
In-Reply-To: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
References: <21d5b3d2-7c27-68a3-2e3e-62a7a983f03c@lma.cnrs-mrs.fr>
Message-ID: <98a27136-99e6-bd8b-870f-0f5e0057d4e4@lma.cnrs-mrs.fr>


Hi all,

It is now released
(in directory "src/inverse_problem_for_model").

When typing "make all" it is compiled automatically with the rest.

Best regards,
Dimitri.

On 06/06/2017 06:26 PM, Dimitri Komatitsch wrote:
> 
> Hi all,
> 
> Vadim Monteiller here in our group has created a nice integrated FWI 
> tool for SPECFEM3D (based on L-BFGS and following 
> http://komatitsch.free.fr/preprints/GJI2_Vadim_2015.pdf Sections 3 and 
> 4). The advantage is that it uses no I/Os to disk, all the inversion 
> routines (computing sensitivity kernels, smoothing them, summing, 
> iterating) are called from a single (parallel) calling program and thus 
> only memory is used.
> We have designed a set of examples, which we will put in BuildBot.
> 
> We will release it before the summer, Vadim is finishing a few final tests.
> 
> Best regards,
> Dimitri.
> 

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