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 ps(ip1jmp1), masse(ip1jmp1, llm) |
REAL, intent(inout):: ps(ip1jmp1) |
27 |
|
REAL masse(ip1jmp1, llm) |
28 |
|
|
29 |
REAL vcovm1(ip1jm, llm), ucovm1(ip1jmp1, llm) |
REAL vcovm1(ip1jm, llm), ucovm1(ip1jmp1, llm) |
30 |
REAL tetam1(ip1jmp1, llm), psm1(ip1jmp1), massem1(ip1jmp1, llm) |
REAL tetam1(ip1jmp1, llm), psm1(ip1jmp1), massem1(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 |