5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, & |
SUBROUTINE caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, & |
8 |
conser, du, dv, dteta, dp, w, pbaru, pbarv, time_0) |
dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, conser) |
9 |
|
|
10 |
! From dyn3d/caldyn.F, v 1.1.1.1 2004/05/19 12:53:06 |
! From dyn3d/caldyn.F, v 1.1.1.1 2004/05/19 12:53:06 |
11 |
! Auteur : P. Le Van |
! Auteur : P. Le Van |
29 |
REAL pkf(ip1jmp1, llm) |
REAL pkf(ip1jmp1, llm) |
30 |
REAL vcont(ip1jm, llm), ucont(ip1jmp1, llm) |
REAL vcont(ip1jm, llm), ucont(ip1jmp1, llm) |
31 |
REAL phi(ip1jmp1, llm), masse(ip1jmp1, llm) |
REAL phi(ip1jmp1, llm), masse(ip1jmp1, llm) |
32 |
REAL dv(ip1jm, llm), du(ip1jmp1, llm) |
REAL dv(ip1jm, llm), dudyn(ip1jmp1, llm) |
33 |
REAL dteta(ip1jmp1, llm) |
REAL dteta(ip1jmp1, llm) |
34 |
real, INTENT(out):: dp(ip1jmp1) |
real, INTENT(out):: dp(ip1jmp1) |
35 |
REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) |
REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) |
60 |
dp = convm(:, 1) / airesurg |
dp = convm(:, 1) / airesurg |
61 |
CALL vitvert(convm, w) |
CALL vitvert(convm, w) |
62 |
CALL tourpot(vcov, ucov, massebxy, vorpot) |
CALL tourpot(vcov, ucov, massebxy, vorpot) |
63 |
CALL dudv1(vorpot, pbaru, pbarv, du, dv) |
CALL dudv1(vorpot, pbaru, pbarv, dudyn, dv) |
64 |
CALL enercin(vcov, ucov, vcont, ucont, ecin) |
CALL enercin(vcov, ucov, vcont, ucont, ecin) |
65 |
CALL bernoui(ip1jmp1, llm, phi, ecin, bern) |
CALL bernoui(ip1jmp1, llm, phi, ecin, bern) |
66 |
CALL dudv2(teta, pkf, bern, du, dv) |
CALL dudv2(teta, pkf, bern, dudyn, dv) |
67 |
|
|
68 |
DO l = 1, llm |
DO l = 1, llm |
69 |
DO ij = 1, ip1jmp1 |
DO ij = 1, ip1jmp1 |
71 |
END DO |
END DO |
72 |
END DO |
END DO |
73 |
|
|
74 |
CALL advect(ang, vcov, teta, w, massebx, masseby, du, dv, dteta, conser) |
CALL advect(ang, vcov, teta, w, massebx, masseby, dudyn, dv, dteta) |
75 |
|
|
76 |
! WARNING probleme de peridocite de dv sur les PC/linux. Pb d'arrondi |
! WARNING probleme de peridocite de dv sur les PC/linux. Pb d'arrondi |
77 |
! probablement. Observe sur le code compile avec pgf90 3.0-1 |
! probablement. Observe sur le code compile avec pgf90 3.0-1 |