19 |
REAL, intent(in):: yu(klon, klev), yv(klon, klev) ! wind, in m s-1 |
REAL, intent(in):: yu(klon, klev), yv(klon, klev) ! wind, in m s-1 |
20 |
REAL, intent(in):: yt(klon, klev) ! temperature, in K |
REAL, intent(in):: yt(klon, klev) ! temperature, in K |
21 |
REAL, intent(in):: yq(klon, klev) |
REAL, intent(in):: yq(klon, klev) |
22 |
|
REAL, intent(in):: ycoefm(:) ! (ngrid) drag coefficient |
23 |
|
|
24 |
REAL, intent(in):: ycoefm(klon) ! drag coefficient |
REAL, intent(inout):: km(:, 2:) ! (knon, 2:klev) |
|
|
|
|
REAL, intent(inout):: km(klon, klev) |
|
25 |
! coefficient de diffusion turbulente de quantité de mouvement (au |
! coefficient de diffusion turbulente de quantité de mouvement (au |
26 |
! bas de chaque couche) (en sortie : la valeur à la fin du pas de |
! bas de chaque couche) (en sortie : la valeur à la fin du pas de |
27 |
! temps), m2 s-1 |
! temps), m2 s-1 |
28 |
|
|
29 |
REAL, intent(inout):: kn(klon, klev) |
REAL, intent(inout):: kn(:, 2:) ! (knon, 2:klev) |
30 |
! coefficient de diffusion turbulente des scalaires (au bas de |
! coefficient de diffusion turbulente des scalaires (au bas de |
31 |
! chaque couche) (en sortie : la valeur à la fin du pas de temps), m2 s-1 |
! chaque couche) (en sortie : la valeur à la fin du pas de temps), m2 s-1 |
32 |
|
|
73 |
enddo |
enddo |
74 |
|
|
75 |
forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1)) |
forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1)) |
76 |
ustar = SQRT(ycoefm(:ngrid) * (yu(:ngrid, 1)**2 + yv(:ngrid, 1)**2)) |
ustar = SQRT(ycoefm * (yu(:ngrid, 1)**2 + yv(:ngrid, 1)**2)) |
77 |
|
|
78 |
! Fin de la partie qui doit être incluse à terme dans clmain |
! Fin de la partie qui doit être incluse à terme dans clmain |
79 |
|
|