Changeset 14757 for NEMO/branches/2021/dev_r14393_HPC-03_Mele_Comm_Cleanup/src/OCE/DYN/dynldf_lap_blp.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/DYN/dynldf_lap_blp.F90 (modified) (5 diffs)

NEMO/branches/2021/dev_r14393_HPC-03_Mele_Comm_Cleanup/src/OCE/DYN/dynldf_lap_blp.F90

@@ -83 +83 @@
          DO jk = 1, jpkm1                                 ! Horizontal slab
             !
-            ! [comm_cleanup] ! DO_2D( 0, 1, 0, 1 ) 
             DO_2D( nn_hls-1, nn_hls, nn_hls-1, nn_hls )
                !                                      ! ahm * e3 * curl  (computed from 1 to jpim1/jpjm1)
@@ -95 +94 @@
             END_2D
             !
-            ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 )                       ! - curl( curl) + grad( div ) 
             DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 )                       ! - curl( curl) + grad( div )
                pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * umask(ji,jj,jk) * (    &    ! * by umask is mandatory for dyn_ldf_blp use
@@ -116 +114 @@
          DO jk = 1, jpkm1                                 ! Horizontal slab
             !
-            ! [comm_cleanup] ! DO_2D( 0, 1, 0, 1 ) 
             DO_2D( nn_hls-1, nn_hls, nn_hls-1, nn_hls )
                !                                      ! shearing stress component (F-point)   NB : ahmf has already been multiplied by fmask
@@ -132 +129 @@
             END_2D
             !
-            ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 ) 
             DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 )
                pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * r1_e1e2u(ji,jj) / e3u(ji,jj,jk,Kmm)                               &
@@ -189 +185 @@
       CALL dyn_ldf_lap( kt, Kbb, Kmm, pu, pv, zulap, zvlap, 1 )   ! rotated laplacian applied to pt (output in zlap,Kbb)
       !
-      IF (nn_hls.eq.1) CALL lbc_lnk( 'dynldf_lap_blp', zulap, 'U', -1.0_wp, zvlap, 'V', -1.0_wp )             ! Lateral boundary conditions
+      IF (nn_hls==1) CALL lbc_lnk( 'dynldf_lap_blp', zulap, 'U', -1.0_wp, zvlap, 'V', -1.0_wp )             ! Lateral boundary conditions
       !
       CALL dyn_ldf_lap( kt, Kbb, Kmm, zulap, zvlap, pu_rhs, pv_rhs, 2 )   ! rotated laplacian applied to zlap (output in pt(:,:,:,:,Krhs))