67 |
REAL, intent(in):: paprs(klon, llm+1) |
REAL, intent(in):: paprs(klon, llm+1) |
68 |
! (pression pour chaque inter-couche, en Pa) |
! (pression pour chaque inter-couche, en Pa) |
69 |
|
|
70 |
REAL pplay(klon, llm) |
REAL, intent(in):: pplay(klon, llm) |
71 |
! (input pression pour le mileu de chaque couche (en Pa)) |
! (input pression pour le mileu de chaque couche (en Pa)) |
72 |
|
|
73 |
REAL pphi(klon, llm) |
REAL pphi(klon, llm) |
1014 |
IF (if_ebil >= 1) THEN |
IF (if_ebil >= 1) THEN |
1015 |
ztit='after dynamic' |
ztit='after dynamic' |
1016 |
CALL diagetpq(airephy, ztit, ip_ebil, 1, 1, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 1, 1, dtime & |
1017 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1018 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1019 |
! Comme les tendances de la physique sont ajoute dans la dynamique, |
! Comme les tendances de la physique sont ajoute dans la dynamique, |
1020 |
! on devrait avoir que la variation d'entalpie par la dynamique |
! on devrait avoir que la variation d'entalpie par la dynamique |
1090 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1091 |
ztit='after reevap' |
ztit='after reevap' |
1092 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 1, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 1, dtime & |
1093 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1094 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1095 |
call diagphy(airephy, ztit, ip_ebil & |
call diagphy(airephy, ztit, ip_ebil & |
1096 |
, zero_v, zero_v, zero_v, zero_v, zero_v & |
, zero_v, zero_v, zero_v, zero_v, zero_v & |
1207 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1208 |
ztit='after clmain' |
ztit='after clmain' |
1209 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
1210 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1211 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1212 |
call diagphy(airephy, ztit, ip_ebil & |
call diagphy(airephy, ztit, ip_ebil & |
1213 |
, zero_v, zero_v, zero_v, zero_v, sens & |
, zero_v, zero_v, zero_v, zero_v, sens & |
1437 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1438 |
ztit='after convect' |
ztit='after convect' |
1439 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
1440 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1441 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1442 |
call diagphy(airephy, ztit, ip_ebil & |
call diagphy(airephy, ztit, ip_ebil & |
1443 |
, zero_v, zero_v, zero_v, zero_v, zero_v & |
, zero_v, zero_v, zero_v, zero_v, zero_v & |
1519 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1520 |
ztit='after dry_adjust' |
ztit='after dry_adjust' |
1521 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
1522 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1523 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1524 |
END IF |
END IF |
1525 |
|
|
1600 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1601 |
ztit='after fisrt' |
ztit='after fisrt' |
1602 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
1603 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1604 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1605 |
call diagphy(airephy, ztit, ip_ebil & |
call diagphy(airephy, ztit, ip_ebil & |
1606 |
, zero_v, zero_v, zero_v, zero_v, zero_v & |
, zero_v, zero_v, zero_v, zero_v, zero_v & |
1688 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1689 |
ztit="after diagcld" |
ztit="after diagcld" |
1690 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
1691 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1692 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1693 |
END IF |
END IF |
1694 |
|
|
1798 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1799 |
ztit='after rad' |
ztit='after rad' |
1800 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
1801 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1802 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1803 |
call diagphy(airephy, ztit, ip_ebil & |
call diagphy(airephy, ztit, ip_ebil & |
1804 |
, topsw, toplw, solsw, sollw, zero_v & |
, topsw, toplw, solsw, sollw, zero_v & |
1920 |
IF (if_ebil >= 2) THEN |
IF (if_ebil >= 2) THEN |
1921 |
ztit='after orography' |
ztit='after orography' |
1922 |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 2, 2, dtime & |
1923 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1924 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1925 |
END IF |
END IF |
1926 |
|
|
1979 |
IF (if_ebil >= 1) THEN |
IF (if_ebil >= 1) THEN |
1980 |
ztit='after physic' |
ztit='after physic' |
1981 |
CALL diagetpq(airephy, ztit, ip_ebil, 1, 1, dtime & |
CALL diagetpq(airephy, ztit, ip_ebil, 1, 1, dtime & |
1982 |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs, pplay & |
, t_seri, q_seri, ql_seri, qs_seri, u_seri, v_seri, paprs & |
1983 |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
, d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) |
1984 |
! Comme les tendances de la physique sont ajoute dans la dynamique, |
! Comme les tendances de la physique sont ajoute dans la dynamique, |
1985 |
! on devrait avoir que la variation d'entalpie par la dynamique |
! on devrait avoir que la variation d'entalpie par la dynamique |