Changeset 4045 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90

Timestamp:

2013-09-25T16:38:24+02:00 (11 years ago)

Author:

clem

Message:

new LIM3

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90 (modified) (2 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90

@@ -35 +35 @@
 #  include "vectopt_loop_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/LIM3 3.4 , UCL - NEMO Consortium (2010)
+   !! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2010)
    !! $Id$
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -137 +137 @@
       END DO                                       ! end of sub-time step loop
 
+      ! -----------------------
+      !!! final step (clem) !!!
+      DO jj = 1, jpjm1                                ! diffusive fluxes in U- and V- direction
+         DO ji = 1 , fs_jpim1   ! vector opt.
+            zflu(ji,jj) = pahu(ji,jj) * e2u(ji,jj) / e1u(ji,jj) * ( ptab(ji+1,jj) - ptab(ji,jj) )
+            zflv(ji,jj) = pahv(ji,jj) * e1v(ji,jj) / e2v(ji,jj) * ( ptab(ji,jj+1) - ptab(ji,jj) )
+         END DO
+      END DO
+      !
+      DO jj= 2, jpjm1                                 ! diffusive trend : divergence of the fluxes
+         DO ji = fs_2 , fs_jpim1   ! vector opt. 
+            zdiv (ji,jj) = (  zflu(ji,jj) - zflu(ji-1,jj  )   &
+                 &            + zflv(ji,jj) - zflv(ji  ,jj-1)  ) / ( e1t (ji,jj) * e2t (ji,jj) )
+            ptab(ji,jj) = ztab0(ji,jj) + 0.5 * ( zdiv(ji,jj) + zdiv0(ji,jj) )
+         END DO
+      END DO
+      CALL lbc_lnk( ptab, 'T', 1. )                   ! lateral boundary condition
+      !!! final step (clem) !!!
+      ! -----------------------
+
       IF(ln_ctl)   THEN
          zrlx(:,:) = ptab(:,:) - ztab0(:,:)