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/tests/WAD/MY_SRC – 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

Location:
NEMO/trunk/tests/WAD/MY_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/tests/WAD/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    7474      !                       !==  grid point position  ==!   (in kilometers) 
    7575      zfact = rn_dx * 1.e-3         ! conversion in km 
    76       DO_2D_11_11 
     76      DO_2D( 1, 1, 1, 1 ) 
    7777         !                       ! longitude 
    7878         plamt(ji,jj) = zfact * (  - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp )  )   

  • NEMO/trunk/tests/WAD/MY_SRC/usrdef_istate.F90

    r12740 r13295  
    176176      ! Apply minimum wetdepth criterion 
    177177      ! 
    178       DO_2D_11_11 
     178      DO_2D( 1, 1, 1, 1 ) 
    179179         IF( ht_0(ji,jj) + pssh(ji,jj) < rn_wdmin1 ) THEN 
    180180            pssh(ji,jj) = ptmask(ji,jj,1)*( rn_wdmin1 - ht_0(ji,jj) ) 

  • NEMO/trunk/tests/WAD/MY_SRC/usrdef_zgr.F90

    r13286 r13295  
    243243      ! at v-point: averaging zht 
    244244      zhv = 0._wp 
    245       DO_2D_00_00 
     245      DO_2D( 0, 0, 0, 0 ) 
    246246         zhv(ji,jj) = 0.5_wp * ( zht(ji,jj) + zht(ji,jj+1) ) 
    247247      END_2D 
    … …  
    280280         ht_0 = zht 
    281281         k_bot(:,:) = jpkm1 * k_top(:,:)  !* bottom ocean = jpk-1 (here use k_top as a land mask) 
    282          DO_2D_11_11 
     282         DO_2D( 1, 1, 1, 1 ) 
    283283            IF( zht(ji,jj) <= -(rn_wdld - rn_wdmin2)) THEN 
    284284               k_bot(ji,jj) = 0 
    … …  
    289289         !                                !* terrain-following coordinate with e3.(k)=cst) 
    290290         !                                !  OVERFLOW case : identical with j-index (T=V, U=F) 
    291          DO_2D_00_00 
     291         DO_2D( 0, 0, 0, 0 ) 
    292292            z1_jpkm1 = 1._wp / REAL( k_bot(ji,jj) - k_top(ji,jj) + 1 , wp) 
    293293            DO jk = 1, jpk 
Note: See TracChangeset for help on using the changeset viewer.