Changeset 11180 for NEMO/branches/UKMO/NEMO_4.0_surge/src/OCE/DYN/dynkeg.F90

Timestamp:

2019-06-25T18:50:27+02:00 (5 years ago)

Author:

clne

Message:

Initial commit of code for 2d (surge) work in NEMO4.
This is aiming to replicate the 3.6 version in branches/UKMO/dev_r5518_Surge_Modelling

File:

• NEMO/branches/UKMO/NEMO_4.0_surge/src/OCE/DYN/dynkeg.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0_surge/src/OCE/DYN/dynkeg.F90

@@ -191 +191 @@
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
-               ua(ji,jj,jk) = ua(ji,jj,jk) - ( zhke(ji+1,jj  ,jk) - zhke(ji,jj,jk) ) / e1u(ji,jj)
-               va(ji,jj,jk) = va(ji,jj,jk) - ( zhke(ji  ,jj+1,jk) - zhke(ji,jj,jk) ) / e2v(ji,jj)
+               ua(ji,jj,jk) = ua(ji,jj,jk) - ( zhke(ji+1,jj  ,jk) - zhke(ji,jj,jk) ) / e1u(ji,jj)* umask(ji,jj,jk)
+               va(ji,jj,jk) = va(ji,jj,jk) - ( zhke(ji  ,jj+1,jk) - zhke(ji,jj,jk) ) / e2v(ji,jj)* vmask(ji,jj,jk)
             END DO 
          END DO