CwatC(I,J,bi,bj)=0.25 _d 0
CwatC(I,J,bi,bj)=SEAICE_waterDrag_south*SQRT(TEMPVAR)
CwatC(I,J,bi,bj)=0.
CwatC(I,J,bi,bj)=SEAICE_waterDrag*SQRT(TEMPVAR)
CwatC(I,J,bi,bj) = CwatC(I,J,bi,bj) * maskC(I,J,kSrf,bi,bj)
O CwatC,
_RL CwatC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)