Heat and momemtum fluxes according to Kondo

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
   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.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
   REALTYPE                  :: cdd,chd,ced
   REALTYPE                  :: tmp
   REALTYPE, parameter       :: eps=1.0e-12