SPECFEM / specfem3d_globe

SPECFEM3D_GLOBE simulates global and regional (continental-scale) seismic wave propagation.
GNU General Public License v3.0
90 stars 96 forks source link

how to change the anisotropic mantle model? #634

Open PhilTsao opened 5 years ago

komatits commented 5 years ago

Dear PhilTsao,

It seems you have manually entered more than 21 anisotropy parameters.

Please also get a recent version of the SPECFEM3D_GLOBE code; there is no AMM_V parameter at the end of the routine in recent versions:

subroutine model_aniso_mantle(r,theta,phi,rho, c11,c12,c13,c14,c15,c16, & c22,c23,c24,c25,c26,c33,c34,c35,c36,c44,c45,c46,c55,c56,c66)

Best regards, Dimitri.

On 1/16/19 5:50 AM, PhilTsao wrote:

In the manual of spacfem3d globe I can see the steep of changing the anisotropic mantle model, but I don't understand the parameter of AMM_V in the command line of call model_aniso_mantle(r,theta,phi,rho, & c11,c12,c13,c14,c15,c16,c22,c23,c24,c25,c26, & c33,c34,c35,c36,c44,c45,c46,c55,c56,c66,AMM_V) when I run the command as call model_aniso_mantle(0.98743,1.065,1.789,3.7471,2236,840,859,0,0,0,840, 2236,859,0,0,0,859,859,2133,0,0,0,0,0,0,652,0,0,0,0,0,0,652,0,0,0,0,0,0,698,AMM_V) it returned bash: Unexpected sign `(' syntax error near please help me! thank you very much.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/geodynamics/specfem3d_globe/issues/634, or mute the thread https://github.com/notifications/unsubscribe-auth/AFjDKcZW8y9UTQgWYtQsYcfgboSfovmRks5vDq-kgaJpZM4aCT0t.

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

PhilTsao commented 5 years ago

THANK YOU~ Dear Komatits

PhilTsao commented 5 years ago

Thank you

发送自 Windows 10 版邮件https://go.microsoft.com/fwlink/?LinkId=550986应用


发件人: Dimitri Komatitsch notifications@github.com 发送时间: Wednesday, January 16, 2019 6:12:20 PM 收件人: geodynamics/specfem3d_globe 抄送: PhilTsao; Author 主题: Re: [geodynamics/specfem3d_globe] how to change the anisotropic mantle model? (#634)

Dear PhilTsao,

It seems you have manually entered more than 21 anisotropy parameters.

Please also get a recent version of the SPECFEM3D_GLOBE code; there is no AMM_V parameter at the end of the routine in recent versions:

subroutine model_aniso_mantle(r,theta,phi,rho, c11,c12,c13,c14,c15,c16, & c22,c23,c24,c25,c26,c33,c34,c35,c36,c44,c45,c46,c55,c56,c66)

Best regards, Dimitri.

On 1/16/19 5:50 AM, PhilTsao wrote:

In the manual of spacfem3d globe I can see the steep of changing the anisotropic mantle model, but I don't understand the parameter of AMM_V in the command line of call model_aniso_mantle(r,theta,phi,rho, & c11,c12,c13,c14,c15,c16,c22,c23,c24,c25,c26, & c33,c34,c35,c36,c44,c45,c46,c55,c56,c66,AMM_V) when I run the command as call model_aniso_mantle(0.98743,1.065,1.789,3.7471,2236,840,859,0,0,0,840, 2236,859,0,0,0,859,859,2133,0,0,0,0,0,0,652,0,0,0,0,0,0,652,0,0,0,0,0,0,698,AMM_V) it returned bash: Unexpected sign `(' syntax error near please help me! thank you very much.

― You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/geodynamics/specfem3d_globe/issues/634, or mute the thread https://github.com/notifications/unsubscribe-auth/AFjDKcZW8y9UTQgWYtQsYcfgboSfovmRks5vDq-kgaJpZM4aCT0t.

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

― You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/geodynamics/specfem3d_globe/issues/634#issuecomment-454724508, or mute the threadhttps://github.com/notifications/unsubscribe-auth/Aftd9Uy_XwABGP74AMPr50j3-FVgkftlks5vDvsEgaJpZM4aCT0t.

PhilTsao commented 5 years ago

sorry~ I try to manually enter the 21 anisotropy parameters in model_aniso_mantle.f90, but the parameters in model_aniso_mantle.f90 so complex. can you teach me step by step? thank you~