169 |
! Calcul des tendances physiques: |
! Calcul des tendances physiques: |
170 |
|
|
171 |
forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps |
forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps |
172 |
CALL exner_hyb(ps, p3d, pks, pk, pkf) |
CALL exner_hyb(ps, p3d, pks, pk) |
173 |
|
|
174 |
rdaym_ini = itau * dtvr / daysec |
rdaym_ini = itau * dtvr / daysec |
175 |
rdayvrai = rdaym_ini + day_ini |
rdayvrai = rdaym_ini + day_ini |
176 |
time = REAL(mod(itau, day_step)) / day_step + time_0 |
time = REAL(mod(itau, day_step)) / day_step + time_0 |
177 |
IF (time > 1.) time = time - 1. |
IF (time > 1.) time = time - 1. |
178 |
|
|
179 |
CALL calfis(rdayvrai, time, ucov, vcov, teta, q, ps, pk, phis, phi, & |
CALL calfis(rdayvrai, time, ucov, vcov, teta, q, pk, phis, phi, w, & |
180 |
w, dufi, dvfi, dtetafi, dqfi, lafin = itau + 1 == itaufin) |
dufi, dvfi, dtetafi, dqfi, lafin = itau + 1 == itaufin) |
181 |
|
|
182 |
! Ajout des tendances physiques: |
! Ajout des tendances physiques: |
183 |
CALL addfi(ucov, vcov, teta, q, dufi, dvfi, dtetafi, dqfi) |
CALL addfi(ucov, vcov, teta, q, dufi, dvfi, dtetafi, dqfi) |