Changeset 15053 for NEMO/trunk/src/OCE/TRA

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.

Location:

NEMO/trunk/src/OCE/TRA

Files:

• eosbn2.F90 (modified) (4 diffs)
• trabbl.F90 (modified) (1 diff)

NEMO/trunk/src/OCE/TRA/eosbn2.F90

@@ -1078 +1078 @@
       z1_T0   = 1._wp/40._wp
       !
-      DO_2D( 1, 1, 1, 1 )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
          !
          zt  = ctmp   (ji,jj) * z1_T0
@@ -1142 +1142 @@
          !
          z1_S0 = 1._wp / 35.16504_wp
-         DO_2D( 1, 1, 1, 1 )
+         DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
             zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 )           ! square root salinity
             ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs &
@@ -1249 +1249 @@
       CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==!
          !
-         DO_3D( 1, 1, 1, 1, 1, jpkm1 )
+         DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 )
             !
             zh  = gdept(ji,jj,jk,Kmm) * r1_Z0                                ! depth
@@ -1308 +1308 @@
       CASE( np_seos )                !==  Vallis (2006) simplified EOS  ==!
          !
-         DO_3D( 1, 1, 1, 1, 1, jpkm1 )
+         DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 )
             zt  = pts(ji,jj,jk,jp_tem) - 10._wp  ! temperature anomaly (t-T0)
             zs = pts (ji,jj,jk,jp_sal) - 35._wp  ! abs. salinity anomaly (s-S0)

NEMO/trunk/src/OCE/TRA/trabbl.F90

@@ -339 +339 @@
       ENDIF
       !                                        !* bottom variables (T, S, alpha, beta, depth, velocity)
-      DO_2D( 1, 1, 1, 1 )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
          ik = mbkt(ji,jj)                             ! bottom T-level index
          zts (ji,jj,jp_tem) = ts(ji,jj,ik,jp_tem,Kbb) ! bottom before T and S