Changeset 1244 for trunk/NEMO/LIM_SRC_3/limrhg.F90

Timestamp:

2009-01-06T12:13:44+01:00 (15 years ago)

Author:

ctlod

Message:

improve the LIM 3.0 stability in using sea surface height, see ticket: #286

File:

• trunk/NEMO/LIM_SRC_3/limrhg.F90 (modified) (4 diffs)

trunk/NEMO/LIM_SRC_3/limrhg.F90

@@ -4 +4 @@
    !!   Ice rheology : sea ice rheology
    !!======================================================================
+   !! History :   -   !  2007-03  (M.A. Morales Maqueda, S. Bouillon) Original code
+   !!            3.0  !  2008-03  (M. Vancoppenolle) LIM3
+   !!             -   !  2008-11  (M. Vancoppenolle, S. Bouillon, Y. Aksenov) add surface tilt in ice rheolohy 
+   !!----------------------------------------------------------------------
 #if defined key_lim3
    !!----------------------------------------------------------------------
@@ -16 +20 @@
    USE dom_oce
    USE dom_ice
+   USE sbc_oce         ! Surface boundary condition: ocean fields
    USE sbc_ice         ! Surface boundary condition: ice fields
    USE ice
@@ -100 +105 @@
       !!                 Vancoppenolle et al., OM08
       !!
-      !! History :
-      !!   1.0  !  07-03  (M.A. Morales Maqueda, S. Bouillon)
-      !!   2.0  !  08-03  M. Vancoppenolle : LIM3
-      !!
       !!-------------------------------------------------------------------
       ! * Arguments
@@ -290 +291 @@
             ! SB On utilise maintenant le gradient de la pente de l'ocean
             ! include it later
-            !    zdsshx =  (ssh_io(ji+1,jj) - ssh_io(ji,jj))/e1u(ji,jj)
-            !    zdsshy =  (ssh_io(ji,jj+1) - ssh_io(ji,jj))/e2v(ji,jj)  
-
-            zdsshx = 0.0
-            zdsshy = 0.0 
+
+            zdsshx =  (ssh_m(ji+1,jj) - ssh_m(ji,jj))/e1u(ji,jj)
+            zdsshy =  (ssh_m(ji,jj+1) - ssh_m(ji,jj))/e2v(ji,jj)
 
             za1ct(ji,jj) = ztagnx - zmass1(ji,jj) * grav * zdsshx