Changeset 4375 for branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2014-01-28T14:55:35+01:00 (10 years ago)

Author:

hliu

Message:

updated gravity filters

File:

• branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (1 diff)

branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -1411 +1411 @@
 !!
       ! HYBRID : 
-      DO jj = 1, jpj
-         DO ji = 1, jpi
-            DO jk = 1, jpkm1
-               IF( scobot(ji,jj) >= fsdept(ji,jj,jk) )   mbathy(ji,jj) = MAX( 2, jk )
-               IF( scobot(ji,jj) == 0._wp            )   mbathy(ji,jj) = 0
-            END DO
-         END DO
-      END DO
+      IF(ln_wad) THEN
+      !! Wetting/Drying case
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               DO jk = 1, jpkm1
+                  IF( scobot(ji,jj) >= fsdept(ji,jj,jk) ) THEN
+                    mbathy(ji,jj) = MAX( 2, jk )
+                  ELSE IF(scobot(ji,jj) <= rn_landele) THEN
+                    mbathy(ji,jj) = 0
+                  ELSE
+                    mbathy(ji,jj) = 2
+                  END IF
+               END DO
+            END DO
+         END DO
+      ELSE
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               DO jk = 1, jpkm1
+                  IF( scobot(ji,jj) >= fsdept(ji,jj,jk) )   mbathy(ji,jj) = MAX( 2, jk )
+                  IF( scobot(ji,jj) == 0._wp            )   mbathy(ji,jj) = 0
+               END DO
+            END DO
+         END DO
+      END IF
       IF( nprint == 1 .AND. lwp ) WRITE(numout,*) ' MIN val mbathy h90 ', MINVAL( mbathy(:,:) ),   &
          &                                                       ' MAX ', MAXVAL( mbathy(:,:) )