39 |
! En entr\'ee : la valeur au d\'ebut du pas de temps ; en sortie : la |
! En entr\'ee : la valeur au d\'ebut du pas de temps ; en sortie : la |
40 |
! valeur \`a la fin du pas de temps. |
! valeur \`a la fin du pas de temps. |
41 |
|
|
42 |
REAL km(:, :) ! (knon, klev + 1) |
REAL, intent(out):: km(:, 2:) ! (knon, 2:klev) |
43 |
! diffusivit\'e turbulente de quantit\'e de mouvement (au bas de |
! diffusivit\'e turbulente de quantit\'e de mouvement (au bas de |
44 |
! chaque couche) (en sortie : la valeur \`a la fin du pas de temps) |
! chaque couche) (en sortie : la valeur \`a la fin du pas de temps) |
45 |
|
|
46 |
REAL kn(:, :) ! (knon, klev + 1) |
REAL, intent(out):: kn(:, 2:) ! (knon, 2:klev) |
47 |
! diffusivit\'e turbulente des scalaires (au bas de chaque couche) |
! diffusivit\'e turbulente des scalaires (au bas de chaque couche) |
48 |
! (en sortie : la valeur \`a la fin du pas de temps) |
! (en sortie : la valeur \`a la fin du pas de temps) |
49 |
|
|
85 |
size(teta, 1), size(ustar), size(q2, 1), size(km, 1), size(kn, 1)], & |
size(teta, 1), size(ustar), size(q2, 1), size(km, 1), size(kn, 1)], & |
86 |
"yamada4 knon") |
"yamada4 knon") |
87 |
call assert(klev == [size(zlev, 2) - 1, size(zlay, 2), size(u, 2), & |
call assert(klev == [size(zlev, 2) - 1, size(zlay, 2), size(u, 2), & |
88 |
size(v, 2), size(teta, 2), size(q2, 2) - 1, size(km, 2) - 1, & |
size(v, 2), size(teta, 2), size(q2, 2) - 1, size(km, 2) + 1, & |
89 |
size(kn, 2) - 1], "yamada4 klev") |
size(kn, 2) + 1], "yamada4 klev") |
90 |
|
|
91 |
ipas = ipas + 1 |
ipas = ipas + 1 |
92 |
|
|