Changeset 9353 for branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

Timestamp:

2018-02-23T16:52:00+01:00 (6 years ago)

Author:

jchanut

Message:

Tilde coordinate update: first draft

File:

• branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90 (modified) (3 diffs)

branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

@@ -70 +70 @@
       !!              - 1D configuration, move Coriolis, u and v at T-point
       !!----------------------------------------------------------------------
-      INTEGER ::   jk          ! dummy loop argument
+      INTEGER ::   jj, jk      ! dummy loop argument
       INTEGER ::   iconf = 0   ! local integers
       !!----------------------------------------------------------------------
@@ -92 +92 @@
       hu_0(:,:) = 0.0_wp                       ! Reference ocean depth at U-points
       hv_0(:,:) = 0.0_wp                       ! Reference ocean depth at V-points
+      hf_0(:,:) = 0.0_wp                       ! Reference ocean depth at F-points
       DO jk = 1, jpk
          ht_0(:,:) = ht_0(:,:) + e3t_0(:,:,jk) * tmask(:,:,jk)
@@ -97 +98 @@
          hv_0(:,:) = hv_0(:,:) + e3v_0(:,:,jk) * vmask(:,:,jk)
       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( hf_0, 'F', 1._wp )
       !
       IF( lk_c1d )           CALL cor_c1d      ! 1D configuration: Coriolis set at T-point