4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE clvent(dtime, u1lay, v1lay, coef, t, ven, paprs, pplay, delp, & |
SUBROUTINE clvent(dtime, u1lay, v1lay, coef, cdrag, t, ven, paprs, pplay, & |
8 |
d_ven, flux_v) |
delp, d_ven, flux_v) |
9 |
|
|
10 |
! Author: Z. X. Li (LMD/CNRS) |
! Author: Z. X. Li (LMD/CNRS) |
11 |
! Date: 1993/08/18 |
! Date: 1993/08/18 |
20 |
REAL, intent(in):: u1lay(:), v1lay(:) ! (knon) |
REAL, intent(in):: u1lay(:), v1lay(:) ! (knon) |
21 |
! vent de la premiere couche (m/s) |
! vent de la premiere couche (m/s) |
22 |
|
|
23 |
REAL, intent(in):: coef(:, :) ! (knon, klev) |
REAL, intent(in):: coef(:, 2:) ! (knon, 2:klev) |
24 |
! Coefficient d'echange (m**2/s) multiplié par le cisaillement du |
! Coefficient d'echange (m**2/s) multiplié par le cisaillement du |
25 |
! vent (dV/dz). La première valeur indique la valeur de Cdrag (sans |
! vent (dV/dz) |
26 |
! unité). |
|
27 |
|
REAL, intent(in):: cdrag(:) ! (knon) sans unité |
28 |
|
|
29 |
REAL, intent(in):: t(:, :) ! (knon, klev) ! temperature (K) |
REAL, intent(in):: t(:, :) ! (knon, klev) ! temperature (K) |
30 |
REAL, intent(in):: ven(:, :) ! (knon, klev) vitesse horizontale (m/s) |
REAL, intent(in):: ven(:, :) ! (knon, klev) vitesse horizontale (m/s) |
54 |
local_ven = ven |
local_ven = ven |
55 |
|
|
56 |
DO i = 1, knon |
DO i = 1, knon |
57 |
zx_coef(i, 1) = coef(i, 1) * (1. + SQRT(u1lay(i)**2 + v1lay(i)**2)) & |
zx_coef(i, 1) = cdrag(i) * (1. + SQRT(u1lay(i)**2 + v1lay(i)**2)) & |
58 |
* pplay(i, 1) / (RD * t(i, 1)) * dtime * RG |
* pplay(i, 1) / (RD * t(i, 1)) * dtime * RG |
59 |
ENDDO |
ENDDO |
60 |
|
|