Changeset 11633

Timestamp:

2019-10-01T21:06:31+02:00 (5 years ago)

Author:

mathiot

Message:

modification for in case a domain is a land domain (MINLOC return 0 as indices for the min) (ticket #2143)

File:

• NEMO/branches/2019/ENHANCE-03_domcfg/src/domngb.F90 (modified) (3 diffs)

NEMO/branches/2019/ENHANCE-03_domcfg/src/domngb.F90

@@ -60 +60 @@
 
       zdist = dist(plon, plat, zglam, zgphi)
+
+      WHERE (zmask(:,:) == 0._wp)
+         zdist(:,:) = HUGE(1._wp)
+      END WHERE
+      zmask = 1._wp
       
       IF( lk_mpp ) THEN  
@@ -74 +79 @@
 
    FUNCTION dist(plonsrc, platsrc, plontrg, plattrg)
-      REAL(wp), INTENT(in) :: plonsrc, platsrc                     ! lat/lon of the source point
+      REAL(wp),                     INTENT(in) :: plonsrc, platsrc ! lat/lon of the source point
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: plontrg, plattrg ! lat/lon of the target (2d in this case)
 
@@ -94 +99 @@
          &      + ( zzs - zzt(:,:) )**2
 
-      dist = ra * SQRT( dist )
+      dist(:,:) = ra * SQRT( dist(:,:) )
 
    END FUNCTION dist