Changeset 789 for trunk/NEMO/OPA_SRC/TRA/trabbl_adv.h90
- Timestamp:
- 2008-01-11T19:04:56+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/trabbl_adv.h90
r719 r789 99 99 ! mbathy= number of w-level, minimum value=1 (cf dommsk.F) 100 100 101 #if defined key_vectopt_loop && ! defined key_mpp_omp101 #if defined key_vectopt_loop 102 102 jj = 1 103 103 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 115 115 zunb(ji,jj) = un(ji,jj,mbku(ji,jj)) 116 116 zvnb(ji,jj) = vn(ji,jj,mbkv(ji,jj)) 117 #if ! defined key_vectopt_loop || defined key_mpp_omp117 #if ! defined key_vectopt_loop 118 118 END DO 119 119 #endif … … 229 229 IF( ln_zps ) THEN ! partial steps correction 230 230 231 # if defined key_vectopt_loop && ! defined key_mpp_omp231 # if defined key_vectopt_loop 232 232 jj = 1 233 233 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 249 249 v_bbl(ji,jj,ikv) = zalphay(ji,jj) * vn(ji,jj,ikv) * ze3v / fse3v(ji,jj,ikv) 250 250 ENDIF 251 # if ! defined key_vectopt_loop || defined key_mpp_omp251 # if ! defined key_vectopt_loop 252 252 END DO 253 253 # endif … … 259 259 ELSE ! if not partial step loop over the whole domain no lbc call 260 260 261 #if defined key_vectopt_loop && ! defined key_mpp_omp261 #if defined key_vectopt_loop 262 262 jj = 1 263 263 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 272 272 v_bbl(ji,jj,ikv) = zalphay(ji,jj) * vn(ji,jj,ikv) 273 273 ENDIF 274 #if ! defined key_vectopt_loop || defined key_mpp_omp274 #if ! defined key_vectopt_loop 275 275 END DO 276 276 #endif … … 284 284 ! ... Second order centered tracer flux at u and v-points 285 285 286 # if defined key_vectopt_loop && ! defined key_mpp_omp286 # if defined key_vectopt_loop 287 287 jj = 1 288 288 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 309 309 zwz(ji,jj) = ( ( zfvj + ABS( zfvj ) ) * zsbb(ji ,jj ) & 310 310 & +( zfvj - ABS( zfvj ) ) * zsbb(ji ,jj+1) ) * 0.5 311 #if ! defined key_vectopt_loop || defined key_mpp_omp311 #if ! defined key_vectopt_loop 312 312 END DO 313 313 #endif 314 314 END DO 315 # if defined key_vectopt_loop && ! defined key_mpp_omp315 # if defined key_vectopt_loop 316 316 jj = 1 317 317 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 331 331 ta(ji,jj,ik) = ta(ji,jj,ik) + zta 332 332 sa(ji,jj,ik) = sa(ji,jj,ik) + zsa 333 #if ! defined key_vectopt_loop || defined key_mpp_omp333 #if ! defined key_vectopt_loop 334 334 END DO 335 335 #endif … … 365 365 IF( ln_zps ) THEN 366 366 367 # if defined key_vectopt_loop && ! defined key_mpp_omp367 # if defined key_vectopt_loop 368 368 jj = 1 369 369 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 383 383 zwu(ji,jj) = zalphax(ji,jj) * e2u(ji,jj) * ze3u 384 384 zwv(ji,jj) = zalphay(ji,jj) * e1v(ji,jj) * ze3v 385 #if ! defined key_vectopt_loop || defined key_mpp_omp385 #if ! defined key_vectopt_loop 386 386 END DO 387 387 #endif … … 390 390 ELSE 391 391 392 # if defined key_vectopt_loop && ! defined key_mpp_omp392 # if defined key_vectopt_loop 393 393 jj = 1 394 394 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 401 401 zwu(ji,jj) = zalphax(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,iku) 402 402 zwv(ji,jj) = zalphay(ji,jj) * e1v(ji,jj) * fse3v(ji,jj,ikv) 403 #if ! defined key_vectopt_loop || defined key_mpp_omp403 #if ! defined key_vectopt_loop 404 404 END DO 405 405 #endif … … 409 409 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 = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 425 425 & ) / zbt 426 426 427 # if ! defined key_vectopt_loop || defined key_mpp_omp427 # if ! defined key_vectopt_loop 428 428 END DO 429 429 # endif
Note: See TracChangeset
for help on using the changeset viewer.