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 14833 – NEMO

Changeset 14833


Ignore:
Timestamp:
2021-05-11T10:50:11+02:00 (3 years ago)
Author:
hadcv
Message:

#2600: Merge in trunk [14820]

Location:
NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE/DYN/divhor.F90

    r14787 r14833  
    8686            &                )                                                             & ! bracket for halo 1 - halo 2 compatibility 
    8787            &              + ( e1v(ji,jj  ) * e3v(ji,jj  ,jk,Kmm) * vv(ji,jj  ,jk,Kmm)     & 
    88             &                - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm)     &  
     88            &                - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm)     & 
    8989            &                )                                                             & ! bracket for halo 1 - halo 2 compatibility 
    9090            &             )  * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 

  • NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE/DYN/dynvor.F90

    r14805 r14833  
    657657         CASE ( 0 )                                   ! original formulation  (masked averaging of e3t divided by 4) 
    658658            DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 
    659                ! NOTE: [halo1-halo2] brackets added to make results independent of nn_hls 
     659               ! round brackets added to fix the order of floating point operations 
     660               ! needed to ensure halo 1 - halo 2 compatibility 
    660661               ze3f = (  (e3t(ji  ,jj+1,jk,Kmm)*tmask(ji  ,jj+1,jk)    & 
    661662                  &    +  e3t(ji+1,jj+1,jk,Kmm)*tmask(ji+1,jj+1,jk))   & 
    … …  
    668669         CASE ( 1 )                                   ! new formulation  (masked averaging of e3t divided by the sum of mask) 
    669670            DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 
    670                ! NOTE: [halo1-halo2] brackets added to make results independent of nn_hls 
     671               ! round brackets added to fix the order of floating point operations 
     672               ! needed to ensure halo 1 - halo 2 compatibility 
    671673               ze3f = (  (e3t(ji  ,jj+1,jk,Kmm)*tmask(ji  ,jj+1,jk)    & 
    672674                  &    +  e3t(ji+1,jj+1,jk,Kmm)*tmask(ji+1,jj+1,jk))   & 
    … …  
    847849         CASE ( np_RVO )                           !* relative vorticity 
    848850            DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 
    849                ! NOTE: [halo1-halo2] brackets added to make results independent of nn_hls 
     851               ! round brackets added to fix the order of floating point operations 
     852               ! needed to ensure halo 1 - halo 2 compatibility 
    850853               zwz(ji,jj,jk) = (  (e2v(ji+1,jj  ) * pv(ji+1,jj  ,jk) - e2v(ji,jj) * pv(ji,jj,jk))    & 
    851854                  &             - (e1u(ji  ,jj+1) * pu(ji  ,jj+1,jk) - e1u(ji,jj) * pu(ji,jj,jk))  ) & 
    … …  
    864867         CASE ( np_CRV )                           !* Coriolis + relative vorticity 
    865868            DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 
    866                ! NOTE: [halo1-halo2] brackets added to make results independent of nn_hls 
     869               ! round brackets added to fix the order of floating point operations 
     870               ! needed to ensure halo 1 - halo 2 compatibility 
    867871               zwz(ji,jj,jk) = (  ff_f(ji,jj) + (  (e2v(ji+1,jj  ) * pv(ji+1,jj  ,jk) - e2v(ji,jj) * pv(ji,jj,jk))    & 
    868872                  &                              - (e1u(ji  ,jj+1) * pu(ji  ,jj+1,jk) - e1u(ji,jj) * pu(ji,jj,jk))  ) & 

  • NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE/TRA/traldf_iso.F90

    r14765 r14833  
    183183            zmskv = wmask(ji,jj,jk) / MAX(   vmask(ji,jj  ,jk-1) + vmask(ji,jj-1,jk)          & 
    184184               &                           + vmask(ji,jj-1,jk-1) + vmask(ji,jj  ,jk) , 1._wp  ) 
     185               ! 
    185186            ! round brackets added to fix the order of floating point operations 
    186187            ! needed to ensure halo 1 - halo 2 compatibility 
    … …  
    188189               &       )                                                           & ! bracket for halo 1 - halo 2 compatibility 
    189190               &       + ( pahu(ji-1,jj,jk-1) + pahu(ji  ,jj,jk)                   & 
    190                &         ) ) * zmsku                                                 ! bracket for halo 1 - halo 2 compatibility 
     191               &         )                                                         & ! bracket for halo 1 - halo 2 compatibility 
     192               &     ) * zmsku 
    191193            zahv_w = ( (  pahv(ji,jj  ,jk-1) + pahv(ji,jj-1,jk)                    & 
    192194               &       )                                                           & ! bracket for halo 1 - halo 2 compatibility 
    193195               &       + ( pahv(ji,jj-1,jk-1) + pahv(ji,jj  ,jk)                   & 
    194                &         ) ) * zmskv                                                 ! bracket for halo 1 - halo 2 compatibility 
     196               &         )                                                         & ! bracket for halo 1 - halo 2 compatibility 
     197               &     ) * zmskv 
    195198               ! 
    196199            ah_wslp2(ji,jj,jk) = zahu_w * wslpi(ji,jj,jk) * wslpi(ji,jj,jk)   & 
    … …  
    208211                  &            + ( ( pahv(ji,jj  ,jk) + pahv(ji,jj  ,jk-1) ) / ( e2v(ji,jj  ) * e2v(ji,jj  ) ) & 
    209212                  &              + ( pahv(ji,jj-1,jk) + pahv(ji,jj-1,jk-1) ) / ( e2v(ji,jj-1) * e2v(ji,jj-1) ) & 
    210                   &              ) )                                                                             ! bracket for halo 1 - halo 2 compatibility 
     213                  &              )                                                                             & ! bracket for halo 1 - halo 2 compatibility 
     214                  &                      ) 
    211215            END_3D 
    212216            ! 

  • NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE/TRA/traldf_triad.F90

    r14787 r14833  
    199199                        &                                      + zah1 * r1_e2v(ji,jj-1) * r1_e2v(ji,jj-1) * vmask(ji,jj-1,jk+kp)   & 
    200200                        &                                      )                                                                   ! bracket for halo 1 - halo 2 compatibility 
    201                ! 
    202201            END_3D 
    203202         END DO 
Note: See TracChangeset for help on using the changeset viewer.