19 |
USE comgeom, ONLY: aire_2d, apoln, apols |
USE comgeom, ONLY: aire_2d, apoln, apols |
20 |
USE disvert_m, ONLY: ap, bp |
USE disvert_m, ONLY: ap, bp |
21 |
USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, & |
USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, & |
22 |
iflag_phys, ok_guide, iecri |
iflag_phys, iecri |
23 |
|
USE conf_guide_m, ONLY: ok_guide |
24 |
USE dimens_m, ONLY: iim, jjm, llm, nqmx |
USE dimens_m, ONLY: iim, jjm, llm, nqmx |
25 |
use dissip_m, only: dissip |
use dissip_m, only: dissip |
26 |
USE dynetat0_m, ONLY: day_ini |
USE dynetat0_m, ONLY: day_ini |
121 |
else |
else |
122 |
! Matsuno |
! Matsuno |
123 |
dt = dtvr |
dt = dtvr |
124 |
if (ok_guide .and. itaufin - itau - 1 > day_step / 4) & |
if (ok_guide) call guide(itau, ucov, vcov, teta, q(:, :, :, 1), ps) |
|
call guide(itau, ucov, vcov, teta, q, ps) |
|
125 |
vcovm1 = vcov |
vcovm1 = vcov |
126 |
ucovm1 = ucov |
ucovm1 = ucov |
127 |
tetam1 = teta |
tetam1 = teta |
128 |
massem1 = masse |
massem1 = masse |
129 |
psm1 = ps |
psm1 = ps |
130 |
finvmaold = masse |
finvmaold = masse |
131 |
CALL filtreg(finvmaold, jjm + 1, llm, - 2, 2, .TRUE.) |
CALL filtreg(finvmaold, direct = .false., intensive = .false.) |
132 |
end if |
end if |
133 |
|
|
134 |
! Calcul des tendances dynamiques: |
! Calcul des tendances dynamiques: |