Changeset 12871 for utils/tools_ticket2457/DOMAINcfg/src/domisf.F90

Timestamp:

2020-05-05T17:57:08+02:00 (9 months ago)

Author:

mathiot

Message:

ticket #2457: add fix suggested in ticket.

File:

• utils/tools_ticket2457/DOMAINcfg/src/domisf.F90 (modified) (2 diffs)

utils/tools_ticket2457/DOMAINcfg/src/domisf.F90

@@ -109 +109 @@
       INTEGER  :: ji, jj                             ! loop indexes
       INTEGER  :: imskjp1, imskjm1, imskip1, imskim1 ! local variable
+      INTEGER  :: ibtest, icompt                     ! local variables
       INTEGER, DIMENSION(jpi,jpj) :: imask           ! isf mask
       !
+      REAL(wp) ::   zdummy       ! dummy variable used for the lnclnk
       REAL(wp) ::   zisfdep_min  ! minimal ice shelf draft allowed
       !!---------------------------------------------------------------------
+      !
+      ! 0.0 fill isolated grid point in the bathymetry
+      ! will be done again later on in zgr_bat_ctl, but need to be done here to adjust risfdep/misfdep respectively
+      icompt = 0
+      DO jj = 2, jpjm1
+         DO ji = 2, jpim1
+            ibtest = MAX(  mbathy(ji-1,jj), mbathy(ji+1,jj),   &
+               &           mbathy(ji,jj-1), mbathy(ji,jj+1)  )
+               IF( ibtest < mbathy(ji,jj) ) THEN
+                  mbathy(ji,jj) = ibtest
+                  bathy(ji,jj)  = gdepw_1d(ibtest+1)
+                  icompt = icompt + 1
+               END IF
+         END DO
+      END DO
+      !
+      ! ensure halo correct 
+      zdummy(:,:) = FLOAT( mbathy(:,:) ) ; CALL lbc_lnk('domisf', zdummy, 'T', 1._wp ) ; mbathy(:,:) = INT( zdummy(:,:) )
+      !
+      IF( lk_mpp )   CALL mpp_sum('domisf', icompt )
+      IF( icompt == 0 ) THEN
+         IF(lwp) WRITE(numout,*)'     no isolated ocean grid points'
+      ELSE
+         IF(lwp) WRITE(numout,*)'    ',icompt,' ocean grid points suppressed'
+      ENDIF
       !
       ! 1.0 set iceshelf to the minimum depth allowed
@@ -220 +247 @@
          WHERE( risfdep(:,:) > 0._wp .AND. risfdep(:,:) >= zdepth )   misfdep(:,:) = jk+1 
       END DO 
-      !
-      ! 2.1 fill isolated grid point in the bathymetry
-      ! will be done again later on in zgr_bat_ctl, but need to be done here to adjust misfdep respectively
-      icompt = 0
-      DO jj = 2, jpjm1
-         DO ji = 2, jpim1
-            ibtest = MAX(  mbathy(ji-1,jj), mbathy(ji+1,jj),   &
-               &           mbathy(ji,jj-1), mbathy(ji,jj+1)  )
-               IF( ibtest < mbathy(ji,jj) ) THEN
-                  mbathy(ji,jj) = ibtest
-                  icompt = icompt + 1
-               END IF
-         END DO
-      END DO
-      !
-      ! ensure halo correct 
-      zdummy(:,:) = FLOAT( mbathy(:,:) ) ; CALL lbc_lnk('domisf', zdummy, 'T', 1._wp ) ; mbathy(:,:) = INT( zdummy(:,:) )
-      !
-      IF( lk_mpp )   CALL mpp_sum('domisf', icompt )
-      IF( icompt == 0 ) THEN
-         IF(lwp) WRITE(numout,*)'     no isolated ocean grid points'
-      ELSE
-         IF(lwp) WRITE(numout,*)'    ',icompt,' ocean grid points suppressed'
-      ENDIF
       !
       ! 2.2 be sure misfdep not below mbathy