Changeset 13295 for NEMO/trunk/tests


Ignore:
Timestamp:
2020-07-10T20:24:21+02:00 (14 months 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
Files:
26 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/BENCH/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    7474      ! Position coordinates (in grid points) 
    7575      !                          ========== 
    76       DO_2D_11_11 
     76      DO_2D( 1, 1, 1, 1 ) 
    7777          
    7878         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
  • NEMO/trunk/tests/BENCH/MY_SRC/usrdef_istate.F90

    r13286 r13295  
    6969      ! we must define z2d as bellow. 
    7070      ! Once we decide to forget trunk compatibility, we must simply define z2d as: 
    71 !!$      DO_2D_00_00 
     71!!$      DO_2D( 0, 0, 0, 0 ) 
    7272!!$         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 
    7373!!$      END_2D 
    7474      igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 
    7575      igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 
    76       DO_2D_00_00 
     76      DO_2D( 0, 0, 0, 0 ) 
    7777         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 
    7878      END_2D 
     
    8181      pssh(:,:) = z2d(:,:)                                                ! +/- 0.05 m 
    8282      ! 
    83       DO_3D_00_00( 1, jpkm1 ) 
     83      DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    8484         zfact = REAL(jk-1,wp) / REAL(jpk-1,wp)   ! 0 to 1 to add a basic stratification 
    8585         ! temperature choosen to lead to ~50% ice at the beginning if rn_thres_sst = 0.5 
  • NEMO/trunk/tests/BENCH/MY_SRC/usrdef_sbc.F90

    r13286 r13295  
    109109      ! we must define z2d as bellow. 
    110110      ! Once we decide to forget trunk compatibility, we must simply define z2d as: 
    111 !!$      DO_2D_00_00 
     111!!$      DO_2D( 0, 0, 0, 0 ) 
    112112!!$         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 
    113113!!$      END_2D 
    114114      igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 
    115115      igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 
    116       DO_2D_00_00 
     116      DO_2D( 0, 0, 0, 0 ) 
    117117         z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 
    118118      END_2D 
  • NEMO/trunk/tests/CANAL/MY_SRC/diawri.F90

    r12740 r13295  
    156156      CALL iom_put(  "sst", ts(:,:,1,jp_tem,Kmm) )    ! surface temperature 
    157157      IF ( iom_use("sbt") ) THEN 
    158          DO_2D_11_11 
     158         DO_2D( 1, 1, 1, 1 ) 
    159159            ikbot = mbkt(ji,jj) 
    160160            z2d(ji,jj) = ts(ji,jj,ikbot,jp_tem,Kmm) 
     
    166166      CALL iom_put(  "sss", ts(:,:,1,jp_sal,Kmm) )    ! surface salinity 
    167167      IF ( iom_use("sbs") ) THEN 
    168          DO_2D_11_11 
     168         DO_2D( 1, 1, 1, 1 ) 
    169169            ikbot = mbkt(ji,jj) 
    170170            z2d(ji,jj) = ts(ji,jj,ikbot,jp_sal,Kmm) 
     
    176176         zztmp = rho0 * 0.25 
    177177         z2d(:,:) = 0._wp 
    178          DO_2D_00_00 
     178         DO_2D( 0, 0, 0, 0 ) 
    179179            zztmp2 = (  ( rCdU_bot(ji+1,jj)+rCdU_bot(ji  ,jj) ) * uu(ji  ,jj,mbku(ji  ,jj),Kmm)  )**2   & 
    180180               &   + (  ( rCdU_bot(ji  ,jj)+rCdU_bot(ji-1,jj) ) * uu(ji-1,jj,mbku(ji-1,jj),Kmm)  )**2   & 
     
    191191      CALL iom_put(  "ssu", uu(:,:,1,Kmm) )            ! surface i-current 
    192192      IF ( iom_use("sbu") ) THEN 
    193          DO_2D_11_11 
     193         DO_2D( 1, 1, 1, 1 ) 
    194194            ikbot = mbku(ji,jj) 
    195195            z2d(ji,jj) = uu(ji,jj,ikbot,Kmm) 
     
    201201      CALL iom_put(  "ssv", vv(:,:,1,Kmm) )            ! surface j-current 
    202202      IF ( iom_use("sbv") ) THEN 
    203          DO_2D_11_11 
     203         DO_2D( 1, 1, 1, 1 ) 
    204204            ikbot = mbkv(ji,jj) 
    205205            z2d(ji,jj) = vv(ji,jj,ikbot,Kmm) 
     
    231231 
    232232      IF ( iom_use("sstgrad") .OR. iom_use("sstgrad2") ) THEN 
    233          DO_2D_00_00 
     233         DO_2D( 0, 0, 0, 0 ) 
    234234            zztmp  = ts(ji,jj,1,jp_tem,Kmm) 
    235235            zztmpx = ( ts(ji+1,jj,1,jp_tem,Kmm) - zztmp ) * r1_e1u(ji,jj) + ( zztmp - ts(ji-1,jj  ,1,jp_tem,Kmm) ) * r1_e1u(ji-1,jj) 
     
    247247      IF( iom_use("heatc") ) THEN 
    248248         z2d(:,:)  = 0._wp  
    249          DO_3D_11_11( 1, jpkm1 ) 
     249         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    250250            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_tem,Kmm) * tmask(ji,jj,jk) 
    251251         END_3D 
     
    255255      IF( iom_use("saltc") ) THEN 
    256256         z2d(:,:)  = 0._wp  
    257          DO_3D_11_11( 1, jpkm1 ) 
     257         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    258258            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) * tmask(ji,jj,jk) 
    259259         END_3D 
     
    263263      IF( iom_use("salt2c") ) THEN 
    264264         z2d(:,:)  = 0._wp  
    265          DO_3D_11_11( 1, jpkm1 ) 
     265         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    266266            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) * tmask(ji,jj,jk) 
    267267         END_3D 
     
    271271      IF ( iom_use("eken") ) THEN 
    272272         z3d(:,:,jpk) = 0._wp  
    273          DO_3D_00_00( 1, jpkm1 ) 
     273         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    274274            zztmp  = 0.25_wp * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 
    275275            z3d(ji,jj,jk) = zztmp * (  uu(ji-1,jj,jk,Kmm)**2 * e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm)   & 
     
    287287         z3d(1,:, : ) = 0._wp 
    288288         z3d(:,1, : ) = 0._wp 
    289          DO_3D_00_00( 1, jpkm1 ) 
     289         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    290290            z3d(ji,jj,jk) = 0.25_wp * ( uu(ji  ,jj,jk,Kmm) * uu(ji  ,jj,jk,Kmm) * e1e2u(ji  ,jj) * e3u(ji  ,jj,jk,Kmm)  & 
    291291               &                      + uu(ji-1,jj,jk,Kmm) * uu(ji-1,jj,jk,Kmm) * e1e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm)  & 
     
    298298 
    299299         z2d(:,:)  = 0._wp  
    300          DO_3D_11_11( 1, jpkm1 ) 
     300         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    301301            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * z3d(ji,jj,jk) * tmask(ji,jj,jk) 
    302302         END_3D 
     
    310310          
    311311         z3d(:,:,jpk) = 0._wp  
    312          DO_3D_00_00( 1, jpkm1 ) 
     312         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    313313            z3d(ji,jj,jk) = (   e2v(ji+1,jj  ) * vv(ji+1,jj  ,jk,Kmm) - e2v(ji,jj) * vv(ji,jj,jk,Kmm)    & 
    314314               &              - e1u(ji  ,jj+1) * uu(ji  ,jj+1,jk,Kmm) + e1u(ji,jj) * uu(ji,jj,jk,Kmm)  ) * r1_e1e2f(ji,jj) 
     
    317317         CALL iom_put( "relvor", z3d )                  ! relative vorticity 
    318318 
    319          DO_3D_11_11( 1, jpkm1 ) 
     319         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    320320            z3d(ji,jj,jk) = ff_f(ji,jj) + z3d(ji,jj,jk)  
    321321         END_3D 
    322322         CALL iom_put( "absvor", z3d )                  ! absolute vorticity 
    323323 
    324          DO_3D_00_00( 1, jpkm1 ) 
     324         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    325325            ze3  = (  e3t(ji,jj+1,jk,Kmm)*tmask(ji,jj+1,jk) + e3t(ji+1,jj+1,jk,Kmm)*tmask(ji+1,jj+1,jk)   & 
    326326               &    + e3t(ji,jj  ,jk,Kmm)*tmask(ji,jj  ,jk) + e3t(ji+1,jj  ,jk,Kmm)*tmask(ji+1,jj  ,jk)  ) 
     
    348348      IF( iom_use("u_heattr") ) THEN 
    349349         z2d(:,:) = 0._wp  
    350          DO_3D_00_00( 1, jpkm1 ) 
     350         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    351351            z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_tem,Kmm) + ts(ji+1,jj,jk,jp_tem,Kmm) ) 
    352352         END_3D 
     
    357357      IF( iom_use("u_salttr") ) THEN 
    358358         z2d(:,:) = 0.e0  
    359          DO_3D_00_00( 1, jpkm1 ) 
     359         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    360360            z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_sal,Kmm) + ts(ji+1,jj,jk,jp_sal,Kmm) ) 
    361361         END_3D 
     
    375375      IF( iom_use("v_heattr") ) THEN 
    376376         z2d(:,:) = 0.e0  
    377          DO_3D_00_00( 1, jpkm1 ) 
     377         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    378378            z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_tem,Kmm) + ts(ji,jj+1,jk,jp_tem,Kmm) ) 
    379379         END_3D 
     
    384384      IF( iom_use("v_salttr") ) THEN 
    385385         z2d(:,:) = 0._wp  
    386          DO_3D_00_00( 1, jpkm1 ) 
     386         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    387387            z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_sal,Kmm) + ts(ji,jj+1,jk,jp_sal,Kmm) ) 
    388388         END_3D 
     
    393393      IF( iom_use("tosmint") ) THEN 
    394394         z2d(:,:) = 0._wp 
    395          DO_3D_00_00( 1, jpkm1 ) 
     395         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    396396            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) *  ts(ji,jj,jk,jp_tem,Kmm) 
    397397         END_3D 
     
    401401      IF( iom_use("somint") ) THEN 
    402402         z2d(:,:)=0._wp 
    403          DO_3D_00_00( 1, jpkm1 ) 
     403         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    404404            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) 
    405405         END_3D 
  • NEMO/trunk/tests/CANAL/MY_SRC/domvvl.F90

    r13286 r13295  
    190190      gdept(:,:,1,Kbb) = 0.5_wp * e3w(:,:,1,Kbb) 
    191191      gdepw(:,:,1,Kbb) = 0.0_wp 
    192       DO_3D_11_11( 2, jpk ) 
     192      DO_3D( 1, 1, 1, 1, 2, jpk ) 
    193193         !    zcoef = tmask - wmask    ! 0 everywhere tmask = wmask, ie everywhere expect at jk = mikt 
    194194         !                             ! 1 everywhere from mbkt to mikt + 1 or 1 (if no isf) 
     
    238238         ENDIF 
    239239         IF ( ln_vvl_zstar_at_eqtor ) THEN   ! use z-star in vicinity of the Equator 
    240             DO_2D_11_11 
     240            DO_2D( 1, 1, 1, 1 ) 
    241241!!gm  case |gphi| >= 6 degrees is useless   initialized just above by default 
    242242               IF( ABS(gphit(ji,jj)) >= 6.) THEN 
     
    407407         zwu(:,:) = 0._wp 
    408408         zwv(:,:) = 0._wp 
    409          DO_3D_10_10( 1, jpkm1 ) 
     409         DO_3D( 1, 0, 1, 0, 1, jpkm1 ) 
    410410            un_td(ji,jj,jk) = rn_ahe3 * umask(ji,jj,jk) * e2_e1u(ji,jj)           & 
    411411               &            * ( tilde_e3t_b(ji,jj,jk) - tilde_e3t_b(ji+1,jj  ,jk) ) 
     
    415415            zwv(ji,jj) = zwv(ji,jj) + vn_td(ji,jj,jk) 
    416416         END_3D 
    417          DO_2D_11_11 
     417         DO_2D( 1, 1, 1, 1 ) 
    418418            un_td(ji,jj,mbku(ji,jj)) = un_td(ji,jj,mbku(ji,jj)) - zwu(ji,jj) 
    419419            vn_td(ji,jj,mbkv(ji,jj)) = vn_td(ji,jj,mbkv(ji,jj)) - zwv(ji,jj) 
    420420         END_2D 
    421          DO_3D_00_00( 1, jpkm1 ) 
     421         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    422422            tilde_e3t_a(ji,jj,jk) = tilde_e3t_a(ji,jj,jk) + (   un_td(ji-1,jj  ,jk) - un_td(ji,jj,jk)    & 
    423423               &                                          +     vn_td(ji  ,jj-1,jk) - vn_td(ji,jj,jk)    & 
     
    647647      gdepw(:,:,1,Kmm) = 0.0_wp 
    648648      gde3w(:,:,1) = gdept(:,:,1,Kmm) - ssh(:,:,Kmm) 
    649       DO_3D_11_11( 2, jpk ) 
     649      DO_3D( 1, 1, 1, 1, 2, jpk ) 
    650650        !    zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk))   ! 0 everywhere tmask = wmask, ie everywhere expect at jk = mikt 
    651651                                                           ! 1 for jk = mikt 
     
    702702         ! 
    703703      CASE( 'U' )                   !* from T- to U-point : hor. surface weighted mean 
    704          DO_3D_10_10( 1, jpk ) 
     704         DO_3D( 1, 0, 1, 0, 1, jpk ) 
    705705            pe3_out(ji,jj,jk) = 0.5_wp * (  umask(ji,jj,jk) * (1.0_wp - zlnwd) + zlnwd ) * r1_e1e2u(ji,jj)   & 
    706706               &                       * (   e1e2t(ji  ,jj) * ( pe3_in(ji  ,jj,jk) - e3t_0(ji  ,jj,jk) )     & 
     
    711711         ! 
    712712      CASE( 'V' )                   !* from T- to V-point : hor. surface weighted mean 
    713          DO_3D_10_10( 1, jpk ) 
     713         DO_3D( 1, 0, 1, 0, 1, jpk ) 
    714714            pe3_out(ji,jj,jk) = 0.5_wp * ( vmask(ji,jj,jk)  * (1.0_wp - zlnwd) + zlnwd ) * r1_e1e2v(ji,jj)   & 
    715715               &                       * (   e1e2t(ji,jj  ) * ( pe3_in(ji,jj  ,jk) - e3t_0(ji,jj  ,jk) )     & 
     
    720720         ! 
    721721      CASE( 'F' )                   !* from U-point to F-point : hor. surface weighted mean 
    722          DO_3D_10_10( 1, jpk ) 
     722         DO_3D( 1, 0, 1, 0, 1, jpk ) 
    723723            pe3_out(ji,jj,jk) = 0.5_wp * (  umask(ji,jj,jk) * umask(ji,jj+1,jk) * (1.0_wp - zlnwd) + zlnwd ) & 
    724724               &                       *    r1_e1e2f(ji,jj)                                                  & 
     
    887887                  ssh(:,:,Kbb) = -ssh_ref 
    888888 
    889                   DO_2D_11_11 
     889                  DO_2D( 1, 1, 1, 1 ) 
    890890                     IF( ht_0(ji,jj)-ssh_ref <  rn_wdmin1 ) THEN ! if total depth is less than min depth 
    891891                        ssh(ji,jj,Kbb) = rn_wdmin1 - (ht_0(ji,jj) ) 
     
    903903               e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) 
    904904 
    905                DO_2D_11_11 
     905               DO_2D( 1, 1, 1, 1 ) 
    906906                  IF ( ht_0(ji,jj) .LE. 0.0 .AND. NINT( ssmask(ji,jj) ) .EQ. 1) THEN 
    907907                     CALL ctl_stop( 'dom_vvl_rst: ht_0 must be positive at potentially wet points' ) 
  • NEMO/trunk/tests/CANAL/MY_SRC/trazdf.F90

    r12740 r13295  
    156156            IF( l_ldfslp ) THEN            ! isoneutral diffusion: add the contribution  
    157157               IF( ln_traldf_msc  ) THEN     ! MSC iso-neutral operator  
    158                   DO_3D_00_00( 2, jpkm1 ) 
     158                  DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 
    159159                     zwt(ji,jj,jk) = zwt(ji,jj,jk) + akz(ji,jj,jk)   
    160160                  END_3D 
    161161               ELSE                          ! standard or triad iso-neutral operator 
    162                   DO_3D_00_00( 2, jpkm1 ) 
     162                  DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 
    163163                     zwt(ji,jj,jk) = zwt(ji,jj,jk) + ah_wslp2(ji,jj,jk) 
    164164                  END_3D 
     
    168168            ! Diagonal, lower (i), upper (s)  (including the bottom boundary condition since avt is masked) 
    169169            IF( ln_zad_Aimp ) THEN         ! Adaptive implicit vertical advection 
    170                DO_3D_00_00( 1, jpkm1 ) 
     170               DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    171171                  zzwi = - p2dt * zwt(ji,jj,jk  ) / e3w(ji,jj,jk  ,Kmm) 
    172172                  zzws = - p2dt * zwt(ji,jj,jk+1) / e3w(ji,jj,jk+1,Kmm) 
     
    177177               END_3D 
    178178            ELSE 
    179                DO_3D_00_00( 1, jpkm1 ) 
     179               DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    180180                  zwi(ji,jj,jk) = - p2dt * zwt(ji,jj,jk  ) / e3w(ji,jj,jk,Kmm) 
    181181                  zws(ji,jj,jk) = - p2dt * zwt(ji,jj,jk+1) / e3w(ji,jj,jk+1,Kmm) 
     
    203203            !   used as a work space array: its value is modified. 
    204204            ! 
    205             DO_2D_00_00 
     205            DO_2D( 0, 0, 0, 0 ) 
    206206               zwt(ji,jj,1) = zwd(ji,jj,1) 
    207207            END_2D 
    208             DO_3D_00_00( 2, jpkm1 ) 
     208            DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 
    209209               zwt(ji,jj,jk) = zwd(ji,jj,jk) - zwi(ji,jj,jk) * zws(ji,jj,jk-1) / zwt(ji,jj,jk-1) 
    210210            END_3D 
     
    212212         ENDIF  
    213213         !          
    214          DO_2D_00_00 
     214         DO_2D( 0, 0, 0, 0 ) 
    215215            pt(ji,jj,1,jn,Kaa) = e3t(ji,jj,1,Kbb) * pt(ji,jj,1,jn,Kbb) + p2dt * e3t(ji,jj,1,Kmm) * pt(ji,jj,1,jn,Krhs) 
    216216         END_2D 
    217          DO_3D_00_00( 2, jpkm1 ) 
     217         DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 
    218218            zrhs = e3t(ji,jj,jk,Kbb) * pt(ji,jj,jk,jn,Kbb) + p2dt * e3t(ji,jj,jk,Kmm) * pt(ji,jj,jk,jn,Krhs)   ! zrhs=right hand side 
    219219            pt(ji,jj,jk,jn,Kaa) = zrhs - zwi(ji,jj,jk) / zwt(ji,jj,jk-1) * pt(ji,jj,jk-1,jn,Kaa) 
    220220         END_3D 
    221221         ! 
    222          DO_2D_00_00 
     222         DO_2D( 0, 0, 0, 0 ) 
    223223            pt(ji,jj,jpkm1,jn,Kaa) = pt(ji,jj,jpkm1,jn,Kaa) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1) 
    224224         END_2D 
    225          DO_3DS_00_00( jpk-2, 1, -1 ) 
     225         DO_3DS( 0, 0, 0, 0, jpk-2, 1, -1 ) 
    226226            pt(ji,jj,jk,jn,Kaa) = ( pt(ji,jj,jk,jn,Kaa) - zws(ji,jj,jk) * pt(ji,jj,jk+1,jn,Kaa) )   & 
    227227               &             / zwt(ji,jj,jk) * tmask(ji,jj,jk) 
  • NEMO/trunk/tests/CANAL/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    9090#endif 
    9191          
    92       DO_2D_11_11          
     92      DO_2D( 1, 1, 1, 1 )          
    9393         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    9494         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
  • NEMO/trunk/tests/CANAL/MY_SRC/usrdef_istate.F90

    r12740 r13295  
    166166         pssh(:,1) = - ff_t(:,1) / grav * pu(:,1,1) * e2t(:,1) 
    167167         DO jl=1, jpnj 
    168             DO_2D_00_00 
     168            DO_2D( 0, 0, 0, 0 ) 
    169169               pssh(ji,jj) = pssh(ji,jj-1) - ff_t(ji,jj) / grav * pu(ji,jj,1) * e2t(ji,jj) 
    170170            END_2D 
     
    183183      CASE(4)    ! geostrophic zonal pulse 
    184184    
    185          DO_2D_11_11 
     185         DO_2D( 1, 1, 1, 1 ) 
    186186            IF ( ABS(glamt(ji,jj)) <= zjetx ) THEN 
    187187               zdu = rn_uzonal 
     
    217217         zP0 = rho0 * zf0 * zumax * zlambda * SQRT(EXP(1._wp)/2._wp) 
    218218         ! 
    219          DO_2D_11_11 
     219         DO_2D( 1, 1, 1, 1 ) 
    220220            zx = glamt(ji,jj) * 1.e3 
    221221            zy = gphit(ji,jj) * 1.e3 
     
    248248         ! velocities: 
    249249         za = 2._wp * zP0 / zlambda**2 
    250          DO_2D_00_00 
     250         DO_2D( 0, 0, 0, 0 ) 
    251251            zx = glamu(ji,jj) * 1.e3 
    252252            zy = gphiu(ji,jj) * 1.e3 
     
    263263         END_2D 
    264264         ! 
    265          DO_2D_00_00 
     265         DO_2D( 0, 0, 0, 0 ) 
    266266            zx = glamv(ji,jj) * 1.e3 
    267267            zy = gphiv(ji,jj) * 1.e3 
  • NEMO/trunk/tests/ICE_ADV1D/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    7878      zphi0 = -REAL( (Nj0glo-2)/2, wp) * 1.e-3 * rn_dy 
    7979 
    80       DO_2D_11_11 
     80      DO_2D( 1, 1, 1, 1 ) 
    8181         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    8282         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
  • NEMO/trunk/tests/ICE_ADV2D/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    9090#endif          
    9191 
    92       DO_2D_11_11 
     92      DO_2D( 1, 1, 1, 1 ) 
    9393         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    9494         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
  • NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    9999#endif          
    100100 
    101       DO_2D_11_11 
     101      DO_2D( 1, 1, 1, 1 ) 
    102102         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    103103         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
  • NEMO/trunk/tests/ISOMIP+/MY_SRC/dtatsd.F90

    r12905 r13295  
    191191         ENDIF 
    192192         ! 
    193          DO_2D_11_11 
     193         DO_2D( 1, 1, 1, 1 ) 
    194194            DO jk = 1, jpk                        ! determines the intepolated T-S profiles at each (i,j) points 
    195195               zl = gdept_0(ji,jj,jk) 
     
    224224         ! 
    225225         IF( ln_zps ) THEN                      ! zps-coordinate (partial steps) interpolation at the last ocean level 
    226             DO_2D_11_11 
     226            DO_2D( 1, 1, 1, 1 ) 
    227227               ik = mbkt(ji,jj)  
    228228               IF( ik > 1 ) THEN 
  • NEMO/trunk/tests/ISOMIP+/MY_SRC/eosbn2.F90

    r12905 r13295  
    243243      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    244244         ! 
    245          DO_3D_11_11( 1, jpkm1 ) 
     245         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    246246            ! 
    247247            zh  = pdep(ji,jj,jk) * r1_Z0                                  ! depth 
     
    279279      CASE( np_seos )                !==  simplified EOS  ==! 
    280280         ! 
    281          DO_3D_11_11( 1, jpkm1 ) 
     281         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    282282            zt  = pts  (ji,jj,jk,jp_tem) - 10._wp 
    283283            zs  = pts  (ji,jj,jk,jp_sal) - 35._wp 
     
    294294      CASE( np_leos )                !==  linear ISOMIP EOS  ==! 
    295295         ! 
    296          DO_3D_11_11( 1, jpkm1 ) 
     296         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    297297            zt  = pts  (ji,jj,jk,jp_tem) - (-1._wp) 
    298298            zs  = pts  (ji,jj,jk,jp_sal) - 34.2_wp 
     
    356356            END DO 
    357357            ! 
    358             DO_3D_11_11( 1, jpkm1 ) 
     358            DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    359359               ! 
    360360               ! compute density (2*nn_sto_eos) times: 
     
    406406         ! Non-stochastic equation of state 
    407407         ELSE 
    408             DO_3D_11_11( 1, jpkm1 ) 
     408            DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    409409               ! 
    410410               zh  = pdep(ji,jj,jk) * r1_Z0                                  ! depth 
     
    444444      CASE( np_seos )                !==  simplified EOS  ==! 
    445445         ! 
    446          DO_3D_11_11( 1, jpkm1 ) 
     446         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    447447            zt  = pts  (ji,jj,jk,jp_tem) - 10._wp 
    448448            zs  = pts  (ji,jj,jk,jp_sal) - 35._wp 
     
    462462      CASE( np_leos )                !==  linear ISOMIP EOS  ==! 
    463463         ! 
    464          DO_3D_11_11( 1, jpkm1 ) 
     464         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    465465            zt  = pts  (ji,jj,jk,jp_tem) - (-1._wp) 
    466466            zs  = pts  (ji,jj,jk,jp_sal) - 34.2_wp 
     
    513513      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    514514         ! 
    515          DO_2D_11_11 
     515         DO_2D( 1, 1, 1, 1 ) 
    516516            ! 
    517517            zh  = pdep(ji,jj) * r1_Z0                                  ! depth 
     
    548548      CASE( np_seos )                !==  simplified EOS  ==! 
    549549         ! 
    550          DO_2D_11_11 
     550         DO_2D( 1, 1, 1, 1 ) 
    551551            ! 
    552552            zt    = pts  (ji,jj,jp_tem)  - 10._wp 
     
    564564      CASE( np_leos )                !==  ISOMIP EOS  ==! 
    565565         ! 
    566          DO_2D_11_11 
     566         DO_2D( 1, 1, 1, 1 ) 
    567567            ! 
    568568            zt    = pts  (ji,jj,jp_tem)  - (-1._wp) 
     
    611611      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    612612         ! 
    613          DO_3D_11_11( 1, jpkm1 ) 
     613         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    614614            ! 
    615615            zh  = gdept(ji,jj,jk,Kmm) * r1_Z0                                ! depth 
     
    664664      CASE( np_seos )                  !==  simplified EOS  ==! 
    665665         ! 
    666          DO_3D_11_11( 1, jpkm1 ) 
     666         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    667667            zt  = pts (ji,jj,jk,jp_tem) - 10._wp   ! pot. temperature anomaly (t-T0) 
    668668            zs  = pts (ji,jj,jk,jp_sal) - 35._wp   ! abs. salinity anomaly (s-S0) 
     
    680680      CASE( np_leos )                  !==  linear ISOMIP EOS  ==! 
    681681         ! 
    682          DO_3D_11_11( 1, jpkm1 ) 
     682         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    683683            zt  = pts (ji,jj,jk,jp_tem) - (-1._wp) 
    684684            zs  = pts (ji,jj,jk,jp_sal) - 34.2_wp   ! abs. salinity anomaly (s-S0) 
     
    734734      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    735735         ! 
    736          DO_2D_11_11 
     736         DO_2D( 1, 1, 1, 1 ) 
    737737            ! 
    738738            zh  = pdep(ji,jj) * r1_Z0                                  ! depth 
     
    787787      CASE( np_seos )                  !==  simplified EOS  ==! 
    788788         ! 
    789          DO_2D_11_11 
     789         DO_2D( 1, 1, 1, 1 ) 
    790790            ! 
    791791            zt    = pts  (ji,jj,jp_tem) - 10._wp   ! pot. temperature anomaly (t-T0) 
     
    803803      CASE( np_leos )                  !==  linear ISOMIP EOS  ==! 
    804804         ! 
    805          DO_2D_11_11 
     805         DO_2D( 1, 1, 1, 1 ) 
    806806            ! 
    807807            zt    = pts  (ji,jj,jp_tem) - (-1._wp)   ! pot. temperature anomaly (t-T0) 
     
    965965      IF( ln_timing )   CALL timing_start('bn2') 
    966966      ! 
    967       DO_3D_11_11( 2, jpkm1 ) 
     967      DO_3D( 1, 1, 1, 1, 2, jpkm1 ) 
    968968         zrw =   ( gdepw(ji,jj,jk  ,Kmm) - gdept(ji,jj,jk,Kmm) )   & 
    969969            &  / ( gdept(ji,jj,jk-1,Kmm) - gdept(ji,jj,jk,Kmm) )  
     
    10131013      z1_T0   = 1._wp/40._wp 
    10141014      ! 
    1015       DO_2D_11_11 
     1015      DO_2D( 1, 1, 1, 1 ) 
    10161016         ! 
    10171017         zt  = ctmp   (ji,jj) * z1_T0 
     
    10661066         ! 
    10671067         z1_S0 = 1._wp / 35.16504_wp 
    1068          DO_2D_11_11 
     1068         DO_2D( 1, 1, 1, 1 ) 
    10691069            zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 )           ! square root salinity 
    10701070            ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 
     
    11731173      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    11741174         ! 
    1175          DO_3D_11_11( 1, jpkm1 ) 
     1175         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    11761176            ! 
    11771177            zh  = gdept(ji,jj,jk,Kmm) * r1_Z0                                ! depth 
     
    12321232      CASE( np_seos )                !==  Vallis (2006) simplified EOS  ==! 
    12331233         ! 
    1234          DO_3D_11_11( 1, jpkm1 ) 
     1234         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    12351235            zt  = pts(ji,jj,jk,jp_tem) - 10._wp  ! temperature anomaly (t-T0) 
    12361236            zs = pts (ji,jj,jk,jp_sal) - 35._wp  ! abs. salinity anomaly (s-S0) 
     
    12481248      CASE( np_leos )                !==  linear ISOMIP EOS  ==! 
    12491249         ! 
    1250          DO_3D_11_11( 1, jpkm1 ) 
     1250         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    12511251            zt  = pts(ji,jj,jk,jp_tem) - (-1._wp)  ! temperature anomaly (t-T0) 
    12521252            zs = pts (ji,jj,jk,jp_sal) - 34.2_wp   ! abs. salinity anomaly (s-S0) 
  • NEMO/trunk/tests/ISOMIP+/MY_SRC/istate.F90

    r12905 r13295  
    105105               ! Apply minimum wetdepth criterion 
    106106               ! 
    107                DO_2D_11_11 
     107               DO_2D( 1, 1, 1, 1 ) 
    108108                  IF( ht_0(ji,jj) + ssh(ji,jj,Kbb)  < rn_wdmin1 ) THEN 
    109109                     ssh(ji,jj,Kbb) = tmask(ji,jj,1)*( rn_wdmin1 - (ht_0(ji,jj)) ) 
     
    159159      ! 
    160160!!gm  the use of umsak & vmask is not necessary below as uu(:,:,:,Kmm), vv(:,:,:,Kmm), uu(:,:,:,Kbb), vv(:,:,:,Kbb) are always masked 
    161       DO_3D_11_11( 1, jpkm1 ) 
     161      DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    162162         uu_b(ji,jj,Kmm) = uu_b(ji,jj,Kmm) + e3u(ji,jj,jk,Kmm) * uu(ji,jj,jk,Kmm) * umask(ji,jj,jk) 
    163163         vv_b(ji,jj,Kmm) = vv_b(ji,jj,Kmm) + e3v(ji,jj,jk,Kmm) * vv(ji,jj,jk,Kmm) * vmask(ji,jj,jk) 
  • NEMO/trunk/tests/ISOMIP+/MY_SRC/tradmp.F90

    r13286 r13295  
    112112      CASE( 0 )                        !*  newtonian damping throughout the water column  *! 
    113113         DO jn = 1, jpts 
    114             DO_3D_00_00( 1, jpkm1 ) 
     114            DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    115115               pts(ji,jj,jk,jn,Krhs) = pts(ji,jj,jk,jn,Krhs)           & 
    116116                  &                  + resto(ji,jj,jk) * ( zts_dta(ji,jj,jk,jn) - pts(ji,jj,jk,jn,Kbb) ) 
     
    119119         ! 
    120120      CASE ( 1 )                       !*  no damping in the turbocline (avt > 5 cm2/s)  *! 
    121          DO_3D_00_00( 1, jpkm1 ) 
     121         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    122122            IF( avt(ji,jj,jk) <= avt_c ) THEN 
    123123               pts(ji,jj,jk,jp_tem,Krhs) = pts(ji,jj,jk,jp_tem,Krhs)   & 
     
    129129         ! 
    130130      CASE ( 2 )                       !*  no damping in the mixed layer   *! 
    131          DO_3D_00_00( 1, jpkm1 ) 
     131         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    132132            IF( gdept(ji,jj,jk,Kmm) >= hmlp (ji,jj) ) THEN 
    133133               pts(ji,jj,jk,jp_tem,Krhs) = pts(ji,jj,jk,jp_tem,Krhs)   & 
  • NEMO/trunk/tests/ISOMIP/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    7777      ! 
    7878      !                       !==  grid point position  ==!   (in degrees) 
    79       DO_2D_11_11 
     79      DO_2D( 1, 1, 1, 1 ) 
    8080         !                       ! longitude   (west coast at lon=0°) 
    8181         plamt(ji,jj) = rn_e1deg * (  - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp )  )   
     
    9191      ! 
    9292      !                       !==  Horizontal scale factors  ==!   (in meters) 
    93       DO_2D_11_11 
     93      DO_2D( 1, 1, 1, 1 ) 
    9494         !                       ! e1   (zonal) 
    9595         pe1t(ji,jj) = ra * rad * COS( rad * pphit(ji,jj) ) * rn_e1deg 
  • NEMO/trunk/tests/ISOMIP/MY_SRC/usrdef_zgr.F90

    r13286 r13295  
    126126         END DO 
    127127         ! top scale factors and depth at T- and W-points 
    128          DO_2D_11_11 
     128         DO_2D( 1, 1, 1, 1 ) 
    129129            ik = k_top(ji,jj) 
    130130            IF ( ik > 2 ) THEN 
     
    145145         END_2D 
    146146         ! bottom scale factors and depth at T- and W-points 
    147          DO_2D_11_11 
     147         DO_2D( 1, 1, 1, 1 ) 
    148148            ik = k_bot(ji,jj) 
    149149            pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) 
     
    158158         pe3u (:,:,:) = pe3t(:,:,:) 
    159159         pe3uw(:,:,:) = pe3w(:,:,:) 
    160          DO_3D_00_00( 1, jpk ) 
     160         DO_3D( 0, 0, 0, 0, 1, jpk ) 
    161161         !                                   ! Computed as the minimum of neighbooring scale factors 
    162162            pe3v (ji,jj,jk) = MIN( pe3t(ji,jj,jk), pe3t(ji,jj+1,jk) ) 
  • NEMO/trunk/tests/LOCK_EXCHANGE/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/OVERFLOW/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/OVERFLOW/MY_SRC/usrdef_zgr.F90

    r13286 r13295  
    184184            pe3vw(:,:,jk) = pe3w_1d (jk) 
    185185         END DO 
    186          DO_2D_11_11 
     186         DO_2D( 1, 1, 1, 1 ) 
    187187            ik = k_bot(ji,jj) 
    188188            pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) 
  • NEMO/trunk/tests/VORTEX/MY_SRC/domvvl.F90

    r13286 r13295  
    190190      gdept(:,:,1,Kbb) = 0.5_wp * e3w(:,:,1,Kbb) 
    191191      gdepw(:,:,1,Kbb) = 0.0_wp 
    192       DO_3D_11_11( 2, jpk ) 
     192      DO_3D( 1, 1, 1, 1, 2, jpk ) 
    193193         !    zcoef = tmask - wmask    ! 0 everywhere tmask = wmask, ie everywhere expect at jk = mikt 
    194194         !                             ! 1 everywhere from mbkt to mikt + 1 or 1 (if no isf) 
     
    238238         ENDIF 
    239239         IF ( ln_vvl_zstar_at_eqtor ) THEN   ! use z-star in vicinity of the Equator 
    240             DO_2D_11_11 
     240            DO_2D( 1, 1, 1, 1 ) 
    241241!!gm  case |gphi| >= 6 degrees is useless   initialized just above by default 
    242242               IF( ABS(gphit(ji,jj)) >= 6.) THEN 
     
    407407         zwu(:,:) = 0._wp 
    408408         zwv(:,:) = 0._wp 
    409          DO_3D_10_10( 1, jpkm1 ) 
     409         DO_3D( 1, 0, 1, 0, 1, jpkm1 ) 
    410410            un_td(ji,jj,jk) = rn_ahe3 * umask(ji,jj,jk) * e2_e1u(ji,jj)           & 
    411411               &            * ( tilde_e3t_b(ji,jj,jk) - tilde_e3t_b(ji+1,jj  ,jk) ) 
     
    415415            zwv(ji,jj) = zwv(ji,jj) + vn_td(ji,jj,jk) 
    416416         END_3D 
    417          DO_2D_11_11 
     417         DO_2D( 1, 1, 1, 1 ) 
    418418            un_td(ji,jj,mbku(ji,jj)) = un_td(ji,jj,mbku(ji,jj)) - zwu(ji,jj) 
    419419            vn_td(ji,jj,mbkv(ji,jj)) = vn_td(ji,jj,mbkv(ji,jj)) - zwv(ji,jj) 
    420420         END_2D 
    421          DO_3D_00_00( 1, jpkm1 ) 
     421         DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    422422            tilde_e3t_a(ji,jj,jk) = tilde_e3t_a(ji,jj,jk) + (   un_td(ji-1,jj  ,jk) - un_td(ji,jj,jk)    & 
    423423               &                                          +     vn_td(ji  ,jj-1,jk) - vn_td(ji,jj,jk)    & 
     
    647647      gdepw(:,:,1,Kmm) = 0.0_wp 
    648648      gde3w(:,:,1) = gdept(:,:,1,Kmm) - ssh(:,:,Kmm) 
    649       DO_3D_11_11( 2, jpk ) 
     649      DO_3D( 1, 1, 1, 1, 2, jpk ) 
    650650        !    zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk))   ! 0 everywhere tmask = wmask, ie everywhere expect at jk = mikt 
    651651                                                           ! 1 for jk = mikt 
     
    702702         ! 
    703703      CASE( 'U' )                   !* from T- to U-point : hor. surface weighted mean 
    704          DO_3D_10_10( 1, jpk ) 
     704         DO_3D( 1, 0, 1, 0, 1, jpk ) 
    705705            pe3_out(ji,jj,jk) = 0.5_wp * (  umask(ji,jj,jk) * (1.0_wp - zlnwd) + zlnwd ) * r1_e1e2u(ji,jj)   & 
    706706               &                       * (   e1e2t(ji  ,jj) * ( pe3_in(ji  ,jj,jk) - e3t_0(ji  ,jj,jk) )     & 
     
    711711         ! 
    712712      CASE( 'V' )                   !* from T- to V-point : hor. surface weighted mean 
    713          DO_3D_10_10( 1, jpk ) 
     713         DO_3D( 1, 0, 1, 0, 1, jpk ) 
    714714            pe3_out(ji,jj,jk) = 0.5_wp * ( vmask(ji,jj,jk)  * (1.0_wp - zlnwd) + zlnwd ) * r1_e1e2v(ji,jj)   & 
    715715               &                       * (   e1e2t(ji,jj  ) * ( pe3_in(ji,jj  ,jk) - e3t_0(ji,jj  ,jk) )     & 
     
    720720         ! 
    721721      CASE( 'F' )                   !* from U-point to F-point : hor. surface weighted mean 
    722          DO_3D_10_10( 1, jpk ) 
     722         DO_3D( 1, 0, 1, 0, 1, jpk ) 
    723723            pe3_out(ji,jj,jk) = 0.5_wp * (  umask(ji,jj,jk) * umask(ji,jj+1,jk) * (1.0_wp - zlnwd) + zlnwd ) & 
    724724               &                       *    r1_e1e2f(ji,jj)                                                  & 
     
    887887                  ssh(:,:,Kbb) = -ssh_ref 
    888888 
    889                   DO_2D_11_11 
     889                  DO_2D( 1, 1, 1, 1 ) 
    890890                     IF( ht_0(ji,jj)-ssh_ref <  rn_wdmin1 ) THEN ! if total depth is less than min depth 
    891891                        ssh(ji,jj,Kbb) = rn_wdmin1 - (ht_0(ji,jj) ) 
     
    903903               e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) 
    904904 
    905                DO_2D_11_11 
     905               DO_2D( 1, 1, 1, 1 ) 
    906906                  IF ( ht_0(ji,jj) .LE. 0.0 .AND. NINT( ssmask(ji,jj) ) .EQ. 1) THEN 
    907907                     CALL ctl_stop( 'dom_vvl_rst: ht_0 must be positive at potentially wet points' ) 
  • NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_hgr.F90

    r13286 r13295  
    9898#endif 
    9999          
    100       DO_2D_11_11 
     100      DO_2D( 1, 1, 1, 1 ) 
    101101         zti = REAL( mig0_oldcmp(ji) - 1, wp )   ! start at i=0 in the global grid without halos 
    102102         ztj = REAL( mjg0_oldcmp(jj) - 1, wp )   ! start at j=0 in the global grid without halos 
  • NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_istate.F90

    r12740 r13295  
    7575      ! Sea level: 
    7676      za = -zP0 * (1._wp-EXP(-zH)) / (grav*(zH-1._wp + EXP(-zH))) 
    77       DO_2D_11_11 
     77      DO_2D( 1, 1, 1, 1 ) 
    7878         zx = glamt(ji,jj) * 1.e3 
    7979         zy = gphit(ji,jj) * 1.e3 
     
    8383      ! 
    8484      ! temperature:          
    85       DO_2D_11_11 
     85      DO_2D( 1, 1, 1, 1 ) 
    8686         zx = glamt(ji,jj) * 1.e3 
    8787         zy = gphit(ji,jj) * 1.e3 
     
    102102      ! velocities: 
    103103      za = 2._wp * zP0 / (zf0 * rho0 * zlambda**2) 
    104       DO_2D_00_00 
     104      DO_2D( 0, 0, 0, 0 ) 
    105105         zx = glamu(ji,jj) * 1.e3 
    106106         zy = gphiu(ji,jj) * 1.e3 
     
    116116      END_2D 
    117117      ! 
    118       DO_2D_00_00 
     118      DO_2D( 0, 0, 0, 0 ) 
    119119         zx = glamv(ji,jj) * 1.e3 
    120120         zy = gphiv(ji,jj) * 1.e3 
  • 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.