Changeset 11633


Ignore:
Timestamp:
2019-10-01T21:06:31+02:00 (13 months 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:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-03_domcfg/src/domngb.F90

    r11201 r11633  
    6060 
    6161      zdist = dist(plon, plat, zglam, zgphi) 
     62 
     63      WHERE (zmask(:,:) == 0._wp) 
     64         zdist(:,:) = HUGE(1._wp) 
     65      END WHERE 
     66      zmask = 1._wp 
    6267       
    6368      IF( lk_mpp ) THEN   
     
    7479 
    7580   FUNCTION dist(plonsrc, platsrc, plontrg, plattrg) 
    76       REAL(wp), INTENT(in) :: plonsrc, platsrc                    ! lat/lon of the source point 
     81      REAL(wp),                     INTENT(in) :: plonsrc, platsrc ! lat/lon of the source point 
    7782      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: plontrg, plattrg ! lat/lon of the target (2d in this case) 
    7883 
     
    9499         &      + ( zzs - zzt(:,:) )**2 
    95100 
    96       dist = ra * SQRT( dist ) 
     101      dist(:,:) = ra * SQRT( dist(:,:) ) 
    97102 
    98103   END FUNCTION dist 
Note: See TracChangeset for help on using the changeset viewer.