Changeset 917 for trunk/NEMO/LIM_SRC_3/limsbc.F90

Timestamp:

2008-05-07T16:50:27+02:00 (16 years ago)

Author:

rblod

Message:

Correct a bug in revision 913 (ice/ocean stress computation)

File:

• trunk/NEMO/LIM_SRC_3/limsbc.F90 (modified) (2 diffs)

trunk/NEMO/LIM_SRC_3/limsbc.F90

@@ -82 +82 @@
       REAL(wp) ::   zfrldu, zfrldv   ! lead fraction at U- & V-points
       REAL(wp) ::   zu_io , zv_io    ! 2 components of the ice-ocean velocity
+      REAL(wp) ::   ztglx , ztgly    ! temporary scalar
       
 #if defined key_coupled    
@@ -279 +280 @@
                zmod  = SQRT( zu_io * zu_io + zv_io * zv_io ) 
                ! quadratic drag formulation
-               ztio_u(ji,jj) = rhoco * zmod * ( cangvg * zu_io - zsang * zv_io ) 
-               ztio_v(ji,jj) = rhoco * zmod * ( cangvg * zv_io + zsang * zu_io ) 
+               ztglx = rhoco * zmod * ( cangvg * zu_io - zsang * zv_io )
+               ztgly = rhoco * zmod * ( cangvg * zv_io + zsang * zu_io )
+               ! IMPORTANT
+               ! these lines are bound to prevent numerical oscillations
+               ! in the ice-ocean stress
+               ! They are physically ill-based. There is a cleaner solution
+               ! to try (remember discussion in Paris Gurvan)
+               ztio_u(ji,jj) = ztglx * exp( - zmod / 0.5 ) 
+               ztio_v(ji,jj) = ztgly * exp( - zmod / 0.5 ) 
                !
             END DO