- Timestamp:
- 2014-04-06T17:28:25+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90
r4381 r4616 105 105 IF ( ln_loglayer.AND.lk_vvl ) THEN ! "log layer" bottom friction coefficient 106 106 107 # if defined key_vectopt_loop108 DO jj = 1, 1109 !CDIR NOVERRCHK110 DO ji = 1, jpij ! vector opt. (forced unrolling)111 # else112 !CDIR NOVERRCHK113 107 DO jj = 1, jpj 114 !CDIR NOVERRCHK115 108 DO ji = 1, jpi 116 # endif117 109 ikbt = mbkt(ji,jj) 118 110 ! JC: possible WAD implementation should modify line below if layers vanish … … 127 119 ENDIF 128 120 129 # if defined key_vectopt_loop130 DO jj = 1, 1131 !CDIR NOVERRCHK132 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling)133 # else134 !CDIR NOVERRCHK135 121 DO jj = 2, jpjm1 136 !CDIR NOVERRCHK137 122 DO ji = 2, jpim1 138 # endif139 123 ikbu = mbku(ji,jj) ! ocean bottom level at u- and v-points 140 124 ikbv = mbkv(ji,jj) ! (deepest ocean u- and v-points) … … 265 249 ! 266 250 IF ( ln_loglayer.AND.(.NOT.lk_vvl) ) THEN ! set "log layer" bottom friction once for all 267 # if defined key_vectopt_loop268 DO jj = 1, 1269 !CDIR NOVERRCHK270 DO ji = 1, jpij ! vector opt. (forced unrolling)271 # else272 !CDIR NOVERRCHK273 251 DO jj = 1, jpj 274 !CDIR NOVERRCHK275 252 DO ji = 1, jpi 276 # endif277 253 ikbt = mbkt(ji,jj) 278 254 ztmp = tmask(ji,jj,ikbt) * ( vkarmn / LOG( 0.5_wp * fse3t_n(ji,jj,ikbt) / rn_bfrz0 ))**2._wp … … 309 285 zmaxbfr = -1.e10_wp ! initialise tracker for maximum of bottom friction coefficient 310 286 ! 311 # if defined key_vectopt_loop312 DO jj = 1, 1313 !CDIR NOVERRCHK314 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling)315 # else316 !CDIR NOVERRCHK317 287 DO jj = 2, jpjm1 318 !CDIR NOVERRCHK319 288 DO ji = 2, jpim1 320 # endif321 289 ikbu = mbku(ji,jj) ! deepest ocean level at u- and v-points 322 290 ikbv = mbkv(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.