Changeset 2772


Ignore:
Timestamp:
2011-06-01T08:31:07+02:00 (9 years ago)
Author:
sga
Message:

NEMO trunk, bugfix to ldfslp.F90: correct jk indices

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90

    r2715 r2772  
    691691            !                          !- bound the slopes: abs(zw.)<= 1/100 and zb..<0 
    692692            !                                kxz max= ah slope max =< e1 e3 /(pi**2 2 dt) 
    693             zbu = MIN(  zbu , -100._wp* ABS( zau ) , -7.e+3_wp/fse3u(ji,jj,ik)* ABS( zau )  ) 
    694             zbv = MIN(  zbv , -100._wp* ABS( zav ) , -7.e+3_wp/fse3v(ji,jj,ik)* ABS( zav )  ) 
     693            zbu = MIN(  zbu , -100._wp* ABS( zau ) , -7.e+3_wp/fse3u(ji,jj,iku)* ABS( zau )  ) 
     694            zbv = MIN(  zbv , -100._wp* ABS( zav ) , -7.e+3_wp/fse3v(ji,jj,ikv)* ABS( zav )  ) 
    695695            !                          !- Slope at u- & v-points (uslpml, vslpml) 
    696             uslpml(ji,jj) = zau / ( zbu - zeps ) * umask(ji,jj,ik) 
    697             vslpml(ji,jj) = zav / ( zbv - zeps ) * vmask(ji,jj,ik) 
     696            uslpml(ji,jj) = zau / ( zbu - zeps ) * umask(ji,jj,iku) 
     697            vslpml(ji,jj) = zav / ( zbv - zeps ) * vmask(ji,jj,ikv) 
    698698            ! 
    699699            !                    !==   i- & j-slopes at w-points just below the Mixed Layer   ==! 
Note: See TracChangeset for help on using the changeset viewer.