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

Timestamp:

2014-02-04T13:12:20+01:00 (10 years ago)

Author:

trackstand2

Message:

Changes to allow jpk to be modified to deepest level within a subdomain. jpkorig holds original value.

File:

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

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

@@ -807 +807 @@
                                     STATUS=ierr)
       IF( ierr == 0)THEN
-         READ(lstr,FMT="(I)",IOSTAT=ierr) lztrim
+         READ(lstr,FMT="(I10)",IOSTAT=ierr) lztrim
          IF(ierr == 0)THEN
             IF (lztrim == 0) domtrim_z = .FALSE.
@@ -817 +817 @@
       ! 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
+      jpkorig = jpk
       IF( domtrim_z ) THEN
          mbkmax(:,:) = MAX(mbkt(:,:)+1, mbku(:,:), mbkv(:,:))
@@ -823 +824 @@
          WRITE(*,*) narea,': ARPDBG: shallowest pt and jpkf = ', &
                     MINVAL(mbkmax(:,:)), jpkf
+         ! Play rather fast and loose and just change the value of jpk
+         ! here...
+         jpk = jpkf
+         jpkm1 = jpk - 1
       ELSE
          WRITE(*,*) narea,': ARPDBG: NOT trimming domain in z'