27 |
real, intent(in):: teta(ip1jmp1, llm) |
real, intent(in):: teta(ip1jmp1, llm) |
28 |
REAL, INTENT(IN):: ps(ip1jmp1), phis(ip1jmp1) |
REAL, INTENT(IN):: ps(ip1jmp1), phis(ip1jmp1) |
29 |
REAL :: vorpot(ip1jm, llm) |
REAL :: vorpot(ip1jm, llm) |
30 |
REAL :: phi(ip1jmp1, llm), bern(ip1jmp1, llm) |
REAL, intent(in):: phi(ip1jmp1, llm) |
31 |
|
real bern(ip1jmp1, llm) |
32 |
REAL :: dp(ip1jmp1) |
REAL :: dp(ip1jmp1) |
33 |
REAL, INTENT (IN):: time_0 |
REAL, INTENT (IN):: time_0 |
34 |
REAL, INTENT (IN):: pk(ip1jmp1, llm) |
REAL, INTENT (IN):: pk(ip1jmp1, llm) |
61 |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
62 |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1) |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1) |
63 |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
64 |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE., 1) |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.) |
65 |
|
|
66 |
! Calcul du moment angulaire |
! Calcul du moment angulaire |
67 |
radsg = rad/g |
radsg = rad/g |