165 |
finvmaold, dtvr, leapf=.false.) |
finvmaold, dtvr, leapf=.false.) |
166 |
end if |
end if |
167 |
|
|
168 |
|
forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps |
169 |
|
CALL exner_hyb(ps, p3d, pks, pk, pkf) |
170 |
|
|
171 |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN |
172 |
! Calcul des tendances physiques: |
! Calcul des tendances physiques: |
173 |
|
|
|
forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps |
|
|
CALL exner_hyb(ps, p3d, pks, pk) |
|
|
|
|
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 |
183 |
CALL addfi(ucov, vcov, teta, q, dufi, dvfi, dtetafi, dqfi) |
CALL addfi(ucov, vcov, teta, q, dufi, dvfi, dtetafi, dqfi) |
184 |
ENDIF |
ENDIF |
185 |
|
|
|
forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps |
|
|
CALL exner_hyb(ps, p3d, pks, pk, pkf) |
|
|
|
|
186 |
IF (MOD(itau + 1, idissip) == 0) THEN |
IF (MOD(itau + 1, idissip) == 0) THEN |
187 |
! Dissipation horizontale et verticale des petites \'echelles |
! Dissipation horizontale et verticale des petites \'echelles |
188 |
|
|