19 |
|
|
20 |
! I. Musat, 01.07.2002 |
! I. Musat, 01.07.2002 |
21 |
|
|
22 |
use coefcdrag_m, only: coefcdrag |
use cdrag_m, only: cdrag |
23 |
use SUPHEC_M, only: RG |
use SUPHEC_M, only: RG |
24 |
|
|
25 |
INTEGER, intent(in):: klon |
INTEGER, intent(in):: klon |
67 |
|
|
68 |
! Richardson at reference level |
! Richardson at reference level |
69 |
|
|
70 |
CALL coefcdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), & |
CALL cdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), gref(:knon), & |
71 |
gref(:knon), psol(:knon), ts(:knon), qsurf(:knon), rugos(:knon), & |
psol(:knon), ts(:knon), qsurf(:knon), rugos(:knon), cdram, cdrah, pref) |
|
cdram, cdrah, pref) |
|
72 |
|
|
73 |
DO i = 1, knon |
DO i = 1, knon |
74 |
delu(i) = ustar(i) / sqrt(cdram(i)) |
delu(i) = ustar(i) / sqrt(cdram(i)) |