Inverse of r_of_s. See there.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | r |
pure elemental function s_of_r(r) result(s) use mod_global_parameters double precision, intent(in) :: r double precision :: s if (log_r0 > zero) then s = dsign(one,r)*dlog(one + dabs(r)/log_r0) else s = dlog(r) end if end function s_of_r