Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rp), | intent(in) | :: | R0 | |||
real(kind=rp), | intent(in) | :: | b | |||
real(kind=rp), | intent(in) | :: | w | |||
real(kind=rp), | intent(in) | :: | dlam | |||
real(kind=rp), | intent(in) | :: | R | |||
real(kind=rp), | intent(in) | :: | Z | |||
real(kind=rp), | intent(in) | :: | T |
function Spong_2D(R0,b,w,dlam,R,Z,T)
REAL(rp), INTENT(IN) :: R0
REAL(rp), INTENT(IN) :: b
REAL(rp), INTENT(IN) :: w
REAL(rp), INTENT(IN) :: dlam
REAL(rp), INTENT(IN) :: R
REAL(rp), INTENT(IN) :: Z
REAL(rp), INTENT(IN) :: T
Real(rp) :: rm
Real(rp) :: lam
REAL(rp) :: Spong_2D
rm=sqrt((R-R0)**2+Z**2)
lam=(sin(deg2rad(T)))**2
Spong_2D=(1-tanh((rm-b)/w))/(1-tanh(-b/w))*exp(-(lam/dlam)**2)
end function Spong_2D