Help : the error when estimate the fault creep

[ThreeIcug]

Dear Sylvain,
Recently, when i use it to estimate my own case for fault creep, it take an error, just like this,

error at line 00203 of source file /home/rkkommu/src/relax/relax-git/src/friction3d.f90

wrong value of (a-b)sigma gives rise to
(a - b) * sigma = 1.000E-03
tau = 1.734E+00
tau_s = 1.571E+00
tau_n = 2.715E-01
tau / ((a - b) sigma) = 1.734E+03

I check the friction3d file , it said if (tau/tauc) .gt. 20 ,then will print the error tips.
I followed your paper of parkfield postseismic deformation, and my total slip similar to 0.5m. So I set the (a-b)sig from 0.2Mpa to 1.4Mpa. Because of using fault displacement in meters and dimension in kilometers , so I multiply 1E-3 and change it to Gpa. The reference rate is from 1 to 50 mm.
Is my understanding wrong? I was wondering if you can help me ,thank you.
Yours,
Jingwei