pbl_profile Subroutine

public subroutine pbl_profile(ps, td2m, zml1, t2m, tml1, u10m, uml1, stress, hf)

                                                               *
                G. WOTAWA, 1995-07-07                          *
                                                               *

                                                               *

DESCRIPTION: CALCULATION OF FRICTION VELOCITY AND SURFACE SENS- * IBLE HEAT FLUX USING THE PROFILE METHOD (BERKOVICZ * AND PRAHM, 1982) * * Output now is surface stress instead of ustar * * *

                                                               *

INPUT: * * * ps surface pressure(Pa) * td2m two metre dew point(K) * zml1 heigth of first model level (m) * t2m two metre temperature (K) * tml1 temperature first model level (K) * u10m ten metre wind speed (ms-1) * uml1 wind speed first model level (ms-1) * *

                                                               *

OUTPUT: * * stress surface stress (i.e., friction velocity (ms-1) squared * multiplied with air density) * hf surface sensible heat flux (Wm-2) * *

ustar friction velocity (ms-1) * maxiter maximum number of iterations *

Wind Speed difference between Model level 1 and 10 m

Monin-Obukhov Theory not applicable --> Set dummy values Potential temperature difference between model level 1 and 10 m

Neutral conditions Successive approximation will not converge Successive approximation successful

Arguments