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

Timestamp:

2019-07-01T12:10:15+02:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-03_domcfg : add management of closed seas in domain cfg by flood filling and lat/lon seed instead of i/j box definition (ticket #2143)

File:

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

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

@@ -11 +11 @@
    !!----------------------------------------------------------------------
    USE dom_oce        ! ocean space and time domain
+   USE phycst
    !
    USE in_out_manager ! I/O manager
@@ -18 +19 @@
    PRIVATE
 
-   PUBLIC   dom_ngb   ! routine called in iom.F90 module
+   PUBLIC   dom_ngb   ! routine called in iom.F90 and domclo.F90 module
+   PUBLIC   dist
 
    !!----------------------------------------------------------------------
@@ -27 +29 @@
 CONTAINS
 
-   SUBROUTINE dom_ngb( plon, plat, kii, kjj, cdgrid, kkk )
+   SUBROUTINE dom_ngb( plon, plat, kii, kjj, rdist, cdgrid, kkk )
       !!----------------------------------------------------------------------
       !!                    ***  ROUTINE dom_ngb  ***
@@ -37 +39 @@
       !!----------------------------------------------------------------------
       REAL(wp)        , INTENT(in   ) ::   plon, plat   ! longitude,latitude of the point
+      REAL(wp)        , INTENT(  out) ::   rdist        ! distance between the located point and the source
       INTEGER         , INTENT(  out) ::   kii, kjj     ! i-,j-index of the closes grid point
       INTEGER         , INTENT(in   ), OPTIONAL :: kkk  ! k-index of the mask level used
@@ -43 +46 @@
       INTEGER :: ik         ! working level
       INTEGER , DIMENSION(2) ::   iloc
-      REAL(wp)               ::   zlon, zmini
       REAL(wp), DIMENSION(jpi,jpj) ::   zglam, zgphi, zmask, zdist
       !!--------------------------------------------------------------------
@@ -57 +59 @@
       END SELECT
 
-      zlon       = MOD( plon       + 720., 360. )                                     ! plon between    0 and 360
-      zglam(:,:) = MOD( zglam(:,:) + 720., 360. )                                     ! glam between    0 and 360
-      IF( zlon > 270. )   zlon = zlon - 360.                                          ! zlon between  -90 and 270
-      IF( zlon <  90. )   WHERE( zglam(:,:) > 180. ) zglam(:,:) = zglam(:,:) - 360.   ! glam between -180 and 180
-      zglam(:,:) = zglam(:,:) - zlon
-
-      zgphi(:,:) = zgphi(:,:) - plat
-      zdist(:,:) = zglam(:,:) * zglam(:,:) + zgphi(:,:) * zgphi(:,:)
+      zdist = dist(plon, plat, zglam, zgphi)
       
       IF( lk_mpp ) THEN  
-         CALL mpp_minloc( 'domngb', zdist(:,:), zmask, zmini, iloc)
+         CALL mpp_minloc( 'domngb', zdist(:,:), zmask, rdist, iloc)
          kii = iloc(1) ; kjj = iloc(2)
       ELSE
          iloc(:) = MINLOC( zdist(:,:), mask = zmask(:,:) == 1.e0 )
+         rdist = MINVAL( zdist(:,:) )
          kii = iloc(1) + nimpp - 1
          kjj = iloc(2) + njmpp - 1
@@ -77 +73 @@
    END SUBROUTINE dom_ngb
 
+   FUNCTION dist(plonsrc, platsrc, plontrg, plattrg)
+      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)
+
+      REAL(wp) :: zxs, zys, zzs
+      REAL(wp), DIMENSION(jpi,jpj) :: zxt, zyt, zzt
+
+      REAL(wp), DIMENSION(jpi,jpj) :: dist ! distance between src and trg
+
+      zxs = COS( rad * platsrc ) * COS( rad * plonsrc )
+      zys = COS( rad * platsrc ) * SIN( rad * plonsrc )
+      zzs = SIN( rad * platsrc )
+
+      zxt = COS( rad * plattrg ) * COS( rad * plontrg )
+      zyt = COS( rad * plattrg ) * SIN( rad * plontrg )
+      zzt = SIN( rad * plattrg )
+
+      dist(:,:) = ( zxs - zxt(:,:) )**2   &
+         &      + ( zys - zyt(:,:) )**2   &
+         &      + ( zzs - zzt(:,:) )**2
+
+      dist = ra * SQRT( dist )
+
+   END FUNCTION dist
+
    !!======================================================================
 END MODULE domngb