_RL advFac, rAdvFac
advFac = 0. _d 0
IF (calcAdvection) advFac = 1. _d 0
rAdvFac = rkSign*advFac
& )*advFac
_RL advFac