- Timestamp:
- 2010-09-29T17:31:36+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/devmercator2010_1/NEMO/OPA_SRC/DYN/dynzdf_imp.F90
r1662 r2132 24 24 USE phycst ! physical constants 25 25 USE in_out_manager ! I/O manager 26 #if defined key_zdfgls 27 USE zdfbfr, ONLY : bfrua, bfrva, wbotu, wbotv ! bottom stresses 28 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 29 #endif 26 30 27 31 IMPLICIT NONE … … 75 79 REAL(wp) :: zzwi ! temporary scalars 76 80 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zwi ! temporary workspace arrays 81 #if defined key_zdfgls 82 INTEGER :: ikbu, ikbv, ikbum1, ikbvm1 83 REAL(wp) :: zcbcu, zcbcv 84 #endif 77 85 !!---------------------------------------------------------------------- 78 86 … … 173 181 END DO 174 182 END DO 183 184 #if defined key_zdfgls 185 ! Save bottom stress for next time step 186 DO jj = 2, jpjm1 187 DO ji = fs_2, fs_jpim1 ! vector opt. 188 ikbu = MIN( mbathy(ji+1,jj ), mbathy(ji,jj) ) 189 ikbum1 = MAX( ikbu-1, 1 ) 190 wbotu(ji,jj) = bfrua(ji,jj) * ua(ji,jj,ikbum1) * umask(ji,jj,ikbum1) 191 END DO 192 END DO 193 CALL lbc_lnk( wbotu(:,:), 'U', -1. ) 194 #endif 175 195 176 196 ! Normalization to obtain the general momentum trend ua … … 272 292 END DO 273 293 294 #if defined key_zdfgls 295 ! Save bottom stress for next time step 296 DO jj = 2, jpjm1 297 DO ji = fs_2, fs_jpim1 ! vector opt. 298 ikbv = MIN( mbathy(ji,jj+1), mbathy(ji,jj) ) 299 ikbvm1 = MAX( ikbv-1, 1 ) 300 wbotv(ji,jj) = bfrva(ji,jj) * va(ji,jj,ikbvm1) * vmask(ji,jj,ikbvm1) 301 END DO 302 END DO 303 CALL lbc_lnk( wbotv(:,:), 'V', -1. ) 304 #endif 305 274 306 ! Normalization to obtain the general momentum trend va 275 307 DO jk = 1, jpkm1
Note: See TracChangeset
for help on using the changeset viewer.