Changeset 12623 for NEMO/branches/2020/r12581_ticket2418/src/OCE/IOM/iom.F90

Timestamp:

2020-03-28T08:38:26+01:00 (4 years ago)

Author:

smasson

Message:

r12581_ticket2418: merge with trunk@12622, see #2418

File:

• NEMO/branches/2020/r12581_ticket2418/src/OCE/IOM/iom.F90 (modified) (2 diffs)

NEMO/branches/2020/r12581_ticket2418/src/OCE/IOM/iom.F90

@@ -111 +111 @@
       CHARACTER(len=lc) :: clname
       INTEGER             :: irefyear, irefmonth, irefday
-      INTEGER           :: ji, jkmin
+      INTEGER           :: ji
       LOGICAL :: llrst_context              ! is context related to restart
       !
@@ -220 +220 @@
           
           ! Add vertical grid bounds
-          jkmin = MIN(2,jpk)  ! in case jpk=1 (i.e. sas2D)
-          zt_bnds(2,:        ) = gdept_1d(:)
-          zt_bnds(1,jkmin:jpk) = gdept_1d(1:jpkm1)
-          zt_bnds(1,1        ) = gdept_1d(1) - e3w_1d(1)
-          zw_bnds(1,:        ) = gdepw_1d(:)
-          zw_bnds(2,1:jpkm1  ) = gdepw_1d(jkmin:jpk)
-          zw_bnds(2,jpk:     ) = gdepw_1d(jpk) + e3t_1d(jpk)
+          zt_bnds(2,:      ) = gdept_1d(:)
+          zt_bnds(1,2:jpk  ) = gdept_1d(1:jpkm1)
+          zt_bnds(1,1      ) = gdept_1d(1) - e3w_1d(1)
+          zw_bnds(1,:      ) = gdepw_1d(:)
+          zw_bnds(2,1:jpkm1) = gdepw_1d(2:jpk)
+          zw_bnds(2,jpk:   ) = gdepw_1d(jpk) + e3t_1d(jpk)
           CALL iom_set_axis_attr(  "deptht", bounds=zw_bnds )
           CALL iom_set_axis_attr(  "depthu", bounds=zw_bnds )