148 |
REAL lwup0(klon, llm + 1), lwup(klon, llm + 1) |
REAL lwup0(klon, llm + 1), lwup(klon, llm + 1) |
149 |
SAVE lwdn0, lwdn, lwup0, lwup |
SAVE lwdn0, lwdn, lwup0, lwup |
150 |
|
|
|
! Amip2 |
|
|
! variables a une pression donnee |
|
|
|
|
|
integer nlevSTD |
|
|
PARAMETER(nlevSTD = 17) |
|
|
|
|
151 |
! prw: precipitable water |
! prw: precipitable water |
152 |
real prw(klon) |
real prw(klon) |
153 |
|
|
156 |
REAL flwp(klon), fiwp(klon) |
REAL flwp(klon), fiwp(klon) |
157 |
REAL flwc(klon, llm), fiwc(klon, llm) |
REAL flwc(klon, llm), fiwc(klon, llm) |
158 |
|
|
|
INTEGER kmax, lmax |
|
|
PARAMETER(kmax = 8, lmax = 8) |
|
|
INTEGER kmaxm1, lmaxm1 |
|
|
PARAMETER(kmaxm1 = kmax - 1, lmaxm1 = lmax - 1) |
|
|
|
|
159 |
! Variables propres a la physique |
! Variables propres a la physique |
160 |
|
|
161 |
INTEGER, save:: radpas |
INTEGER, save:: radpas |
320 |
REAL, PARAMETER:: t_coup = 234. |
REAL, PARAMETER:: t_coup = 234. |
321 |
REAL zphi(klon, llm) |
REAL zphi(klon, llm) |
322 |
|
|
323 |
! cf. AM Variables pour la CLA (hbtm2) |
! cf. Anne Mathieu variables pour la couche limite atmosphérique (hbtm) |
324 |
|
|
325 |
REAL, SAVE:: pblh(klon, nbsrf) ! Hauteur de couche limite |
REAL, SAVE:: pblh(klon, nbsrf) ! Hauteur de couche limite |
326 |
REAL, SAVE:: plcl(klon, nbsrf) ! Niveau de condensation de la CLA |
REAL, SAVE:: plcl(klon, nbsrf) ! Niveau de condensation de la CLA |
332 |
REAL, SAVE:: trmb1(klon, nbsrf) ! deep_cape |
REAL, SAVE:: trmb1(klon, nbsrf) ! deep_cape |
333 |
REAL, SAVE:: trmb2(klon, nbsrf) ! inhibition |
REAL, SAVE:: trmb2(klon, nbsrf) ! inhibition |
334 |
REAL, SAVE:: trmb3(klon, nbsrf) ! Point Omega |
REAL, SAVE:: trmb3(klon, nbsrf) ! Point Omega |
335 |
! Grdeurs de sorties |
! Grandeurs de sorties |
336 |
REAL s_pblh(klon), s_lcl(klon), s_capCL(klon) |
REAL s_pblh(klon), s_lcl(klon), s_capCL(klon) |
337 |
REAL s_oliqCL(klon), s_cteiCL(klon), s_pblt(klon) |
REAL s_oliqCL(klon), s_cteiCL(klon), s_pblt(klon) |
338 |
REAL s_therm(klon), s_trmb1(klon), s_trmb2(klon) |
REAL s_therm(klon), s_trmb1(klon), s_trmb2(klon) |
370 |
REAL prfl(klon, llm + 1), psfl(klon, llm + 1) |
REAL prfl(klon, llm + 1), psfl(klon, llm + 1) |
371 |
|
|
372 |
INTEGER, save:: ibas_con(klon), itop_con(klon) |
INTEGER, save:: ibas_con(klon), itop_con(klon) |
373 |
|
real ema_pct(klon) ! Emanuel pressure at cloud top, in Pa |
374 |
|
|
375 |
REAL rain_con(klon), rain_lsc(klon) |
REAL rain_con(klon), rain_lsc(klon) |
376 |
REAL, save:: snow_con(klon) |
REAL, save:: snow_con(klon) ! neige (mm / s) |
377 |
real snow_lsc(klon) |
real snow_lsc(klon) |
378 |
REAL d_ts(klon, nbsrf) |
REAL d_ts(klon, nbsrf) |
379 |
|
|
487 |
real zmasse(klon, llm) |
real zmasse(klon, llm) |
488 |
! (column-density of mass of air in a cell, in kg m-2) |
! (column-density of mass of air in a cell, in kg m-2) |
489 |
|
|
|
real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2 |
|
490 |
integer, save:: ncid_startphy, itau_phy |
integer, save:: ncid_startphy, itau_phy |
491 |
|
|
492 |
namelist /physiq_nml/ ok_journe, ok_mensuel, ok_instan, fact_cldcon, & |
namelist /physiq_nml/ ok_journe, ok_mensuel, ok_instan, fact_cldcon, & |
850 |
mp = 0. |
mp = 0. |
851 |
phi = 0. |
phi = 0. |
852 |
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, & |
853 |
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, & |
854 |
ibas_con, itop_con, upwd, dnwd, dnwd0, Ma, cape, iflagctrl, & |
itop_con, upwd, dnwd, dnwd0, Ma, cape, iflagctrl, qcondc, wd, & |
855 |
qcondc, wd, pmflxr, da, phi, mp) |
pmflxr, da, phi, mp) |
856 |
|
snow_con = 0. |
857 |
clwcon0 = qcondc |
clwcon0 = qcondc |
858 |
mfu = upwd + dnwd |
mfu = upwd + dnwd |
859 |
IF (.NOT. ok_gust) wd = 0. |
IF (.NOT. ok_gust) wd = 0. |
870 |
call clouds_gno(klon, llm, q_seri, zqsat, clwcon0, ptconv, ratqsc, & |
call clouds_gno(klon, llm, q_seri, zqsat, clwcon0, ptconv, ratqsc, & |
871 |
rnebcon0) |
rnebcon0) |
872 |
|
|
873 |
|
forall (i = 1:klon) ema_pct(i) = paprs(i,itop_con(i) + 1) |
874 |
mfd = 0. |
mfd = 0. |
875 |
pen_u = 0. |
pen_u = 0. |
876 |
pen_d = 0. |
pen_d = 0. |
1386 |
! Ecriture des sorties |
! Ecriture des sorties |
1387 |
|
|
1388 |
use dimens_m, only: iim, jjm |
use dimens_m, only: iim, jjm |
1389 |
|
use gr_fi_ecrit_m, only: gr_fi_ecrit |
1390 |
USE histsync_m, ONLY: histsync |
USE histsync_m, ONLY: histsync |
1391 |
USE histwrite_m, ONLY: histwrite |
USE histwrite_m, ONLY: histwrite |
1392 |
|
|
1578 |
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) |
1579 |
CALL histwrite(nid_ins, "s_trmb3", itau_w, zx_tmp_2d) |
CALL histwrite(nid_ins, "s_trmb3", itau_w, zx_tmp_2d) |
1580 |
|
|
1581 |
|
if (conv_emanuel) then |
1582 |
|
CALL gr_fi_ecrit(1, klon, iim, jjm + 1, ema_pct, zx_tmp_2d) |
1583 |
|
CALL histwrite(nid_ins, "ptop", itau_w, zx_tmp_2d) |
1584 |
|
end if |
1585 |
|
|
1586 |
! Champs 3D: |
! Champs 3D: |
1587 |
|
|
1588 |
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) |