libf/phylmd/ustarhb.f

!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ustarhb.F,v 1.1 2004/06/22 11:45:35 lmdzadmin Exp $
!
      SUBROUTINE ustarhb(knon,u,v,cd_m, ustar)
      use dimens_m
      use dimphy
      use SUPHEC_M
      use yoethf_m
      use fcttre
      IMPLICIT none
c======================================================================
c Laurent Li (LMD/CNRS), le 30 septembre 1998
c Couche limite non-locale. Adaptation du code du CCM3.
c Code non teste, donc a ne pas utiliser.
c======================================================================
c Nonlocal scheme that determines eddy diffusivities based on a
c diagnosed boundary layer height and a turbulent velocity scale.
c Also countergradient effects for heat and moisture are included.
c
c For more information, see Holtslag, A.A.M., and B.A. Boville, 1993:
c Local versus nonlocal boundary-layer diffusion in a global climate
c model. J. of Climate, vol. 6, 1825-1842.
c======================================================================
c
c Arguments:
c
      INTEGER knon ! nombre de points a calculer
      REAL u(klon,klev) ! vitesse U (m/s)
      REAL v(klon,klev) ! vitesse V (m/s)
      REAL cd_m(klon) ! coefficient de friction au sol pour vitesse
      REAL ustar(klon)
c
      INTEGER i, k
      REAL zxt, zxq, zxu, zxv, zxmod, taux, tauy
      REAL zx_alf1, zx_alf2 ! parametres pour extrapolation
      LOGICAL unssrf(klon)  ! unstb pbl w/lvls within srf pbl lyr
      LOGICAL unsout(klon)  ! unstb pbl w/lvls in outer pbl lyr
      LOGICAL check(klon)   ! True=>chk if Richardson no.>critcal
c
      DO i = 1, knon
        zx_alf1 = 1.0
        zx_alf2 = 1.0 - zx_alf1
        zxu = u(i,1)*zx_alf1+u(i,2)*zx_alf2
        zxv = v(i,1)*zx_alf1+v(i,2)*zx_alf2
        zxmod = 1.0+SQRT(zxu**2+zxv**2)
        taux = zxu *zxmod*cd_m(i)
        tauy = zxv *zxmod*cd_m(i)
        ustar(i) = SQRT(taux**2+tauy**2)
c       print*,'Ust ',zxu,zxmod,taux,ustar(i)
      ENDDO
c
      return
      end
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ustarhb.F,v 1.1 2004/06/22 11:45:35 lmdzadmin Exp $
3			!
4			SUBROUTINE ustarhb(knon,u,v,cd_m, ustar)
5			use dimens_m
6			use dimphy
7	guez	38	use SUPHEC_M
8			use yoethf_m
9	guez	3	use fcttre
10			IMPLICIT none
11			c======================================================================
12			c Laurent Li (LMD/CNRS), le 30 septembre 1998
13			c Couche limite non-locale. Adaptation du code du CCM3.
14			c Code non teste, donc a ne pas utiliser.
15			c======================================================================
16			c Nonlocal scheme that determines eddy diffusivities based on a
17			c diagnosed boundary layer height and a turbulent velocity scale.
18			c Also countergradient effects for heat and moisture are included.
19			c
20			c For more information, see Holtslag, A.A.M., and B.A. Boville, 1993:
21			c Local versus nonlocal boundary-layer diffusion in a global climate
22			c model. J. of Climate, vol. 6, 1825-1842.
23			c======================================================================
24			c
25			c Arguments:
26			c
27			INTEGER knon ! nombre de points a calculer
28			REAL u(klon,klev) ! vitesse U (m/s)
29			REAL v(klon,klev) ! vitesse V (m/s)
30			REAL cd_m(klon) ! coefficient de friction au sol pour vitesse
31			REAL ustar(klon)
32			c
33			INTEGER i, k
34			REAL zxt, zxq, zxu, zxv, zxmod, taux, tauy
35			REAL zx_alf1, zx_alf2 ! parametres pour extrapolation
36			LOGICAL unssrf(klon) ! unstb pbl w/lvls within srf pbl lyr
37			LOGICAL unsout(klon) ! unstb pbl w/lvls in outer pbl lyr
38			LOGICAL check(klon) ! True=>chk if Richardson no.>critcal
39			c
40			DO i = 1, knon
41			zx_alf1 = 1.0
42			zx_alf2 = 1.0 - zx_alf1
43			zxu = u(i,1)zx_alf1+u(i,2)zx_alf2
44			zxv = v(i,1)zx_alf1+v(i,2)zx_alf2
45			zxmod = 1.0+SQRT(zxu2+zxv2)
46			taux = zxu zxmodcd_m(i)
47			tauy = zxv zxmodcd_m(i)
48			ustar(i) = SQRT(taux2+tauy2)
49			c print*,'Ust ',zxu,zxmod,taux,ustar(i)
50			ENDDO
51			c
52			return
53			end