Changeset 2344


Ignore:
Timestamp:
2010-10-29T15:25:31+02:00 (10 years ago)
Author:
gm
Message:

v3.3beta: Suppress obsolete key_mpp_omp

File:
1 edited

Legend:

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

    r2287 r2344  
    99   !!   NEMO     0.5  ! 2002-10  (G. Madec)  Free form, F90 
    1010   !!            1.0  ! 2005-10  (A. Beckmann)  correction for s-coordinates 
     11   !!            3.3  ! 2006-10  (C. Harris, G. Nurser)  add ldf_slp_grif (Griffies operator) 
    1112   !!---------------------------------------------------------------------- 
    1213#if   defined key_ldfslp   ||   defined key_esopa 
     
    3435   PUBLIC   ldf_slp         ! routine called by step.F90 
    3536   PUBLIC   ldf_slp_init    ! routine called by opa.F90 
    36    PUBLIC ldf_slp_grif   !              " 
     37   PUBLIC   ldf_slp_grif    ! routine called by step.F90 
    3738 
    3839   LOGICAL , PUBLIC, PARAMETER              ::   lk_ldfslp = .TRUE.   !: slopes flag 
     
    5758   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    5859   !! $Id$ 
    59    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     60   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6061   !!---------------------------------------------------------------------- 
    61  
    6262CONTAINS 
    6363 
     
    348348   END SUBROUTINE ldf_slp 
    349349 
     350 
    350351   SUBROUTINE ldf_slp_grif ( kt ) 
    351352     !!---------------------------------------------------------------------- 
     
    364365     !! ** Action : - alpha, beta 
    365366     !!               wslp2 squared slope of neutral surfaces at w-points. 
     367     !!---------------------------------------------------------------------- 
     368     USE oce, zdit  => ua   ! use ua as workspace 
     369     USE oce, zdis  => va   ! use va as workspace 
     370     USE oce, zdjt  => ta   ! use ta as workspace 
     371     USE oce, zdjs  => sa   ! use sa as workspace 
    366372     !! 
    367      !! History : 
    368      !!   9.0  !  06-10  (C. Harris)  New subroutine 
    369      !!---------------------------------------------------------------------- 
    370      !! * Modules used 
    371      USE oce            , zdit  => ua,  &  ! use ua as workspace 
    372           zdis  => va,  &  ! use va as workspace 
    373           zdjt  => ta,  &  ! use ta as workspace 
    374           zdjs  => sa      ! use sa as workspace 
    375      !! * Arguments 
    376373     INTEGER, INTENT( in ) ::   kt         ! ocean time-step index 
    377  
    378      !! * Local declarations 
     374     !! 
    379375     INTEGER  ::   ji, jj, jk, ip, jp, kp  ! dummy loop indices 
    380      INTEGER  ::   iku, ikv                ! temporary integer 
    381      REAL(wp) ::   & 
    382           zt, zs, zh, zt2, zsp5, zp1t1,   &  ! temporary scalars 
    383           zdenr, zrhotmp, zdndt, zdddt,   &  !     "        " 
    384           zdnds, zddds, znum, zden,       &  !     "        " 
    385           zslope, za_sxe, zslopec, zdsloper,&  !     "        " 
    386           zfact, zepsln, zatempw,zatempu,zatempv, &   !     "        " 
    387           ze1ur,ze2vr,ze3wr,zdxt,zdxs,zdyt,zdys,zdzt,zdzs,zvolf,& 
    388           zr_slpmax,zdxrho,zdyrho,zabs_dzrho 
    389      REAL(wp), DIMENSION(jpi,jpj,jpk,0:1,0:1) ::   & 
    390           zsx,zsy 
    391      REAL(wp), DIMENSION(jpi,jpj,0:1,0:1) ::   & 
    392           zsx_ml_base,zsy_ml_base 
    393      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
    394           zdkt,zdks 
    395      REAL(wp), DIMENSION(jpi,jpj) ::   & 
    396           zr_ml_basew 
     376     INTEGER  ::   iku, ikv                ! local integer 
     377     REAL(wp) ::   zt, zs, zh, zt2, zsp5, zp1t1           ! local scalars 
     378     REAL(wp) ::     zdenr, zrhotmp, zdndt, zdddt         !   -      - 
     379     REAL(wp) ::     zdnds, zddds, znum, zden             !   -      - 
     380     REAL(wp) ::     zslope, za_sxe, zslopec, zdsloper    !   -      - 
     381     REAL(wp) ::     zfact, zepsln, zatempw,zatempu,zatempv    !   -      - 
     382     REAL(wp) ::     ze1ur, ze2vr, ze3wr, zdxt, zdxs, zdyt, zdys, zdzt, zdzs, zvolf  
     383     REAL(wp) ::     zr_slpmax, zdxrho, zdyrho, zabs_dzrho 
     384     REAL(wp), DIMENSION(jpi,jpj,jpk,0:1,0:1) ::   zsx,zsy 
     385     REAL(wp), DIMENSION(jpi,jpj    ,0:1,0:1) ::   zsx_ml_base, zsy_ml_base 
     386     REAL(wp), DIMENSION(jpi,jpj,jpk)         ::   zdkt, zdks 
     387     REAL(wp), DIMENSION(jpi,jpj) ::   zr_ml_basew 
    397388     !!---------------------------------------------------------------------- 
    398389 
     
    444435                      &                               + 0.380374e-04 ) * zh) 
    445436 
    446               ENDDO 
    447            ENDDO 
    448         ENDDO 
     437              END DO 
     438           END DO 
     439        END DO 
    449440 
    450441     CASE ( 1 ) 
     
    521512 
    522513     IF( ln_zps ) THEN      ! partial steps correction at the last level 
    523 # if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    524      jj = 1 
     514# if defined key_vectopt_loop 
     515     DO jj = 1, 1 
    525516        DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
    526517# else 
     
    537528                 zdjt (ji,jj,ikv) = gtsv(ji,jj,jp_tem) 
    538529                 zdjs (ji,jj,ikv) = gtsv(ji,jj,jp_sal) 
    539 # if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    540530              END DO 
    541 # endif 
    542531           END DO 
    543532        ENDIF 
     
    733722        DO jp=0,1 
    734723           DO kp=0,1 
    735  
    736724              DO jk = 1, jpkm1 
    737  
    738725                 DO jj = 1, jpjm1 
    739  
    740726                    DO ji = 1, fs_jpim1   ! vector opt. 
    741727                       ! k index of uppermost point(s) of triad is jk+kp-1 
     
    744730                       zfact = 1 - 1/(1 + (jk+kp-1)/nmln(ji,jj+jp)) 
    745731                       zsy(ji,jj+jp,jk,1-jp,kp) = zfact*zsy(ji,jj+jp,jk,1-jp,kp) + & 
    746                             & (1.0_wp-zfact)*(fsdepw(ji,jj+jp,jk+kp)*zr_ml_basew(ji,jj+jp))*zsy_ml_base(ji,jj+jp,1-jp,kp)  
     732                           & (1.0_wp-zfact)*(fsdepw(ji,jj+jp,jk+kp)*zr_ml_basew(ji,jj+jp))*zsy_ml_base(ji,jj+jp,1-jp,kp)  
    747733                    END DO 
    748734 
     
    756742        DO jp=0,1 
    757743           DO kp=0,1 
    758  
    759744              DO jk = 1, jpkm1 
    760  
    761745                 DO jj = 1, jpjm1 
    762  
    763746                    DO ji = 1, fs_jpim1   ! vector opt. 
    764747 
     
    786769                    END DO 
    787770                 END DO 
    788  
    789771              END DO 
    790772           END DO 
    791773        END DO 
    792774 
    793         tfw(:,:,1)=0.0 
    794         sfw(:,:,1)=0.0 
    795         wslp2(:,:,1)=0.0 
     775        tfw(:,:,1)=0.e0 
     776        sfw(:,:,1)=0.e0 
     777        wslp2(:,:,1)=0.e0 
    796778 
    797779        CALL lbc_lnk( wslp2, 'W', 1. ) 
     
    808790        CALL lbc_lnk( psix_eiv, 'U', -1. ) 
    809791        CALL lbc_lnk( psiy_eiv, 'V', -1. ) 
    810  
    811  
    812       END SUBROUTINE ldf_slp_grif 
     792      ! 
     793   END SUBROUTINE ldf_slp_grif 
    813794 
    814795 
Note: See TracChangeset for help on using the changeset viewer.