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 1693 – NEMO

Changeset 1693


Ignore:
Timestamp:
2009-10-30T16:33:58+01:00 (14 years ago)
Author:
smasson
Message:

slight bug when sending surface current in coupled mode, see ticket:575

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

    r1558 r1693  
    11631163               DO jj = 2, jpjm1 
    11641164                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1165                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)   
     1165                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)   
    11661166                     zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)   
    11671167                     zitx1(ji,jj) = 0.25 * ( u_ice(ji+1,jj+1) + u_ice(ji,jj+1)                     & 
    … …  
    11741174               DO jj = 2, jpjm1 
    11751175                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1176                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)   
     1176                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)   
    11771177                     zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)   
    11781178                     zitx1(ji,jj) = 0.25 * ( u_ice(ji-1,jj-1) + u_ice(ji,jj-1)                     & 
    … …  
    11981198               DO jj = 2, jpjm1 
    11991199                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1200                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)   &    
     1200                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)   &    
    12011201                        &         + 0.25 * ( u_ice(ji+1,jj+1) + u_ice(ji,jj+1)                     & 
    12021202                        &                  + u_ice(ji+1,jj  ) + u_ice(ji,jj  )  ) *  fr_i(ji,jj) 
    … …  
    12091209               DO jj = 2, jpjm1 
    12101210                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    1211                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj-1,1) ) * zfr_l(ji,jj)   &    
     1211                     zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)   &    
    12121212                        &         + 0.25 * ( u_ice(ji-1,jj-1) + u_ice(ji,jj-1)                     & 
    12131213                        &                  + u_ice(ji-1,jj  ) + u_ice(ji,jj  )  ) *  fr_i(ji,jj) 
Note: See TracChangeset for help on using the changeset viewer.