Changeset 888 for trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90
- Timestamp:
- 2008-04-11T19:05:03+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90
r719 r888 21 21 USE dom_oce ! ocean space and time domain 22 22 USE zdf_oce ! ocean vertical physics 23 USE sbc_oce ! surface boundary condition: ocean 23 24 USE phycst ! physical constants 24 USE taumod ! surface stress25 25 USE eosbn2 ! equation of state 26 USE ocesbc ! thermohaline fluxes27 26 USE zdfddm ! double diffusion mixing 28 27 USE in_out_manager ! I/O manager … … 148 147 !!---------------------------------------------------------------------- 149 148 !! OPA 9.0 , LOCEAN-IPSL (2005) 150 !! $ Header$149 !! $Id$ 151 150 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 152 151 !!---------------------------------------------------------------------- … … 460 459 zBosol(ji,jj) = grav * zthermal * qsr(ji,jj) 461 460 ! Non radiative surface buoyancy force 462 zBo (ji,jj) = grav * zthermal * ( qt(ji,jj) - qsr(ji,jj)) - grav * zhalin * emp(ji,jj)461 zBo (ji,jj) = grav * zthermal * qns(ji,jj) - grav * zhalin * emp(ji,jj) 463 462 ! Surface Temperature flux for non-local term 464 wt0(ji,jj) = - qt(ji,jj)* ro0cpr * tmask(ji,jj,1)463 wt0(ji,jj) = - ( qsr(ji,jj) + qns(ji,jj) )* ro0cpr * tmask(ji,jj,1) 465 464 ! Surface salinity flux for non-local term 466 465 ws0(ji,jj) = - ( emp(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1) … … 476 475 zrhos = rhop(ji,jj,1) + zflageos * rau0 * ( 1. - tmask(ji,jj,1) ) 477 476 ! Friction velocity (zustar), at T-point : LMD94 eq. 2 478 ztx = 0.5 * ( taux(ji,jj) + taux(ji - 1, jj ) )479 zty = 0.5 * ( tauy(ji,jj) + tauy(ji , jj - 1) )477 ztx = 0.5 * ( utau(ji,jj) + utau(ji - 1, jj ) ) 478 zty = 0.5 * ( vtau(ji,jj) + vtau(ji , jj - 1) ) 480 479 ztau = SQRT( ztx * ztx + zty * zty ) 481 480 zustar(ji,jj) = SQRT( ztau / ( zrhos + epsln ) )
Note: See TracChangeset
for help on using the changeset viewer.