5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE coefkz(nsrf, paprs, pplay, ksta, ksta_ter, ts, rugos, u, v, t, & |
SUBROUTINE coefkz(nsrf, paprs, pplay, ksta, ksta_ter, ts, rugos, u, v, t, & |
8 |
q, qsurf, coefm, coefh) |
q, qsurf, coefm, coefh, ycdragm, ycdragh) |
9 |
|
|
10 |
! Authors: F. Hourdin, M. Forichon, Z. X. Li (LMD/CNRS) |
! Authors: F. Hourdin, M. Forichon, Z. X. Li (LMD/CNRS) |
11 |
! Date: September 22nd, 1993 |
! Date: September 22nd, 1993 |
35 |
REAL, intent(in):: t(:, :) ! (klon, klev) temperature (K) |
REAL, intent(in):: t(:, :) ! (klon, klev) temperature (K) |
36 |
real, intent(in):: q(:, :) ! (klon, klev) vapeur d'eau (kg/kg) |
real, intent(in):: q(:, :) ! (klon, klev) vapeur d'eau (kg/kg) |
37 |
real, intent(in):: qsurf(:) ! (knon) |
real, intent(in):: qsurf(:) ! (knon) |
38 |
REAL, intent(out):: coefm(:, :) ! (knon, klev) coefficient, vitesse |
REAL, intent(out):: coefm(:, 2:) ! (knon, 2:klev) coefficient, vitesse |
39 |
|
|
40 |
real, intent(out):: coefh(:, :) ! (knon, klev) |
real, intent(out):: coefh(:, 2:) ! (knon, 2:klev) |
41 |
! coefficient, chaleur et humidité |
! coefficient, chaleur et humidité |
42 |
|
|
43 |
|
real, intent(out):: ycdragm(:), ycdragh(:) ! (knon) |
44 |
|
|
45 |
! Local: |
! Local: |
46 |
|
|
47 |
INTEGER knon ! nombre de points a traiter |
INTEGER knon ! nombre de points a traiter |
134 |
z1(i) = zgeop(i, 1) |
z1(i) = zgeop(i, 1) |
135 |
ENDDO |
ENDDO |
136 |
|
|
137 |
CALL clcdrag(nsrf, u1, v1, t1, q1, z1, ts, qsurf, rugos, coefm(:, 1), & |
CALL clcdrag(nsrf, u1, v1, t1, q1, z1, ts, qsurf, rugos, ycdragm, ycdragh) |
|
coefh(:, 1)) |
|
138 |
|
|
139 |
! Calculer les coefficients turbulents dans l'atmosphere |
! Calculer les coefficients turbulents dans l'atmosphere |
140 |
|
|