Changeset 2388 for trunk

Timestamp:

2010-11-15T16:24:52+01:00 (13 years ago)

Author:

smasson

Message:

ticket #757, LIM-2 - bulk CLIO bug in surface ice stress

File:

• trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

@@ -509 +509 @@
          END DO
       CASE( 'B' )                          ! B-grid ice dynamics
+         ! Change from wind speed to wind stress over OCEAN (cao is used)
+         zcoef = 0.5 * cai / cao 
          ! stress from ocean U- and V-points to ice U,V point
 !CDIR COLLAPSE
          DO jj = 2, jpj
             DO ji = 2, jpi   ! B grid : no vector opt.
-               p_taui(ji,jj) = 0.5 * ( utau(ji-1,jj  ) + utau(ji-1,jj-1) )
-               p_tauj(ji,jj) = 0.5 * ( vtau(ji  ,jj-1) + vtau(ji-1,jj-1) )
+               p_taui(ji,jj) = zcoef * ( utau(ji-1,jj  ) + utau(ji-1,jj-1) )
+               p_tauj(ji,jj) = zcoef * ( vtau(ji  ,jj-1) + vtau(ji-1,jj-1) )
             END DO
          END DO