Changeset 15053 for NEMO/trunk/src/OCE/ISF/isftbl.F90

Timestamp:

2021-06-24T17:39:38+02:00 (3 years ago)

Author:

clem

Message:

nn_hls=2: solve a number of floating points because some arrays are not defined on the first halo. Looks like this is not the end.

File:

• NEMO/trunk/src/OCE/ISF/isftbl.F90 (modified) (3 diffs)

NEMO/trunk/src/OCE/ISF/isftbl.F90

@@ -78 +78 @@
          ! compute tbl property at T point
          pvarout(1,:) = 0._wp
-         DO_2D( 0, 1, 1, 1 )
+         DO_2D( nn_hls-1, nn_hls, nn_hls, nn_hls )
             pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji-1,jj))
          END_2D
@@ -99 +99 @@
          ! pvarout is an averaging of wet point
          pvarout(:,1) = 0._wp
-         DO_2D( 1, 1, 0, 1 )
+         DO_2D( nn_hls, nn_hls, nn_hls-1, nn_hls )
             pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji,jj-1))
          END_2D
@@ -138 +138 @@
       !
       ! compute tbl top.bottom level and thickness
-      DO_2D( 1, 1, 1, 1 )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
          !
          ! tbl top/bottom indices initialisation