Changeset 11002 for NEMO/branches/2019/fix_vvl_ticket1791/src/OCE/DOM/domain.F90

Timestamp:

2019-05-20T14:53:58+02:00 (5 years ago)

Author:

jchanut

Message:

vvl interp, #1791: 1) Add resting depth at F-points. 2) Rewrite dom_vvl_interpol with 2 aditionnal methods for horizontal interpolation at U-V-F points (nmet=0,1,2 in subroutine header). Set nmet=1, i.e. interpolate interfaces instead of scale factors anomaly. 3) Changes in the rest of the code: consider input array in dom_vvl_interpol at T point for F scale factors interpolation.

File:

• NEMO/branches/2019/fix_vvl_ticket1791/src/OCE/DOM/domain.F90 (modified) (2 diffs)

NEMO/branches/2019/fix_vvl_ticket1791/src/OCE/DOM/domain.F90

@@ -150 +150 @@
       hu_0(:,:) = 0._wp
       hv_0(:,:) = 0._wp
+      hf_0(:,:) = 0._wp
       DO jk = 1, jpk
          ht_0(:,:) = ht_0(:,:) + e3t_0(:,:,jk) * tmask(:,:,jk)
@@ -156 +157 @@
       END DO
       !
+      DO jj = 1, jpjm1
+         DO jk = 1, jpk
+            hf_0(:,jj) = hf_0(:,jj) + e3f_0(:,jj,jk) * umask(:,jj,jk) * umask(:,jj+1,jk)
+         END DO
+      END DO
+      CALL lbc_lnk( 'dom_init' , hf_0, 'F', 1._wp )
       !           !==  time varying part of coordinate system  ==!
       !