4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE coefkz2(nsrf, paprs, pplay, t, pcfm, pcfh) |
SUBROUTINE coefkz2(nsrf, paprs, pplay, t, coefm, coefh) |
8 |
|
|
9 |
! J'introduit un peu de diffusion sauf dans les endroits o\`u une |
! J'introduit un peu de diffusion sauf dans les endroits o\`u une |
10 |
! forte inversion est pr\'esente. On peut dire que la diffusion |
! forte inversion est pr\'esente. On peut dire que la diffusion |
24 |
|
|
25 |
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
26 |
|
|
27 |
REAL, intent(out):: pcfm(:, 2:) ! (knon, 2:klev) coefficient vitesse |
REAL, intent(out):: coefm(:, 2:) ! (knon, 2:klev) coefficient vitesse |
28 |
|
|
29 |
REAL, intent(out):: pcfh(:, 2:) ! (knon, 2:klev) |
REAL, intent(out):: coefh(:, 2:) ! (knon, 2:klev) |
30 |
! coefficient chaleur et humidite) |
! coefficient chaleur et humidite) |
31 |
|
|
32 |
! Local: |
! Local: |
56 |
! Initialiser les sorties |
! Initialiser les sorties |
57 |
DO k = 2, klev |
DO k = 2, klev |
58 |
DO i = 1, knon |
DO i = 1, knon |
59 |
pcfm(i, k) = 0.0 |
coefm(i, k) = 0.0 |
60 |
pcfh(i, k) = 0.0 |
coefh(i, k) = 0.0 |
61 |
ENDDO |
ENDDO |
62 |
ENDDO |
ENDDO |
63 |
|
|
88 |
((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN |
((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN |
89 |
zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) & |
zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) & |
90 |
/(paprs(i, 2)-paprs(i, klev+1))))**2 |
/(paprs(i, 2)-paprs(i, klev+1))))**2 |
91 |
pcfm(i, k)= zl2(i)* kstable |
coefm(i, k)= zl2(i)* kstable |
92 |
pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different |
coefh(i, k) = coefm(i, k) /prandtl ! h et m different |
93 |
ENDIF |
ENDIF |
94 |
ENDDO |
ENDDO |
95 |
ENDDO |
ENDDO |