Changeset 2938 for branches/2011/dev_r2802_NOCL_Smagorinsky

Timestamp:

2011-10-16T10:02:39+02:00 (13 years ago)

Author:

hliu

Message:

change the order of multiply diffusivity, smagorinsky+bilaplacian not recommended

File:

• branches/2011/dev_r2802_NOCL_Smagorinsky/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90 (modified) (2 diffs)

branches/2011/dev_r2802_NOCL_Smagorinsky/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90

@@ -131 +131 @@
                DO ji = fs_2, fs_jpim1   ! vector opt.
                   zbtr = 1.0 / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )
-                  zlt(ji,jj) = fsahtt(ji,jj,jk) * zbtr * (   ztu(ji,jj,jk) - ztu(ji-1,jj,jk)   &
+                  zlt(ji,jj) = ((1-nkaht_smag)*fsahtt(ji,jj,jk)+nkaht_smag) * zbtr *           &
+                     &                                   (   ztu(ji,jj,jk) - ztu(ji-1,jj,jk)   &
                      &                                     + ztv(ji,jj,jk) - ztv(ji,jj-1,jk)   )
                END DO
             END DO
+
+
+
             CALL lbc_lnk( zlt, 'T', 1. )     ! Lateral boundary conditions (unchanged sgn)
 
@@ -151 +155 @@
                   ztra = zbtr * (  ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk)  )
                   ! add it to the general tracer trends
-                  pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + ztra
+                  pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + ztra*(nkaht_smag*fsahtt(ji,jj,jk)+1-nkaht_smag)
                END DO
             END DO