Heat and momemtum fluxes according to Kondo (Source File: kondo.F90)

INTERFACE:

    subroutine kondo(sst,airt,u10,v10,precip,evap,taux,tauy,qe,qh)

DESCRIPTION:

Based on the model sea surface temperature, the wind vector at 10 m height, the air pressure at 2 m, the dry air temperature and the air pressure at 2 m, and the relative humidity (either directly given or recalculated from the wet bulb or the dew point temperature), this routine first computes the transfer coefficients for the surface momentum flux vector, $(\tau_x^s,\tau_y^s)$ ($c_{dd}$), the latent heat flux, $Q_e$, ($c_{ed}$) and the sensible heat flux, $Q_h$, ($c_{hd}$) heat flux according to the Kondo (1975) bulk formulae. Afterwards, these fluxes are calculated according to the following formulae:

$$\begin{array}{rcl} \tau_x^s &=& c_{dd} \rho_a W_x W \\ \\ \ta... ...q_a) \\ \\ Q_h &=& c_{hd} C_{pa} \rho_a W (T_w-T_a) \end{array}$$ (229)

with the air density $\rho_a$, the wind speed at 10 m, $W$, the $x$- and the $y$-component of the wind velocity vector, $W_x$ and $W_y$, respectively, the specific evaporation heat of sea water, $L$, the specific saturation humidity, $q_s$, the actual specific humidity $q_a$, the specific heat capacity of air at constant pressure, $C_{pa}$, the sea surface temperature, $T_w$ and the dry air temperature, $T_a$.

USES:

    use airsea_variables, only: kelvin,const06,rgas,rho_0
    use airsea_variables, only: qs,qa,rhoa
    use airsea_variables, only: cpa,cpw
    use airsea, only: rain_impact,calc_evaporation
    IMPLICIT NONE

INPUT PARAMETERS:

    REALTYPE, intent(in)                :: sst,airt,u10,v10,precip

INPUT/OUTPUT PARAMETERS:

    REALTYPE, intent(inout)             :: evap

OUTPUT PARAMETERS:

    REALTYPE, intent(out)               :: taux,tauy,qe,qh

REVISION HISTORY:

    Original author(s): Hans Burchard and Karsten Bolding
    $Log: kondo.F90,v $
    Revision 1.2  2007-12-21 12:38:03  kb
    added precip/evap to kondo + cleaned
    Revision 1.1  2007-09-25 10:06:10  kbk
    modularized the airsea module - added Fairall method

LOCAL VARIABLES:

    REALTYPE                  :: w,L
    REALTYPE                  :: s,s0
    REALTYPE                  :: ae_d,be_d,pe_d
    REALTYPE                  :: ae_h,be_h,ce_h,pe_h
    REALTYPE                  :: ae_e,be_e,ce_e,pe_e
    REALTYPE                  :: x,x1,x2,x3
    REALTYPE                  :: ta,ta_k,tw,tw_k
    REALTYPE                  :: cdd,chd,ced
    REALTYPE                  :: tmp,rainfall,cd_rain
    REALTYPE, parameter       :: eps=1.0e-12