diff --git a/src/OCE/SBC/sbcmod.F90 b/src/OCE/SBC/sbcmod.F90 index 30ee55c7c0dbb66d89024718ff621af8b2c4f73d..8881e2f4aff9d366301bb136bc3d40ea2e2c8d5f 100644 --- a/src/OCE/SBC/sbcmod.F90 +++ b/src/OCE/SBC/sbcmod.F90 @@ -444,8 +444,8 @@ CONTAINS ! IF( ln_wave .AND. ln_tauoc ) THEN ! Wave stress reduction DO_2D( 0, 0, 0, 0) - utau(ji,jj) = utau(ji,jj) * ( tauoc_wave(ji,jj) + tauoc_wave(ji-1,jj) ) * 0.5_wp - vtau(ji,jj) = vtau(ji,jj) * ( tauoc_wave(ji,jj) + tauoc_wave(ji,jj-1) ) * 0.5_wp + utau(ji,jj) = utau(ji,jj) * ( tauoc_wave(ji,jj) + tauoc_wave(ji+1,jj) ) * 0.5_wp + vtau(ji,jj) = vtau(ji,jj) * ( tauoc_wave(ji,jj) + tauoc_wave(ji,jj+1) ) * 0.5_wp END_2D ! CALL lbc_lnk( 'sbcwave', utau, 'U', -1. )