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 13295 for NEMO/trunk/src/SWE/dynkeg.F90 – NEMO

Ignore:
Timestamp:
2020-07-10T20:24:21+02:00 (4 years ago)
Author:
acc
Message:

Replace do-loop macros in the trunk with alternative forms with greater flexibility for extra halo applications. This alters a lot of routines but does not change any behaviour or results. do_loop_substitute.h90 is greatly simplified by this change. SETTE results are identical to those with the previous revision

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/SWE/dynkeg.F90

    r12983 r13295  
    103103!!an45 to be ADDED : que cas C2 - "wet points only" il suffit de x2 le terme quadratic a la coast (nn_dynkeg_adv = 2) 
    104104      CASE ( nkeg_C2_wpo )                          !--  Standard scheme  --! 
    105          DO_3D_01_01( 1, jpkm1 ) 
     105         DO_3D( 0, 1, 0, 1, 1, jpkm1 ) 
    106106            zu =  (   puu(ji-1,jj  ,jk,Kmm) * puu(ji-1,jj  ,jk,Kmm)   & 
    107107               &    + puu(ji  ,jj  ,jk,Kmm) * puu(ji  ,jj  ,jk,Kmm)   ) * ( 2._wp - umask(ji-1,jj,jk) * umask(ji,jj,jk) ) 
     
    113113      ! 
    114114      CASE ( nkeg_C2 )                          !--  Standard scheme  --! 
    115          DO_3D_01_01( 1, jpkm1 ) 
     115         DO_3D( 0, 1, 0, 1, 1, jpkm1 ) 
    116116            zu =    puu(ji-1,jj  ,jk,Kmm) * puu(ji-1,jj  ,jk,Kmm)   & 
    117117               &  + puu(ji  ,jj  ,jk,Kmm) * puu(ji  ,jj  ,jk,Kmm) 
     
    121121         END_3D 
    122122      CASE ( nkeg_HW )                          !--  Hollingsworth scheme  --! 
    123          DO_3D_00_00( 1, jpkm1 ) 
     123         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    124124            zu = 8._wp * ( puu(ji-1,jj  ,jk,Kmm) * puu(ji-1,jj  ,jk,Kmm)    & 
    125125               &         + puu(ji  ,jj  ,jk,Kmm) * puu(ji  ,jj  ,jk,Kmm) )  & 
     
    137137      END SELECT  
    138138      ! 
    139       DO_3D_00_00( 1, jpkm1 ) 
     139      DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    140140         puu(ji,jj,jk,Krhs) = puu(ji,jj,jk,Krhs) - ( zhke(ji+1,jj  ,jk) - zhke(ji,jj,jk) ) / e1u(ji,jj) 
    141141         pvv(ji,jj,jk,Krhs) = pvv(ji,jj,jk,Krhs) - ( zhke(ji  ,jj+1,jk) - zhke(ji,jj,jk) ) / e2v(ji,jj) 
Note: See TracChangeset for help on using the changeset viewer.