42 |
REAL massebx(ip1jmp1, llm), masseby((iim + 1) * jjm, llm) |
REAL massebx(ip1jmp1, llm), masseby((iim + 1) * jjm, llm) |
43 |
REAL vorpot(iim + 1, jjm, llm) |
REAL vorpot(iim + 1, jjm, llm) |
44 |
real ecin(iim + 1, jjm + 1, llm), convm(ip1jmp1, llm) |
real ecin(iim + 1, jjm + 1, llm), convm(ip1jmp1, llm) |
|
REAL bern(iim + 1, jjm + 1, llm) |
|
45 |
REAL massebxy(iim + 1, jjm, llm), dp(ip1jmp1) |
REAL massebxy(iim + 1, jjm, llm), dp(ip1jmp1) |
46 |
INTEGER l |
INTEGER l |
47 |
|
|
60 |
CALL vitvert(convm, w) |
CALL vitvert(convm, w) |
61 |
CALL tourpot(vcov, ucov, massebxy, vorpot) |
CALL tourpot(vcov, ucov, massebxy, vorpot) |
62 |
CALL enercin(vcov, ucov, vcont, ucont, ecin) |
CALL enercin(vcov, ucov, vcont, ucont, ecin) |
63 |
CALL bernoui(phi, ecin, bern) |
CALL sortvarc(ucov, teta, ps, masse, pk, phis, vorpot, phi, & |
64 |
CALL sortvarc(ucov, teta, ps, masse, pk, phis, vorpot, phi, bern, dp, & |
bernoui(phi, ecin), dp, resetvarc = .true.) |
|
resetvarc = .true.) |
|
65 |
|
|
66 |
END SUBROUTINE caldyn0 |
END SUBROUTINE caldyn0 |
67 |
|
|