57 |
enddo |
enddo |
58 |
enddo |
enddo |
59 |
|
|
60 |
zdt=dtime/float(nsplit_thermals) |
zdt=dtime/real(nsplit_thermals) |
61 |
do isplit = 1, nsplit_thermals |
do isplit = 1, nsplit_thermals |
62 |
CALL thermcell(klon, klev, zdt, pplay, paprs, pphi, u_seri, v_seri, & |
CALL thermcell(klon, klev, zdt, pplay, paprs, pphi, u_seri, v_seri, & |
63 |
t_seri, q_seri, d_u_the, d_v_the, d_t_the, d_q_the, zfm_therm, & |
t_seri, q_seri, d_u_the, d_v_the, d_t_the, d_q_the, zfm_therm, & |
65 |
tho_thermals) |
tho_thermals) |
66 |
|
|
67 |
! transformation de la derivee en tendance |
! transformation de la derivee en tendance |
68 |
d_t_the=d_t_the*dtime/float(nsplit_thermals) |
d_t_the=d_t_the*dtime/real(nsplit_thermals) |
69 |
d_u_the=d_u_the*dtime/float(nsplit_thermals) |
d_u_the=d_u_the*dtime/real(nsplit_thermals) |
70 |
d_v_the=d_v_the*dtime/float(nsplit_thermals) |
d_v_the=d_v_the*dtime/real(nsplit_thermals) |
71 |
d_q_the=d_q_the*dtime/float(nsplit_thermals) |
d_q_the=d_q_the*dtime/real(nsplit_thermals) |
72 |
fm_therm=fm_therm +zfm_therm/float(nsplit_thermals) |
fm_therm=fm_therm +zfm_therm/real(nsplit_thermals) |
73 |
entr_therm=entr_therm +zentr_therm/float(nsplit_thermals) |
entr_therm=entr_therm +zentr_therm/real(nsplit_thermals) |
74 |
fm_therm(:, klev+1)=0. |
fm_therm(:, klev+1)=0. |
75 |
|
|
76 |
! accumulation de la tendance |
! accumulation de la tendance |