IF ( ABS(Rj).LT.smallNo .OR.
& ABS(Rjp).LT.smallNo ) THEN
thetaM=(Rjp+smallNo)/(smallNo+Rj)
_RL smallNo
smallNo = 1.0D-20
& ABS(Rjm).LT.smallNo ) THEN
thetaP=(Rjm+smallNo)/(smallNo+Rj)