Changeset 1569

Timestamp:

2009-07-31T18:37:36+02:00 (15 years ago)

Author:

ctlod

Message:

correction of a bug in the Coriolis term in the rheology of LIM2, see ticket: #504

File:

• trunk/NEMO/LIM_SRC_2/limrhg_2.F90 (modified) (1 diff)

trunk/NEMO/LIM_SRC_2/limrhg_2.F90

@@ -173 +173 @@
             zmass(ji,jj) = ( zmasst(ji,jj  ) * wght(ji,jj,2,2) + zmasst(ji-1,jj  ) * wght(ji,jj,1,2)   &
                &           + zmasst(ji,jj-1) * wght(ji,jj,2,1) + zmasst(ji-1,jj-1) * wght(ji,jj,1,1) ) * zusw
-            zcorl(ji,jj) = zmass(ji,jj) * fcor(ji,jj)
+            zcorl(ji,jj) = zmass(ji,jj) &
+               &           *( fcor(ji,jj  ) * wght(ji,jj,2,2) + fcor(ji-1,jj  )*wght(ji,jj,1,2)   &
+               &            + fcor(ji,jj-1) * wght(ji,jj,2,1) + fcor(ji-1,jj-1)*wght(ji,jj,1,1) ) * zusw
 
             ! Wind stress.