I’m currently trying to model spontaneous earthquakes with rate-state friction in subduction zones. I read the section 7.10.9 of the v2.2.1 manual and ran stepo6 successfully. I noticed the reference slip rate used in step06 is 1.0cm/year.
However, in most of the literature I have read on the modeling of spontaneous earthquake in subduction zones, people tend to use 1.0e-6m/s as reference slip rate. Therefore, I tried to change the reference slip rate in step06 to 1.0e-6m/s(Figure 1). And I also changed the state variable to 5.0e+5*s. I didn’t make any other changes, but eventually the non linear solver fails to converge. I wonder what the problem is.