Changeset 1684

Timestamp:

2009-10-28T11:27:35+01:00 (15 years ago)

Author:

rblod

Message:

Correct a bug in LIM3 when averaging ocean surface velocity, see ticket #538

File:

• trunk/NEMO/LIM_SRC_3/limsbc.F90 (modified) (1 diff)

trunk/NEMO/LIM_SRC_3/limsbc.F90

@@ -204 +204 @@
                zat_v = ( at_i(ji,jj) + at_i(ji,jj+1) ) * 0.5 
                !
-               zu_ico = u_ice(ji,jj) - 0.5 * ( ub(ji,jj,1) - ssu_mb(ji,jj) )   ! ice-oce velocity using ub and ssu_mb
-               zv_ico = v_ice(ji,jj) - 0.5 * ( vb(ji,jj,1) - ssv_mb(ji,jj) )
+               zu_ico = u_ice(ji,jj) - 0.5 * ( ub(ji,jj,1) + ssu_mb(ji,jj) )   ! ice-oce velocity using ub and ssu_mb
+               zv_ico = v_ice(ji,jj) - 0.5 * ( vb(ji,jj,1) + ssv_mb(ji,jj) )
                !                                        ! quadratic drag formulation with rotation
 !!gm still an error in the rotation, but usually the angle is zero (zsang=0, cangvg=1)