source: LMDZ6/branches/Amaury_dev/libf/phylmd/ustarhb.F90 @ 5220

Last change on this file since 5220 was 5105, checked in by abarral, 4 months ago

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.6 KB
Line 
1
2! $Header$
3
4SUBROUTINE ustarhb(klon, klev, knon, u, v, cd_m, ustar)
5  IMPLICIT NONE
6  ! ======================================================================
7  ! Laurent Li (LMD/CNRS), le 30 septembre 1998
8  ! Couche limite non-locale. Adaptation du code du CCM3.
9  ! Code non teste, donc a ne pas utiliser.
10  ! ======================================================================
11  ! Nonlocal scheme that determines eddy diffusivities based on a
12  ! diagnosed boundary layer height and a turbulent velocity scale.
13  ! Also countergradient effects for heat and moisture are included.
14
15  ! For more information, see Holtslag, A.A.M., and B.A. Boville, 1993:
16  ! Local versus nonlocal boundary-layer diffusion in a global climate
17  ! model. J. of Climate, vol. 6, 1825-1842.
18  ! ======================================================================
19
20  ! Arguments:
21
22  INTEGER, INTENT(IN) :: klon, klev, knon ! nombre de points a calculer
23  REAL, DIMENSION(klon, klev), INTENT(IN) :: u,v ! vent horizontal (m/s)
24  REAL, DIMENSION(klon), INTENT(IN) :: cd_m ! coefficient de friction au sol pour vitesse
25  REAL, DIMENSION(klon), INTENT(OUT) :: ustar
26
27  INTEGER :: i, k
28  REAL :: zxt, zxq, zxu, zxv, zxmod, taux, tauy
29  REAL :: zx_alf1, zx_alf2 ! parametres pour extrapolation
30
31  DO i = 1, knon
32    zx_alf1 = 1.0
33    zx_alf2 = 1.0 - zx_alf1
34    zxu = u(i, 1)*zx_alf1 + u(i, 2)*zx_alf2
35    zxv = v(i, 1)*zx_alf1 + v(i, 2)*zx_alf2
36    zxmod = 1.0 + sqrt(zxu**2+zxv**2)
37    taux = zxu*zxmod*cd_m(i)
38    tauy = zxv*zxmod*cd_m(i)
39    ustar(i) = sqrt(taux**2+tauy**2)
40  END DO
41
42
43END SUBROUTINE ustarhb
Note: See TracBrowser for help on using the repository browser.