Changeset 1082


Ignore:
Timestamp:
2008-06-05T16:04:42+02:00 (12 years ago)
Author:
rblod
Message:

Bugs correction in KPP, see tickets #69 and #101

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

    r899 r1082  
    350350                  ze3sqr = 1. / ( fse3w(ji,jj,jk) * fse3w(ji,jj,jk) ) 
    351351                  ! Square of vertical shear  at interfaces 
    352                   zsh2   = 0.5 * 0.5 * ( zdku2 + zdkv2 ) * ze3sqr 
     352                  zsh2   = 0.5 * ( zdku2 + zdkv2 ) * ze3sqr 
    353353                  zri    = MAX( rn2(ji,jj,jk), 0. ) / ( zsh2 + epsln )  
    354354#if defined key_c1d 
     
    461461            zBosol(ji,jj) = grav * zthermal * qsr(ji,jj) 
    462462            ! Non radiative surface buoyancy force 
    463             zBo   (ji,jj) = grav * zthermal * qns(ji,jj) -  grav * zhalin * emp(ji,jj) 
     463            zBo   (ji,jj) = grav * zthermal * qns(ji,jj) -  grav * zhalin * emps(ji,jj) 
    464464            ! Surface Temperature flux for non-local term 
    465465            wt0(ji,jj) = - ( qsr(ji,jj) + qns(ji,jj) )* ro0cpr * tmask(ji,jj,1) 
    466466            ! Surface salinity flux for non-local term 
    467             ws0(ji,jj) = - ( emp(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1) 
     467            ws0(ji,jj) = - ( emps(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1) 
    468468         ENDDO 
    469469      ENDDO 
Note: See TracChangeset for help on using the changeset viewer.