Changeset 10272 for NEMO/trunk/src/OCE/DYN/dynspg_ts.F90

Timestamp:

2018-11-02T12:15:39+01:00 (6 years ago)

Author:

deazer

Message:

Bug fix for wd case with implicit bottom friction.
(Currently only done for ln_dynadv_vec)
Trivial, just move lines so that the are done in all advection cases

see Ticket:

https://forge.ipsl.jussieu.fr/nemo/ticket/2149

File:

• NEMO/trunk/src/OCE/DYN/dynspg_ts.F90 (modified) (2 diffs)

NEMO/trunk/src/OCE/DYN/dynspg_ts.F90

@@ -1114 +1114 @@
                             &   ) * ssvmask(ji,jj)
  
-!jth implicit bottom friction:
-                  IF ( ll_wd ) THEN ! revert to explicit for bit comparison tests in non wad runs
-                     ua_e(ji,jj) =  ua_e(ji,jj) /(1.0 -   rdtbt * zCdU_u(ji,jj) * hur_e(ji,jj))
-                     va_e(ji,jj) =  va_e(ji,jj) /(1.0 -   rdtbt * zCdU_v(ji,jj) * hvr_e(ji,jj))
-                  ENDIF
-
                END DO
             END DO
@@ -1144 +1138 @@
                             &               +    hv_n(ji,jj)  * zv_frc(ji,jj) ) &
                             &   ) * zhvra
+               END DO
+            END DO
+         ENDIF
+!jth implicit bottom friction:
+         IF ( ll_wd ) THEN ! revert to explicit for bit comparison tests in non wad runs
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                     ua_e(ji,jj) =  ua_e(ji,jj) /(1.0 -   rdtbt * zCdU_u(ji,jj) * hur_e(ji,jj))
+                     va_e(ji,jj) =  va_e(ji,jj) /(1.0 -   rdtbt * zCdU_v(ji,jj) * hvr_e(ji,jj))
                END DO
             END DO