Changeset 11048 for NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdylib.F90

Timestamp:

2019-05-23T18:36:06+02:00 (5 years ago)

Author:

girrmann

Message:

dev_r10984_HPC-13 : Step 1, boundary is now detected all over the local domain, this does not change the result. Improve bdy treatment for bdy_rnf in bdytra.F90, this changes the result when keyword runoff is specified in namelist

File:

• NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdylib.F90 (modified) (3 diffs)

NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdylib.F90

@@ -56 +56 @@
             ii = idx%nbi(ib,igrd) 
             ij = idx%nbj(ib,igrd)
+            IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove
             zwgt = idx%nbw(ib,igrd)
             pta(ii,ij,ik) = ( pta(ii,ij,ik) + zwgt * (dta(ib,ik) - pta(ii,ij,ik) ) ) * tmask(ii,ij,ik)
@@ -83 +84 @@
          ii = idx%nbi(ib,igrd)
          ij = idx%nbj(ib,igrd)
+         IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove
          DO ik = 1, jpkm1
             pta(ii,ij,ik) = dta(ib,ik) * tmask(ii,ij,ik)
@@ -466 +468 @@
          ii = idx%nbi(ib,igrd)
          ij = idx%nbj(ib,igrd)
-         !
+         IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove
          SELECT CASE( idx%ntreat(ib,igrd) )   ! select free ocean neighbours
             !     o