20 |
|
|
21 |
! Arguments: |
! Arguments: |
22 |
|
|
23 |
REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm), teta(ip1jmp1, llm) |
REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) |
24 |
|
real, intent(inout):: teta(ip1jmp1, llm) |
25 |
REAL q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nq) |
REAL q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nq) |
26 |
REAL, intent(inout):: ps(ip1jmp1) |
REAL, intent(inout):: ps(ip1jmp1) |
27 |
REAL masse(ip1jmp1, llm) |
REAL masse(ip1jmp1, llm) |
114 |
|
|
115 |
DO ij = 1, ip1jmp1 |
DO ij = 1, ip1jmp1 |
116 |
hscr(ij) = teta(ij, l) |
hscr(ij) = teta(ij, l) |
117 |
teta(ij, l) = tetam1(ij, l)*massem1(ij, l)/masse(ij, l) + & |
teta(ij, l) = tetam1(ij, l) * massem1(ij, l) / masse(ij, l) & |
118 |
dt*dteta(ij, l)/masse(ij, l) |
+ dt * dteta(ij, l) / masse(ij, l) |
119 |
END DO |
END DO |
120 |
|
|
121 |
! Calcul de la valeur moyenne, unique aux poles pour teta |
! Calcul de la valeur moyenne, unique aux poles pour teta |