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/OFF/dtadyn.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/OFF/dtadyn.F90

    r13286 r13295  
    371371        gdepw(:,:,1,Kmm) = 0.0_wp 
    372372 
    373         DO_3D_11_11( 2, jpk ) 
     373        DO_3D( 1, 1, 1, 1, 2, jpk ) 
    374374          !    zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk))   ! 0 everywhere 
    375375          !    tmask = wmask, ie everywhere expect at jk = mikt 
     
    396396         ! 
    397397         nk_rnf(:,:) = 0                               ! set the number of level over which river runoffs are applied 
    398          DO_2D_11_11 
     398         DO_2D( 1, 1, 1, 1 ) 
    399399            IF( h_rnf(ji,jj) > 0._wp ) THEN 
    400400               jk = 2 
     
    410410         END_2D 
    411411!!st pourquoi on n'utilise pas le gde3w ici plutôt que de faire une boucle ?  
    412          DO_2D_11_11 
     412         DO_2D( 1, 1, 1, 1 ) 
    413413            h_rnf(ji,jj) = 0._wp 
    414414            DO jk = 1, nk_rnf(ji,jj) 
     
    605605      gdepw(:,:,1,Kmm) = 0.0_wp 
    606606      ! 
    607       DO_3D_11_11( 2, jpk ) 
     607      DO_3D( 1, 1, 1, 1, 2, jpk ) 
    608608         zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk)) 
    609609         gdepw(ji,jj,jk,Kmm) = gdepw(ji,jj,jk-1,Kmm) + e3t(ji,jj,jk-1,Kmm) 
     
    686686      ! 
    687687!!st code dupliqué même remarque que plus haut pourquoi ne pas utiliser gdepw ? 
    688       DO_2D_11_11 
     688      DO_2D( 1, 1, 1, 1 ) 
    689689         h_rnf(ji,jj) = 0._wp 
    690690         DO jk = 1, nk_rnf(ji,jj)                           ! recalculates h_rnf to be the depth in metres 
Note: See TracChangeset for help on using the changeset viewer.