Changeset 13288 for NEMO/trunk/src/OCE/DYN/dynhpg.F90

Timestamp:

2020-07-09T18:23:18+02:00 (4 years ago)

Author:

smasson

Message:

trunk: bugfix following #2474

File:

• NEMO/trunk/src/OCE/DYN/dynhpg.F90 (modified) (1 diff)

NEMO/trunk/src/OCE/DYN/dynhpg.F90

@@ -961 +961 @@
       ! Preparing vertical density profile "zrhh(:,:,:)" for hybrid-sco coordinate
       DO_2D_11_11
-       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(ji,jj,jkk  ), gde3w(ji,jj,jkk-1),   &
-                &                      gde3w(ji,jj,jkk-2), rhd    (ji,jj,jkk-1), rhd(ji,jj,jkk-2))
+                &                      gde3w(ji,jj,jkk-2), zrhh (ji,jj,jkk-1), zrhh(ji,jj,jkk-2))
           END DO
        ENDIF