60 |
use readsulfate_m, only: readsulfate |
use readsulfate_m, only: readsulfate |
61 |
use readsulfate_preind_m, only: readsulfate_preind |
use readsulfate_preind_m, only: readsulfate_preind |
62 |
use yoegwd, only: sugwd |
use yoegwd, only: sugwd |
63 |
USE suphec_m, ONLY: ra, rcpd, retv, rg, rlvtt, romega, rsigma, rtt |
USE suphec_m, ONLY: rcpd, retv, rg, rlvtt, romega, rsigma, rtt |
64 |
USE temps, ONLY: itau_phy |
USE temps, ONLY: itau_phy |
65 |
|
use transp_m, only: transp |
66 |
use unit_nml_m, only: unit_nml |
use unit_nml_m, only: unit_nml |
67 |
USE ymds2ju_m, ONLY: ymds2ju |
USE ymds2ju_m, ONLY: ymds2ju |
68 |
USE yoethf_m, ONLY: r2es, rvtmp2 |
USE yoethf_m, ONLY: r2es, rvtmp2 |
290 |
! Declaration des procedures appelees |
! Declaration des procedures appelees |
291 |
|
|
292 |
EXTERNAL nuage ! calculer les proprietes radiatives |
EXTERNAL nuage ! calculer les proprietes radiatives |
|
EXTERNAL transp ! transport total de l'eau et de l'energie |
|
293 |
|
|
294 |
! Variables locales |
! Variables locales |
295 |
|
|
523 |
|
|
524 |
IF (if_ebil >= 1) zero_v = 0. |
IF (if_ebil >= 1) zero_v = 0. |
525 |
IF (nqmx < 2) CALL abort_gcm('physiq', & |
IF (nqmx < 2) CALL abort_gcm('physiq', & |
526 |
'eaux vapeur et liquide sont indispensables', 1) |
'eaux vapeur et liquide sont indispensables') |
527 |
|
|
528 |
test_firstcal: IF (firstcal) THEN |
test_firstcal: IF (firstcal) THEN |
529 |
! initialiser |
! initialiser |
1303 |
ENDDO |
ENDDO |
1304 |
ENDDO |
ENDDO |
1305 |
|
|
1306 |
CALL aaam_bud(ra, rg, romega, rlat, rlon, pphis, zustrdr, zustrli, & |
CALL aaam_bud(rg, romega, rlat, rlon, pphis, zustrdr, zustrli, zustrph, & |
1307 |
zustrph, zvstrdr, zvstrli, zvstrph, paprs, u, v, aam, torsfc) |
zvstrdr, zvstrli, zvstrph, paprs, u, v, aam, torsfc) |
1308 |
|
|
1309 |
IF (if_ebil >= 2) CALL diagetpq(airephy, 'after orography', ip_ebil, 2, & |
IF (if_ebil >= 2) CALL diagetpq(airephy, 'after orography', ip_ebil, 2, & |
1310 |
2, dtphys, t_seri, q_seri, ql_seri, u_seri, v_seri, paprs, d_h_vcol, & |
2, dtphys, t_seri, q_seri, ql_seri, u_seri, v_seri, paprs, d_h_vcol, & |
1321 |
pctsrf, frac_impa, frac_nucl, pphis, airephy, dtphys, itap) |
pctsrf, frac_impa, frac_nucl, pphis, airephy, dtphys, itap) |
1322 |
|
|
1323 |
! Calculer le transport de l'eau et de l'energie (diagnostique) |
! Calculer le transport de l'eau et de l'energie (diagnostique) |
1324 |
CALL transp(paprs, zxtsol, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, & |
CALL transp(paprs, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, ue, uq) |
|
ue, uq) |
|
1325 |
|
|
1326 |
! diag. bilKP |
! diag. bilKP |
1327 |
|
|