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.
lbc_lnk_multi_generic.h90 in NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_DYN_optimization/src/OCE/LBC – NEMO

source: NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_DYN_optimization/src/OCE/LBC/lbc_lnk_multi_generic.h90 @ 11380

Last change on this file since 11380 was 11380, checked in by girrmann, 5 years ago

dev_r10984_HPC-13 : adding extra halos in dyn_spg_ts is now possible, only works with a single halo when used with tide or bdy, see #2308
  • Property svn:keywords set to Id
  • Property svn:mime-type set to text/x-fortran
File size: 7.1 KB
Line 
1#if defined DIM_2d
2#   define ARRAY_TYPE(i,j,k,l)   REAL(wp), DIMENSION(i,j)
3#   define PTR_TYPE              TYPE(PTR_2D)
4#   define PTR_ptab              pt2d
5#endif
6#if defined DIM_3d
7#   define ARRAY_TYPE(i,j,k,l)   REAL(wp), DIMENSION(i,j,k)
8#   define PTR_TYPE              TYPE(PTR_3D)
9#   define PTR_ptab              pt3d
10#endif
11#if defined DIM_4d
12#   define ARRAY_TYPE(i,j,k,l)   REAL(wp), DIMENSION(i,j,k,l)
13#   define PTR_TYPE              TYPE(PTR_4D)
14#   define PTR_ptab              pt4d
15#endif
16
17   SUBROUTINE ROUTINE_MULTI( cdname                                                                              &
18      &                    , pt1 , cdna1 , psgn1 , pt2 , cdna2 , psgn2 , pt3 , cdna3 , psgn3 , pt4 , cdna4 , psgn4  &
19      &                    , pt5 , cdna5 , psgn5 , pt6 , cdna6 , psgn6 , pt7 , cdna7 , psgn7 , pt8 , cdna8 , psgn8  &
20      &                    , pt9 , cdna9 , psgn9 , pt10, cdna10, psgn10, pt11, cdna11, psgn11, pt12, cdna12, psgn12 &
21      &                    , pt13, cdna13, psgn13, pt14, cdna14, psgn14, pt15, cdna15, psgn15, pt16, cdna16, psgn16 &
22      &                    , pt17, cdna17, psgn17, pt18, cdna18, psgn18, pt19, cdna19, psgn19, pt20, cdna20, psgn20 &
23      &                    , kfillmode, pfillval, ldsend, ldrecv, khlcom )
24      !!---------------------------------------------------------------------
25      CHARACTER(len=*)   ,                   INTENT(in   ) :: cdname  ! name of the calling subroutine
26      ARRAY_TYPE(:,:,:,:)          , TARGET, INTENT(inout) :: pt1     ! arrays on which the lbc is applied
27      ARRAY_TYPE(:,:,:,:), OPTIONAL, TARGET, INTENT(inout) :: pt2  , pt3  , pt4  , pt5  , pt6  , pt7  , pt8  , pt9  , pt10  , pt11
28      ARRAY_TYPE(:,:,:,:), OPTIONAL, TARGET, INTENT(inout) :: pt12  , pt13  , pt14  , pt15  , pt16  , pt17  , pt18  , pt19  , pt20
29      CHARACTER(len=1)                     , INTENT(in   ) :: cdna1   ! nature of pt2D. array grid-points
30      CHARACTER(len=1)   , OPTIONAL        , INTENT(in   ) :: cdna2, cdna3, cdna4, cdna5, cdna6, cdna7, cdna8, cdna9, cdna10, cdna11
31      CHARACTER(len=1)   , OPTIONAL        , INTENT(in   ) :: cdna12, cdna13, cdna14, cdna15, cdna16, cdna17, cdna18, cdna19, cdna20
32      REAL(wp)                             , INTENT(in   ) :: psgn1   ! sign used across the north fold
33      REAL(wp)           , OPTIONAL        , INTENT(in   ) :: psgn2, psgn3, psgn4, psgn5, psgn6, psgn7, psgn8, psgn9, psgn10, psgn11
34      REAL(wp)           , OPTIONAL        , INTENT(in   ) :: psgn12, psgn13, psgn14, psgn15, psgn16, psgn17, psgn18, psgn19, psgn20
35      INTEGER            , OPTIONAL        , INTENT(in   ) :: kfillmode   ! filling method for halo over land (default = constant)
36      REAL(wp)           , OPTIONAL        , INTENT(in   ) :: pfillval    ! background value (used at closed boundaries)
37      LOGICAL, DIMENSION(4), OPTIONAL      , INTENT(in   ) :: ldsend, ldrecv   ! indicate how communications are to be carried out
38      INTEGER            , OPTIONAL        , INTENT(in   ) :: khlcom           ! number of ranks and rows to be communicated
39      !!
40      INTEGER                          ::   kfld        ! number of elements that will be attributed
41      PTR_TYPE         , DIMENSION(20) ::   ptab_ptr    ! pointer array
42      CHARACTER(len=1) , DIMENSION(20) ::   cdna_ptr    ! nature of ptab_ptr grid-points
43      REAL(wp)         , DIMENSION(20) ::   psgn_ptr    ! sign used across the north fold boundary
44      !!---------------------------------------------------------------------
45      !
46      kfld = 0          ! initial array of pointer size
47      !
48      !                 ! Load the first array
49      CALL ROUTINE_LOAD( pt1, cdna1, psgn1, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
50      !
51      !                 ! Look if more arrays are added
52      IF( PRESENT(psgn2 ) )   CALL ROUTINE_LOAD( pt2 , cdna2 , psgn2 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
53      IF( PRESENT(psgn3 ) )   CALL ROUTINE_LOAD( pt3 , cdna3 , psgn3 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
54      IF( PRESENT(psgn4 ) )   CALL ROUTINE_LOAD( pt4 , cdna4 , psgn4 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
55      IF( PRESENT(psgn5 ) )   CALL ROUTINE_LOAD( pt5 , cdna5 , psgn5 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
56      IF( PRESENT(psgn6 ) )   CALL ROUTINE_LOAD( pt6 , cdna6 , psgn6 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
57      IF( PRESENT(psgn7 ) )   CALL ROUTINE_LOAD( pt7 , cdna7 , psgn7 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
58      IF( PRESENT(psgn8 ) )   CALL ROUTINE_LOAD( pt8 , cdna8 , psgn8 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
59      IF( PRESENT(psgn9 ) )   CALL ROUTINE_LOAD( pt9 , cdna9 , psgn9 , ptab_ptr, cdna_ptr, psgn_ptr, kfld )
60      IF( PRESENT(psgn10) )   CALL ROUTINE_LOAD( pt10, cdna10, psgn10, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
61      IF( PRESENT(psgn11) )   CALL ROUTINE_LOAD( pt11, cdna11, psgn11, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
62      IF( PRESENT(psgn12) )   CALL ROUTINE_LOAD( pt12, cdna12, psgn12, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
63      IF( PRESENT(psgn13) )   CALL ROUTINE_LOAD( pt13, cdna13, psgn13, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
64      IF( PRESENT(psgn14) )   CALL ROUTINE_LOAD( pt14, cdna14, psgn14, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
65      IF( PRESENT(psgn15) )   CALL ROUTINE_LOAD( pt15, cdna15, psgn15, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
66      IF( PRESENT(psgn16) )   CALL ROUTINE_LOAD( pt16, cdna16, psgn16, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
67      IF( PRESENT(psgn17) )   CALL ROUTINE_LOAD( pt17, cdna17, psgn17, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
68      IF( PRESENT(psgn18) )   CALL ROUTINE_LOAD( pt18, cdna18, psgn18, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
69      IF( PRESENT(psgn19) )   CALL ROUTINE_LOAD( pt19, cdna19, psgn19, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
70      IF( PRESENT(psgn20) )   CALL ROUTINE_LOAD( pt20, cdna20, psgn20, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
71      !
72      CALL lbc_lnk_ptr    ( cdname,               ptab_ptr, cdna_ptr, psgn_ptr, kfld, kfillmode, pfillval, ldsend, ldrecv, khlcom )
73      !
74   END SUBROUTINE ROUTINE_MULTI
75
76
77   SUBROUTINE ROUTINE_LOAD( ptab, cdna, psgn, ptab_ptr, cdna_ptr, psgn_ptr, kfld )
78      !!---------------------------------------------------------------------
79      ARRAY_TYPE(:,:,:,:)   , TARGET, INTENT(inout) ::   ptab       ! arrays on which the lbc is applied
80      CHARACTER(len=1)              , INTENT(in   ) ::   cdna       ! nature of pt2d array grid-points
81      REAL(wp)                      , INTENT(in   ) ::   psgn       ! sign used across the north fold boundary
82      PTR_TYPE        , DIMENSION(:), INTENT(inout) ::   ptab_ptr   ! array of pointers
83      CHARACTER(len=1), DIMENSION(:), INTENT(inout) ::   cdna_ptr   ! nature of pt2d_array array grid-points
84      REAL(wp)        , DIMENSION(:), INTENT(inout) ::   psgn_ptr   ! sign used across the north fold boundary
85      INTEGER                       , INTENT(inout) ::   kfld       ! number of elements that has been attributed
86      !!---------------------------------------------------------------------
87      !
88      kfld                    =  kfld + 1
89      ptab_ptr(kfld)%PTR_ptab => ptab
90      cdna_ptr(kfld)          =  cdna
91      psgn_ptr(kfld)          =  psgn
92      !
93   END SUBROUTINE ROUTINE_LOAD
94
95#undef ARRAY_TYPE
96#undef PTR_TYPE
97#undef PTR_ptab
Note: See TracBrowser for help on using the repository browser.