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 9176 for branches/UKMO/dev_r5518_GO6_package_OMP/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90 – NEMO

Ignore:
Timestamp:
2018-01-04T13:30:03+01:00 (6 years ago)
Author:
andmirek
Message:

#2001: OMP directives

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_OMP/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

    r6793 r9176  
    231231      CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    232232         ! 
     233!$OMP PARALLEL DO PRIVATE(zh, zt, zs, ztm, zn3, zn2, zn1, zn0, zn) 
    233234         DO jk = 1, jpkm1 
    234235            DO jj = 1, jpj 
     
    271272      CASE( 1 )                !==  simplified EOS  ==! 
    272273         ! 
     274!$OMP PARALLEL DO PRIVATE(zt, zs, zh, ztm, zn) 
    273275         DO jk = 1, jpkm1 
    274276            DO jj = 1, jpj 
     
    393395         ! Non-stochastic equation of state 
    394396         ELSE 
     397!$OMP PARALLEL DO PRIVATE(zh, zt, zs, ztm, zn3, zn2, zn1, zn0, zn0, zn) 
    395398            DO jk = 1, jpkm1 
    396399               DO jj = 1, jpj 
     
    435438      CASE( 1 )                !==  simplified EOS  ==! 
    436439         ! 
     440!$OMP PARALLEL DO PRIVATE(zt, zs, zh, ztm, zn) 
    437441         DO jk = 1, jpkm1 
    438442            DO jj = 1, jpj 
     
    493497      CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    494498         ! 
     499!$OMP PARALLEL DO PRIVATE(zh, zt, zs, zn3, zn2, zn1, zn0, zn) 
    495500         DO jj = 1, jpjm1 
    496501            DO ji = 1, fs_jpim1   ! vector opt. 
     
    532537      CASE( 1 )                !==  simplified EOS  ==! 
    533538         ! 
     539!$OMP PARALLEL DO PRIVATE(zt, zs, zh, zn) 
    534540         DO jj = 1, jpjm1 
    535541            DO ji = 1, fs_jpim1   ! vector opt. 
     
    583589      CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    584590         ! 
     591!$OMP PARALLEL DO PRIVATE(zh, zt, zs, ztm, zn3, zn2, zn1, zn0, zn) 
    585592         DO jk = 1, jpkm1 
    586593            DO jj = 1, jpj 
     
    640647      CASE( 1 )                  !==  simplified EOS  ==! 
    641648         ! 
     649!$OMP PARALLEL DO PRIVATE(zt, zs, zh, ztm, zn) 
    642650         DO jk = 1, jpkm1 
    643651            DO jj = 1, jpj 
     
    697705      CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    698706         ! 
     707!$OMP PARALLEL DO PRIVATE(zh, zt, zs, zn3, zn2, zn1, zn0, zn) 
    699708         DO jj = 1, jpjm1 
    700709            DO ji = 1, fs_jpim1   ! vector opt. 
     
    755764      CASE( 1 )                  !==  simplified EOS  ==! 
    756765         ! 
     766!$OMP PARALLEL DO PRIVATE(zt, zs, zh, zn) 
    757767         DO jj = 1, jpjm1 
    758768            DO ji = 1, fs_jpim1   ! vector opt. 
     
    910920      IF( nn_timing == 1 ) CALL timing_start('bn2') 
    911921      ! 
     922!$OMP PARALLEL DO PRIVATE(zrw, zaw, zbw) 
    912923      DO jk = 2, jpkm1           ! interior points only (2=< jk =< jpkm1 ) 
    913924         DO jj = 1, jpj          ! surface and bottom value set to zero one for all in istate.F90 
     
    962973      z1_T0   = 1._wp/40._wp 
    963974      ! 
     975!$OMP PARALLEL DO PRIVATE(zt, zs,ztm, zn, zd) 
    964976      DO jj = 1, jpj 
    965977         DO ji = 1, jpi 
     
    10161028      CASE ( -1, 1 )                !==  CT,SA (TEOS-10 formulation) ==! 
    10171029         ! 
     1030!$OMP PARALLEL DO PRIVATE(zs) 
    10181031         DO jj = 1, jpj 
    10191032            DO ji = 1, jpi 
     
    10231036            END DO 
    10241037         END DO 
    1025          ptf(:,:) = ptf(:,:) * psal(:,:) 
     1038!$OMP PARALLEL DO 
     1039         DO jj = 1, jpj 
     1040            DO ji = 1, jpi 
     1041               ptf(ji,jj) = ptf(ji,jj) * psal(ji,jj) 
     1042            END DO 
     1043         END DO 
    10261044         ! 
    10271045         IF( PRESENT( pdep ) )   ptf(:,:) = ptf(:,:) - 7.53e-4 * pdep(:,:) 
     
    10291047      CASE ( 0 )                     !==  PT,SP (UNESCO formulation)  ==! 
    10301048         ! 
    1031          ptf(:,:) = ( - 0.0575_wp + 1.710523e-3_wp * SQRT( psal(:,:) )   & 
    1032             &                     - 2.154996e-4_wp *       psal(:,:)   ) * psal(:,:) 
     1049!$OMP PARALLEL DO 
     1050         DO jj = 1, jpj 
     1051            DO ji = 1, jpi 
     1052               ptf(ji,jj) = ( - 0.0575_wp + 1.710523e-3_wp * SQRT( psal(ji,jj) )   & 
     1053            &                     - 2.154996e-4_wp *       psal(ji,jj)   ) * psal(ji,jj) 
     1054            END DO 
     1055         END DO 
    10331056            ! 
    10341057         IF( PRESENT( pdep ) )   ptf(:,:) = ptf(:,:) - 7.53e-4 * pdep(:,:) 
     
    11251148      CASE( -1, 0 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    11261149         ! 
     1150!$OMP PARALLEL DO PRIVATE(zh, zt, zs, ztm, zn2, zn1, zn0, zn) 
    11271151         DO jk = 1, jpkm1 
    11281152            DO jj = 1, jpj 
     
    11881212      CASE( 1 )                !==  Vallis (2006) simplified EOS  ==! 
    11891213         ! 
     1214!$OMP PARALLEL DO PRIVATE(zt, zs, zh, ztm, zn) 
    11901215         DO jk = 1, jpkm1 
    11911216            DO jj = 1, jpj 
Note: See TracChangeset for help on using the changeset viewer.