trunk/phylmd/ustarhb.f90

module ustarhb_m

  IMPLICIT none

contains

  SUBROUTINE ustarhb(knon, u, v, cd_m, ustar)

    ! From LMDZ4/libf/phylmd/ustarhb.F, version 1.1 2004/06/22 11:45:35

    ! Laurent Li (LMD/CNRS), le 30 septembre 1998
    ! Couche limite non-locale. Adaptation du code du CCM3.
    ! Code non teste, donc a ne pas utiliser.

    ! Nonlocal scheme that determines eddy diffusivities based on a
    ! diagnosed boundary layer height and a turbulent velocity scale.
    ! Also countergradient effects for heat and moisture are included.

    ! For more information, see Holtslag, A.A.M., and B.A. Boville, 1993:
    ! Local versus nonlocal boundary-layer diffusion in a global climate
    ! model. J. of Climate, vol. 6, 1825-1842.

    USE dimphy, ONLY: klev, klon

    ! Arguments:

    INTEGER knon ! nombre de points a calculer
    REAL u(klon, klev) ! vitesse U (m/s)
    REAL v(klon, klev) ! vitesse V (m/s)

    REAL, intent(in):: cd_m(:) 
    ! (knon) coefficient de friction au sol pour vitesse

    REAL ustar(klon)

    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

    !---------------------------------------------------------------

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

  end SUBROUTINE ustarhb

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