Changeset 10673 for NEMO/branches/2019/fix_ticket2238_solution1/src/OCE/ICB/icbutl.F90

Timestamp:

2019-02-13T14:27:25+01:00 (5 years ago)

Author:

mathiot

Message:

creation of branch for ticket #2238 suggestion 1

Location:

NEMO/branches/2019/fix_ticket2238_solution1

Files:

• . (copied) (copied from NEMO/trunk)
• src/OCE/ICB/icbutl.F90 (modified) (1 diff)

NEMO/branches/2019/fix_ticket2238_solution1/src/OCE/ICB/icbutl.F90

@@ -223 +223 @@
       !
       ! find position in this processor. Prevent near edge problems (see #1389)
-      !
-      IF    ( ii < mig( 1 ) ) THEN   ;   ii = 1
-      ELSEIF( ii > mig(jpi) ) THEN   ;   ii = jpi
-      ELSE                           ;   ii = mi1(ii)
-      ENDIF
-      IF    ( ij < mjg( 1 ) ) THEN   ;   ij = 1
-      ELSEIF( ij > mjg(jpj) ) THEN   ;   ij = jpj
-      ELSE                           ;   ij  = mj1(ij)
-      ENDIF
-      !
-      IF( ii == jpi )   ii = ii-1      
-      IF( ij == jpj )   ij = ij-1
+      ! (PM) will be useless if extra halo is used in NEMO
+      !
+      IF    ( ii <= mig(1)-1 ) THEN   ;   ii = 0
+      ELSEIF( ii  > mig(jpi) ) THEN   ;   ii = jpi
+      ELSE                            ;   ii = mi1(ii)
+      ENDIF
+      IF    ( ij <= mjg(1)-1 ) THEN   ;   ij = 0
+      ELSEIF( ij  > mjg(jpj) ) THEN   ;   ij = jpj
+      ELSE                            ;   ij = mj1(ij)
+      ENDIF
+      !
       !
       icb_utl_bilin_h = ( pfld(ii,ij  ) * (1.-zi) + pfld(ii+1,ij  ) * zi ) * (1.-zj)   &