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 |
|
REAL, intent(in):: ycoefm(klon) ! drag coefficient |
|
23 |
|
|
24 |
REAL, intent(inout):: km(klon, 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 |
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 |
|
|