5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE physiq(lafin, rdayvrai, time, dtphys, paprs, play, pphi, pphis, & |
SUBROUTINE physiq(lafin, rdayvrai, time, dtphys, paprs, play, pphi, pphis, & |
8 |
u, v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps, dudyn, PVteta) |
u, v, t, qx, omega, d_u, d_v, d_t, d_qx, d_ps, dudyn) |
9 |
|
|
10 |
! From phylmd/physiq.F, version 1.22 2006/02/20 09:38:28 |
! From phylmd/physiq.F, version 1.22 2006/02/20 09:38:28 |
11 |
! (subversion revision 678) |
! (subversion revision 678) |
107 |
INTEGER nbteta |
INTEGER nbteta |
108 |
PARAMETER(nbteta = 3) |
PARAMETER(nbteta = 3) |
109 |
|
|
|
REAL PVteta(klon, nbteta) |
|
|
! (output vorticite potentielle a des thetas constantes) |
|
|
|
|
110 |
LOGICAL ok_gust ! pour activer l'effet des gust sur flux surface |
LOGICAL ok_gust ! pour activer l'effet des gust sur flux surface |
111 |
PARAMETER (ok_gust = .FALSE.) |
PARAMETER (ok_gust = .FALSE.) |
112 |
|
|
1252 |
ratqs = ratqss |
ratqs = ratqss |
1253 |
endif |
endif |
1254 |
|
|
|
! Processus de condensation ŕ grande echelle et processus de |
|
|
! précipitation : |
|
1255 |
CALL fisrtilp(dtphys, paprs, play, t_seri, q_seri, ptconv, ratqs, & |
CALL fisrtilp(dtphys, paprs, play, t_seri, q_seri, ptconv, ratqs, & |
1256 |
d_t_lsc, d_q_lsc, d_ql_lsc, rneb, cldliq, rain_lsc, snow_lsc, & |
d_t_lsc, d_q_lsc, d_ql_lsc, rneb, cldliq, rain_lsc, snow_lsc, & |
1257 |
pfrac_impa, pfrac_nucl, pfrac_1nucl, frac_impa, frac_nucl, prfl, & |
pfrac_impa, pfrac_nucl, pfrac_1nucl, frac_impa, frac_nucl, prfl, & |
1553 |
|
|
1554 |
! Calcul des tendances traceurs |
! Calcul des tendances traceurs |
1555 |
call phytrac(rnpb, itap, lmt_pas, julien, time, firstcal, lafin, nqmx-2, & |
call phytrac(rnpb, itap, lmt_pas, julien, time, firstcal, lafin, nqmx-2, & |
1556 |
dtphys, u, t, paprs, play, mfu, mfd, pen_u, pde_u, pen_d, pde_d, & |
dtphys, u, t, paprs, play, mfu, mfd, pde_u, pen_d, ycoefh, fm_therm, & |
1557 |
ycoefh, fm_therm, entr_therm, yu1, yv1, ftsol, pctsrf, frac_impa, & |
entr_therm, yu1, yv1, ftsol, pctsrf, frac_impa, frac_nucl, pphis, & |
1558 |
frac_nucl, pphis, albsol, rhcl, cldfra, rneb, diafra, cldliq, & |
albsol, rhcl, cldfra, rneb, diafra, cldliq, pmflxr, pmflxs, prfl, & |
1559 |
pmflxr, pmflxs, prfl, psfl, da, phi, mp, upwd, dnwd, tr_seri, zmasse) |
psfl, da, phi, mp, upwd, dnwd, tr_seri, zmasse) |
1560 |
|
|
1561 |
IF (offline) THEN |
IF (offline) call phystokenc(dtphys, rlon, rlat, t, mfu, mfd, pen_u, & |
1562 |
call phystokenc(dtphys, rlon, rlat, t, mfu, mfd, pen_u, pde_u, & |
pde_u, pen_d, pde_d, fm_therm, entr_therm, ycoefh, yu1, yv1, ftsol, & |
1563 |
pen_d, pde_d, fm_therm, entr_therm, ycoefh, yu1, yv1, ftsol, & |
pctsrf, frac_impa, frac_nucl, pphis, airephy, dtphys, itap) |
|
pctsrf, frac_impa, frac_nucl, pphis, airephy, dtphys, itap) |
|
|
ENDIF |
|
1564 |
|
|
1565 |
! Calculer le transport de l'eau et de l'energie (diagnostique) |
! Calculer le transport de l'eau et de l'energie (diagnostique) |
1566 |
CALL transp(paprs, zxtsol, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, & |
CALL transp(paprs, zxtsol, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, & |