Changeset 1711 for trunk/NEMO/LIM_SRC_2

Timestamp:

2009-11-04T17:01:10+01:00 (14 years ago)

Author:

rblod

Message:

Add tilt off sea-surface in LIM2 rheology, see ticket #587

File:

• trunk/NEMO/LIM_SRC_2/limrhg_2.F90 (modified) (3 diffs)

trunk/NEMO/LIM_SRC_2/limrhg_2.F90

@@ -18 +18 @@
    !!----------------------------------------------------------------------
    USE par_oce        ! ocean parameter
+   USE dom_oce        ! ocean space and time domain
+   USE sbc_oce        ! surface boundary condition: ocean variables
    USE sbc_ice        ! surface boundary condition: ice variables
    USE dom_ice_2      ! domaine: ice variables
@@ -63 +65 @@
       INTEGER ::   iter, jter          ! temporary integers
       CHARACTER (len=50) ::   charout
-      REAL(wp) ::   ze11  , ze12  , ze22  , ze21       ! temporary scalars
-      REAL(wp) ::   zt11  , zt12  , zt21  , zt22       !    "         "
-      REAL(wp) ::   zvis11, zvis21, zvis12, zvis22     !    "         "
-      REAL(wp) ::   zgphsx, ztagnx, zunw, zur, zusw    !    "         "
-      REAL(wp) ::   zgphsy, ztagny, zvnw, zvr          !    "         "
+      REAL(wp) ::   ze11  , ze12  , ze22  , ze21               ! temporary scalars
+      REAL(wp) ::   zt11  , zt12  , zt21  , zt22               !    "         "
+      REAL(wp) ::   zvis11, zvis21, zvis12, zvis22             !    "         "
+      REAL(wp) ::   zgphsx, ztagnx, zgsshx, zunw, zur, zusw    !    "         "
+      REAL(wp) ::   zgphsy, ztagny, zgsshy, zvnw, zvr          !    "         "
       REAL(wp) ::   zresm,  za, zac, zmod
       REAL(wp) ::   zmpzas, zstms, zindu, zusdtp, zmassdt, zcorlal
@@ -195 +197 @@
                &     + ( alambd(ji,jj,1,1,1,2) - alambd(ji,jj,1,2,1,2) ) * zpresh(ji-1,jj  )
 
+            ! Gradient of the sea surface height
+            zgsshx =  (   (ssh_m(ji  ,jj  ) - ssh_m(ji-1,jj  ))/e1u(ji-1,jj  )   &
+               &       +  (ssh_m(ji  ,jj-1) - ssh_m(ji-1,jj-1))/e1u(ji-1,jj-1)   ) * 0.5
+            zgsshy =  (   (ssh_m(ji  ,jj  ) - ssh_m(ji  ,jj-1))/e2v(ji  ,jj-1)   &
+               &       +  (ssh_m(ji-1,jj  ) - ssh_m(ji-1,jj-1))/e2v(ji-1,jj-1)   ) * 0.5
+
             ! Computation of the velocity field taking into account the ice-ice interaction.                                 
             ! Terms that are independent of the ice velocity field.
-            za1ct(ji,jj) = ztagnx - zcorl(ji,jj) * v_oce(ji,jj) - zgphsx
-            za2ct(ji,jj) = ztagny + zcorl(ji,jj) * u_oce(ji,jj) - zgphsy
+            za1ct(ji,jj) = ztagnx - zmass(ji,jj) * grav * zgsshx - zgphsx
+            za2ct(ji,jj) = ztagny - zmass(ji,jj) * grav * zgsshy - zgphsy
          END DO
       END DO