Changeset 5014 for branches/2015/dev_r4826_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

Timestamp:

2015-01-07T19:03:53+01:00 (6 years ago)

Author:

hliu

Message:

upload the modifications for W/D based on r:4826

File:

• branches/2015/dev_r4826_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90 (modified) (3 diffs)

branches/2015/dev_r4826_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

@@ -145 +145 @@
       fse3t_a(:,:,jpk) = e3t_0(:,:,jpk)
 
+      IF(ln_wd) THEN
+        DO jj = 1, jpj
+          DO ji = 1, jpi
+            IF(e3t_0(ji,jj,1) < 0._wp) fse3t_a(ji,jj,1:2) = 0.5_wp * rn_wdmin1 
+          ENDDO
+        ENDDO
+      END IF
+
       ! Reconstruction of all vertical scale factors at now and before time steps
       ! =============================================================================
@@ -793 +801 @@
       CHARACTER(len=*), INTENT(in) ::   cdrw   ! "READ"/"WRITE" flag
       !! * Local declarations
-      INTEGER ::   jk
+      INTEGER ::   ji, jj, jk
       INTEGER ::   id1, id2, id3, id4, id5     ! local integers
       !!----------------------------------------------------------------------
@@ -867 +875 @@
             fse3t_n(:,:,:) = e3t_0(:,:,:)
             sshn(:,:) = 0.0_wp
+
+            IF(ln_wd) THEN
+              DO jj = 1, jpj
+                DO ji = 1, jpi
+                  IF(e3t_0(ji,jj,1) < 0._wp) THEN
+                    fse3t_b(ji,jj,1:2) = 0.5_wp * rn_wdmin1 
+                    fse3t_n(ji,jj,1:2) = 0.5_wp * rn_wdmin1 
+                    sshb(ji,jj) = rn_wdmin1 - bathy(ji,jj)
+                    sshn(ji,jj) = rn_wdmin1 - bathy(ji,jj)
+                  ENDIF
+                ENDDO
+              ENDDO
+            END IF
+
             IF( ln_vvl_ztilde .OR. ln_vvl_layer) THEN
                tilde_e3t_b(:,:,:) = 0.0_wp