Changeset 1197


Ignore:
Timestamp:
2008-09-22T10:03:39+02:00 (12 years ago)
Author:
cetlod
Message:

remove useless CPP key, key_mpp_omp in passive tracers transport routines, see ticket:49

Location:
trunk/NEMO/TOP_SRC/TRP
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/TRP/trcbbc.F90

    r1193 r1197  
    7575      INTEGER, INTENT( in ) ::   kt                         ! ocean time-step index 
    7676 
    77 #if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     77#if defined key_vectopt_loop  
    7878      INTEGER ::   ji, jn                                   ! dummy loop indices 
    7979#else 
     
    106106 
    107107 
    108 #if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     108#if defined key_vectopt_loop 
    109109            DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
    110110               tra(ji,1,nbotlevt(ji,1),jn) = tra(ji,1,nbotlevt(ji,1),jn) + qgh_trd(ji,1) 
     
    220220      CASE ( 1:2 )                !  geothermal heat flux 
    221221 
    222 #if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     222#if defined key_vectopt_loop 
    223223         DO ji = 1, jpij   ! vector opt. (forced unrolling) 
    224224            qgh_trd(ji,1) = ro0cpr * qgh_trd(ji,1) / fse3t(ji,1,nbotlevt(ji,1) ) 
  • trunk/NEMO/TOP_SRC/TRP/trcbbl.F90

    r1175 r1197  
    146146      ! mbathy= number of w-level, minimum value=1 (cf dommsk.F) 
    147147 
    148 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     148#  if defined key_vectopt_loop 
    149149      jj = 1 
    150150      DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    157157            zsnb(ji,jj) = sn(ji,jj,ik) * tmask(ji,jj,1) 
    158158            zdep(ji,jj) = fsdept(ji,jj,ik)                ! depth of the ocean bottom T-level 
    159 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     159#  if ! defined key_vectopt_loop 
    160160         END DO 
    161161#  endif 
     
    164164      IF( ln_zps ) THEN      ! partial steps correction 
    165165 
    166 #   if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     166#   if defined key_vectopt_loop 
    167167         jj = 1 
    168168         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    179179               zahu(ji,jj) = atrcbbl * e2u(ji,jj) * ze3u / e1u(ji,jj) * umask(ji,jj,1) 
    180180               zahv(ji,jj) = atrcbbl * e1v(ji,jj) * ze3v / e2v(ji,jj) * vmask(ji,jj,1) 
    181 #   if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     181#   if ! defined key_vectopt_loop 
    182182            END DO 
    183183#   endif 
    184184         END DO 
    185185      ELSE                  ! z-coordinate - full steps or s-coordinate 
    186 #   if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     186#   if defined key_vectopt_loop 
    187187         jj = 1 
    188188         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    195195               zahu(ji,jj) = atrcbbl * e2u(ji,jj) * fse3u(ji,jj,iku) / e1u(ji,jj) * umask(ji,jj,1) 
    196196               zahv(ji,jj) = atrcbbl * e1v(ji,jj) * fse3v(ji,jj,ikv) / e2v(ji,jj) * vmask(ji,jj,1) 
    197 #   if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     197#   if ! defined key_vectopt_loop 
    198198            END DO 
    199199#   endif 
     
    220220      
    221221         ! first derivative (gradient)          
    222 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     222#  if defined key_vectopt_loop 
    223223         jj = 1 
    224224         DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    229229               ik = mbkt(ji,jj)  
    230230               ztrb(ji,jj) = trb(ji,jj,ik,jn) * tmask(ji,jj,1) 
    231 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    232             END DO 
    233 #  endif 
    234          END DO 
    235  
    236 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     231#  if ! defined key_vectopt_loop 
     232            END DO 
     233#  endif 
     234         END DO 
     235 
     236#  if defined key_vectopt_loop 
    237237         jj = 1 
    238238         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    243243               zkx(ji,jj) = bblx(ji,jj) * zahu(ji,jj) * ( ztrb(ji+1,jj) - ztrb(ji,jj) ) 
    244244               zky(ji,jj) = bbly(ji,jj) * zahv(ji,jj) * ( ztrb(ji,jj+1) - ztrb(ji,jj) ) 
    245 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     245#  if ! defined key_vectopt_loop 
    246246            END DO 
    247247#  endif 
     
    261261      CASE ( 0 )                 ! Jackett and McDougall (1994) formulation 
    262262          
    263 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     263#  if defined key_vectopt_loop 
    264264      jj = 1 
    265265      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    280280            zsign = SIGN( 0.5, - zgdrho * ( zdep(ji+1,jj) - zdep(ji,jj) ) ) 
    281281            zki(ji,jj) = ( 0.5 - zsign ) * zahu(ji,jj) 
    282 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    283          END DO 
    284 #  endif 
    285       END DO 
    286  
    287 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     282#  if ! defined key_vectopt_loop 
     283         END DO 
     284#  endif 
     285      END DO 
     286 
     287#  if defined key_vectopt_loop 
    288288      jj = 1 
    289289      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    304304            zsign = SIGN( 0.5, -zgdrho * ( zdep(ji,jj+1) - zdep(ji,jj) ) ) 
    305305            zkj(ji,jj) = ( 0.5 - zsign ) * zahv(ji,jj) 
    306 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     306#  if ! defined key_vectopt_loop 
    307307         END DO 
    308308#  endif 
     
    311311      CASE ( 1 )                 ! Linear formulation function of temperature only 
    312312 
    313 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     313#  if defined key_vectopt_loop 
    314314      jj = 1 
    315315      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    323323            zsign = SIGN( 0.5, - zgdrho * ( zdep(ji+1,jj) - zdep(ji,jj) ) ) 
    324324            zki(ji,jj) = ( 0.5 - zsign ) * zahu(ji,jj) 
    325 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    326          END DO 
    327 #  endif 
    328       END DO 
    329  
    330 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     325#  if ! defined key_vectopt_loop 
     326         END DO 
     327#  endif 
     328      END DO 
     329 
     330#  if defined key_vectopt_loop 
    331331      jj = 1 
    332332      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    341341            zkj(ji,jj) = ( 0.5 - zsign ) * zahv(ji,jj) 
    342342 
    343 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     343#  if ! defined key_vectopt_loop 
    344344         END DO 
    345345#  endif 
     
    385385 
    386386         ! first derivative (gradient) 
    387 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     387#  if defined key_vectopt_loop 
    388388         jj = 1 
    389389         DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    394394               ik = mbkt(ji,jj) 
    395395               ztrb(ji,jj) = trb(ji,jj,ik,jn) * tmask(ji,jj,1) 
    396 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    397             END DO 
    398 #  endif 
    399          END DO 
    400 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     396#  if ! defined key_vectopt_loop 
     397            END DO 
     398#  endif 
     399         END DO 
     400#  if defined key_vectopt_loop 
    401401         jj = 1 
    402402         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    407407               zkx(ji,jj) = zki(ji,jj) * ( ztrb(ji+1,jj) - ztrb(ji,jj) ) 
    408408               zky(ji,jj) = zkj(ji,jj) * ( ztrb(ji,jj+1) - ztrb(ji,jj) ) 
    409 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     409#  if ! defined key_vectopt_loop 
    410410            END DO 
    411411#  endif 
     
    449449          
    450450         ! second derivative (divergence) and add to the general tracer trend 
    451 #  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     451#  if defined key_vectopt_loop 
    452452         jj = 1 
    453453         DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    461461                  &    + zky(ji,jj) - zky(ji  ,jj-1)  ) * zbtr 
    462462               tra(ji,jj,ik,jn) = tra(ji,jj,ik,jn) + ztra 
    463 #  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     463#  if ! defined key_vectopt_loop 
    464464            END DO 
    465465#  endif 
  • trunk/NEMO/TOP_SRC/TRP/trcbbl_adv.h90

    r1152 r1197  
    110110      ! mbathy= number of w-level, minimum value=1 (cf dommsk.F) 
    111111 
    112 #if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     112#if defined key_vectopt_loop 
    113113      jj = 1 
    114114      DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    124124            zunb(ji,jj) = un(ji,jj,mbku(ji,jj))  
    125125            zvnb(ji,jj) = vn(ji,jj,mbkv(ji,jj))  
    126 #if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     126#if ! defined key_vectopt_loop 
    127127         END DO 
    128128#endif 
     
    260260      ! partial steps correction    
    261261       
    262 #if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     262#if defined key_vectopt_loop 
    263263       jj = 1 
    264264       DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    280280                v_trc_bbl(ji,jj,ikv) = zalphay(ji,jj) * vn(ji,jj,ikv) * ze3v / fse3v(ji,jj,ikv)        
    281281             ENDIF 
    282 #if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     282#if ! defined key_vectopt_loop 
    283283          END DO 
    284284#endif 
     
    291291       ! if not partial step loop over the whole domain no lbc call 
    292292 
    293 #if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     293#if defined key_vectopt_loop 
    294294      jj = 1 
    295295      DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    304304               v_trc_bbl(ji,jj,ikv) = zalphay(ji,jj) * vn(ji,jj,ikv)        
    305305            ENDIF 
    306 #if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     306#if ! defined key_vectopt_loop 
    307307         END DO 
    308308#endif 
     
    313313       DO jn = 1, jptra 
    314314 
    315 #if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     315#if defined key_vectopt_loop 
    316316          jj = 1 
    317317          DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    322322                ik = mbkt(ji,jj)                               ! index of the bottom ocean T-level 
    323323                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_omp 
     324#if ! defined key_vectopt_loop 
    325325             END DO 
    326326#endif 
     
    333333          ! ... Second order centered tracer flux at u and v-points 
    334334        
    335 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     335# if defined key_vectopt_loop 
    336336          jj = 1 
    337337          DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    349349                zwy(ji,jj) = ( ( zfvj + ABS( zfvj ) ) * ztrb(ji  ,jj  )   & 
    350350                   &          +( zfvj - ABS( zfvj ) ) * ztrb(ji  ,jj+1) ) * 0.5 
    351 #if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     351#if ! defined key_vectopt_loop 
    352352             END DO 
    353353#endif 
    354354          END DO 
    355355      
    356 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     356# if defined key_vectopt_loop 
    357357          jj = 1 
    358358          DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    369369                ! add it to the general tracer trends 
    370370                tra(ji,jj,ik,jn) = tra(ji,jj,ik,jn) + ztra 
    371 #if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     371#if ! defined key_vectopt_loop 
    372372             END DO 
    373373#endif 
     
    409409       IF( ln_zps ) THEN  
    410410      
    411 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     411# if defined key_vectopt_loop 
    412412          jj = 1 
    413413          DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    427427                zwu(ji,jj) = zalphax(ji,jj) * e2u(ji,jj) * ze3u   
    428428                zwv(ji,jj) = zalphay(ji,jj) * e1v(ji,jj) * ze3v 
    429 #if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     429#if ! defined key_vectopt_loop 
    430430             END DO 
    431431#endif 
     
    434434       ELSE 
    435435 
    436 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     436# if defined key_vectopt_loop 
    437437          jj = 1 
    438438          DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    445445                zwu(ji,jj) = zalphax(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,iku)  
    446446                zwv(ji,jj) = zalphay(ji,jj) * e1v(ji,jj) * fse3v(ji,jj,ikv)  
    447 #if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     447#if ! defined key_vectopt_loop  
    448448             END DO 
    449449#endif 
     
    451451       ENDIF 
    452452 
    453 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     453# if defined key_vectopt_loop 
    454454       jj = 1 
    455455       DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    467467                &   ) / zbt 
    468468 
    469 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     469# if ! defined key_vectopt_loop 
    470470          END DO 
    471471# endif 
  • trunk/NEMO/TOP_SRC/TRP/zpshde_trc.F90

    r1152 r1197  
    125125      DO jn = 1, jptra 
    126126         ! Interpolation of passive tracers at the last ocean level 
    127 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     127# if defined key_vectopt_loop 
    128128         jj = 1 
    129129         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolled) 
     
    170170                  pgtrv(ji,jj,jn) = vmask(ji,jj,1) * ( ptra(ji,jj+1,ikv,jn) - ztraj(ji,jj,jn) ) 
    171171               ENDIF 
    172 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     172# if ! defined key_vectopt_loop 
    173173            END DO 
    174174# endif 
Note: See TracChangeset for help on using the changeset viewer.