Changeset 14757 for NEMO/branches/2021/dev_r14393_HPC-03_Mele_Comm_Cleanup/src/OCE/TRA/traldf_iso.F90

Timestamp:

2021-04-27T17:33:44+02:00 (3 years ago)

Author:

francesca

Message:

Fortran 77 '.EQ.' operator replacement in conditional statements; [comm_cleanup] tags removal - ticket #2607

File:

• NEMO/branches/2021/dev_r14393_HPC-03_Mele_Comm_Cleanup/src/OCE/TRA/traldf_iso.F90 (modified) (17 diffs)

NEMO/branches/2021/dev_r14393_HPC-03_Mele_Comm_Cleanup/src/OCE/TRA/traldf_iso.F90

@@ -147 +147 @@
          ENDIF
          !
-         ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpk )
          DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk ) 
             akz     (ji,jj,jk) = 0._wp
@@ -173 +172 @@
       IF( kpass == 1 ) THEN                  !==  first pass only  ==!
          !
-         ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )
          DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
             !
@@ -196 +194 @@
          !
          IF( ln_traldf_msc ) THEN                ! stabilizing vertical diffusivity coefficient
-            ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )
             DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
                ! round brackets added to fix the order of floating point operations
@@ -210 +207 @@
             !
             IF( ln_traldf_blp ) THEN                ! bilaplacian operator
-               ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )
                DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
                   akz(ji,jj,jk) = 16._wp   &
@@ -219 +215 @@
                END_3D
             ELSEIF( ln_traldf_lap ) THEN              ! laplacian operator
-               ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )
                DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
                   ze3w_2 = e3w(ji,jj,jk,Kmm) * e3w(ji,jj,jk,Kmm)
@@ -228 +223 @@
            !
          ELSE                                    ! 33 flux set to zero with akz=ah_wslp2 ==>> computed in full implicit
-            ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpk )
             DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk ) 
                akz(ji,jj,jk) = ah_wslp2(ji,jj,jk)
@@ -248 +242 @@
 
          ! Horizontal tracer gradient
-         ! [comm_cleanup] ! DO_3D( 1, 0, 1, 0, 1, jpkm1 )
          DO_3D( nn_hls, nn_hls-1, nn_hls, nn_hls-1, 1, jpkm1 ) 
             zdit(ji,jj,jk) = ( pt(ji+1,jj  ,jk,jn) - pt(ji,jj,jk,jn) ) * umask(ji,jj,jk)
@@ -254 +247 @@
          END_3D
          IF( ln_zps ) THEN      ! botton and surface ocean correction of the horizontal gradient
-            ! [comm_cleanup] ! DO_2D( 1, 0, 1, 0 )           ! bottom correction (partial bottom cell)
                DO_2D( nn_hls, nn_hls-1, nn_hls, nn_hls-1 )            ! bottom correction (partial bottom cell)
                zdit(ji,jj,mbku(ji,jj)) = pgu(ji,jj,jn)
@@ -260 +252 @@
             END_2D
             IF( ln_isfcav ) THEN      ! first wet level beneath a cavity
-               ! [comm_cleanup] ! DO_2D( 1, 0, 1, 0 )
                DO_2D( nn_hls, nn_hls-1, nn_hls, nn_hls-1 ) 
                   IF( miku(ji,jj) > 1 )   zdit(ji,jj,miku(ji,jj)) = pgui(ji,jj,jn)
@@ -274 +265 @@
          DO jk = 1, jpkm1                                 ! Horizontal slab
             !
-            ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 )
             DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
                !                             !== Vertical tracer gradient
@@ -284 +274 @@
             END_2D
             !
-            ! [comm_cleanup] ! DO_2D( 1, 0, 1, 0 )           !==  Horizontal fluxes
             DO_2D( nn_hls, nn_hls-1, nn_hls, nn_hls-1 )           !==  Horizontal fluxes
                zabe1 = pahu(ji,jj,jk) * e2_e1u(ji,jj) * e3u(ji,jj,jk,Kmm)
@@ -314 +303 @@
             END_2D
             !
-            ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 )           !== horizontal divergence and add to pta
             DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 )           !== horizontal divergence and add to pta
                ! round brackets added to fix the order of floating point operations
@@ -336 +324 @@
          ztfw(:,:, 1 ) = 0._wp      ;      ztfw(:,:,jpk) = 0._wp
 
-         ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )    ! interior (2=<jk=<jpk-1)
          DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 )    ! interior (2=<jk=<jpk-1)
             !
@@ -367 +354 @@
          !                                !==  add the vertical 33 flux  ==!
          IF( ln_traldf_lap ) THEN               ! laplacian case: eddy coef = ah_wslp2 - akz
-            ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )
             DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
                ztfw(ji,jj,jk) = ztfw(ji,jj,jk) + e1e2t(ji,jj) / e3w(ji,jj,jk,Kmm) * wmask(ji,jj,jk)   &
@@ -377 +363 @@
             SELECT CASE( kpass )
             CASE(  1  )                            ! 1st pass : eddy coef = ah_wslp2
-               ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )
                DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
                   ztfw(ji,jj,jk) =   &
@@ -384 +369 @@
                END_3D
             CASE(  2  )                         ! 2nd pass : eddy flux = ah_wslp2 and akz applied on pt  and pt2 gradients, resp.
-               ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 )
                DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
                   ztfw(ji,jj,jk) = ztfw(ji,jj,jk) + e1e2t(ji,jj) / e3w(ji,jj,jk,Kmm) * wmask(ji,jj,jk)                  &
@@ -393 +377 @@
          ENDIF
          !
-         ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpkm1 )    !==  Divergence of vertical fluxes added to pta  ==!
          DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpkm1 )    !==  Divergence of vertical fluxes added to pta  ==! 
             pt_rhs(ji,jj,jk,jn) = pt_rhs(ji,jj,jk,jn) + zsign * (  ztfw (ji,jj,jk) - ztfw(ji,jj,jk+1)  ) * r1_e1e2t(ji,jj)   &