IF ( ABS(Rj)*thetaMax .LE. ABS(Rjp) ) THEN
thetaM=SIGN(thetaMax,Rjp*Rj)
_RL thetaMax
PARAMETER( thetaMax = 1.D+20 )
IF ( ABS(Rj)*thetaMax .LE. ABS(Rjm) ) THEN
thetaP=SIGN(thetaMax,Rjm*Rj)