Changeset 1326


Ignore:
Timestamp:
2009-02-20T11:08:16+01:00 (12 years ago)
Author:
cetlod
Message:

suppression of useless CPP key key_mpp_omp, see ticket:350

Location:
trunk/NEMO/OFF_SRC
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OFF_SRC/LDF/ldfeiv.F90

    r1265 r1326  
    3737 
    3838CONTAINS 
    39  
    40 # if defined key_mpp_omp 
    41    !!---------------------------------------------------------------------- 
    42    !!   'key_mpp_omp' :                  OpenMP /  NEC autotasking (j-slab) 
    43    !!---------------------------------------------------------------------- 
    44  
    45    SUBROUTINE ldf_eiv( kt ) 
    46       !!---------------------------------------------------------------------- 
    47       !!                  ***  ROUTINE ldf_eiv  *** 
    48       !! 
    49       !! ** Purpose :   Compute the eddy induced velocity coefficient from the 
    50       !!      growth rate of baroclinic instability. 
    51       !! 
    52       !! ** Method : Specific to the offline model. Computes the horizontal 
    53       !!             values from the vertical value 
    54       !! 
    55       !! History : 
    56       !!   9.0  !  06-03  (O. Aumont)  Free form, F90 
    57       !!---------------------------------------------------------------------- 
    58       !! * Arguments 
    59       INTEGER, INTENT( in ) ::   kt     ! ocean time-step inedx 
    60  
    61       !! * Local declarations 
    62       INTEGER ::   ji, jj, jk           ! dummy loop indices 
    63       !!---------------------------------------------------------------------- 
    64  
    65       IF( kt == nit000 ) THEN 
    66          IF(lwp) WRITE(numout,*) 
    67          IF(lwp) WRITE(numout,*) 'ldf_eiv : eddy induced velocity coefficients' 
    68          IF(lwp) WRITE(numout,*) '~~~~~~~   NEC autotasking / OpenMP : j-slab' 
    69       ENDIF 
    70  
    71       ! Average the diffusive coefficient at u- v- points 
    72       DO jj = 2, jpjm1 
    73          DO ji = fs_2, fs_jpim1   ! vector opt. 
    74             aeiu(ji,jj) = .5 * (aeiw(ji,jj) + aeiw(ji+1,jj  )) 
    75             aeiv(ji,jj) = .5 * (aeiw(ji,jj) + aeiw(ji  ,jj+1)) 
    76          END DO 
    77       END DO 
    78       !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,synchro,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    79  
    80       ! lateral boundary condition on aeiu, aeiv 
    81       CALL lbc_lnk( aeiu, 'U', 1. ) 
    82       CALL lbc_lnk( aeiv, 'V', 1. ) 
    83  
    84    END SUBROUTINE ldf_eiv 
    85  
    86 # else 
    87    !!---------------------------------------------------------------------- 
    88    !!   Default key                                             k-j-i loops 
    89    !!---------------------------------------------------------------------- 
    9039 
    9140   SUBROUTINE ldf_eiv( kt ) 
     
    12978   END SUBROUTINE ldf_eiv 
    13079 
    131 # endif 
    13280 
    13381#else 
  • trunk/NEMO/OFF_SRC/LDF/ldfslp.F90

    r1156 r1326  
    140140 
    141141      IF( ln_zps ) THEN      ! partial steps correction at the bottom ocean level (zps_hde routine) 
    142 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     142# if defined key_vectopt_loop 
    143143         jj = 1 
    144144         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    152152               zgru(ji,jj,iku) = gru(ji,jj)  
    153153               zgrv(ji,jj,ikv) = grv(ji,jj)                
    154 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     154# if ! defined key_vectopt_loop  
    155155            END DO 
    156156# endif 
     
    491491      ! mask for mixed layer 
    492492      DO jk = 1, jpk 
    493 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     493# if defined key_vectopt_loop 
    494494         jj = 1 
    495495         DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    505505                  omlmask(ji,jj,jk) = 0.e0 
    506506               ENDIF 
    507 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     507# if ! defined key_vectopt_loop  
    508508            END DO 
    509509# endif 
     
    523523      zwy(:,jpj) = 0.e0 
    524524      zwy(jpi,:) = 0.e0 
    525 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     525# if defined key_vectopt_loop  
    526526      jj = 1 
    527527      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    536536               &             * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) )   & 
    537537               &             / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 
    538 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     538# if ! defined key_vectopt_loop  
    539539         END DO 
    540540# endif 
     
    544544 
    545545      ! Slope at u points 
    546 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     546# if defined key_vectopt_loop  
    547547      jj = 1 
    548548      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    561561            ! uslpml 
    562562            uslpml (ji,jj) = zau / ( zbu - zeps ) * umask (ji,jj,ik) 
    563 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     563# if ! defined key_vectopt_loop 
    564564         END DO 
    565565# endif 
     
    573573      zwy ( :, jpj) = 0.e0 
    574574      zwy ( jpi, :) = 0.e0 
    575 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     575# if defined key_vectopt_loop 
    576576      jj = 1 
    577577      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    585585               &             * ( pn2(ji,jj,ik) + pn2(ji,jj,ik+1) )   & 
    586586               &             / MAX( tmask(ji,jj,ik) + tmask (ji,jj,ik+1), 1. ) 
    587 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     587# if ! defined key_vectopt_loop 
    588588         END DO 
    589589# endif 
     
    594594 
    595595      ! Slope at v points 
    596 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     596# if defined key_vectopt_loop 
    597597      jj = 1 
    598598      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    611611            ! vslpml 
    612612            vslpml (ji,jj) = zav / ( zbv - zeps ) * vmask (ji,jj,ik) 
    613 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     613# if ! defined key_vectopt_loop 
    614614         END DO 
    615615# endif 
     
    625625      ! Local vertical density gradient evaluated from N^2 
    626626      ! zwy = d/dz(prd)= - mk ( prd ) / grav * pn2 -- at w point 
    627 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     627# if defined key_vectopt_loop 
    628628      jj = 1 
    629629      DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    637637            zwy (ji,jj) = zm05g * pn2 (ji,jj,ik) *     & 
    638638               &             ( prd (ji,jj,ik) + prd (ji,jj,ikm1) + 2. ) 
    639 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     639# if ! defined key_vectopt_loop 
    640640         END DO 
    641641# endif 
     
    643643 
    644644      ! Slope at w point 
    645 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     645# if defined key_vectopt_loop 
    646646      jj = 1 
    647647      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    673673            wslpiml (ji,jj) = zai / ( zbi - zeps) * tmask (ji,jj,ik) 
    674674            wslpjml (ji,jj) = zaj / ( zbj - zeps) * tmask (ji,jj,ik) 
    675 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     675# if ! defined key_vectopt_loop 
    676676         END DO 
    677677# endif 
  • trunk/NEMO/OFF_SRC/TRA/zpshde.F90

    r1152 r1326  
    130130 
    131131      ! Interpolation of T and S at the last ocean level 
    132 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     132# if defined key_vectopt_loop 
    133133         jj = 1 
    134134         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolled) 
     
    193193               pgsv(ji,jj) = vmask(ji,jj,1) * ( psal(ji,jj+1,ikv) - zsj(ji,jj) ) 
    194194            ENDIF 
    195 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     195# if ! defined key_vectopt_loop 
    196196         END DO 
    197197# endif 
     
    205205 
    206206      ! Gradient of density at the last level  
    207 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
     207# if defined key_vectopt_loop  
    208208         jj = 1 
    209209         DO ji = 1, jpij-jpi   ! vector opt. (forced unrolled) 
     
    226226               pgrv(ji,jj) = vmask(ji,jj,1) * ( prd(ji,jj+1,ikv) - zrj(ji,jj) ) 
    227227            ENDIF 
    228 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     228# if ! defined key_vectopt_loop 
    229229         END DO 
    230230# endif 
  • trunk/NEMO/OFF_SRC/ZDF/zdfmxl.F90

    r1152 r1326  
    4444 
    4545CONTAINS 
    46  
    47 # if defined key_mpp_omp 
    48    !!---------------------------------------------------------------------- 
    49    !!   'key_mpp_omp'                                   j-k-i loop (j-slab) 
    50    !!---------------------------------------------------------------------- 
    51  
    52    SUBROUTINE zdf_mxl( kt ) 
    53       !!---------------------------------------------------------------------- 
    54       !!                    ***  ROUTINE zdfmxl  *** 
    55       !!                    
    56       !! ** Purpose :   Compute the turbocline depth and the mixed layer depth 
    57       !!      with a density criteria. 
    58       !! 
    59       !! ** Method  :   The turbocline depth is the depth at which the vertical  
    60       !!      eddy diffusivity coefficient (resulting from the vertical physics 
    61       !!      alone, not the isopycnal part, see trazdf.F) fall below a given 
    62       !!      value defined locally (avt_c here taken equal to 5 cm/s2) 
    63       !! 
    64       !! ** Action  : 
    65       !! 
    66       !!---------------------------------------------------------------------- 
    67       !! * Arguments 
    68       INTEGER, INTENT( in ) ::   kt         ! ocean time-step index 
    69  
    70       !! * Local declarations 
    71       INTEGER ::   ji, jj, jk     ! dummy loop indices 
    72       INTEGER ::   ik             ! temporary integer 
    73       INTEGER, DIMENSION(jpi,jpj) ::   & 
    74          imld                     ! temporary workspace 
    75       !!---------------------------------------------------------------------- 
    76  
    77       IF( kt == nit000 ) THEN 
    78          IF(lwp) WRITE(numout,*) 
    79          IF(lwp) WRITE(numout,*) 'zdf_mxl : mixed layer depth, j-k-i loops' 
    80          IF(lwp) WRITE(numout,*) '~~~~~~~' 
    81       ENDIF 
    82  
    83       !                                                ! =============== 
    84       DO jj = 1, jpj                                   !  Vertical slab 
    85          !                                             ! =============== 
    86  
    87          ! 1. Turbocline depth 
    88          ! ------------------- 
    89          ! last w-level at which avt<avt_c (starting from the bottom jk=jpk) 
    90          ! (since avt(.,.,jpk)=0, we have jpk=< imld =< 2 ) 
    91          DO jk = jpk, 2, -1 
    92             DO ji = 1, jpi 
    93                IF( avt(ji,jj,jk) < avt_c ) imld(ji,jj) = jk  
    94             END DO 
    95          END DO 
    96  
    97          ! Turbocline depth and sub-turbocline temperature 
    98          DO ji = 1, jpi 
    99             ik = imld(ji,jj) 
    100             hmld (ji,jj) = fsdepw(ji,jj,ik) * tmask(ji,jj,1) 
    101          END DO 
    102  
    103 !!gm idea 
    104 !!    
    105 !!gm     DO jk = jpk, 2, -1 
    106 !!gm        DO ji = 1, jpi 
    107 !!gm           IF( avt(ji,jj,jk) < avt_c ) hmld(ji,jj) = fsdepw(ji,jj,jk) * tmask(ji,jj,1) 
    108 !!gm        END DO 
    109 !!gm     END DO 
    110 !!gm 
    111  
    112          ! 2. Mixed layer depth 
    113          ! -------------------- 
    114          ! Initialization to the number of w ocean point mbathy 
    115          nmln(:,jj) = mbathy(:,jj) 
    116  
    117          ! Last w-level at which rhop>=rho surf+rho_c (starting from jpk-1) 
    118          ! (rhop defined at t-point, thus jk-1 for w-level just above) 
    119          DO jk = jpkm1, 2, -1 
    120             DO ji = 1, jpi 
    121                IF( rhop(ji,jj,jk) > rhop(ji,jj,1) + rho_c )   nmln(ji,jj) = jk 
    122             END DO 
    123          END DO 
    124  
    125          ! Mixed layer depth 
    126          DO ji = 1, jpi 
    127             ik = nmln(ji,jj) 
    128             hmlp (ji,jj) = fsdepw(ji,jj,ik) * tmask(ji,jj,1) 
    129             hmlpt(ji,jj) = fsdept(ji,jj,ik-1) 
    130          END DO 
    131          !                                             ! =============== 
    132       END DO                                           !   End of slab 
    133       !                                                ! =============== 
    134  
    135       IF(ln_ctl)   CALL prt_ctl( tab2d_1=REAL(nmln,wp), clinfo1=' nmln : ', tab2d_2=hmld, clinfo2=' hmld : ', ovlap=1 ) 
    136  
    137    END SUBROUTINE zdf_mxl 
    138  
    139 # else 
    140    !!---------------------------------------------------------------------- 
    141    !!   Default option :                                         k-j-i loop 
    142    !!---------------------------------------------------------------------- 
    14346 
    14447   SUBROUTINE zdf_mxl( kt ) 
     
    237140 
    238141   END SUBROUTINE zdf_mxl 
    239 #endif 
    240142 
    241143   !!====================================================================== 
  • trunk/NEMO/OFF_SRC/eosbn2.F90

    r1152 r1326  
    487487         DO jj = 1, jpjm1 
    488488!CDIR NOVERRCHK 
    489 #if defined key_mpp_omp 
    490             DO ji = 1, jpim1 
    491 #else 
    492489            DO ji = 1, fs_jpim1   ! vector opt. 
    493 #endif 
    494490               zws(ji,jj) = SQRT( ABS( psal(ji,jj) ) ) 
    495491            END DO 
     
    499495         DO jj = 1, jpjm1                                 ! Horizontal slab 
    500496            !                                             ! =============== 
    501 #if defined key_mpp_omp 
    502             DO ji = 1, jpim1 
    503 #else 
    504497            DO ji = 1, fs_jpim1   ! vector opt. 
    505 #endif 
    506  
    507498               zmask = tmask(ji,jj,1)      ! land/sea bottom mask = surf. mask 
    508499 
     
    554545         DO jj = 1, jpjm1                                 ! Horizontal slab 
    555546            !                                             ! =============== 
    556 #if defined key_mpp_omp 
    557             DO ji = 1, jpim1 
    558 #else 
    559547            DO ji = 1, fs_jpim1   ! vector opt. 
    560 #endif 
    561548               prd(ji,jj) = ( 0.0285 - ralpha * ptem(ji,jj) ) * tmask(ji,jj,1) 
    562549            END DO 
     
    571558         DO jj = 1, jpjm1                                 ! Horizontal slab 
    572559            !                                             ! =============== 
    573 #if defined key_mpp_omp 
    574             DO ji = 1, jpim1 
    575 #else 
    576560            DO ji = 1, fs_jpim1   ! vector opt. 
    577 #endif 
    578561               prd(ji,jj) = ( rbeta * psal(ji,jj) - ralpha * ptem(ji,jj) ) * tmask(ji,jj,1)  
    579562            END DO 
  • trunk/NEMO/OFF_SRC/opa.F90

    r1325 r1326  
    201201      CALL istate_init                      ! ocean initial state (Dynamics and tracers) 
    202202 
    203 !!!!      CALL trc_ini                           ! Passive tracers 
    204  
    205203      CALL day_init                          ! Calendar 
    206204 
Note: See TracChangeset for help on using the changeset viewer.