Changeset 13095 for NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DYN/dynhpg.F90

Timestamp:

2020-06-11T14:08:55+02:00 (4 years ago)

Author:

jamesharle

Message:

Update indexing in hpg_prj (see ticket #2474)

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DYN/dynhpg.F90 (modified) (1 diff)

NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DYN/dynhpg.F90

@@ -1027 +1027 @@
       DO jj = 1, jpj
         DO ji = 1, jpi
-          jk = mbkt(ji,jj)+1
-          IF(     jk <=  0   ) THEN   ;   zrhh(ji,jj,    :   ) = 0._wp
-          ELSEIF( jk ==  1   ) THEN   ;   zrhh(ji,jj,jk+1:jpk) = rhd(ji,jj,jk)
+          jk = mbkt(ji,jj)
+          IF(     jk <=  1   ) THEN   ;   zrhh(ji,jj,    :   ) = 0._wp
+          ELSEIF( jk ==  2   ) THEN   ;   zrhh(ji,jj,jk+1:jpk) = rhd(ji,jj,jk)
           ELSEIF( jk < jpkm1 ) THEN
              DO jkk = jk+1, jpk
                 zrhh(ji,jj,jkk) = interp1(gde3w_n(ji,jj,jkk  ), gde3w_n(ji,jj,jkk-1),   &
-                   &                      gde3w_n(ji,jj,jkk-2), rhd    (ji,jj,jkk-1), rhd(ji,jj,jkk-2))
+                   &                      gde3w_n(ji,jj,jkk-2), zrhh    (ji,jj,jkk-1), zrhh(ji,jj,jkk-2))
              END DO
           ENDIF