Changeset 11201 for NEMO/branches/2019/ENHANCE-03_domcfg/src/dommsk.F90

Timestamp:

2019-07-01T12:10:15+02:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-03_domcfg : add management of closed seas in domain cfg by flood filling and lat/lon seed instead of i/j box definition (ticket #2143)

File:

• NEMO/branches/2019/ENHANCE-03_domcfg/src/dommsk.F90 (modified) (6 diffs)

NEMO/branches/2019/ENHANCE-03_domcfg/src/dommsk.F90

@@ -24 +24 @@
    !!----------------------------------------------------------------------
    USE dom_oce        ! ocean space and time domain
+   USE domwri         ! domain: write the meshmask file
    USE bdy_oce        ! open boundary
    !
@@ -60 +61 @@
       !!      and ko_bot, the indices of the fist and last ocean t-levels which 
       !!      are either defined in usrdef_zgr or read in zgr_read.
-      !!                The velocity masks (umask, vmask, wmask, wumask, wvmask) 
+      !!                The velocity masks (umask, vmask) 
       !!      are deduced from a product of the two neighboring tmask.
       !!                The vorticity mask (fmask) is deduced from tmask taking
@@ -75 +76 @@
       !!                due to cyclic or North Fold boundaries as well as MPP halos.
       !!
-      !! ** Action :   tmask, umask, vmask, wmask, wumask, wvmask : land/ocean mask 
+      !! ** Action :   tmask, umask, vmask, wmask : land/ocean mask 
       !!                         at t-, u-, v- w, wu-, and wv-points (=0. or 1.)
       !!               fmask   : land/ocean mask at f-point (=0., or =1., or 
@@ -188 +189 @@
       !-----------------------------------------
       wmask (:,:,1) = tmask(:,:,1)     ! surface
-      wumask(:,:,1) = umask(:,:,1)
-      wvmask(:,:,1) = vmask(:,:,1)
       DO jk = 2, jpk                   ! interior values
          wmask (:,:,jk) = tmask(:,:,jk) * tmask(:,:,jk-1)
-         wumask(:,:,jk) = umask(:,:,jk) * umask(:,:,jk-1)   
-         wvmask(:,:,jk) = vmask(:,:,jk) * vmask(:,:,jk-1)
       END DO
 
@@ -202 +199 @@
       ssumask(:,:) = MAXVAL( umask(:,:,:), DIM=3 )
       ssvmask(:,:) = MAXVAL( vmask(:,:,:), DIM=3 )
-
 
       ! Interior domain mask  (used for global sum)
@@ -289 +285 @@
       ENDIF
       !
+      ! write out mesh mask
+      IF ( nn_msh > 0 ) CALL dom_wri
+      !
    END SUBROUTINE dom_msk