Changes from branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90 at r8627 to trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90 at r4990

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90 (modified) (3 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90

@@ -17 +17 @@
    USE oce             ! ocean dynamics and tracers
    USE dom_oce         ! ocean space and time domain
-   USE phycst          ! physical constants
    USE ldfdyn_oce      ! ocean dynamics: lateral physics
    USE zdf_oce         ! ocean vertical physics
    !
    USE in_out_manager  ! I/O manager
-   USE iom             ! I/O library
    USE timing          ! Timing
 
@@ -64 +62 @@
       INTEGER  ::   ji, jj, jk             ! dummy loop indices
       REAL(wp) ::   zua, zva, ze2u, ze1v   ! local scalars
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   z2d   ! 2D workspace 
       !!----------------------------------------------------------------------
       !
@@ -96 +93 @@
       END DO                                           !   End of slab
       !                                                ! ===============
-      
-      IF( iom_use('dispkexyfo') ) THEN   ! ocean Kinetic Energy dissipation per unit area
-         !                               ! due to lateral friction (xy=lateral) 
-         !   see NEMO_book appendix C, §C.7.2 (N.B. here averaged at t-points)
-         !   local dissipation of KE at t-point due to laplacian operator is given by :
-         !      - ahmt hdivb**2 - mi( mj(ahmf rotb**2 e1f*e2f*e3t) ) / (e1e2t*e3t)
-         !
-         ALLOCATE( z2d(jpi,jpj) )
-         z2d(:,:) = 0._wp
-         DO jk = 1, jpkm1
-            DO jj = 2, jpjm1
-               DO ji = 2, jpim1
-                  z2d(ji,jj) = z2d(ji,jj)          -   (                                                         &
-                     &   hdivb(ji,jj,jk)**2 * fsahmt(ji,jj,jk) * fse3t_n(ji,jj,jk)                               &
-                     & + 0.25_wp * (                                                                             &
-                     &   rotb (ji  ,jj  ,jk)**2 * fsahmf(ji  ,jj  ,jk) * e12f(ji  ,jj  ) * fse3f(ji  ,jj  ,jk)   &
-                     & + rotb (ji-1,jj  ,jk)**2 * fsahmf(ji-1,jj  ,jk) * e12f(ji-1,jj  ) * fse3f(ji-1,jj  ,jk)   &
-                     & + rotb (ji  ,jj-1,jk)**2 * fsahmf(ji  ,jj-1,jk) * e12f(ji  ,jj-1) * fse3f(ji  ,jj-1,jk)   &
-                     & + rotb (ji-1,jj-1,jk)**2 * fsahmf(ji-1,jj-1,jk) * e12f(ji-1,jj-1) * fse3f(ji-1,jj-1,jk)   &
-                     &             ) * r1_e12t(ji,jj)  ) * tmask(ji,jj,jk)
-               END DO
-            END DO
-         END DO
-         z2d(:,:) = rau0 * z2d(:,:) 
-         CALL lbc_lnk( z2d,'T', 1. )
-         CALL iom_put( 'dispkexyfo', z2d )
-         DEALLOCATE( z2d )
-      ENDIF
-
       IF( nn_timing == 1 )  CALL timing_stop('dyn_ldf_lap')
       !