Changeset 1197 for trunk/NEMO/TOP_SRC/TRP/trcbbl_adv.h90
- Timestamp:
- 2008-09-22T10:03:39+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/TRP/trcbbl_adv.h90
r1152 r1197 110 110 ! mbathy= number of w-level, minimum value=1 (cf dommsk.F) 111 111 112 #if defined key_vectopt_loop && ! defined key_mpp_omp112 #if defined key_vectopt_loop 113 113 jj = 1 114 114 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 124 124 zunb(ji,jj) = un(ji,jj,mbku(ji,jj)) 125 125 zvnb(ji,jj) = vn(ji,jj,mbkv(ji,jj)) 126 #if ! defined key_vectopt_loop || defined key_mpp_omp126 #if ! defined key_vectopt_loop 127 127 END DO 128 128 #endif … … 260 260 ! partial steps correction 261 261 262 #if defined key_vectopt_loop && ! defined key_mpp_omp262 #if defined key_vectopt_loop 263 263 jj = 1 264 264 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 280 280 v_trc_bbl(ji,jj,ikv) = zalphay(ji,jj) * vn(ji,jj,ikv) * ze3v / fse3v(ji,jj,ikv) 281 281 ENDIF 282 #if ! defined key_vectopt_loop || defined key_mpp_omp282 #if ! defined key_vectopt_loop 283 283 END DO 284 284 #endif … … 291 291 ! if not partial step loop over the whole domain no lbc call 292 292 293 #if defined key_vectopt_loop && ! defined key_mpp_omp293 #if defined key_vectopt_loop 294 294 jj = 1 295 295 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 304 304 v_trc_bbl(ji,jj,ikv) = zalphay(ji,jj) * vn(ji,jj,ikv) 305 305 ENDIF 306 #if ! defined key_vectopt_loop || defined key_mpp_omp306 #if ! defined key_vectopt_loop 307 307 END DO 308 308 #endif … … 313 313 DO jn = 1, jptra 314 314 315 #if defined key_vectopt_loop && ! defined key_mpp_omp315 #if defined key_vectopt_loop 316 316 jj = 1 317 317 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 322 322 ik = mbkt(ji,jj) ! index of the bottom ocean T-level 323 323 ztrb(ji,jj) = trb(ji,jj,ik,jn) * tmask(ji,jj,1) ! masked now T at the ocean bottom 324 #if ! defined key_vectopt_loop || defined key_mpp_omp324 #if ! defined key_vectopt_loop 325 325 END DO 326 326 #endif … … 333 333 ! ... Second order centered tracer flux at u and v-points 334 334 335 # if defined key_vectopt_loop && ! defined key_mpp_omp335 # if defined key_vectopt_loop 336 336 jj = 1 337 337 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 349 349 zwy(ji,jj) = ( ( zfvj + ABS( zfvj ) ) * ztrb(ji ,jj ) & 350 350 & +( zfvj - ABS( zfvj ) ) * ztrb(ji ,jj+1) ) * 0.5 351 #if ! defined key_vectopt_loop || defined key_mpp_omp351 #if ! defined key_vectopt_loop 352 352 END DO 353 353 #endif 354 354 END DO 355 355 356 # if defined key_vectopt_loop && ! defined key_mpp_omp356 # if defined key_vectopt_loop 357 357 jj = 1 358 358 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 369 369 ! add it to the general tracer trends 370 370 tra(ji,jj,ik,jn) = tra(ji,jj,ik,jn) + ztra 371 #if ! defined key_vectopt_loop || defined key_mpp_omp371 #if ! defined key_vectopt_loop 372 372 END DO 373 373 #endif … … 409 409 IF( ln_zps ) THEN 410 410 411 # if defined key_vectopt_loop && ! defined key_mpp_omp411 # if defined key_vectopt_loop 412 412 jj = 1 413 413 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 427 427 zwu(ji,jj) = zalphax(ji,jj) * e2u(ji,jj) * ze3u 428 428 zwv(ji,jj) = zalphay(ji,jj) * e1v(ji,jj) * ze3v 429 #if ! defined key_vectopt_loop || defined key_mpp_omp429 #if ! defined key_vectopt_loop 430 430 END DO 431 431 #endif … … 434 434 ELSE 435 435 436 # if defined key_vectopt_loop && ! defined key_mpp_omp436 # if defined key_vectopt_loop 437 437 jj = 1 438 438 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 445 445 zwu(ji,jj) = zalphax(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,iku) 446 446 zwv(ji,jj) = zalphay(ji,jj) * e1v(ji,jj) * fse3v(ji,jj,ikv) 447 #if ! defined key_vectopt_loop || defined key_mpp_omp447 #if ! defined key_vectopt_loop 448 448 END DO 449 449 #endif … … 451 451 ENDIF 452 452 453 # if defined key_vectopt_loop && ! defined key_mpp_omp453 # if defined key_vectopt_loop 454 454 jj = 1 455 455 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 467 467 & ) / zbt 468 468 469 # if ! defined key_vectopt_loop || defined key_mpp_omp469 # if ! defined key_vectopt_loop 470 470 END DO 471 471 # endif
Note: See TracChangeset
for help on using the changeset viewer.