Changeset 15795

Timestamp:

2022-04-23T10:33:25+02:00 (2 years ago)

Author:

smasson

Message:

bugfix in ar5, see ​https://forge.nemo-ocean.eu/nemo/nemo/-/issues/39

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diaar5.F90 (modified) (2 diffs)

NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diaar5.F90

@@ -377 +377 @@
       !!----------------------------------------------------------------------
       INTEGER  ::   inum
-      INTEGER  ::   ik, idep
+      INTEGER  ::   ik
       INTEGER  ::   ji, jj, jk  ! dummy loop indices
       REAL(wp) ::   zztmp  
@@ -404 +404 @@
             DO jj = 1, jpj               ! interpolation of salinity at the last ocean level (i.e. the partial step)
                DO ji = 1, jpi
-                  idep = tmask(ji,jj,jk) * e3t_0(ji,jj,jk)
-                  zvol0 (ji,jj) = zvol0 (ji,jj) +  idep * e1e2t(ji,jj)
-                  thick0(ji,jj) = thick0(ji,jj) +  idep    
+                  zztmp = tmask(ji,jj,jk) * e3t_0(ji,jj,jk)
+                  zvol0 (ji,jj) = zvol0 (ji,jj) + zztmp * e1e2t(ji,jj)
+                  thick0(ji,jj) = thick0(ji,jj) + zztmp    
                END DO
             END DO