Changeset 811 for branches/dev_001_SBC/NEMO/OPA_SRC/LDF/ldfslp.F90
- Timestamp:
- 2008-02-07T17:00:12+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_SBC/NEMO/OPA_SRC/LDF/ldfslp.F90
r699 r811 141 141 142 142 IF( ln_zps ) THEN ! partial steps correction at the bottom ocean level (zps_hde routine) 143 # if defined key_vectopt_loop && ! defined key_mpp_omp143 # if defined key_vectopt_loop 144 144 jj = 1 145 145 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 153 153 zgru(ji,jj,iku) = gru(ji,jj) 154 154 zgrv(ji,jj,ikv) = grv(ji,jj) 155 # if ! defined key_vectopt_loop || defined key_mpp_omp155 # if ! defined key_vectopt_loop 156 156 END DO 157 157 # endif … … 492 492 ! mask for mixed layer 493 493 DO jk = 1, jpk 494 # if defined key_vectopt_loop && ! defined key_mpp_omp494 # if defined key_vectopt_loop 495 495 jj = 1 496 496 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 506 506 omlmask(ji,jj,jk) = 0.e0 507 507 ENDIF 508 # if ! defined key_vectopt_loop || defined key_mpp_omp508 # if ! defined key_vectopt_loop 509 509 END DO 510 510 # endif … … 524 524 zwy(:,jpj) = 0.e0 525 525 zwy(jpi,:) = 0.e0 526 # if defined key_vectopt_loop && ! defined key_mpp_omp526 # if defined key_vectopt_loop 527 527 jj = 1 528 528 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 537 537 & * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) ) & 538 538 & / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 539 # if ! defined key_vectopt_loop || defined key_mpp_omp539 # if ! defined key_vectopt_loop 540 540 END DO 541 541 # endif … … 545 545 546 546 ! Slope at u points 547 # if defined key_vectopt_loop && ! defined key_mpp_omp547 # if defined key_vectopt_loop 548 548 jj = 1 549 549 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 562 562 ! uslpml 563 563 uslpml (ji,jj) = zau / ( zbu - zeps ) * umask (ji,jj,ik) 564 # if ! defined key_vectopt_loop || defined key_mpp_omp564 # if ! defined key_vectopt_loop 565 565 END DO 566 566 # endif … … 574 574 zwy ( :, jpj) = 0.e0 575 575 zwy ( jpi, :) = 0.e0 576 # if defined key_vectopt_loop && ! defined key_mpp_omp576 # if defined key_vectopt_loop 577 577 jj = 1 578 578 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) … … 586 586 & * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) ) & 587 587 & / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 588 # if ! defined key_vectopt_loop || defined key_mpp_omp588 # if ! defined key_vectopt_loop 589 589 END DO 590 590 # endif … … 595 595 596 596 ! Slope at v points 597 # if defined key_vectopt_loop && ! defined key_mpp_omp597 # if defined key_vectopt_loop 598 598 jj = 1 599 599 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 612 612 ! vslpml 613 613 vslpml (ji,jj) = zav / ( zbv - zeps ) * vmask (ji,jj,ik) 614 # if ! defined key_vectopt_loop || defined key_mpp_omp614 # if ! defined key_vectopt_loop 615 615 END DO 616 616 # endif … … 626 626 ! Local vertical density gradient evaluated from N^2 627 627 ! zwy = d/dz(prd)= - mk ( prd ) / grav * pn2 -- at w point 628 # if defined key_vectopt_loop && ! defined key_mpp_omp628 # if defined key_vectopt_loop 629 629 jj = 1 630 630 DO ji = 1, jpij ! vector opt. (forced unrolling) … … 638 638 zwy (ji,jj) = zm05g * pn2 (ji,jj,ik) * & 639 639 & ( prd (ji,jj,ik) + prd (ji,jj,ikm1) + 2. ) 640 # if ! defined key_vectopt_loop || defined key_mpp_omp640 # if ! defined key_vectopt_loop 641 641 END DO 642 642 # endif … … 644 644 645 645 ! Slope at w point 646 # if defined key_vectopt_loop && ! defined key_mpp_omp646 # if defined key_vectopt_loop 647 647 jj = 1 648 648 DO ji = jpi+2, jpij-jpi-1 ! vector opt. (forced unrolling) … … 674 674 wslpiml (ji,jj) = zai / ( zbi - zeps) * tmask (ji,jj,ik) 675 675 wslpjml (ji,jj) = zaj / ( zbj - zeps) * tmask (ji,jj,ik) 676 # if ! defined key_vectopt_loop || defined key_mpp_omp676 # if ! defined key_vectopt_loop 677 677 END DO 678 678 # endif
Note: See TracChangeset
for help on using the changeset viewer.