Changeset 13295 for NEMO/trunk/src/SWE/dynldf_lap_blp.F90

Timestamp:

2020-07-10T20:24:21+02:00 (4 years ago)

Author:

acc

Message:

Replace do-loop macros in the trunk with alternative forms with greater flexibility for extra halo applications. This alters a lot of routines but does not change any behaviour or results. do_loop_substitute.h90 is greatly simplified by this change. SETTE results are identical to those with the previous revision

File:

• NEMO/trunk/src/SWE/dynldf_lap_blp.F90 (modified) (6 diffs)

NEMO/trunk/src/SWE/dynldf_lap_blp.F90

@@ -97 +97 @@
             DO jk = 1, jpkm1                                 ! Horizontal slab
                !                                             
-               DO_2D_01_01
+               DO_2D( 0, 1, 0, 1 )
                !                                      ! ahm * e3 * curl  (computed from 1 to jpim1/jpjm1)
 !!gm open question here : e3f  at before or now ?    probably now... 
@@ -112 +112 @@
                END_2D
                !
-               DO_2D_00_00
+               DO_2D( 0, 0, 0, 0 )
                   pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * (                                             &
                      &              - ( zcur(ji  ,jj) - zcur(ji,jj-1) ) * r1_e2u(ji,jj) / e3u(ji,jj,jk,Kmm)   &
@@ -128 +128 @@
             DO jk = 1, jpkm1                                 ! Horizontal slab
                !
-               DO_2D_01_01
+               DO_2D( 0, 1, 0, 1 )
                   !                                      ! shearing stress component (F-point)   NB : ahmf has already been multiplied by fmask
                   zshe(ji-1,jj-1) = ahmf(ji-1,jj-1,jk)                                                              &
@@ -143 +143 @@
                END_2D
                !
-               DO_2D_00_00
+               DO_2D( 0, 0, 0, 0 )
                   pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * r1_e1e2u(ji,jj) / e3u(ji,jj,jk,Kmm)                               &
                      &    * (   (   zten(ji+1,jj  ) * e2t(ji+1,jj  )*e2t(ji+1,jj  ) * e3t(ji+1,jj  ,jk,Kmm)                       &
@@ -164 +164 @@
             DO jk = 1, jpkm1                                 ! Horizontal slab
                !
-               DO_2D_01_01
+               DO_2D( 0, 1, 0, 1 )
                   !                                      ! shearing stress component (F-point)   NB : ahmf has already been multiplied by fmask
                   zshe(ji-1,jj-1) = ahmf(ji-1,jj-1,jk)                                           &
@@ -175 +175 @@
                END_2D
                !
-               DO_2D_00_00
+               DO_2D( 0, 0, 0, 0 )
                   pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * r1_e1e2u(ji,jj) / e3u(ji,jj,jk,Kmm)   &
                      &    * (   zten(ji+1,jj  ) * e2t(ji+1,jj  ) * e3t(ji+1,jj  ,jk,Kmm)              &