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

[ThreeIcug]

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]

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]

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.