9 |
|
|
10 |
contains |
contains |
11 |
|
|
12 |
SUBROUTINE physiq(firstcal, lafin, rdayvrai, gmtime, pdtphys, paprs, & |
SUBROUTINE physiq(lafin, rdayvrai, gmtime, pdtphys, paprs, & |
13 |
pplay, pphi, pphis, u, v, t, qx, omega, d_u, d_v, & |
pplay, pphi, pphis, u, v, t, qx, omega, d_u, d_v, & |
14 |
d_t, d_qx, d_ps, dudyn, PVteta) |
d_t, d_qx, d_ps, dudyn, PVteta) |
15 |
|
|
16 |
! From phylmd/physiq.F, v 1.22 2006/02/20 09:38:28 |
! From phylmd/physiq.F, version 1.22 2006/02/20 09:38:28 |
17 |
|
|
18 |
! Author : Z.X. Li (LMD/CNRS), date: 1993/08/18 |
! Author : Z.X. Li (LMD/CNRS), date: 1993/08/18 |
19 |
|
|
71 |
|
|
72 |
REAL, intent(in):: gmtime ! heure de la journée en fraction de jour |
REAL, intent(in):: gmtime ! heure de la journée en fraction de jour |
73 |
REAL, intent(in):: pdtphys ! pas d'integration pour la physique (seconde) |
REAL, intent(in):: pdtphys ! pas d'integration pour la physique (seconde) |
|
LOGICAL, intent(in):: firstcal ! first call to "calfis" |
|
74 |
logical, intent(in):: lafin ! dernier passage |
logical, intent(in):: lafin ! dernier passage |
75 |
|
|
76 |
REAL, intent(in):: paprs(klon, llm+1) |
REAL, intent(in):: paprs(klon, llm+1) |
85 |
REAL pphis(klon) ! input geopotentiel du sol |
REAL pphis(klon) ! input geopotentiel du sol |
86 |
|
|
87 |
REAL u(klon, llm) ! input vitesse dans la direction X (de O a E) en m/s |
REAL u(klon, llm) ! input vitesse dans la direction X (de O a E) en m/s |
88 |
REAL v(klon, llm) ! input vitesse Y (de S a N) en m/s |
REAL, intent(in):: v(klon, llm) ! vitesse Y (de S a N) en m/s |
89 |
REAL t(klon, llm) ! input temperature (K) |
REAL t(klon, llm) ! input temperature (K) |
90 |
|
|
91 |
REAL, intent(in):: qx(klon, llm, nqmx) |
REAL, intent(in):: qx(klon, llm, nqmx) |
98 |
REAL d_qx(klon, llm, nqmx) ! output tendance physique de "qx" (kg/kg/s) |
REAL d_qx(klon, llm, nqmx) ! output tendance physique de "qx" (kg/kg/s) |
99 |
REAL d_ps(klon) ! output tendance physique de la pression au sol |
REAL d_ps(klon) ! output tendance physique de la pression au sol |
100 |
|
|
101 |
|
LOGICAL:: firstcal = .true. |
102 |
|
|
103 |
INTEGER nbteta |
INTEGER nbteta |
104 |
PARAMETER(nbteta=3) |
PARAMETER(nbteta=3) |
105 |
|
|
1802 |
END IF |
END IF |
1803 |
|
|
1804 |
! Calcul des tendances traceurs |
! Calcul des tendances traceurs |
1805 |
call phytrac(rnpb, itap, lmt_pas, julien, gmtime, firstcal, lafin, nqmx-2, & |
call phytrac(rnpb, itap, lmt_pas, julien, gmtime, firstcal, lafin, & |
1806 |
pdtphys, u, v, t, paprs, pplay, pmfu, pmfd, pen_u, pde_u, pen_d, & |
nqmx-2, pdtphys, u, t, paprs, pplay, pmfu, pmfd, pen_u, pde_u, & |
1807 |
pde_d, ycoefh, fm_therm, entr_therm, yu1, yv1, ftsol, pctsrf, & |
pen_d, pde_d, ycoefh, fm_therm, entr_therm, yu1, yv1, ftsol, pctsrf, & |
1808 |
frac_impa, frac_nucl, pphis, pphi, albsol, rhcl, cldfra, rneb, & |
frac_impa, frac_nucl, pphis, pphi, albsol, rhcl, cldfra, rneb, & |
1809 |
diafra, cldliq, pmflxr, pmflxs, prfl, psfl, da, phi, mp, upwd, dnwd, & |
diafra, cldliq, pmflxr, pmflxs, prfl, psfl, da, phi, mp, upwd, dnwd, & |
1810 |
tr_seri, zmasse) |
tr_seri, zmasse) |
1811 |
|
|
1914 |
t_ancien, q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0) |
t_ancien, q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0) |
1915 |
ENDIF |
ENDIF |
1916 |
|
|
1917 |
|
firstcal = .FALSE. |
1918 |
|
|
1919 |
contains |
contains |
1920 |
|
|
1921 |
subroutine write_histday |
subroutine write_histday |