Changeset 1708 for trunk/NEMO/OPA_SRC/TRD/trdicp.F90

Timestamp:

2009-11-04T14:24:34+01:00 (15 years ago)

Author:

rblod

Message:

Stability for explicit bottom friction, see ticket #588

File:

• trunk/NEMO/OPA_SRC/TRD/trdicp.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/TRD/trdicp.F90

@@ -117 +117 @@
                   umo(ktrd) = umo(ktrd) + ptrd2dx(ji,jj) * e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,1)
                   vmo(ktrd) = vmo(ktrd) + ptrd2dy(ji,jj) * e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,1)
-               END DO
-            END DO
-            !
-         CASE( jpdyn_trd_bfr )         ! bottom friction fluxes
-            DO jj = 1, jpj
-               DO ji = 1, jpi
-                  umo(ktrd) = umo(ktrd) + ptrd2dx(ji,jj)
-                  vmo(ktrd) = vmo(ktrd) + ptrd2dy(ji,jj)
                END DO
             END DO
@@ -520 +512 @@
             WRITE (numout,9530) hke(jpicpd_swf) / tvolt
             WRITE (numout,9531) hke(jpicpd_dat) / tvolt
-            WRITE (numout,9532)
-            WRITE (numout,9533)   &
+            WRITE (numout,9532) hke(jpicpd_bfr) / tvolt
+            WRITE (numout,9533)
+            WRITE (numout,9534)   &
             &     (  hke(jpicpd_hpg) + hke(jpicpd_keg) + hke(jpicpd_rvo) + hke(jpicpd_pvo) + hke(jpicpd_ldf)   &
             &      + hke(jpicpd_had) + hke(jpicpd_zad) + hke(jpicpd_zdf) + hke(jpicpd_spg) + hke(jpicpd_dat)   &
-            &      + hke(jpicpd_swf) ) / tvolt
+            &      + hke(jpicpd_swf) + hke(jpicpd_bfr) ) / tvolt
          ENDIF
 
@@ -539 +532 @@
  9530    FORMAT(' surface wind forcing      u2= ', e20.13)
  9531    FORMAT(' dampimg term              u2= ', e20.13)
- 9532    FORMAT(' --------------------------------------------------')
- 9533    FORMAT(' total trend               u2= ', e20.13)
+ 9532    FORMAT(' bottom flux               u2= ', e20.13)
+ 9533    FORMAT(' --------------------------------------------------')
+ 9534    FORMAT(' total trend               u2= ', e20.13)
 
          IF(lwp) THEN