trunk/phylmd/ustarhb.f

module ustarhb_m

  IMPLICIT none

contains

  pure function ustarhb(u, v, cd_m)

    ! 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 test\'e, donc \`a ne pas utiliser.

    ! Non-local 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.

    REAL, intent(in):: u(:), v(:) ! wind in first layer (m/s)
    REAL, intent(in):: cd_m(:) ! coefficient de friction au sol pour vitesse
    REAL ustarhb(size(u))

    ! Local:
    INTEGER i
    REAL wind

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

    DO i = 1, size(u)
       wind = SQRT(u(i)**2 + v(i)**2)
       ustarhb(i) = wind * (1 + wind) * cd_m(i)
    ENDDO

  end function ustarhb

end module ustarhb_m
1	guez	62	module ustarhb_m
2
3			IMPLICIT none
4
5			contains
6
7	guez	227	pure function ustarhb(u, v, cd_m)
8	guez	62
9	guez	227	! From LMDZ4/libf/phylmd/ustarhb.F, version 1.1, 2004/06/22 11:45:35
10	guez	62
11			! Laurent Li (LMD/CNRS), le 30 septembre 1998
12			! Couche limite non-locale. Adaptation du code du CCM3.
13	guez	227	! Code non test\'e, donc \`a ne pas utiliser.
14	guez	62
15	guez	238	! Non-local scheme that determines eddy diffusivities based on a
16	guez	62	! diagnosed boundary layer height and a turbulent velocity scale.
17			! Also countergradient effects for heat and moisture are included.
18
19	guez	227	! For more information, see Holtslag, A.A.M. and B.A. Boville, 1993:
20	guez	62	! Local versus nonlocal boundary-layer diffusion in a global climate
21			! model. J. of Climate, vol. 6, 1825-1842.
22
23	guez	227	REAL, intent(in):: u(:), v(:) ! wind in first layer (m/s)
24			REAL, intent(in):: cd_m(:) ! coefficient de friction au sol pour vitesse
25			REAL ustarhb(size(u))
26	guez	62
27	guez	227	! Local:
28	guez	105	INTEGER i
29	guez	245	REAL wind
30	guez	62
31			!---------------------------------------------------------------
32
33	guez	227	DO i = 1, size(u)
34	guez	245	wind = SQRT(u(i)2 + v(i)2)
35			ustarhb(i) = wind * (1 + wind) * cd_m(i)
36	guez	62	ENDDO
37
38	guez	227	end function ustarhb
39	guez	62
40			end module ustarhb_m