_RL taux_onIce_cntr, tauy_onIce_cntr, uvel_cntr, vvel_cntr
tauy_onIce_cntr=HALF*
rhs_y= -tauy_onIce_cntr +mIceCor*uvel_cntr