Changeset 15104 for NEMO/trunk/src/ABL

Timestamp:

2021-07-07T16:36:00+02:00 (3 years ago)

Author:

clem

Message:

nn_hls=2: repair some loops.

File:

• NEMO/trunk/src/ABL/ablmod.F90 (modified) (10 diffs)

NEMO/trunk/src/ABL/ablmod.F90

@@ -128 +128 @@
       !! needed for surface boundary condition of TKE
       !! pwndm contains | U10m - U_oce | (see blk_oce_1 in sbcblk)
-      DO_2D( 1, 1, 1, 1 )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
          zzoce         = pCd_du    (ji,jj) * pwndm    (ji,jj)
 #if defined key_si3
@@ -250 +250 @@
             !
             ! Advance u_abl & v_abl to time n+1
-            DO_2D( 1, 1, 1, 1 )
+            DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
                zcff = ( fft_abl(ji,jj) * rDt_abl )*( fft_abl(ji,jj) * rDt_abl )  ! (f dt)**2
 
@@ -488 +488 @@
          DO jk = 2, jpka    ! outer loop
          !-------------
-            DO_2D( 0, 1, 0, 1 )
+            DO_2D( nn_hls-1, nn_hls, nn_hls-1, nn_hls )
                zcff1 = pblh( ji, jj )
                zsig  = ght_abl(jk) / MAX( jp_pblh_min,  MIN(  jp_pblh_max, zcff1  ) )
@@ -511 +511 @@
       DO jk = 2, jpka    ! outer loop
       !-------------
-         DO_2D( 1, 1, 1, 1 )
+         DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
             zcff1 = pblh( ji, jj )
             zsig  = ght_abl(jk) / MAX( jp_pblh_min,  MIN(  jp_pblh_max, zcff1  ) )
@@ -585 +585 @@
       !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
       !
-      DO_2D( 1, 1, 1, 1 )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
          ztemp          =  tq_abl( ji, jj, 2, nt_a, jp_ta )
          zhumi          =  tq_abl( ji, jj, 2, nt_a, jp_qa )
@@ -595 +595 @@
       END_2D
 
-      DO_2D( 0, 1, 0, 1 )
+      DO_2D( nn_hls-1, nn_hls, nn_hls-1, nn_hls )
          zwnd_i(ji,jj) = u_abl(ji  ,jj,2,nt_a) - 0.5_wp * ( pssu(ji  ,jj) + pssu(ji-1,jj) )
          zwnd_j(ji,jj) = v_abl(ji,jj  ,2,nt_a) - 0.5_wp * ( pssv(ji,jj  ) + pssv(ji,jj-1) )
@@ -603 +603 @@
       !
       ! ... scalar wind ( = | U10m - U_oce | ) at T-point (masked)
-      DO_2D( 1, 1, 1, 1 )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
          zcff          = SQRT(  zwnd_i(ji,jj) * zwnd_i(ji,jj)   &
             &                 + zwnd_j(ji,jj) * zwnd_j(ji,jj) )   ! * msk_abl(ji,jj)
@@ -1265 +1265 @@
       smth_b = 1._wp / 4._wp
       !
-      DO_2D( 1, 0, 1, 1 )
+      DO_2D( nn_hls, nn_hls-1, nn_hls, nn_hls )
          zumsk = msk(ji,jj) * msk(ji+1,jj)
          zdX ( ji, jj ) = ( pvar2d( ji+1,jj ) - pvar2d( ji  ,jj ) ) * zumsk
       END_2D
 
-      DO_2D( 1, 1, 1, 0 )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls-1 )
          zvmsk = msk(ji,jj) * msk(ji,jj+1)
          zdY ( ji, jj ) = ( pvar2d( ji, jj+1 ) - pvar2d( ji  ,jj ) ) * zvmsk
       END_2D
 
-      DO_2D( 0, 0, 1, 0 )
+      DO_2D( nn_hls-1, nn_hls-1, nn_hls, nn_hls-1 )
          zFY ( ji, jj  ) =   zdY ( ji, jj   )                        &
             & +  smth_a*  ( (zdX ( ji, jj+1 ) - zdX( ji-1, jj+1 ))   &
@@ -1281 +1281 @@
       END_2D
 
-      DO_2D( 1, 0, 0, 0 )
+      DO_2D( nn_hls, nn_hls-1, nn_hls-1, nn_hls-1 )
          zFX( ji, jj  ) =    zdX( ji, jj   )                         &
            &    + smth_a*(  (zdY( ji+1, jj ) - zdY( ji+1, jj-1))     &
@@ -1287 +1287 @@
       END_2D
 
-      DO_2D( 0, 0, 0, 0 )
+      DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 )
          pvar2d( ji  ,jj ) = pvar2d( ji  ,jj )              &
   &         + msk(ji,jj) * smth_b * (                       &