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 14215 for NEMO/trunk/src/OCE/SBC/sbcice_cice.F90 – NEMO

Ignore:
Timestamp:
2020-12-18T14:49:22+01:00 (4 years ago)
Author:
acc
Message:

trunk changes to swap the order of arguments to the DO LOOP macros. These changes result in a more natural i-j-k ordering as explained in #2595. SETTE is passed before and after these changes and results are unchanged. This fixes #2595

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/SBC/sbcice_cice.F90

    r14053 r14215  
    312312! x comp of wind stress (CI_1) 
    313313! U point to F point 
    314          DO_2D( 1, 0, 1, 1 ) 
     314         DO_2D( 1, 1, 1, 0 ) 
    315315            ztmp(ji,jj) = 0.5 * (  fr_iu(ji,jj) * utau(ji,jj)      & 
    316316                                 + fr_iu(ji,jj+1) * utau(ji,jj+1) ) * fmask(ji,jj,1) 
     
    320320! y comp of wind stress (CI_2) 
    321321! V point to F point 
    322          DO_2D( 1, 1, 1, 0 ) 
     322         DO_2D( 1, 0, 1, 1 ) 
    323323            ztmp(ji,jj) = 0.5 * (  fr_iv(ji,jj) * vtau(ji,jj)      & 
    324324                                 + fr_iv(ji+1,jj) * vtau(ji+1,jj) ) * fmask(ji,jj,1) 
     
    437437! x comp and y comp of surface ocean current 
    438438! U point to F point 
    439       DO_2D( 1, 0, 1, 1 ) 
     439      DO_2D( 1, 1, 1, 0 ) 
    440440         ztmp(ji,jj)=0.5*(ssu_m(ji,jj)+ssu_m(ji,jj+1))*fmask(ji,jj,1) 
    441441      END_2D 
     
    443443 
    444444! V point to F point 
    445       DO_2D( 1, 1, 1, 0 ) 
     445      DO_2D( 1, 0, 1, 1 ) 
    446446         ztmp(ji,jj)=0.5*(ssv_m(ji,jj)+ssv_m(ji+1,jj))*fmask(ji,jj,1) 
    447447      END_2D 
     
    513513! F point to V point 
    514514 
    515       DO_2D( 1, 0, 0, 0 ) 
     515      DO_2D( 0, 0, 1, 0 ) 
    516516         ss_iov(ji,jj) = 0.5 * ( ztmp1(ji-1,jj) + ztmp1(ji,jj) ) * vmask(ji,jj,1) 
    517517      END_2D 
Note: See TracChangeset for help on using the changeset viewer.