- Timestamp:
- 2017-12-13T14:57:33+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/LBC/mpp_lnk_generic.h90
r8882 r9012 72 72 ipf = F_SIZE(ptab) ! 5th - use in "multi" case (array of pointers) 73 73 ! 74 ALLOCATE( zt3ns(jpi, jprecj,ipk,ipl,ipf,2), zt3sn(jpi,jprecj,ipk,ipl,ipf,2), &75 & zt3ew(jpj, jpreci,ipk,ipl,ipf,2), zt3we(jpj,jpreci,ipk,ipl,ipf,2) )74 ALLOCATE( zt3ns(jpi,nn_hls,ipk,ipl,ipf,2), zt3sn(jpi,nn_hls,ipk,ipl,ipf,2), & 75 & zt3ew(jpj,nn_hls,ipk,ipl,ipf,2), zt3we(jpj,nn_hls,ipk,ipl,ipf,2) ) 76 76 ! 77 77 ll_Iperio = nbondi == 2 .AND. (nperio == 1 .OR. nperio == 4 .OR. nperio == 6) … … 116 116 ARRAY_IN(jpi,:,:,:,jf) = ARRAY_IN( 2 ,:,:,:,jf) 117 117 ELSE !* closed 118 IF( .NOT. NAT_IN(jf) == 'F' ) ARRAY_IN( 1 : jpreci,:,:,:,jf) = zland ! east except F-point119 ARRAY_IN(nlci- jpreci+1:jpi ,:,:,:,jf) = zland ! west118 IF( .NOT. NAT_IN(jf) == 'F' ) ARRAY_IN( 1 :nn_hls,:,:,:,jf) = zland ! east except F-point 119 ARRAY_IN(nlci-nn_hls+1:jpi ,:,:,:,jf) = zland ! west 120 120 ENDIF 121 121 ! ! North-South boundaries … … 124 124 ARRAY_IN(:,jpj,:,:,jf) = ARRAY_IN(:, 2 ,:,:,jf) 125 125 ELSE !* closed 126 IF( .NOT. NAT_IN(jf) == 'F' ) ARRAY_IN(:, 1 : jprecj,:,:,jf) = zland ! south except F-point127 ARRAY_IN(:,nlcj- jprecj+1:jpj ,:,:,jf) = zland ! north126 IF( .NOT. NAT_IN(jf) == 'F' ) ARRAY_IN(:, 1 :nn_hls,:,:,jf) = zland ! south except F-point 127 ARRAY_IN(:,nlcj-nn_hls+1:jpj ,:,:,jf) = zland ! north 128 128 ENDIF 129 129 END DO … … 142 142 DO jl = 1, ipl 143 143 DO jk = 1, ipk 144 DO jh = 1, jpreci145 zt3ew(:,jh,jk,jl,jf,1) = ARRAY_IN( jpreci+jh,:,jk,jl,jf)144 DO jh = 1, nn_hls 145 zt3ew(:,jh,jk,jl,jf,1) = ARRAY_IN(nn_hls+jh,:,jk,jl,jf) 146 146 zt3we(:,jh,jk,jl,jf,1) = ARRAY_IN(iihom +jh,:,jk,jl,jf) 147 147 END DO … … 152 152 ! 153 153 ! ! Migrations 154 imigr = jpreci* jpj * ipk * ipl * ipf154 imigr = nn_hls * jpj * ipk * ipl * ipf 155 155 ! 156 156 SELECT CASE ( nbondi ) … … 173 173 ! 174 174 ! ! Write Dirichlet lateral conditions 175 iihom = nlci- jpreci175 iihom = nlci-nn_hls 176 176 ! 177 177 SELECT CASE ( nbondi ) … … 180 180 DO jl = 1, ipl 181 181 DO jk = 1, ipk 182 DO jh = 1, jpreci182 DO jh = 1, nn_hls 183 183 ARRAY_IN(iihom+jh,:,jk,jl,jf) = zt3ew(:,jh,jk,jl,jf,2) 184 184 END DO … … 190 190 DO jl = 1, ipl 191 191 DO jk = 1, ipk 192 DO jh = 1, jpreci192 DO jh = 1, nn_hls 193 193 ARRAY_IN(jh ,:,jk,jl,jf) = zt3we(:,jh,jk,jl,jf,2) 194 194 ARRAY_IN(iihom+jh,:,jk,jl,jf) = zt3ew(:,jh,jk,jl,jf,2) … … 201 201 DO jl = 1, ipl 202 202 DO jk = 1, ipk 203 DO jh = 1, jpreci203 DO jh = 1, nn_hls 204 204 ARRAY_IN(jh ,:,jk,jl,jf) = zt3we(:,jh,jk,jl,jf,2) 205 205 END DO … … 218 218 DO jl = 1, ipl 219 219 DO jk = 1, ipk 220 DO jh = 1, jprecj220 DO jh = 1, nn_hls 221 221 zt3sn(:,jh,jk,jl,jf,1) = ARRAY_IN(:,ijhom +jh,jk,jl,jf) 222 zt3ns(:,jh,jk,jl,jf,1) = ARRAY_IN(:, jprecj+jh,jk,jl,jf)222 zt3ns(:,jh,jk,jl,jf,1) = ARRAY_IN(:,nn_hls+jh,jk,jl,jf) 223 223 END DO 224 224 END DO … … 228 228 ! 229 229 ! ! Migrations 230 imigr = jprecj* jpi * ipk * ipl * ipf230 imigr = nn_hls * jpi * ipk * ipl * ipf 231 231 ! 232 232 SELECT CASE ( nbondj ) … … 249 249 ! 250 250 ! ! Write Dirichlet lateral conditions 251 ijhom = nlcj- jprecj251 ijhom = nlcj-nn_hls 252 252 ! 253 253 SELECT CASE ( nbondj ) … … 256 256 DO jl = 1, ipl 257 257 DO jk = 1, ipk 258 DO jh = 1, jprecj258 DO jh = 1, nn_hls 259 259 ARRAY_IN(:,ijhom+jh,jk,jl,jf) = zt3ns(:,jh,jk,jl,jf,2) 260 260 END DO … … 266 266 DO jl = 1, ipl 267 267 DO jk = 1, ipk 268 DO jh = 1, jprecj268 DO jh = 1, nn_hls 269 269 ARRAY_IN(:, jh,jk,jl,jf) = zt3sn(:,jh,jk,jl,jf,2) 270 270 ARRAY_IN(:,ijhom+jh,jk,jl,jf) = zt3ns(:,jh,jk,jl,jf,2) … … 277 277 DO jl = 1, ipl 278 278 DO jk = 1, ipk 279 DO jh = 1, jprecj279 DO jh = 1, nn_hls 280 280 ARRAY_IN(:,jh,jk,jl,jf) = zt3sn(:,jh,jk,jl,jf,2) 281 281 END DO
Note: See TracChangeset
for help on using the changeset viewer.