Changeset 2309

Timestamp:

2010-10-25T16:42:27+02:00 (13 years ago)

Author:

cbricaud

Message:

save bottom stress for next time step for zdfgls with explicit vertical diffusion scheme

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90 (modified) (2 diffs)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90

@@ -46 +46 @@
 
    !! * Substitutions
+#  include "vectopt_loop_substitute.h90"
 #  include "domzgr_substitute.h90"
    !!----------------------------------------------------------------------
@@ -119 +120 @@
       ENDIF
       !
+
+      !zdfgls case with explicit vertical diffusion scheme: save bottom stress for next time step
+
+#if defined key_zdfgls && defined key_dynspg_exp
+      DO jj = 2, jpjm1
+         DO ji = fs_2, fs_jpim1   ! vector opt.
+            ! 
+            ikbu   = MIN( mbathy(ji+1,jj  ), mbathy(ji,jj) )
+            ikbum1 = MAX( ikbu-1, 1 )
+            wbotu(ji,jj) = bfrua(ji,jj) * un(ji,jj,ikbum1) * umask(ji,jj,ikbum1)
+            !
+            ikbv   = MIN( mbathy(ji,jj+1), mbathy(ji,jj) )
+            ikbvm1 = MAX( ikbv-1, 1 )
+            wbotv(ji,jj) = bfrva(ji,jj) * vn(ji,jj,ikbvm1) * vmask(ji,jj,ikbvm1)
+            !
+         END DO
+      END DO
+      CALL lbc_lnk( wbotu(:,:), 'U', -1. )
+      CALL lbc_lnk( wbotv(:,:), 'V', -1. )
+#endif
+
    END SUBROUTINE zdf_bfr