381 |
REAL prfl(klon, llm + 1), psfl(klon, llm + 1) |
REAL prfl(klon, llm + 1), psfl(klon, llm + 1) |
382 |
|
|
383 |
INTEGER, save:: ibas_con(klon), itop_con(klon) |
INTEGER, save:: ibas_con(klon), itop_con(klon) |
384 |
|
real ema_pct(klon) ! Emanuel pressure at cloud top, in Pa |
385 |
|
|
386 |
REAL rain_con(klon), rain_lsc(klon) |
REAL rain_con(klon), rain_lsc(klon) |
387 |
REAL, save:: snow_con(klon) |
REAL, save:: snow_con(klon) ! neige (mm / s) |
388 |
real snow_lsc(klon) |
real snow_lsc(klon) |
389 |
REAL d_ts(klon, nbsrf) |
REAL d_ts(klon, nbsrf) |
390 |
|
|
862 |
mp = 0. |
mp = 0. |
863 |
phi = 0. |
phi = 0. |
864 |
CALL concvl(dtphys, paprs, play, t_seri, q_seri, u_seri, v_seri, sig1, & |
CALL concvl(dtphys, paprs, play, t_seri, q_seri, u_seri, v_seri, sig1, & |
865 |
w01, d_t_con, d_q_con, d_u_con, d_v_con, rain_con, snow_con, & |
w01, d_t_con, d_q_con, d_u_con, d_v_con, rain_con, ibas_con, & |
866 |
ibas_con, itop_con, upwd, dnwd, dnwd0, Ma, cape, iflagctrl, & |
itop_con, upwd, dnwd, dnwd0, Ma, cape, iflagctrl, qcondc, wd, & |
867 |
qcondc, wd, pmflxr, da, phi, mp) |
pmflxr, da, phi, mp) |
868 |
|
snow_con = 0. |
869 |
clwcon0 = qcondc |
clwcon0 = qcondc |
870 |
mfu = upwd + dnwd |
mfu = upwd + dnwd |
871 |
IF (.NOT. ok_gust) wd = 0. |
IF (.NOT. ok_gust) wd = 0. |
882 |
call clouds_gno(klon, llm, q_seri, zqsat, clwcon0, ptconv, ratqsc, & |
call clouds_gno(klon, llm, q_seri, zqsat, clwcon0, ptconv, ratqsc, & |
883 |
rnebcon0) |
rnebcon0) |
884 |
|
|
885 |
|
forall (i = 1:klon) ema_pct(i) = paprs(i,itop_con(i) + 1) |
886 |
mfd = 0. |
mfd = 0. |
887 |
pen_u = 0. |
pen_u = 0. |
888 |
pen_d = 0. |
pen_d = 0. |
1589 |
CALL gr_fi_ecrit(1, klon, iim, jjm + 1, s_trmb3, zx_tmp_2d) |
CALL gr_fi_ecrit(1, klon, iim, jjm + 1, s_trmb3, zx_tmp_2d) |
1590 |
CALL histwrite(nid_ins, "s_trmb3", itau_w, zx_tmp_2d) |
CALL histwrite(nid_ins, "s_trmb3", itau_w, zx_tmp_2d) |
1591 |
|
|
1592 |
|
if (conv_emanuel) then |
1593 |
|
CALL gr_fi_ecrit(1, klon, iim, jjm + 1, ema_pct, zx_tmp_2d) |
1594 |
|
CALL histwrite(nid_ins, "ptop", itau_w, zx_tmp_2d) |
1595 |
|
end if |
1596 |
|
|
1597 |
! Champs 3D: |
! Champs 3D: |
1598 |
|
|
1599 |
CALL gr_fi_ecrit(llm, klon, iim, jjm + 1, t_seri, zx_tmp_3d) |
CALL gr_fi_ecrit(llm, klon, iim, jjm + 1, t_seri, zx_tmp_3d) |