16 |
! Arguments: |
! Arguments: |
17 |
|
|
18 |
REAL, INTENT (IN) :: ucov(ip1jmp1, llm) |
REAL, INTENT (IN) :: ucov(ip1jmp1, llm) |
19 |
REAL :: teta(ip1jmp1, llm), masse(ip1jmp1, llm) |
REAL, INTENT(IN):: teta(ip1jmp1, llm) |
20 |
|
real masse(ip1jmp1, llm) |
21 |
REAL, INTENT (IN) :: ps(ip1jmp1) |
REAL, INTENT (IN) :: ps(ip1jmp1) |
22 |
REAL, INTENT (IN) :: phis(ip1jmp1) |
REAL, INTENT (IN) :: phis(ip1jmp1) |
23 |
REAL :: vorpot(ip1jm, llm) |
REAL :: vorpot(ip1jm, llm) |
47 |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
48 |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt / (iim * jjp1)) |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt / (iim * jjp1)) |
49 |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
50 |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE., 1) |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.) |
51 |
|
|
52 |
! Calcul du moment angulaire |
! Calcul du moment angulaire |
53 |
|
|