229 |
|
|
230 |
forall(l= 1: llm) v(:, l) = pack(zvfi(:, :, l), dyn_phy) |
forall(l= 1: llm) v(:, l) = pack(zvfi(:, :, l), dyn_phy) |
231 |
|
|
232 |
!IM calcul PV a teta=350, 380, 405K |
! Compute potential vorticity at theta = 350, 380 and 405 K: |
233 |
CALL PVtheta(klon, llm, ucov, vcov, teta, t, play, paprs, ntetaSTD, & |
CALL PVtheta(klon, llm, ucov, vcov, teta, t, play, paprs, ntetaSTD, & |
234 |
rtetaSTD, PVteta) |
rtetaSTD, PVteta) |
235 |
|
|
236 |
! Appel de la physique : |
! Appel de la physique : |
237 |
CALL physiq(lafin, rdayvrai, time, dtphys, paprs, play, pphi, pphis, u, & |
CALL physiq(lafin, rdayvrai, time, dtphys, paprs, play, pphi, pphis, u, & |
238 |
v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps, dudyn, PVteta) |
v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps, dudyn) |
239 |
|
|
240 |
! transformation des tendances physiques en tendances dynamiques: |
! transformation des tendances physiques en tendances dynamiques: |
241 |
|
|