Ignore:
Timestamp:
2007-10-16T16:59:56+02:00 (14 years ago)
Author:
ctlod
Message:

get back to the nemo_v2_3 version for trunk

File:
1 edited

Legend:

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

    • Property svn:keywords changed from Id to Author Date Id Revision
    r714 r719  
    2121   USE dom_oce         ! ocean space and time domain 
    2222   USE zdf_oce         ! ocean vertical physics 
    23    USE sbc_oce         ! surface boundary condition: ocean 
    2423   USE phycst          ! physical constants 
     24   USE taumod          ! surface stress 
    2525   USE eosbn2          ! equation of state 
     26   USE ocesbc          ! thermohaline fluxes 
    2627   USE zdfddm          ! double diffusion mixing 
    2728   USE in_out_manager  ! I/O manager 
     
    147148   !!---------------------------------------------------------------------- 
    148149   !!   OPA 9.0 , LOCEAN-IPSL   (2005) 
    149    !! $Id$ 
     150   !! $Header$  
    150151   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    151152   !!---------------------------------------------------------------------- 
     
    459460            zBosol(ji,jj) = grav * zthermal * qsr(ji,jj) 
    460461            ! Non radiative surface buoyancy force 
    461             zBo   (ji,jj) = grav * zthermal * qns(ji,jj) -  grav * zhalin * emp(ji,jj) 
     462            zBo   (ji,jj) = grav * zthermal * ( qt(ji,jj) - qsr(ji,jj) ) -  grav * zhalin * emp(ji,jj) 
    462463            ! Surface Temperature flux for non-local term 
    463             wt0(ji,jj) = - ( qsr(ji,jj) + qns(ji,jj) )* ro0cpr * tmask(ji,jj,1) 
     464            wt0(ji,jj) = - qt(ji,jj) * ro0cpr * tmask(ji,jj,1) 
    464465            ! Surface salinity flux for non-local term 
    465466            ws0(ji,jj) = - ( emp(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1) 
     
    475476            zrhos         = rhop(ji,jj,1) + zflageos * rau0 * ( 1. - tmask(ji,jj,1) )   
    476477            ! Friction velocity (zustar), at T-point : LMD94 eq. 2 
    477             ztx           = 0.5 * ( utau(ji,jj) + utau(ji - 1, jj    ) ) 
    478             zty           = 0.5 * ( vtau(ji,jj) + vtau(ji    , jj - 1) ) 
     478            ztx           = 0.5 * ( taux(ji,jj) + taux(ji - 1, jj    ) ) 
     479            zty           = 0.5 * ( tauy(ji,jj) + tauy(ji    , jj - 1) ) 
    479480            ztau          = SQRT( ztx * ztx + zty * zty ) 
    480481            zustar(ji,jj) = SQRT( ztau / ( zrhos +  epsln ) ) 
Note: See TracChangeset for help on using the changeset viewer.