106 |
call assert(shape(ucov) == (/iim + 1, jjm + 1, llm/), "leapfrog") |
call assert(shape(ucov) == (/iim + 1, jjm + 1, llm/), "leapfrog") |
107 |
|
|
108 |
itaufin = nday * day_step |
itaufin = nday * day_step |
109 |
! "day_step" is a multiple of "iperiod", therefore "itaufin" is one too |
! "day_step" is a multiple of "iperiod", therefore so is "itaufin". |
110 |
|
|
111 |
dq = 0. |
dq = 0. |
112 |
|
|
129 |
massem1 = masse |
massem1 = masse |
130 |
psm1 = ps |
psm1 = ps |
131 |
finvmaold = masse |
finvmaold = masse |
132 |
CALL filtreg(finvmaold, jjm + 1, llm, - 2, 2, .TRUE., 1) |
CALL filtreg(finvmaold, jjm + 1, llm, - 2, 2, .TRUE.) |
133 |
end if |
end if |
134 |
|
|
135 |
! Calcul des tendances dynamiques: |
! Calcul des tendances dynamiques: |