Changeset 14260

Timestamp:

2021-01-05T14:35:55+01:00 (4 years ago)

Author:

smueller

Message:

Two bugfixes and complete initialisation of local arrays in internal subroutines of module subroutine zdf_osm (ticket #2353)

File:

• NEMO/branches/2020/dev_r14122_ENHANCE-14_smueller_OSMOSIS_streamlining/src/OCE/ZDF/zdfosm.F90 (modified) (4 diffs)

NEMO/branches/2020/dev_r14122_ENHANCE-14_smueller_OSMOSIS_streamlining/src/OCE/ZDF/zdfosm.F90

@@ -1681 +1681 @@
              zpe_mle_ref = 0._wp
              DO jk = ibld(ji,jj), mld_prof(ji,jj)
+               zthermal = rab_n(ji,jj,jk,jp_tem)
+               zbeta    = rab_n(ji,jj,jk,jp_sal)
                zbuoy = grav * ( zthermal * ts(ji,jj,jk,jp_tem,Kmm) - zbeta * ts(ji,jj,jk,jp_sal,Kmm) )
                zpe_mle_layer = zpe_mle_layer + zbuoy * gdepw(ji,jj,jk,Kmm) * e3w(ji,jj,jk,Kmm)
@@ -1715 +1717 @@
                   lpyc(ji,jj) = .FALSE.
                 ELSE
-                   lpyc = .TRUE.
+                   lpyc(ji,jj) = .TRUE.
                 ENDIF
               ELSE
@@ -2373 +2375 @@
       ! calculate horizontal gradients at u & v points
 
+      zmld_midu(:,:)   = 0.0_wp
+      ztsm_midu(:,:,:) = 0.0_wp
       DO_2D( 0, 0, 1, 0 )
          zdtdx(ji,jj) = ( ztm(ji+1,jj) - ztm( ji,jj) )  * umask(ji,jj,1) / e1u(ji,jj)
@@ -2381 +2385 @@
       END_2D
 
+      zmld_midv(:,:)   = 0.0_wp
+      ztsm_midv(:,:,:) = 0.0_wp
       DO_2D( 1, 0, 0, 0 )
          zdtdy(ji,jj) = ( ztm(ji,jj+1) - ztm( ji,jj) ) * vmask(ji,jj,1) / e1v(ji,jj)