IF ( ABS(y_in(j+1)).LT.yPole .AND. ABS(y_in(j)).GT.yPole ) THEN
y_in(j) = -yPole
y_in(j-1) = -2.*yPole - y_in(1)
IF ( ABS(y_in(j+1)).GT.yPole .AND. ABS(y_in(j)).GT.yPole ) THEN
IF ( ABS(y_in(j-1)).LT.yPole .AND. ABS(y_in(j)).GT.yPole ) THEN
y_in(j) = yPole
y_in(j+1) = 2.*yPole - y_in(j-1)
IF ( poleSymmetry .AND. ABS(y_in(j)).GT.yPole ) THEN
IF ( ABS(y_in(j)).EQ.yPole ) THEN
_RL yPole, symSign, poleValue
PARAMETER ( yPole = 90. )
_RL yPole, symSign, poleU, poleV