Rnondim = -depthZ/(LrhoW(i,j)*Smod)
Rnondim = 1.
IF ( Rnondim.GE.1. _d 0 ) THEN
f2 = op5*( 1. _d 0 + SIN( fpi*(Rnondim-op5) ))
Rnondim = -depthZ/(LrhoS(i,j)*Smod)
_RL f1,Smod,f2,Rnondim
Rnondim= -depthZ/(Lrho(i,j)*Smod)