Ignore:
Timestamp:
2008-04-11T19:05:03+02:00 (13 years ago)
Author:
ctlod
Message:

merge dev_001_SBC branche with the trunk to include the New Surface Module package, see ticket: #113

File:
1 edited

Legend:

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

    r719 r888  
    2121   USE dom_oce         ! ocean space and time domain 
    2222   USE zdf_oce         ! ocean vertical physics 
     23   USE sbc_oce         ! surface boundary condition: ocean 
    2324   USE phycst          ! physical constants 
    24    USE taumod          ! surface stress 
    2525   USE eosbn2          ! equation of state 
    26    USE ocesbc          ! thermohaline fluxes 
    2726   USE zdfddm          ! double diffusion mixing 
    2827   USE in_out_manager  ! I/O manager 
     
    148147   !!---------------------------------------------------------------------- 
    149148   !!   OPA 9.0 , LOCEAN-IPSL   (2005) 
    150    !! $Header$  
     149   !! $Id$ 
    151150   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    152151   !!---------------------------------------------------------------------- 
     
    460459            zBosol(ji,jj) = grav * zthermal * qsr(ji,jj) 
    461460            ! 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) 
    463462            ! 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) 
    465464            ! Surface salinity flux for non-local term 
    466465            ws0(ji,jj) = - ( emp(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1) 
     
    476475            zrhos         = rhop(ji,jj,1) + zflageos * rau0 * ( 1. - tmask(ji,jj,1) )   
    477476            ! 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) ) 
    480479            ztau          = SQRT( ztx * ztx + zty * zty ) 
    481480            zustar(ji,jj) = SQRT( ztau / ( zrhos +  epsln ) ) 
Note: See TracChangeset for help on using the changeset viewer.