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{displaymath}\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}\end{displaymath} (231)

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
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

Karsten Bolding 2012-12-28