Changeset 1326 for trunk/NEMO/OFF_SRC/LDF/ldfslp.F90
- Timestamp:
- 2009-02-20T11:08:16+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OFF_SRC/LDF/ldfslp.F90
r1156 r1326 140 140 141 141 IF( ln_zps ) THEN ! partial steps correction at the bottom ocean level (zps_hde routine) 142 # if defined key_vectopt_loop && ! defined key_mpp_omp142 # if defined key_vectopt_loop 143 143 jj = 1 144 144 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 152 152 zgru(ji,jj,iku) = gru(ji,jj) 153 153 zgrv(ji,jj,ikv) = grv(ji,jj) 154 # if ! defined key_vectopt_loop || defined key_mpp_omp154 # if ! defined key_vectopt_loop 155 155 END DO 156 156 # endif … … 491 491 ! mask for mixed layer 492 492 DO jk = 1, jpk 493 # if defined key_vectopt_loop && ! defined key_mpp_omp493 # if defined key_vectopt_loop 494 494 jj = 1 495 495 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 505 505 omlmask(ji,jj,jk) = 0.e0 506 506 ENDIF 507 # if ! defined key_vectopt_loop || defined key_mpp_omp507 # if ! defined key_vectopt_loop 508 508 END DO 509 509 # endif … … 523 523 zwy(:,jpj) = 0.e0 524 524 zwy(jpi,:) = 0.e0 525 # if defined key_vectopt_loop && ! defined key_mpp_omp525 # if defined key_vectopt_loop 526 526 jj = 1 527 527 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 536 536 & * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) ) & 537 537 & / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 538 # if ! defined key_vectopt_loop || defined key_mpp_omp538 # if ! defined key_vectopt_loop 539 539 END DO 540 540 # endif … … 544 544 545 545 ! Slope at u points 546 # if defined key_vectopt_loop && ! defined key_mpp_omp546 # if defined key_vectopt_loop 547 547 jj = 1 548 548 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 561 561 ! uslpml 562 562 uslpml (ji,jj) = zau / ( zbu - zeps ) * umask (ji,jj,ik) 563 # if ! defined key_vectopt_loop || defined key_mpp_omp563 # if ! defined key_vectopt_loop 564 564 END DO 565 565 # endif … … 573 573 zwy ( :, jpj) = 0.e0 574 574 zwy ( jpi, :) = 0.e0 575 # if defined key_vectopt_loop && ! defined key_mpp_omp575 # if defined key_vectopt_loop 576 576 jj = 1 577 577 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 585 585 & * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) ) & 586 586 & / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 587 # if ! defined key_vectopt_loop || defined key_mpp_omp587 # if ! defined key_vectopt_loop 588 588 END DO 589 589 # endif … … 594 594 595 595 ! Slope at v points 596 # if defined key_vectopt_loop && ! defined key_mpp_omp596 # if defined key_vectopt_loop 597 597 jj = 1 598 598 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 611 611 ! vslpml 612 612 vslpml (ji,jj) = zav / ( zbv - zeps ) * vmask (ji,jj,ik) 613 # if ! defined key_vectopt_loop || defined key_mpp_omp613 # if ! defined key_vectopt_loop 614 614 END DO 615 615 # endif … … 625 625 ! Local vertical density gradient evaluated from N^2 626 626 ! zwy = d/dz(prd)= - mk ( prd ) / grav * pn2 -- at w point 627 # if defined key_vectopt_loop && ! defined key_mpp_omp627 # if defined key_vectopt_loop 628 628 jj = 1 629 629 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 637 637 zwy (ji,jj) = zm05g * pn2 (ji,jj,ik) * & 638 638 & ( prd (ji,jj,ik) + prd (ji,jj,ikm1) + 2. ) 639 # if ! defined key_vectopt_loop || defined key_mpp_omp639 # if ! defined key_vectopt_loop 640 640 END DO 641 641 # endif … … 643 643 644 644 ! Slope at w point 645 # if defined key_vectopt_loop && ! defined key_mpp_omp645 # if defined key_vectopt_loop 646 646 jj = 1 647 647 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 673 673 wslpiml (ji,jj) = zai / ( zbi - zeps) * tmask (ji,jj,ik) 674 674 wslpjml (ji,jj) = zaj / ( zbj - zeps) * tmask (ji,jj,ik) 675 # if ! defined key_vectopt_loop || defined key_mpp_omp675 # if ! defined key_vectopt_loop 676 676 END DO 677 677 # endif
Note: See TracChangeset
for help on using the changeset viewer.