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