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/ICE/icethd.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/ICE/icethd.F90

    r13226 r13295  
    109109         zu_io(:,:) = u_ice(:,:) - ssu_m(:,:) 
    110110         zv_io(:,:) = v_ice(:,:) - ssv_m(:,:) 
    111          DO_2D_00_00 
     111         DO_2D( 0, 0, 0, 0 ) 
    112112            zfric(ji,jj) = rn_cio * ( 0.5_wp *  & 
    113113               &                    (  zu_io(ji,jj) * zu_io(ji,jj) + zu_io(ji-1,jj) * zu_io(ji-1,jj)   & 
     
    115115         END_2D 
    116116      ELSE      !  if no ice dynamics => transmit directly the atmospheric stress to the ocean 
    117          DO_2D_00_00 
     117         DO_2D( 0, 0, 0, 0 ) 
    118118            zfric(ji,jj) = r1_rho0 * SQRT( 0.5_wp *  & 
    119119               &                         (  utau(ji,jj) * utau(ji,jj) + utau(ji-1,jj) * utau(ji-1,jj)   & 
     
    126126      ! Partial computation of forcing for the thermodynamic sea ice model 
    127127      !--------------------------------------------------------------------! 
    128       DO_2D_11_11 
     128      DO_2D( 1, 1, 1, 1 ) 
    129129         rswitch  = tmask(ji,jj,1) * MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi10 ) ) ! 0 if no ice 
    130130         ! 
     
    196196         ! select ice covered grid points 
    197197         npti = 0 ; nptidx(:) = 0 
    198          DO_2D_11_11 
     198         DO_2D( 1, 1, 1, 1 ) 
    199199            IF ( a_i(ji,jj,jl) > epsi10 ) THEN      
    200200               npti         = npti  + 1 
Note: See TracChangeset for help on using the changeset viewer.