[3186] | 1 | MODULE wrk_nemo |
---|
[3134] | 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE wrk_nemo *** |
---|
| 4 | !! NEMO work space: define and allocate work-space arrays used in |
---|
| 5 | !! all components of NEMO |
---|
| 6 | !!====================================================================== |
---|
| 7 | !! History : 4.0 ! 2011-01 (A Porter) Original code |
---|
| 8 | !!---------------------------------------------------------------------- |
---|
| 9 | |
---|
| 10 | !!---------------------------------------------------------------------- |
---|
| 11 | !! wrk_alloc : get work space arrays |
---|
| 12 | !! wrk_dealloc : release work space arrays |
---|
| 13 | !! |
---|
| 14 | !! 1d arrays: |
---|
| 15 | !! REAL(wp), POINTER, DIMENSION(:) :: arr1, arr2, ... arr10 |
---|
| 16 | !! or |
---|
| 17 | !! INTEGER, POINTER, DIMENSION(:) :: arr1, arr2, ... arr10 |
---|
| 18 | !! ... |
---|
[3144] | 19 | !! CALL wrk_alloc( nx, arr1, arr2, ... arr10, kistart = kistart ) |
---|
[3134] | 20 | !! ... |
---|
[3144] | 21 | !! CALL wrk_dealloc( nx, arr1, arr2, ... arr10, kistart = kistart) |
---|
[3134] | 22 | !! with: |
---|
| 23 | !! - arr*: 1d arrays. real or (not and) integer |
---|
| 24 | !! - nx: size of the 1d arr* arrays |
---|
| 25 | !! - arr2, ..., arr10: optional parameters |
---|
[3144] | 26 | !! - kistart: optional parameter to lower bound of the 1st dimension (default = 1) |
---|
[3134] | 27 | !! |
---|
| 28 | !! 2d arrays: |
---|
| 29 | !! REAL(wp), POINTER, DIMENSION(:,:) :: arr1, arr2, ... arr10 |
---|
| 30 | !! or |
---|
| 31 | !! INTEGER, POINTER, DIMENSION(:,:) :: arr1, arr2, ... arr10 |
---|
| 32 | !! ... |
---|
[3144] | 33 | !! CALL wrk_alloc( nx, ny, arr1, arr2, ... arr10, kistart = kistart, kjstart = kjstart ) |
---|
[3134] | 34 | !! ... |
---|
[3144] | 35 | !! CALL wrk_dealloc( nx, ny, arr1, arr2, ... arr10, kistart = kistart, kjstart = kjstart ) |
---|
[3134] | 36 | !! with: |
---|
| 37 | !! - arr* 2d arrays. real or (not and) integer |
---|
| 38 | !! - nx, ny: size of the 2d arr* arrays |
---|
| 39 | !! - arr2, ..., arr10: optional parameters |
---|
[3144] | 40 | !! - kistart, kjstart: optional parameters to lower bound of the 1st/2nd dimension (default = 1) |
---|
[3134] | 41 | !! |
---|
| 42 | !! 3d arrays: |
---|
| 43 | !! REAL(wp), POINTER, DIMENSION(:,:,:) :: arr1, arr2, ... arr10 |
---|
| 44 | !! or |
---|
| 45 | !! INTEGER, POINTER, DIMENSION(:,:,:) :: arr1, arr2, ... arr10 |
---|
| 46 | !! ... |
---|
[3144] | 47 | !! CALL wrk_alloc( nx, ny, nz, arr1, arr2, ... arr10, kistart = kistart, kjstart = kjstart, kkstart = kkstart ) |
---|
[3134] | 48 | !! ... |
---|
[3144] | 49 | !! CALL wrk_dealloc( nx, ny, nz, arr1, arr2, ... arr10, kistart = kistart, kjstart = kjstart, kkstart = kkstart ) |
---|
[3134] | 50 | !! with: |
---|
| 51 | !! - arr* 3d arrays. real or (not and) integer |
---|
| 52 | !! - nx, ny, nz: size of the 3d arr* arrays |
---|
| 53 | !! - arr2, ..., arr10: optional parameters |
---|
[3144] | 54 | !! - kistart, kjstart, kkstart: optional parameters to lower bound of the 1st/2nd/3rd dimension (default = 1) |
---|
[3134] | 55 | !! |
---|
| 56 | !! 4d arrays: |
---|
| 57 | !! REAL(wp), POINTER, DIMENSION(:,:,:,:) :: arr1, arr2, ... arr10 |
---|
| 58 | !! or |
---|
| 59 | !! INTEGER, POINTER, DIMENSION(:,:,:,:) :: arr1, arr2, ... arr10 |
---|
| 60 | !! ... |
---|
[3144] | 61 | !! CALL wrk_alloc( nx, ny, nz, nl, arr1, arr2, ... arr10, & |
---|
| 62 | !! & kistart = kistart, kjstart = kjstart, kkstart = kkstart, klstart = klstart ) |
---|
[3134] | 63 | !! ... |
---|
[3144] | 64 | !! CALL wrk_dealloc( nx, ny, nz, nl, arr1, arr2, ... arr10, & |
---|
| 65 | !! & kistart = kistart, kjstart = kjstart, kkstart = kkstart, klstart = klstart ) |
---|
[3134] | 66 | !! with: |
---|
| 67 | !! - arr* 3d arrays. real or (not and) integer |
---|
| 68 | !! - nx, ny, nz, nl: size of the 4d arr* arrays |
---|
| 69 | !! - arr2, ..., arr10: optional parameters |
---|
[3144] | 70 | !! - kistart, kjstart, kkstart, klstart: optional parameters to lower bound of the 1st/2nd/3rd/4th dimension (default = 1) |
---|
[3134] | 71 | !! |
---|
| 72 | !!---------------------------------------------------------------------- |
---|
| 73 | USE par_oce ! ocean parameters |
---|
| 74 | |
---|
| 75 | IMPLICIT NONE |
---|
| 76 | PRIVATE |
---|
| 77 | |
---|
| 78 | PUBLIC wrk_alloc, wrk_dealloc, wrk_list |
---|
| 79 | |
---|
| 80 | INTERFACE wrk_alloc |
---|
| 81 | MODULE PROCEDURE wrk_alloc_1dr, wrk_alloc_2dr, wrk_alloc_3dr, wrk_alloc_4dr, & |
---|
| 82 | & wrk_alloc_1di, wrk_alloc_2di, wrk_alloc_3di, wrk_alloc_4di |
---|
| 83 | END INTERFACE |
---|
| 84 | |
---|
| 85 | INTERFACE wrk_dealloc |
---|
| 86 | MODULE PROCEDURE wrk_dealloc_1dr, wrk_dealloc_2dr, wrk_dealloc_3dr, wrk_dealloc_4dr, & |
---|
| 87 | & wrk_dealloc_1di, wrk_dealloc_2di, wrk_dealloc_3di, wrk_dealloc_4di |
---|
| 88 | END INTERFACE |
---|
| 89 | |
---|
| 90 | |
---|
| 91 | INTEGER, PARAMETER :: jparray = 1000 |
---|
| 92 | INTEGER, PARAMETER :: jpmaxdim = 4 |
---|
| 93 | |
---|
| 94 | INTEGER, PARAMETER :: jpnotdefined = 0 |
---|
| 95 | INTEGER, PARAMETER :: jpinteger = 1 |
---|
| 96 | INTEGER, PARAMETER :: jpreal = 2 |
---|
| 97 | |
---|
| 98 | TYPE leaf |
---|
| 99 | LOGICAL :: in_use |
---|
| 100 | INTEGER :: indic |
---|
| 101 | INTEGER , DIMENSION(:) , POINTER :: iwrk1d => NULL() |
---|
| 102 | INTEGER , DIMENSION(:,:) , POINTER :: iwrk2d => NULL() |
---|
| 103 | INTEGER , DIMENSION(:,:,:) , POINTER :: iwrk3d => NULL() |
---|
| 104 | INTEGER , DIMENSION(:,:,:,:), POINTER :: iwrk4d => NULL() |
---|
| 105 | REAL(wp), DIMENSION(:) , POINTER :: zwrk1d => NULL() |
---|
| 106 | REAL(wp), DIMENSION(:,:) , POINTER :: zwrk2d => NULL() |
---|
| 107 | REAL(wp), DIMENSION(:,:,:) , POINTER :: zwrk3d => NULL() |
---|
| 108 | REAL(wp), DIMENSION(:,:,:,:), POINTER :: zwrk4d => NULL() |
---|
| 109 | TYPE (leaf), POINTER :: next => NULL() |
---|
| 110 | TYPE (leaf), POINTER :: prev => NULL() |
---|
| 111 | END TYPE leaf |
---|
| 112 | |
---|
| 113 | TYPE branch |
---|
| 114 | INTEGER :: itype |
---|
[3144] | 115 | INTEGER, DIMENSION(jpmaxdim) :: ishape, istart |
---|
[3134] | 116 | TYPE(leaf), POINTER :: start => NULL() |
---|
| 117 | TYPE(leaf), POINTER :: current => NULL() |
---|
| 118 | END TYPE branch |
---|
| 119 | |
---|
[3150] | 120 | TYPE(branch), SAVE, DIMENSION(jparray) :: tree |
---|
[3134] | 121 | |
---|
| 122 | LOGICAL :: linit = .FALSE. |
---|
| 123 | !!---------------------------------------------------------------------- |
---|
| 124 | !! NEMO/OPA 4.0 , NEMO Consortium (2011) |
---|
[5239] | 125 | !! $Id$ |
---|
[3134] | 126 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
| 127 | !!---------------------------------------------------------------------- |
---|
| 128 | CONTAINS |
---|
| 129 | |
---|
| 130 | SUBROUTINE wrk_list |
---|
| 131 | ! to list 3d arrays in use, to be duplicated for all cases |
---|
| 132 | WRITE(*,*) 'Arrays in use :' |
---|
| 133 | ! CALL listage(tree_3d(1)%s_wrk_3d_start) |
---|
| 134 | WRITE(*,*) '' |
---|
| 135 | |
---|
| 136 | END SUBROUTINE wrk_list |
---|
| 137 | |
---|
| 138 | |
---|
| 139 | RECURSIVE SUBROUTINE listage(ptr) |
---|
| 140 | |
---|
| 141 | TYPE(leaf), POINTER, INTENT(in) :: ptr |
---|
| 142 | ! |
---|
| 143 | IF( ASSOCIATED(ptr%next) ) CALL listage(ptr%next) |
---|
| 144 | WRITE(*,*) ptr%in_use, ptr%indic |
---|
| 145 | |
---|
| 146 | END SUBROUTINE listage |
---|
| 147 | |
---|
| 148 | |
---|
[3144] | 149 | SUBROUTINE wrk_alloc_1dr( kidim, p1d01, p1d02, p1d03, p1d04, p1d05, p1d06, p1d07, p1d08, p1d09, p1d10, kistart ) |
---|
[3134] | 150 | INTEGER , INTENT(in ) :: kidim ! dimensions size |
---|
| 151 | REAL(wp), POINTER, DIMENSION(:), INTENT(inout) :: p1d01 |
---|
| 152 | REAL(wp), POINTER, DIMENSION(:), INTENT(inout), OPTIONAL :: p1d02,p1d03,p1d04,p1d05,p1d06,p1d07,p1d08,p1d09,p1d10 |
---|
[3144] | 153 | INTEGER , INTENT(in ), OPTIONAL :: kistart |
---|
[3134] | 154 | ! |
---|
[3144] | 155 | CALL wrk_alloc_xd( kidim, 0, 0, 0, kistart, 1, 1, 1, & |
---|
| 156 | & p1d01 = p1d01, p1d02 = p1d02, p1d03 = p1d03, p1d04 = p1d04, p1d05 = p1d05, & |
---|
| 157 | & p1d06 = p1d06, p1d07 = p1d07, p1d08 = p1d08, p1d09 = p1d09, p1d10 = p1d10 ) |
---|
[3134] | 158 | ! |
---|
| 159 | END SUBROUTINE wrk_alloc_1dr |
---|
| 160 | |
---|
| 161 | |
---|
[3144] | 162 | SUBROUTINE wrk_alloc_1di( kidim, k1d01, k1d02, k1d03, k1d04, k1d05, k1d06, k1d07, k1d08, k1d09, k1d10, kistart ) |
---|
[3134] | 163 | INTEGER , INTENT(in ) :: kidim ! dimensions size |
---|
| 164 | INTEGER , POINTER, DIMENSION(:), INTENT(inout) :: k1d01 |
---|
| 165 | INTEGER , POINTER, DIMENSION(:), INTENT(inout), OPTIONAL :: k1d02,k1d03,k1d04,k1d05,k1d06,k1d07,k1d08,k1d09,k1d10 |
---|
[3144] | 166 | INTEGER , INTENT(in ), OPTIONAL :: kistart |
---|
[3134] | 167 | ! |
---|
[3144] | 168 | CALL wrk_alloc_xd( kidim, 0, 0, 0, kistart, 1, 1, 1, & |
---|
| 169 | & k1d01 = k1d01, k1d02 = k1d02, k1d03 = k1d03, k1d04 = k1d04, k1d05 = k1d05, & |
---|
| 170 | & k1d06 = k1d06, k1d07 = k1d07, k1d08 = k1d08, k1d09 = k1d09, k1d10 = k1d10 ) |
---|
[3134] | 171 | ! |
---|
| 172 | END SUBROUTINE wrk_alloc_1di |
---|
| 173 | |
---|
| 174 | |
---|
[3144] | 175 | SUBROUTINE wrk_alloc_2dr( kidim, kjdim, p2d01, p2d02, p2d03, p2d04, p2d05, p2d06, p2d07, p2d08, p2d09, p2d10, kistart, kjstart ) |
---|
[3134] | 176 | INTEGER , INTENT(in ) :: kidim, kjdim ! dimensions size |
---|
| 177 | REAL(wp), POINTER, DIMENSION(:,:), INTENT(inout) :: p2d01 |
---|
| 178 | REAL(wp), POINTER, DIMENSION(:,:), INTENT(inout), OPTIONAL :: p2d02,p2d03,p2d04,p2d05,p2d06,p2d07,p2d08,p2d09,p2d10 |
---|
[3144] | 179 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart |
---|
[3134] | 180 | ! |
---|
[3144] | 181 | CALL wrk_alloc_xd( kidim, kjdim, 0, 0, kistart, kjstart, 1, 1, & |
---|
| 182 | & p2d01 = p2d01, p2d02 = p2d02, p2d03 = p2d03, p2d04 = p2d04, p2d05 = p2d05, & |
---|
| 183 | & p2d06 = p2d06, p2d07 = p2d07, p2d08 = p2d08, p2d09 = p2d09, p2d10 = p2d10 ) |
---|
[3134] | 184 | ! |
---|
| 185 | END SUBROUTINE wrk_alloc_2dr |
---|
| 186 | |
---|
| 187 | |
---|
[3144] | 188 | SUBROUTINE wrk_alloc_2di( kidim, kjdim, k2d01, k2d02, k2d03, k2d04, k2d05, k2d06, k2d07, k2d08, k2d09, k2d10, kistart, kjstart ) |
---|
[3134] | 189 | INTEGER , INTENT(in ) :: kidim, kjdim ! dimensions size |
---|
| 190 | INTEGER , POINTER, DIMENSION(:,:), INTENT(inout) :: k2d01 |
---|
| 191 | INTEGER , POINTER, DIMENSION(:,:), INTENT(inout), OPTIONAL :: k2d02,k2d03,k2d04,k2d05,k2d06,k2d07,k2d08,k2d09,k2d10 |
---|
[3144] | 192 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart |
---|
[3134] | 193 | ! |
---|
[3144] | 194 | CALL wrk_alloc_xd( kidim, kjdim, 0, 0, kistart, kjstart, 1, 1, & |
---|
| 195 | & k2d01 = k2d01, k2d02 = k2d02, k2d03 = k2d03, k2d04 = k2d04, k2d05 = k2d05, & |
---|
| 196 | & k2d06 = k2d06, k2d07 = k2d07, k2d08 = k2d08, k2d09 = k2d09, k2d10 = k2d10 ) |
---|
[3134] | 197 | ! |
---|
| 198 | END SUBROUTINE wrk_alloc_2di |
---|
| 199 | |
---|
| 200 | |
---|
[3144] | 201 | SUBROUTINE wrk_alloc_3dr( kidim, kjdim, kkdim, p3d01, p3d02, p3d03, p3d04, p3d05, p3d06, p3d07, p3d08, p3d09, p3d10, & |
---|
| 202 | & kistart, kjstart, kkstart ) |
---|
[3134] | 203 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim ! dimensions size |
---|
| 204 | REAL(wp), POINTER, DIMENSION(:,:,:), INTENT(inout) :: p3d01 |
---|
| 205 | REAL(wp), POINTER, DIMENSION(:,:,:), INTENT(inout), OPTIONAL :: p3d02,p3d03,p3d04,p3d05,p3d06,p3d07,p3d08,p3d09,p3d10 |
---|
[3144] | 206 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart |
---|
[3134] | 207 | ! |
---|
[3144] | 208 | CALL wrk_alloc_xd( kidim, kjdim, kkdim, 0, kistart, kjstart, kkstart, 1, & |
---|
| 209 | & p3d01 = p3d01, p3d02 = p3d02, p3d03 = p3d03, p3d04 = p3d04, p3d05 = p3d05, & |
---|
| 210 | & p3d06 = p3d06, p3d07 = p3d07, p3d08 = p3d08, p3d09 = p3d09, p3d10 = p3d10 ) |
---|
[3134] | 211 | ! |
---|
| 212 | END SUBROUTINE wrk_alloc_3dr |
---|
| 213 | |
---|
| 214 | |
---|
[3144] | 215 | SUBROUTINE wrk_alloc_3di( kidim, kjdim, kkdim, k3d01, k3d02, k3d03, k3d04, k3d05, k3d06, k3d07, k3d08, k3d09, k3d10, & |
---|
| 216 | & kistart, kjstart, kkstart ) |
---|
[3134] | 217 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim ! dimensions size |
---|
| 218 | INTEGER , POINTER, DIMENSION(:,:,:), INTENT(inout) :: k3d01 |
---|
| 219 | INTEGER , POINTER, DIMENSION(:,:,:), INTENT(inout), OPTIONAL :: k3d02,k3d03,k3d04,k3d05,k3d06,k3d07,k3d08,k3d09,k3d10 |
---|
[3144] | 220 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart |
---|
[3134] | 221 | ! |
---|
[3144] | 222 | CALL wrk_alloc_xd( kidim, kjdim, kkdim, 0, kistart, kjstart, kkstart, 1, & |
---|
| 223 | & k3d01 = k3d01, k3d02 = k3d02, k3d03 = k3d03, k3d04 = k3d04, k3d05 = k3d05, & |
---|
| 224 | & k3d06 = k3d06, k3d07 = k3d07, k3d08 = k3d08, k3d09 = k3d09, k3d10 = k3d10 ) |
---|
[3134] | 225 | ! |
---|
| 226 | END SUBROUTINE wrk_alloc_3di |
---|
| 227 | |
---|
| 228 | |
---|
[3144] | 229 | SUBROUTINE wrk_alloc_4dr( kidim, kjdim, kkdim, kldim, p4d01, p4d02, p4d03, p4d04, p4d05, p4d06, p4d07, p4d08, p4d09, p4d10, & |
---|
| 230 | & kistart, kjstart, kkstart, klstart ) |
---|
[3134] | 231 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim, kldim ! dimensions size |
---|
| 232 | REAL(wp), POINTER, DIMENSION(:,:,:,:), INTENT(inout) :: p4d01 |
---|
| 233 | REAL(wp), POINTER, DIMENSION(:,:,:,:), INTENT(inout), OPTIONAL :: p4d02,p4d03,p4d04,p4d05,p4d06,p4d07,p4d08,p4d09,p4d10 |
---|
[3144] | 234 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart, klstart |
---|
[3134] | 235 | ! |
---|
[3144] | 236 | CALL wrk_alloc_xd( kidim, kjdim, kkdim, kldim, kistart, kjstart, kkstart, klstart, & |
---|
| 237 | & p4d01 = p4d01, p4d02 = p4d02, p4d03 = p4d03, p4d04 = p4d04, p4d05 = p4d05, & |
---|
| 238 | & p4d06 = p4d06, p4d07 = p4d07, p4d08 = p4d08, p4d09 = p4d09, p4d10 = p4d10 ) |
---|
[3134] | 239 | ! |
---|
| 240 | END SUBROUTINE wrk_alloc_4dr |
---|
| 241 | |
---|
| 242 | |
---|
[3144] | 243 | SUBROUTINE wrk_alloc_4di( kidim, kjdim, kkdim, kldim, k4d01, k4d02, k4d03, k4d04, k4d05, k4d06, k4d07, k4d08, k4d09, k4d10, & |
---|
| 244 | & kistart, kjstart, kkstart, klstart ) |
---|
[3134] | 245 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim, kldim ! dimensions size |
---|
| 246 | INTEGER , POINTER, DIMENSION(:,:,:,:), INTENT(inout) :: k4d01 |
---|
| 247 | INTEGER , POINTER, DIMENSION(:,:,:,:), INTENT(inout), OPTIONAL :: k4d02,k4d03,k4d04,k4d05,k4d06,k4d07,k4d08,k4d09,k4d10 |
---|
[3144] | 248 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart, klstart |
---|
[3134] | 249 | ! |
---|
[3144] | 250 | CALL wrk_alloc_xd( kidim, kjdim, kkdim, kldim, kistart, kjstart, kkstart, klstart, & |
---|
| 251 | & k4d01 = k4d01, k4d02 = k4d02, k4d03 = k4d03, k4d04 = k4d04, k4d05 = k4d05, & |
---|
| 252 | & k4d06 = k4d06, k4d07 = k4d07, k4d08 = k4d08, k4d09 = k4d09, k4d10 = k4d10 ) |
---|
[3134] | 253 | ! |
---|
| 254 | END SUBROUTINE wrk_alloc_4di |
---|
| 255 | |
---|
| 256 | |
---|
[3144] | 257 | SUBROUTINE wrk_dealloc_1dr( kidim, p1d01, p1d02, p1d03, p1d04, p1d05, p1d06, p1d07, p1d08, p1d09, p1d10, kistart ) |
---|
[3134] | 258 | INTEGER , INTENT(in ) :: kidim ! dimensions size |
---|
| 259 | REAL(wp), POINTER, DIMENSION(:), INTENT(inout) :: p1d01 |
---|
| 260 | REAL(wp), POINTER, DIMENSION(:), INTENT(inout), OPTIONAL :: p1d02,p1d03,p1d04,p1d05,p1d06,p1d07,p1d08,p1d09,p1d10 |
---|
[3144] | 261 | INTEGER , INTENT(in ), OPTIONAL :: kistart |
---|
[3134] | 262 | ! |
---|
| 263 | INTEGER :: icnt, jn |
---|
| 264 | icnt = 1 + COUNT( (/ PRESENT(p1d02),PRESENT(p1d03),PRESENT(p1d04),PRESENT(p1d05), & |
---|
| 265 | & PRESENT(p1d06),PRESENT(p1d07),PRESENT(p1d08),PRESENT(p1d09),PRESENT(p1d10) /) ) |
---|
[3144] | 266 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpreal, kidim, 0, 0, 0, kistart, 1, 1, 1) ; END DO |
---|
[3134] | 267 | ! |
---|
| 268 | END SUBROUTINE wrk_dealloc_1dr |
---|
| 269 | |
---|
| 270 | |
---|
[3144] | 271 | SUBROUTINE wrk_dealloc_1di( kidim, k1d01, k1d02, k1d03, k1d04, k1d05, k1d06, k1d07, k1d08, k1d09, k1d10, kistart ) |
---|
[3134] | 272 | INTEGER , INTENT(in ) :: kidim ! dimensions size |
---|
| 273 | INTEGER , POINTER, DIMENSION(:), INTENT(inout) :: k1d01 |
---|
| 274 | INTEGER , POINTER, DIMENSION(:), INTENT(inout), OPTIONAL :: k1d02,k1d03,k1d04,k1d05,k1d06,k1d07,k1d08,k1d09,k1d10 |
---|
[3144] | 275 | INTEGER , INTENT(in ), OPTIONAL :: kistart |
---|
[3134] | 276 | ! |
---|
| 277 | INTEGER :: icnt, jn |
---|
| 278 | icnt = 1 + COUNT( (/ PRESENT(k1d02),PRESENT(k1d03),PRESENT(k1d04),PRESENT(k1d05), & |
---|
| 279 | & PRESENT(k1d06),PRESENT(k1d07),PRESENT(k1d08),PRESENT(k1d09),PRESENT(k1d10) /) ) |
---|
[3144] | 280 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpinteger, kidim, 0, 0, 0, kistart, 1, 1, 1 ) ; END DO |
---|
[3134] | 281 | ! |
---|
| 282 | END SUBROUTINE wrk_dealloc_1di |
---|
| 283 | |
---|
| 284 | |
---|
[3144] | 285 | SUBROUTINE wrk_dealloc_2dr( kidim, kjdim, p2d01, p2d02, p2d03, p2d04, p2d05, p2d06, p2d07, p2d08, p2d09, p2d10, kistart,kjstart ) |
---|
[3134] | 286 | INTEGER , INTENT(in ) :: kidim, kjdim ! dimensions size |
---|
| 287 | REAL(wp), POINTER, DIMENSION(:,:), INTENT(inout) :: p2d01 |
---|
| 288 | REAL(wp), POINTER, DIMENSION(:,:), INTENT(inout), OPTIONAL :: p2d02,p2d03,p2d04,p2d05,p2d06,p2d07,p2d08,p2d09,p2d10 |
---|
[3144] | 289 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart |
---|
[3134] | 290 | ! |
---|
| 291 | INTEGER :: icnt, jn |
---|
| 292 | icnt = 1 + COUNT( (/ PRESENT(p2d02),PRESENT(p2d03),PRESENT(p2d04),PRESENT(p2d05), & |
---|
| 293 | & PRESENT(p2d06),PRESENT(p2d07),PRESENT(p2d08),PRESENT(p2d09),PRESENT(p2d10) /) ) |
---|
[3144] | 294 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpreal, kidim, kjdim, 0, 0, kistart, kjstart, 1, 1 ) ; END DO |
---|
[3134] | 295 | ! |
---|
| 296 | END SUBROUTINE wrk_dealloc_2dr |
---|
| 297 | |
---|
| 298 | |
---|
[3144] | 299 | SUBROUTINE wrk_dealloc_2di( kidim, kjdim, k2d01, k2d02, k2d03, k2d04, k2d05, k2d06, k2d07, k2d08, k2d09, k2d10, kistart,kjstart ) |
---|
[3134] | 300 | INTEGER , INTENT(in ) :: kidim, kjdim ! dimensions size |
---|
| 301 | INTEGER , POINTER, DIMENSION(:,:), INTENT(inout) :: k2d01 |
---|
| 302 | INTEGER , POINTER, DIMENSION(:,:), INTENT(inout), OPTIONAL :: k2d02,k2d03,k2d04,k2d05,k2d06,k2d07,k2d08,k2d09,k2d10 |
---|
[3144] | 303 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart |
---|
[3134] | 304 | ! |
---|
| 305 | INTEGER :: icnt, jn |
---|
| 306 | icnt = 1 + COUNT( (/ PRESENT(k2d02),PRESENT(k2d03),PRESENT(k2d04),PRESENT(k2d05), & |
---|
| 307 | & PRESENT(k2d06),PRESENT(k2d07),PRESENT(k2d08),PRESENT(k2d09),PRESENT(k2d10) /) ) |
---|
[3144] | 308 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpinteger, kidim, kjdim, 0, 0, kistart, kjstart, 1, 1 ) ; END DO |
---|
[3134] | 309 | ! |
---|
| 310 | END SUBROUTINE wrk_dealloc_2di |
---|
| 311 | |
---|
| 312 | |
---|
[3144] | 313 | SUBROUTINE wrk_dealloc_3dr( kidim, kjdim, kkdim, p3d01, p3d02, p3d03, p3d04, p3d05, p3d06, p3d07, p3d08, p3d09, p3d10, & |
---|
| 314 | & kistart, kjstart, kkstart ) |
---|
[3134] | 315 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim ! dimensions size |
---|
| 316 | REAL(wp), POINTER, DIMENSION(:,:,:), INTENT(inout) :: p3d01 |
---|
| 317 | REAL(wp), POINTER, DIMENSION(:,:,:), INTENT(inout), OPTIONAL :: p3d02,p3d03,p3d04,p3d05,p3d06,p3d07,p3d08,p3d09,p3d10 |
---|
[3144] | 318 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart |
---|
[3134] | 319 | ! |
---|
| 320 | INTEGER :: icnt, jn |
---|
| 321 | icnt = 1 + COUNT( (/ PRESENT(p3d02),PRESENT(p3d03),PRESENT(p3d04),PRESENT(p3d05), & |
---|
| 322 | & PRESENT(p3d06),PRESENT(p3d07),PRESENT(p3d08),PRESENT(p3d09),PRESENT(p3d10) /) ) |
---|
[3144] | 323 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpreal, kidim, kjdim, kkdim, 0, kistart, kjstart, kkstart, 1 ) ; END DO |
---|
[3134] | 324 | ! |
---|
| 325 | END SUBROUTINE wrk_dealloc_3dr |
---|
| 326 | |
---|
| 327 | |
---|
[3144] | 328 | SUBROUTINE wrk_dealloc_3di( kidim, kjdim, kkdim, k3d01, k3d02, k3d03, k3d04, k3d05, k3d06, k3d07, k3d08, k3d09, k3d10, & |
---|
| 329 | & kistart, kjstart, kkstart ) |
---|
[3134] | 330 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim ! dimensions size |
---|
| 331 | INTEGER , POINTER, DIMENSION(:,:,:), INTENT(inout) :: k3d01 |
---|
| 332 | INTEGER , POINTER, DIMENSION(:,:,:), INTENT(inout), OPTIONAL :: k3d02,k3d03,k3d04,k3d05,k3d06,k3d07,k3d08,k3d09,k3d10 |
---|
[3144] | 333 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart |
---|
[3134] | 334 | ! |
---|
| 335 | INTEGER :: icnt, jn |
---|
| 336 | icnt = 1 + COUNT( (/ PRESENT(k3d02),PRESENT(k3d03),PRESENT(k3d04),PRESENT(k3d05), & |
---|
| 337 | & PRESENT(k3d06),PRESENT(k3d07),PRESENT(k3d08),PRESENT(k3d09),PRESENT(k3d10) /) ) |
---|
[3144] | 338 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpinteger, kidim, kjdim, kkdim, 0, kistart, kjstart, kkstart, 1 ) ; END DO |
---|
[3134] | 339 | ! |
---|
| 340 | END SUBROUTINE wrk_dealloc_3di |
---|
| 341 | |
---|
| 342 | |
---|
[3144] | 343 | SUBROUTINE wrk_dealloc_4dr( kidim, kjdim, kkdim, kldim, p4d01, p4d02, p4d03, p4d04, p4d05, p4d06, p4d07, p4d08, p4d09, p4d10, & |
---|
| 344 | & kistart, kjstart, kkstart, klstart ) |
---|
[3134] | 345 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim, kldim ! dimensions size |
---|
| 346 | REAL(wp), POINTER, DIMENSION(:,:,:,:), INTENT(inout) :: p4d01 |
---|
| 347 | REAL(wp), POINTER, DIMENSION(:,:,:,:), INTENT(inout), OPTIONAL :: p4d02,p4d03,p4d04,p4d05,p4d06,p4d07,p4d08,p4d09,p4d10 |
---|
[3144] | 348 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart, klstart |
---|
[3134] | 349 | ! |
---|
| 350 | INTEGER :: icnt, jn |
---|
| 351 | icnt = 1 + COUNT( (/ PRESENT(p4d02),PRESENT(p4d03),PRESENT(p4d04),PRESENT(p4d05), & |
---|
| 352 | & PRESENT(p4d06),PRESENT(p4d07),PRESENT(p4d08),PRESENT(p4d09),PRESENT(p4d10) /) ) |
---|
[3144] | 353 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpreal, kidim, kjdim, kkdim, kldim, kistart, kjstart, kkstart, klstart ) ; END DO |
---|
[3134] | 354 | ! |
---|
| 355 | END SUBROUTINE wrk_dealloc_4dr |
---|
| 356 | |
---|
| 357 | |
---|
[3144] | 358 | SUBROUTINE wrk_dealloc_4di( kidim, kjdim, kkdim, kldim, k4d01, k4d02, k4d03, k4d04, k4d05, k4d06, k4d07, k4d08, k4d09, k4d10, & |
---|
| 359 | & kistart, kjstart, kkstart, klstart ) |
---|
[3134] | 360 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim, kldim ! dimensions size |
---|
| 361 | INTEGER , POINTER, DIMENSION(:,:,:,:), INTENT(inout) :: k4d01 |
---|
| 362 | INTEGER , POINTER, DIMENSION(:,:,:,:), INTENT(inout), OPTIONAL :: k4d02,k4d03,k4d04,k4d05,k4d06,k4d07,k4d08,k4d09,k4d10 |
---|
[3144] | 363 | INTEGER , INTENT(in ), OPTIONAL :: kistart, kjstart, kkstart, klstart |
---|
[3134] | 364 | ! |
---|
| 365 | INTEGER :: icnt, jn |
---|
| 366 | icnt = 1 + COUNT( (/ PRESENT(k4d02),PRESENT(k4d03),PRESENT(k4d04),PRESENT(k4d05), & |
---|
| 367 | & PRESENT(k4d06),PRESENT(k4d07),PRESENT(k4d08),PRESENT(k4d09),PRESENT(k4d10) /) ) |
---|
[3144] | 368 | DO jn = 1, icnt ; CALL wrk_deallocbase( jpinteger, kidim, kjdim, kkdim, kldim, kistart, kjstart, kkstart, klstart ) ; END DO |
---|
[3134] | 369 | ! |
---|
| 370 | END SUBROUTINE wrk_dealloc_4di |
---|
| 371 | |
---|
| 372 | |
---|
| 373 | SUBROUTINE wrk_alloc_xd( kidim, kjdim, kkdim, kldim, & |
---|
[3144] | 374 | & kisrt, kjsrt, kksrt, klsrt, & |
---|
[3134] | 375 | & k1d01, k1d02, k1d03, k1d04, k1d05, k1d06, k1d07, k1d08, k1d09, k1d10, & |
---|
| 376 | & k2d01, k2d02, k2d03, k2d04, k2d05, k2d06, k2d07, k2d08, k2d09, k2d10, & |
---|
| 377 | & k3d01, k3d02, k3d03, k3d04, k3d05, k3d06, k3d07, k3d08, k3d09, k3d10, & |
---|
| 378 | & k4d01, k4d02, k4d03, k4d04, k4d05, k4d06, k4d07, k4d08, k4d09, k4d10, & |
---|
| 379 | & p1d01, p1d02, p1d03, p1d04, p1d05, p1d06, p1d07, p1d08, p1d09, p1d10, & |
---|
| 380 | & p2d01, p2d02, p2d03, p2d04, p2d05, p2d06, p2d07, p2d08, p2d09, p2d10, & |
---|
| 381 | & p3d01, p3d02, p3d03, p3d04, p3d05, p3d06, p3d07, p3d08, p3d09, p3d10, & |
---|
| 382 | & p4d01, p4d02, p4d03, p4d04, p4d05, p4d06, p4d07, p4d08, p4d09, p4d10 ) |
---|
| 383 | INTEGER ,INTENT(in ) :: kidim, kjdim, kkdim, kldim ! dimensions size |
---|
[3144] | 384 | INTEGER ,INTENT(in ),OPTIONAL:: kisrt, kjsrt, kksrt, klsrt |
---|
[3134] | 385 | INTEGER , POINTER, DIMENSION(: ),INTENT(inout),OPTIONAL:: k1d01,k1d02,k1d03,k1d04,k1d05,k1d06,k1d07,k1d08,k1d09,k1d10 |
---|
| 386 | INTEGER , POINTER, DIMENSION(:,: ),INTENT(inout),OPTIONAL:: k2d01,k2d02,k2d03,k2d04,k2d05,k2d06,k2d07,k2d08,k2d09,k2d10 |
---|
| 387 | INTEGER , POINTER, DIMENSION(:,:,: ),INTENT(inout),OPTIONAL:: k3d01,k3d02,k3d03,k3d04,k3d05,k3d06,k3d07,k3d08,k3d09,k3d10 |
---|
| 388 | INTEGER , POINTER, DIMENSION(:,:,:,:),INTENT(inout),OPTIONAL:: k4d01,k4d02,k4d03,k4d04,k4d05,k4d06,k4d07,k4d08,k4d09,k4d10 |
---|
| 389 | REAL(wp), POINTER, DIMENSION(: ),INTENT(inout),OPTIONAL:: p1d01,p1d02,p1d03,p1d04,p1d05,p1d06,p1d07,p1d08,p1d09,p1d10 |
---|
| 390 | REAL(wp), POINTER, DIMENSION(:,: ),INTENT(inout),OPTIONAL:: p2d01,p2d02,p2d03,p2d04,p2d05,p2d06,p2d07,p2d08,p2d09,p2d10 |
---|
| 391 | REAL(wp), POINTER, DIMENSION(:,:,: ),INTENT(inout),OPTIONAL:: p3d01,p3d02,p3d03,p3d04,p3d05,p3d06,p3d07,p3d08,p3d09,p3d10 |
---|
| 392 | REAL(wp), POINTER, DIMENSION(:,:,:,:),INTENT(inout),OPTIONAL:: p4d01,p4d02,p4d03,p4d04,p4d05,p4d06,p4d07,p4d08,p4d09,p4d10 |
---|
| 393 | ! |
---|
| 394 | LOGICAL :: llpres |
---|
[3144] | 395 | INTEGER :: jn, iisrt, ijsrt, iksrt, ilsrt |
---|
[3134] | 396 | ! |
---|
| 397 | IF( .NOT. linit ) THEN |
---|
| 398 | tree(:)%itype = jpnotdefined |
---|
[3144] | 399 | DO jn = 1, jparray ; tree(jn)%ishape(:) = 0 ; tree(jn)%istart(:) = 0 ; END DO |
---|
[3134] | 400 | linit = .TRUE. |
---|
| 401 | ENDIF |
---|
| 402 | |
---|
[3144] | 403 | IF( PRESENT(kisrt) ) THEN ; iisrt = kisrt ; ELSE ; iisrt = 1 ; ENDIF |
---|
| 404 | IF( PRESENT(kjsrt) ) THEN ; ijsrt = kjsrt ; ELSE ; ijsrt = 1 ; ENDIF |
---|
| 405 | IF( PRESENT(kksrt) ) THEN ; iksrt = kksrt ; ELSE ; iksrt = 1 ; ENDIF |
---|
| 406 | IF( PRESENT(klsrt) ) THEN ; ilsrt = klsrt ; ELSE ; ilsrt = 1 ; ENDIF |
---|
| 407 | |
---|
[3134] | 408 | llpres = PRESENT(k1d01) .OR. PRESENT(k2d01) .OR. PRESENT(k3d01) .OR. PRESENT(k4d01) & |
---|
| 409 | & .OR. PRESENT(p1d01) .OR. PRESENT(p2d01) .OR. PRESENT(p3d01) .OR. PRESENT(p4d01) |
---|
[3144] | 410 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 411 | & k1d01, k2d01, k3d01, k4d01, p1d01, p2d01, p3d01, p4d01 ) |
---|
[3134] | 412 | llpres = PRESENT(k1d02) .OR. PRESENT(k2d02) .OR. PRESENT(k3d02) .OR. PRESENT(k4d02) & |
---|
| 413 | & .OR. PRESENT(p1d02) .OR. PRESENT(p2d02) .OR. PRESENT(p3d02) .OR. PRESENT(p4d02) |
---|
[3144] | 414 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 415 | & k1d02, k2d02, k3d02, k4d02, p1d02, p2d02, p3d02, p4d02 ) |
---|
[3134] | 416 | llpres = PRESENT(k1d03) .OR. PRESENT(k2d03) .OR. PRESENT(k3d03) .OR. PRESENT(k4d03) & |
---|
| 417 | & .OR. PRESENT(p1d03) .OR. PRESENT(p2d03) .OR. PRESENT(p3d03) .OR. PRESENT(p4d03) |
---|
[3144] | 418 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 419 | & k1d03, k2d03, k3d03, k4d03, p1d03, p2d03, p3d03, p4d03 ) |
---|
[3134] | 420 | llpres = PRESENT(k1d04) .OR. PRESENT(k2d04) .OR. PRESENT(k3d04) .OR. PRESENT(k4d04) & |
---|
| 421 | & .OR. PRESENT(p1d04) .OR. PRESENT(p2d04) .OR. PRESENT(p3d04) .OR. PRESENT(p4d04) |
---|
[3144] | 422 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 423 | & k1d04, k2d04, k3d04, k4d04, p1d04, p2d04, p3d04, p4d04 ) |
---|
[3134] | 424 | llpres = PRESENT(k1d05) .OR. PRESENT(k2d05) .OR. PRESENT(k3d05) .OR. PRESENT(k4d05) & |
---|
| 425 | & .OR. PRESENT(p1d05) .OR. PRESENT(p2d05) .OR. PRESENT(p3d05) .OR. PRESENT(p4d05) |
---|
[3144] | 426 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 427 | & k1d05, k2d05, k3d05, k4d05, p1d05, p2d05, p3d05, p4d05 ) |
---|
[3134] | 428 | llpres = PRESENT(k1d06) .OR. PRESENT(k2d06) .OR. PRESENT(k3d06) .OR. PRESENT(k4d06) & |
---|
| 429 | & .OR. PRESENT(p1d06) .OR. PRESENT(p2d06) .OR. PRESENT(p3d06) .OR. PRESENT(p4d06) |
---|
[3144] | 430 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 431 | & k1d06, k2d06, k3d06, k4d06, p1d06, p2d06, p3d06, p4d06 ) |
---|
[3134] | 432 | llpres = PRESENT(k1d07) .OR. PRESENT(k2d07) .OR. PRESENT(k3d07) .OR. PRESENT(k4d07) & |
---|
| 433 | & .OR. PRESENT(p1d07) .OR. PRESENT(p2d07) .OR. PRESENT(p3d07) .OR. PRESENT(p4d07) |
---|
[3144] | 434 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 435 | & k1d07, k2d07, k3d07, k4d07, p1d07, p2d07, p3d07, p4d07 ) |
---|
[3134] | 436 | llpres = PRESENT(k1d08) .OR. PRESENT(k2d08) .OR. PRESENT(k3d08) .OR. PRESENT(k4d08) & |
---|
| 437 | & .OR. PRESENT(p1d08) .OR. PRESENT(p2d08) .OR. PRESENT(p3d08) .OR. PRESENT(p4d08) |
---|
[3144] | 438 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 439 | & k1d08, k2d08, k3d08, k4d08, p1d08, p2d08, p3d08, p4d08 ) |
---|
[3134] | 440 | llpres = PRESENT(k1d09) .OR. PRESENT(k2d09) .OR. PRESENT(k3d09) .OR. PRESENT(k4d09) & |
---|
| 441 | & .OR. PRESENT(p1d09) .OR. PRESENT(p2d09) .OR. PRESENT(p3d09) .OR. PRESENT(p4d09) |
---|
[3144] | 442 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 443 | & k1d09, k2d09, k3d09, k4d09, p1d09, p2d09, p3d09, p4d09 ) |
---|
[3134] | 444 | llpres = PRESENT(k1d10) .OR. PRESENT(k2d10) .OR. PRESENT(k3d10) .OR. PRESENT(k4d10) & |
---|
| 445 | & .OR. PRESENT(p1d10) .OR. PRESENT(p2d10) .OR. PRESENT(p3d10) .OR. PRESENT(p4d10) |
---|
[3144] | 446 | IF( llpres ) CALL wrk_allocbase( kidim, kjdim, kkdim, kldim, iisrt, ijsrt, iksrt, ilsrt, & |
---|
| 447 | & k1d10, k2d10, k3d10, k4d10, p1d10, p2d10, p3d10, p4d10 ) |
---|
[3134] | 448 | |
---|
| 449 | END SUBROUTINE wrk_alloc_xd |
---|
| 450 | |
---|
| 451 | |
---|
[3144] | 452 | SUBROUTINE wrk_allocbase( kidim , kjdim , kkdim , kldim , kisrt , kjsrt , kksrt , klsrt , & |
---|
| 453 | & kwrk1d, kwrk2d, kwrk3d, kwrk4d, pwrk1d, pwrk2d, pwrk3d, pwrk4d ) |
---|
[3134] | 454 | INTEGER , INTENT(in ) :: kidim, kjdim, kkdim, kldim |
---|
[3144] | 455 | INTEGER , INTENT(in ) :: kisrt, kjsrt, kksrt, klsrt |
---|
[3134] | 456 | INTEGER , POINTER, DIMENSION(:) , INTENT(inout), OPTIONAL :: kwrk1d |
---|
| 457 | INTEGER , POINTER, DIMENSION(:,:) , INTENT(inout), OPTIONAL :: kwrk2d |
---|
| 458 | INTEGER , POINTER, DIMENSION(:,:,:) , INTENT(inout), OPTIONAL :: kwrk3d |
---|
| 459 | INTEGER , POINTER, DIMENSION(:,:,:,:), INTENT(inout), OPTIONAL :: kwrk4d |
---|
| 460 | REAL(wp), POINTER, DIMENSION(:) , INTENT(inout), OPTIONAL :: pwrk1d |
---|
| 461 | REAL(wp), POINTER, DIMENSION(:,:) , INTENT(inout), OPTIONAL :: pwrk2d |
---|
| 462 | REAL(wp), POINTER, DIMENSION(:,:,:) , INTENT(inout), OPTIONAL :: pwrk3d |
---|
| 463 | REAL(wp), POINTER, DIMENSION(:,:,:,:), INTENT(inout), OPTIONAL :: pwrk4d |
---|
| 464 | ! |
---|
[3144] | 465 | INTEGER, DIMENSION(jpmaxdim) :: ishape, isrt, iend |
---|
[3134] | 466 | INTEGER :: itype |
---|
| 467 | INTEGER :: ii |
---|
| 468 | |
---|
| 469 | ! define the shape to be given to the work array |
---|
| 470 | ishape(:) = (/ kidim, kjdim, kkdim, kldim /) |
---|
[3144] | 471 | ! define the starting index of the dimension shape to be given to the work array |
---|
| 472 | isrt (:) = (/ kisrt, kjsrt, kksrt, klsrt /) |
---|
| 473 | iend (:) = ishape(:) + isrt(:) - 1 |
---|
[3134] | 474 | |
---|
| 475 | ! is it integer or real array? |
---|
| 476 | IF( PRESENT(kwrk1d) .OR. PRESENT(kwrk2d) .OR. PRESENT(kwrk3d) .OR. PRESENT(kwrk4d) ) itype = jpinteger |
---|
| 477 | IF( PRESENT(pwrk1d) .OR. PRESENT(pwrk2d) .OR. PRESENT(pwrk3d) .OR. PRESENT(pwrk4d) ) itype = jpreal |
---|
| 478 | |
---|
[3144] | 479 | ! find the branch with the matching shape, staring index and type or get the first "free" branch |
---|
[3134] | 480 | ii = 1 |
---|
[3144] | 481 | DO WHILE( ( ANY( tree(ii)%ishape /= ishape ) .OR. ANY( tree(ii)%istart /= isrt ) .OR. tree(ii)%itype /= itype ) & |
---|
| 482 | & .AND. SUM( tree(ii)%ishape ) /= 0 ) |
---|
[3134] | 483 | ii = ii + 1 |
---|
| 484 | IF (ii > jparray) STOP ! increase the value of jparray (should not be needed as already very big!) |
---|
| 485 | END DO |
---|
| 486 | |
---|
| 487 | IF( SUM( tree(ii)%ishape ) == 0 ) THEN ! create a new branch |
---|
| 488 | tree(ii)%itype = itype ! define the type of this branch |
---|
[3144] | 489 | tree(ii)%ishape(:) = ishape(:) ! define the shape of this branch |
---|
| 490 | tree(ii)%istart(:) = isrt(:) ! define the lower bounds of this branch |
---|
[3134] | 491 | ALLOCATE( tree(ii)%start ) ! allocate its start |
---|
| 492 | ALLOCATE( tree(ii)%current) ! allocate the current leaf (the first leaf) |
---|
| 493 | |
---|
| 494 | tree(ii)%start%in_use = .FALSE. ! Never use the start as work array |
---|
| 495 | tree(ii)%start%indic = 0 |
---|
| 496 | tree(ii)%start%prev => NULL() ! nothing before the start |
---|
| 497 | tree(ii)%start%next => tree(ii)%current ! first leaf link to the start |
---|
| 498 | |
---|
| 499 | tree(ii)%current%in_use = .FALSE. ! first leaf is not yet used |
---|
| 500 | tree(ii)%current%indic = 1 ! first leaf |
---|
| 501 | tree(ii)%current%prev => tree(ii)%start ! previous leaf is the start |
---|
| 502 | tree(ii)%current%next => NULL() ! next leaf is not yet defined |
---|
| 503 | ! allocate the array of the first leaf |
---|
[3144] | 504 | IF( PRESENT(kwrk1d) ) ALLOCATE( tree(ii)%current%iwrk1d(isrt(1):iend(1) ) ) |
---|
| 505 | IF( PRESENT(kwrk2d) ) ALLOCATE( tree(ii)%current%iwrk2d(isrt(1):iend(1),isrt(2):iend(2) ) ) |
---|
| 506 | IF( PRESENT(kwrk3d) ) ALLOCATE( tree(ii)%current%iwrk3d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3) ) ) |
---|
| 507 | IF( PRESENT(kwrk4d) ) ALLOCATE( tree(ii)%current%iwrk4d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3),isrt(4):iend(4)) ) |
---|
| 508 | IF( PRESENT(pwrk1d) ) ALLOCATE( tree(ii)%current%zwrk1d(isrt(1):iend(1) ) ) |
---|
| 509 | IF( PRESENT(pwrk2d) ) ALLOCATE( tree(ii)%current%zwrk2d(isrt(1):iend(1),isrt(2):iend(2) ) ) |
---|
| 510 | IF( PRESENT(pwrk3d) ) ALLOCATE( tree(ii)%current%zwrk3d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3) ) ) |
---|
| 511 | IF( PRESENT(pwrk4d) ) ALLOCATE( tree(ii)%current%zwrk4d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3),isrt(4):iend(4)) ) |
---|
[3134] | 512 | |
---|
| 513 | ELSE IF( .NOT. ASSOCIATED(tree(ii)%current%next) ) THEN ! all leafs used -> define a new one |
---|
| 514 | ALLOCATE( tree(ii)%current%next ) ! allocate the new leaf |
---|
| 515 | tree(ii)%current%next%in_use = .FALSE. ! this leaf is not yet used |
---|
| 516 | tree(ii)%current%next%indic = tree(ii)%current%indic + 1 ! number of this leaf |
---|
| 517 | tree(ii)%current%next%prev => tree(ii)%current ! previous leaf of the new leaf is the current leaf |
---|
| 518 | tree(ii)%current%next%next => NULL() ! next leaf is not yet defined |
---|
| 519 | |
---|
| 520 | tree(ii)%current => tree(ii)%current%next ! the current leaf becomes the new one |
---|
| 521 | |
---|
| 522 | ! allocate the array of the new leaf |
---|
[3144] | 523 | IF( PRESENT(kwrk1d) ) ALLOCATE( tree(ii)%current%iwrk1d(isrt(1):iend(1) ) ) |
---|
| 524 | IF( PRESENT(kwrk2d) ) ALLOCATE( tree(ii)%current%iwrk2d(isrt(1):iend(1),isrt(2):iend(2) ) ) |
---|
| 525 | IF( PRESENT(kwrk3d) ) ALLOCATE( tree(ii)%current%iwrk3d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3) ) ) |
---|
| 526 | IF( PRESENT(kwrk4d) ) ALLOCATE( tree(ii)%current%iwrk4d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3),isrt(4):iend(4)) ) |
---|
| 527 | IF( PRESENT(pwrk1d) ) ALLOCATE( tree(ii)%current%zwrk1d(isrt(1):iend(1) ) ) |
---|
| 528 | IF( PRESENT(pwrk2d) ) ALLOCATE( tree(ii)%current%zwrk2d(isrt(1):iend(1),isrt(2):iend(2) ) ) |
---|
| 529 | IF( PRESENT(pwrk3d) ) ALLOCATE( tree(ii)%current%zwrk3d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3) ) ) |
---|
| 530 | IF( PRESENT(pwrk4d) ) ALLOCATE( tree(ii)%current%zwrk4d(isrt(1):iend(1),isrt(2):iend(2),isrt(3):iend(3),isrt(4):iend(4)) ) |
---|
[3134] | 531 | |
---|
| 532 | ELSE |
---|
| 533 | tree(ii)%current => tree(ii)%current%next ! the current leaf becomes the next one |
---|
| 534 | ENDIF |
---|
| 535 | ! |
---|
| 536 | ! use the array of the current leaf as a work array |
---|
| 537 | IF( PRESENT(kwrk1d) ) kwrk1d => tree(ii)%current%iwrk1d |
---|
| 538 | IF( PRESENT(kwrk2d) ) kwrk2d => tree(ii)%current%iwrk2d |
---|
| 539 | IF( PRESENT(kwrk3d) ) kwrk3d => tree(ii)%current%iwrk3d |
---|
| 540 | IF( PRESENT(kwrk4d) ) kwrk4d => tree(ii)%current%iwrk4d |
---|
| 541 | IF( PRESENT(pwrk1d) ) pwrk1d => tree(ii)%current%zwrk1d |
---|
| 542 | IF( PRESENT(pwrk2d) ) pwrk2d => tree(ii)%current%zwrk2d |
---|
| 543 | IF( PRESENT(pwrk3d) ) pwrk3d => tree(ii)%current%zwrk3d |
---|
| 544 | IF( PRESENT(pwrk4d) ) pwrk4d => tree(ii)%current%zwrk4d |
---|
| 545 | tree(ii)%current%in_use = .TRUE. ! this leaf is now used |
---|
| 546 | ! |
---|
| 547 | END SUBROUTINE wrk_allocbase |
---|
| 548 | |
---|
| 549 | |
---|
[3144] | 550 | SUBROUTINE wrk_deallocbase( ktype, kidim, kjdim, kkdim, kldim, kisrt, kjsrt, kksrt, klsrt ) |
---|
| 551 | INTEGER, INTENT(in ) :: ktype |
---|
| 552 | INTEGER, INTENT(in ) :: kidim, kjdim, kkdim, kldim |
---|
| 553 | INTEGER, INTENT(in ), OPTIONAL :: kisrt, kjsrt, kksrt, klsrt |
---|
[3134] | 554 | ! |
---|
[3144] | 555 | INTEGER, DIMENSION(jpmaxdim) :: ishape, istart |
---|
[3134] | 556 | INTEGER :: ii |
---|
| 557 | |
---|
[3144] | 558 | ishape(:) = (/ kidim, kjdim, kkdim, kldim /) |
---|
| 559 | IF( PRESENT(kisrt) ) THEN ; istart(1) = kisrt ; ELSE ; istart(1) = 1 ; ENDIF |
---|
| 560 | IF( PRESENT(kjsrt) ) THEN ; istart(2) = kjsrt ; ELSE ; istart(2) = 1 ; ENDIF |
---|
| 561 | IF( PRESENT(kksrt) ) THEN ; istart(3) = kksrt ; ELSE ; istart(3) = 1 ; ENDIF |
---|
| 562 | IF( PRESENT(klsrt) ) THEN ; istart(4) = klsrt ; ELSE ; istart(4) = 1 ; ENDIF |
---|
| 563 | |
---|
[3134] | 564 | ! find the branch with the matcing shape and type or get the first "free" branch |
---|
| 565 | ii = 1 |
---|
[3144] | 566 | DO WHILE( ANY( tree(ii)%ishape /= ishape ) .OR. ANY( tree(ii)%istart /= istart ) .OR. tree(ii)%itype /= ktype ) |
---|
[3134] | 567 | ii = ii + 1 |
---|
| 568 | END DO |
---|
| 569 | ! |
---|
| 570 | tree(ii)%current%in_use = .FALSE. ! current leaf is no more used |
---|
| 571 | tree(ii)%current => tree(ii)%current%prev ! move back toward previous leaf |
---|
| 572 | ! |
---|
| 573 | END SUBROUTINE wrk_deallocbase |
---|
| 574 | |
---|
| 575 | |
---|
| 576 | SUBROUTINE wrk_stop(cmsg) |
---|
| 577 | !!---------------------------------------------------------------------- |
---|
| 578 | !! *** ROUTINE wrk_stop *** |
---|
| 579 | !! ** Purpose : to act as local alternative to ctl_stop. |
---|
| 580 | !! Avoids dependency on in_out_manager module. |
---|
| 581 | !!---------------------------------------------------------------------- |
---|
| 582 | CHARACTER(LEN=*), INTENT(in) :: cmsg |
---|
| 583 | !!---------------------------------------------------------------------- |
---|
| 584 | ! |
---|
| 585 | ! WRITE(kumout, cform_err2) |
---|
| 586 | WRITE(*,*) TRIM(cmsg) |
---|
| 587 | ! ARPDBG - would like to CALL mppstop here to force a stop but that |
---|
| 588 | ! introduces a dependency on lib_mpp. Could CALL mpi_abort() directly |
---|
| 589 | ! but that's fairly brutal. Better to rely on CALLing routine to |
---|
| 590 | ! deal with the error passed back from the wrk_X routine? |
---|
| 591 | !CALL mppstop |
---|
| 592 | ! |
---|
| 593 | END SUBROUTINE wrk_stop |
---|
| 594 | |
---|
| 595 | !!===================================================================== |
---|
[3186] | 596 | END MODULE wrk_nemo |
---|