New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
crslbclnk.F90 in branches/UKMO/r8395_mix-lyr_diag/NEMOGCM/NEMO/OPA_SRC/CRS – NEMO

source: branches/UKMO/r8395_mix-lyr_diag/NEMOGCM/NEMO/OPA_SRC/CRS/crslbclnk.F90 @ 11290

Last change on this file since 11290 was 11290, checked in by jcastill, 5 years ago

Remove svn keywords

File size: 6.2 KB
Line 
1MODULE crslbclnk
2   !!======================================================================
3   !!                       ***  MODULE  crslbclnk  ***
4   !!               A temporary solution for lbclnk for coarsened grid.
5   !! Ocean        : lateral boundary conditions for grid coarsening
6   !!=====================================================================
7   !! History :   ! 2012-06  (J. Simeon, G. Madec, C. Ethe, C. Calone)     Original code
8   !!----------------------------------------------------------------------
9   USE par_kind, ONLY: wp
10   USE dom_oce
11   USE crs
12   !
13   USE lbclnk
14   USE in_out_manager
15   
16   INTERFACE crs_lbc_lnk
17      MODULE PROCEDURE crs_lbc_lnk_3d, crs_lbc_lnk_3d_gather, crs_lbc_lnk_2d
18   END INTERFACE
19   
20   PUBLIC crs_lbc_lnk
21   
22   !!----------------------------------------------------------------------
23   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
24   !! $Id$
25   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
26   !!----------------------------------------------------------------------
27CONTAINS
28
29   SUBROUTINE crs_lbc_lnk_3d( pt3d1, cd_type1, psgn, cd_mpp, pval )
30      !!---------------------------------------------------------------------
31      !!                  ***  SUBROUTINE crs_lbc_lnk  ***
32      !!
33      !! ** Purpose :   set lateral boundary conditions for coarsened grid
34      !!
35      !! ** Method  :   Swap domain indices from full to coarse domain
36      !!                before arguments are passed directly to lbc_lnk.
37      !!                Upon exiting, switch back to full domain indices.
38      !!----------------------------------------------------------------------
39      CHARACTER(len=1)                        , INTENT(in   ) ::   cd_type1 ! grid type
40      REAL(wp)                                , INTENT(in   ) ::   psgn     ! control of the sign
41      REAL(wp), DIMENSION(jpi_crs,jpj_crs,jpk), INTENT(inout) ::   pt3d1    ! 3D array on which the lbc is applied
42      REAL(wp)                      , OPTIONAL, INTENT(in   ) ::   pval     ! valeur sur les halo
43      CHARACTER(len=3)              , OPTIONAL, INTENT(in   ) ::   cd_mpp   ! MPP only (here do nothing)
44      !
45      LOGICAL  ::   ll_grid_crs
46      REAL(wp) ::   zval   ! valeur sur les halo
47      !!----------------------------------------------------------------------
48      !
49      ll_grid_crs = ( jpi == jpi_crs )
50      !
51      IF( PRESENT(pval) ) THEN  ;  zval = pval
52      ELSE                      ;  zval = 0._wp
53      ENDIF
54      !
55      IF( .NOT.ll_grid_crs )   CALL dom_grid_crs   ! Save the parent grid information  & Switch to coarse grid domain
56      !
57      IF( PRESENT( cd_mpp ) ) THEN ; CALL lbc_lnk( pt3d1, cd_type1, psgn, cd_mpp, pval=zval  )
58      ELSE                         ; CALL lbc_lnk( pt3d1, cd_type1, psgn, pval=zval  )
59      ENDIF
60      !
61      IF( .NOT.ll_grid_crs )   CALL dom_grid_glo   ! Return to parent grid domain
62      !
63   END SUBROUTINE crs_lbc_lnk_3d
64   
65   
66   SUBROUTINE crs_lbc_lnk_3d_gather( pt3d1, cd_type1, pt3d2, cd_type2, psgn )
67      !!---------------------------------------------------------------------
68      !!                  ***  SUBROUTINE crs_lbc_lnk  ***
69      !!
70      !! ** Purpose :   set lateral boundary conditions for coarsened grid
71      !!
72      !! ** Method  :   Swap domain indices from full to coarse domain
73      !!                before arguments are passed directly to lbc_lnk.
74      !!                Upon exiting, switch back to full domain indices.
75      !!----------------------------------------------------------------------
76      CHARACTER(len=1)                        , INTENT(in   ) ::   cd_type1, cd_type2 ! grid type
77      REAL(wp)                                , INTENT(in   ) ::   psgn               ! control of the sign
78      REAL(wp), DIMENSION(jpi_crs,jpj_crs,jpk), INTENT(inout) ::   pt3d1   , pt3d2    ! 3D array on which the lbc is applied
79      !
80      LOGICAL ::   ll_grid_crs
81      !!----------------------------------------------------------------------
82      !
83      ll_grid_crs = ( jpi == jpi_crs )
84      !
85      IF( .NOT.ll_grid_crs )   CALL dom_grid_crs   ! Save the parent grid information  & Switch to coarse grid domain
86      !
87      CALL lbc_lnk( pt3d1, cd_type1, pt3d2, cd_type2, psgn  )
88      !
89      IF( .NOT.ll_grid_crs )   CALL dom_grid_glo   ! Return to parent grid domain
90      !
91   END SUBROUTINE crs_lbc_lnk_3d_gather
92
93   
94   
95   SUBROUTINE crs_lbc_lnk_2d(pt2d, cd_type, psgn, cd_mpp, pval)
96      !!---------------------------------------------------------------------
97      !!                  ***  SUBROUTINE crs_lbc_lnk  ***
98      !!
99      !! ** Purpose :   set lateral boundary conditions for coarsened grid
100      !!
101      !! ** Method  :   Swap domain indices from full to coarse domain
102      !!                before arguments are passed directly to lbc_lnk.
103      !!                Upon exiting, switch back to full domain indices.
104      !!----------------------------------------------------------------------
105      CHARACTER(len=1)                    , INTENT(in   ) ::   cd_type  ! grid type
106      REAL(wp)                            , INTENT(in   ) ::   psgn     ! control of the sign
107      REAL(wp), DIMENSION(jpi_crs,jpj_crs), INTENT(inout) ::   pt2d     ! 2D array on which the lbc is applied
108      REAL(wp)                  , OPTIONAL, INTENT(in   ) ::   pval     ! valeur sur les halo
109      CHARACTER(len=3)          , OPTIONAL, INTENT(in   ) ::   cd_mpp   ! MPP only (here do nothing)
110      !     
111      LOGICAL  ::   ll_grid_crs
112      REAL(wp) ::   zval     ! valeur sur les halo
113      !!----------------------------------------------------------------------
114      !
115      ll_grid_crs = ( jpi == jpi_crs )
116      !
117      IF( PRESENT(pval) ) THEN  ;  zval = pval
118      ELSE                      ;  zval = 0._wp
119      ENDIF
120      !
121      IF( .NOT.ll_grid_crs )   CALL dom_grid_crs   ! Save the parent grid information  & Switch to coarse grid domain
122      !
123      IF( PRESENT( cd_mpp ) ) THEN ; CALL lbc_lnk( pt2d, cd_type, psgn, cd_mpp, pval=zval  )
124      ELSE                         ; CALL lbc_lnk( pt2d, cd_type, psgn, pval=zval  )
125      ENDIF
126      !
127      IF( .NOT.ll_grid_crs )   CALL dom_grid_glo   ! Return to parent grid domain
128      !
129   END SUBROUTINE crs_lbc_lnk_2d
130
131   !!======================================================================
132END MODULE crslbclnk
Note: See TracBrowser for help on using the repository browser.