search

geodynamics / relax

Relax implements a semi-analytic Fourier-domain solver and equivalent body forces to compute quasi-static relaxation of stress perturbation.
GNU General Public License v3.0
25 stars 18 forks source link

chi-chi example can't get the same result of figure3a #3

Open ThreeIcug opened 2 years ago

ThreeIcug commented 2 years ago

Dear author, The Relax is a good method to research the postseismic. And the name of Relax is full of wisdom. Recently ,i use the chi-chi earthquake example (Rousset et al,2012)to study the code. I want to get the the same result of figure3a. so i use the coseismic.sh and change the # observation depth (dispalcement amd stress) 0 15 Then ,I use the grdmath method which is told in 6.4 chapter to calculate the norm of the deviatoric stress. However, I get result like figure, the value is smaller than the paper. I don't change other parameters. I was wonder if you can help me.

sbarbot commented 2 years ago

Figure 3a shows the norm of the coseismic stress change at 15 km depth. This can be reproduced with the input file "coseismic.sh" modified as you did. If the results are different, perhaps by a constant factor, that may be due to the definition of the norm. I have seen common definitions that vary by a factor of sqrt(2), for example. Please check the details of you calculation.

ThreeIcug commented 2 years ago

Thanks you for you help. I have reviewed a lot of literature, and I think the method in 6.4 chapter, which is used to calculate norm of the deviatoric stress (second invariant), is right or is normalised. Because that paper donot introduce the detail calculation formula, so i think that will wait for the help of the author. Thank you.