source: LMDZ4/trunk/libf/phylmd/ustarhb.F @ 1550

Last change on this file since 1550 was 776, checked in by Laurent Fairhead, 17 years ago

Suite du merge entre la version et la HEAD: quelques modifications
de Yann sur le

LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.9 KB
RevLine 
[541]1!
2! $Header$
3!
4      SUBROUTINE ustarhb(knon,u,v,cd_m, ustar)
[766]5      use dimphy
[541]6      IMPLICIT none
7c======================================================================
8c Laurent Li (LMD/CNRS), le 30 septembre 1998
9c Couche limite non-locale. Adaptation du code du CCM3.
10c Code non teste, donc a ne pas utiliser.
11c======================================================================
12c Nonlocal scheme that determines eddy diffusivities based on a
13c diagnosed boundary layer height and a turbulent velocity scale.
14c Also countergradient effects for heat and moisture are included.
15c
16c For more information, see Holtslag, A.A.M., and B.A. Boville, 1993:
17c Local versus nonlocal boundary-layer diffusion in a global climate
18c model. J. of Climate, vol. 6, 1825-1842.
19c======================================================================
[766]20cym#include "dimensions.h"
21cym#include "dimphy.h"
[541]22#include "YOMCST.h"
23c
24c Arguments:
25c
26      INTEGER knon ! nombre de points a calculer
27      REAL u(klon,klev) ! vitesse U (m/s)
28      REAL v(klon,klev) ! vitesse V (m/s)
29      REAL cd_m(klon) ! coefficient de friction au sol pour vitesse
30      REAL ustar(klon)
31c
32      INTEGER i, k
33      REAL zxt, zxq, zxu, zxv, zxmod, taux, tauy
34      REAL zx_alf1, zx_alf2 ! parametres pour extrapolation
35      LOGICAL unssrf(klon)  ! unstb pbl w/lvls within srf pbl lyr
36      LOGICAL unsout(klon)  ! unstb pbl w/lvls in outer pbl lyr
37      LOGICAL check(klon)   ! True=>chk if Richardson no.>critcal
38c
39#include "YOETHF.h"
40#include "FCTTRE.h"
41      DO i = 1, knon
42        zx_alf1 = 1.0
43        zx_alf2 = 1.0 - zx_alf1
44        zxu = u(i,1)*zx_alf1+u(i,2)*zx_alf2
45        zxv = v(i,1)*zx_alf1+v(i,2)*zx_alf2
46        zxmod = 1.0+SQRT(zxu**2+zxv**2)
47        taux = zxu *zxmod*cd_m(i)
48        tauy = zxv *zxmod*cd_m(i)
49        ustar(i) = SQRT(taux**2+tauy**2)
50c       print*,'Ust ',zxu,zxmod,taux,ustar(i)
51      ENDDO
52c
53      return
54      end
Note: See TracBrowser for help on using the repository browser.