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	!
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	use SUPHEC_M
8	use yoethf_m
9	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