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/OCE/DYN/wet_dry.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/OCE/DYN/wet_dry.F90

    r13237 r13295  
    174174      ! 
    175175      wdmask(:,:) = 1._wp 
    176       DO_2D_01_01 
     176      DO_2D( 0, 1, 0, 1 ) 
    177177         ! 
    178178         IF( tmask(ji,jj,1)        < 0.5_wp )   CYCLE    ! we don't care about land cells 
     
    198198      wdramp(:,:) = min((ht_0(:,:) + psshb1(:,:) - rn_wdmin1)/(rn_wdmin0 - rn_wdmin1),1.0_wp) 
    199199      !jth assume don't need a lbc_lnk here 
    200       DO_2D_10_10 
     200      DO_2D( 1, 0, 1, 0 ) 
    201201         wdrampu(ji,jj) = MIN( wdramp(ji,jj) , wdramp(ji+1,jj) ) 
    202202         wdrampv(ji,jj) = MIN( wdramp(ji,jj) , wdramp(ji,jj+1) ) 
     
    211211         jflag = 0     ! flag indicating if any further iterations are needed 
    212212         ! 
    213          DO_2D_01_01 
     213         DO_2D( 0, 1, 0, 1 ) 
    214214            IF( tmask(ji, jj, 1) < 0.5_wp )   CYCLE  
    215215            IF( ht_0(ji,jj)      > zdepwd )   CYCLE 
     
    307307      zwdlmtv(:,:) = 1._wp 
    308308      ! 
    309       DO_2D_01_01 
     309      DO_2D( 0, 1, 0, 1 ) 
    310310         ! 
    311311         IF( tmask(ji, jj, 1 ) < 0.5_wp) CYCLE   ! we don't care about land cells 
     
    333333         jflag = 0     ! flag indicating if any further iterations are needed 
    334334         ! 
    335          DO_2D_01_01 
     335         DO_2D( 0, 1, 0, 1 ) 
    336336            ! 
    337337            IF( tmask(ji, jj, 1 ) < 0.5_wp )   CYCLE  
Note: See TracChangeset for help on using the changeset viewer.