bond_ms_set Subroutine

private subroutine bond_ms_set(params, lp, rmax)

Setter for Marko-Siggia bond.

   U = [-(1/2)*rtilde^2 + 0.25/(1-rtilde)^2 + 0.25*rtilde]*(rmax/lp), r < rmax
   where rtilde = r/rmax.

User-set parameters:

  • params(1) = lp (persistence length)
  • params(2) = rmax

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(inout), dimension(:):: params
real(kind=rp), intent(in), optional :: lp
real(kind=rp), intent(in), optional :: rmax

Contents

None