Skip to content

add flash heating into slip law#1

Open
Egao0206 wants to merge 14 commits intoHuihuiweng:thermpresfrom
Egao0206:thermpres
Open

add flash heating into slip law#1
Egao0206 wants to merge 14 commits intoHuihuiweng:thermpresfrom
Egao0206:thermpres

Conversation

@Egao0206
Copy link

Add regularized RSF (slip law with flash heating) in "bc_dynflt_rsf.f90", where the case is "kind = 5" . However, Compiled successfully, verified failed.

end type rsf_type

public :: rsf_type, rsf_read, rsf_init, rsf_mu, rsf_solver, rsf_qs_solver, rsf_timestep
public :: rsf_type, rsf_read, rsf_init, rsf_mu, rsf_solver, rsf_qs_solver, rsf_timestep, State_to_theta,Theta_to_state
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

There is no need to make these two function public, as they are only called within this subroutine.

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Fixed, removed from public

! as in Ampuero and Ben-Zion (2008)
! 2 = logarithmic rate-and-state with aging state law
! 3 = logarithmic rate-and-state with slip state law
! 4 = Regularized rate-and-state with slip state law, flash heating(Vw as parameter)
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

4 is the case with V-shape rate-and-state. You shall use 5 here.

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Fixed, updated case 4 and 5

endif
enddo
Q_ss = f%a * log((2 * f%Vstar / abs(v)) * sinh(mu_ss / f%a))
dQ = abs(v)/f%Dc * (Q - Q_ss)
Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Correct: dQ = abs(v) / f%Dc * (Q_ss - Q)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants