Changeset 4405 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2014-02-04T13:11:58+01:00 (10 years ago)

Author:

trackstand2

Message:

Add option to switch trimming of sub-domains in z via NEMO_TRIM_DOMZ env var

File:

• branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (2 diffs)

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -772 +772 @@
       USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
       USE wrk_nemo, ONLY:   zmbk => wrk_2d_1
-      USE par_oce,  ONLY: jpkf, jpkfm1
+      USE par_oce,  ONLY:   jpkf, jpkfm1
       !!
       INTEGER ::   ji, jj   ! dummy loop indices
+      LOGICAL :: domtrim_z
+      CHARACTER(LEN=8) :: lstr            ! Local string for reading env. var.
+      INTEGER          :: lztrim, ierr    ! Local int for      "      "    "
       !!----------------------------------------------------------------------
       !
@@ -799 +802 @@
       zmbk(:,:) = REAL( mbkv(:,:), wp )   ;   CALL lbc_lnk(zmbk,'V',1.)   ;   mbkv  (:,:) = MAX( INT( zmbk(:,:) ), 1 )
       !
-      ! Compute and store the deepest bottom level of any grid-type at each grid point
-      ! For use in removing data below ocean floor from compute loops
-      mbkmax(:,:) = MAX(mbkt(:,:)+1, mbku(:,:), mbkv(:,:))
-      jpkf = MAXVAL( mbkmax(:,:) )
+
+      domtrim_z = .TRUE.
+      CALL GET_ENVIRONMENT_VARIABLE(NAME='NEMO_TRIM_DOMZ', VALUE=lstr, &
+                                    STATUS=ierr)
+      IF( ierr == 0)THEN
+         READ(lstr,FMT="(I)",IOSTAT=ierr) lztrim
+         IF(ierr == 0)THEN
+            IF (lztrim == 0) domtrim_z = .FALSE.
+         ELSE
+            CALL ctl_warn('zgr_bot_level: failed to parse value of NEMO_TRIM_DOMZ environment variable: '//TRIM(lstr))
+         END IF
+      END IF
+
+      ! Compute and store the deepest bottom level of any grid-type at each grid
+      ! point. For use in removing data below ocean floor from compute loops
+      IF( domtrim_z ) THEN
+         mbkmax(:,:) = MAX(mbkt(:,:)+1, mbku(:,:), mbkv(:,:))
+!         jpkf = MIN(jpk, 1 + MAXVAL( mbkmax(:,:) ) )
+         jpkf = MAXVAL( mbkmax(:,:) )
+         WRITE(*,*) narea,': ARPDBG: shallowest pt and jpkf = ', &
+                    MINVAL(mbkmax(:,:)), jpkf
+      ELSE
+         WRITE(*,*) narea,': ARPDBG: NOT trimming domain in z'
+         jpkf = jpk
+      END IF
       jpkfm1 = jpkf - 1
-      WRITE(*,*) narea,': ARPDBG: shallowest pt and jpkf = ', MINVAL(mbkmax(:,:)), jpkf
 
       !