Changeset 3837 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/LBC/lbclnk.F90

Timestamp:

2013-03-12T15:55:32+01:00 (11 years ago)

Author:

trackstand2

Message:

Merge of finiss

File:

• branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/LBC/lbclnk.F90 (modified) (7 diffs)

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/LBC/lbclnk.F90

@@ -18 +18 @@
    USE exchmod          ! Comms for irregular domain decomposition
 
+   ! This is important - it determines which set of comms routines are
+   ! called when lbc_lnk() is invoked.
    INTERFACE lbc_lnk
 #if  defined key_mpp_rkpart
@@ -104 +106 @@
 
 
-   SUBROUTINE lbc_lnk_3d( pt3d, cd_type, psgn, cd_mpp, pval )
+   SUBROUTINE lbc_lnk_3d( pt3d, cd_type, psgn, cd_mpp, pval, lzero )
+      USE lib_mpp, ONLY: ctl_stop
       !!---------------------------------------------------------------------
       !!                  ***  ROUTINE lbc_lnk_3d  ***
@@ -124 +127 @@
       CHARACTER(len=3)                , INTENT(in   ), OPTIONAL ::   cd_mpp    ! MPP only (here do nothing)
       REAL(wp)                        , INTENT(in   ), OPTIONAL ::   pval      ! background value (for closed boundaries)
+      LOGICAL                         , INTENT(in   ), OPTIONAL ::   lzero     ! Whether to zero halos on closed boundaries
+
       !!
       REAL(wp) ::   zland
@@ -132 +137 @@
       ENDIF
 
+      IF( PRESENT( lzero ) )THEN
+         CALL ctl_stop('STOP','lbc_lnk_3d: IMPLEMENT lzero option!')
+      ENDIF
 
       IF( PRESENT( cd_mpp ) ) THEN
@@ -197 +205 @@
 
 
-   SUBROUTINE lbc_lnk_2d( pt2d, cd_type, psgn, cd_mpp, pval )
+   SUBROUTINE lbc_lnk_2d( pt2d, cd_type, psgn, cd_mpp, pval, lzero )
+      USE lib_mpp, ONLY: ctl_stop
       !!---------------------------------------------------------------------
       !!                 ***  ROUTINE lbc_lnk_2d  ***
@@ -214 +223 @@
       CHARACTER(len=3)            , INTENT(in   ), OPTIONAL ::   cd_mpp    ! MPP only (here do nothing)
       REAL(wp)                    , INTENT(in   ), OPTIONAL ::   pval      ! background value (for closed boundaries)
+      LOGICAL                     , INTENT(in   ), OPTIONAL ::   lzero    ! Whether to zero halos on closed boundaries
       !!
       REAL(wp) ::   zland
@@ -220 +230 @@
       IF( PRESENT( pval ) ) THEN   ;   zland = pval      ! set land value (zero by default)
       ELSE                         ;   zland = 0.e0
+      ENDIF
+
+      IF( PRESENT( lzero ) )THEN
+         CALL ctl_stop('STOP','lbc_lnk_2d: IMPLEMENT lzero option!')
       ENDIF