New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 811 for branches/dev_001_SBC/NEMO/OPA_SRC/LDF/ldfslp.F90 – NEMO

Ignore:
Timestamp:
2008-02-07T17:00:12+01:00 (16 years ago)
Author:
ctlod
Message:

dev_001_SBC: merge with the trunk last changesets: #780, 782, 783, 784, 785, 788, 789, 793, 794

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/dev_001_SBC/NEMO/OPA_SRC/LDF/ldfslp.F90

    r699 r811  
    141141 
    142142      IF( ln_zps ) THEN      ! partial steps correction at the bottom ocean level (zps_hde routine) 
    143 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     143# if defined key_vectopt_loop   
    144144         jj = 1 
    145145         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    153153               zgru(ji,jj,iku) = gru(ji,jj)  
    154154               zgrv(ji,jj,ikv) = grv(ji,jj)                
    155 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     155# if ! defined key_vectopt_loop  
    156156            END DO 
    157157# endif 
     
    492492      ! mask for mixed layer 
    493493      DO jk = 1, jpk 
    494 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     494# if defined key_vectopt_loop 
    495495         jj = 1 
    496496         DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    506506                  omlmask(ji,jj,jk) = 0.e0 
    507507               ENDIF 
    508 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     508# if ! defined key_vectopt_loop 
    509509            END DO 
    510510# endif 
     
    524524      zwy(:,jpj) = 0.e0 
    525525      zwy(jpi,:) = 0.e0 
    526 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     526# if defined key_vectopt_loop 
    527527      jj = 1 
    528528      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    537537               &             * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) )   & 
    538538               &             / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 
    539 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     539# if ! defined key_vectopt_loop 
    540540         END DO 
    541541# endif 
     
    545545 
    546546      ! Slope at u points 
    547 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     547# if defined key_vectopt_loop 
    548548      jj = 1 
    549549      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    562562            ! uslpml 
    563563            uslpml (ji,jj) = zau / ( zbu - zeps ) * umask (ji,jj,ik) 
    564 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     564# if ! defined key_vectopt_loop 
    565565         END DO 
    566566# endif 
     
    574574      zwy ( :, jpj) = 0.e0 
    575575      zwy ( jpi, :) = 0.e0 
    576 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     576# if defined key_vectopt_loop 
    577577      jj = 1 
    578578      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    586586               &             * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) )   & 
    587587               &             / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 
    588 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     588# if ! defined key_vectopt_loop 
    589589         END DO 
    590590# endif 
     
    595595 
    596596      ! Slope at v points 
    597 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     597# if defined key_vectopt_loop 
    598598      jj = 1 
    599599      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    612612            ! vslpml 
    613613            vslpml (ji,jj) = zav / ( zbv - zeps ) * vmask (ji,jj,ik) 
    614 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     614# if ! defined key_vectopt_loop 
    615615         END DO 
    616616# endif 
     
    626626      ! Local vertical density gradient evaluated from N^2 
    627627      ! zwy = d/dz(prd)= - mk ( prd ) / grav * pn2 -- at w point 
    628 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     628# if defined key_vectopt_loop 
    629629      jj = 1 
    630630      DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    638638            zwy (ji,jj) = zm05g * pn2 (ji,jj,ik) *     & 
    639639               &             ( prd (ji,jj,ik) + prd (ji,jj,ikm1) + 2. ) 
    640 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     640# if ! defined key_vectopt_loop 
    641641         END DO 
    642642# endif 
     
    644644 
    645645      ! Slope at w point 
    646 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     646# if defined key_vectopt_loop 
    647647      jj = 1 
    648648      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    674674            wslpiml (ji,jj) = zai / ( zbi - zeps) * tmask (ji,jj,ik) 
    675675            wslpjml (ji,jj) = zaj / ( zbj - zeps) * tmask (ji,jj,ik) 
    676 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     676# if ! defined key_vectopt_loop 
    677677         END DO 
    678678# endif 
Note: See TracChangeset for help on using the changeset viewer.