6 |
|
|
7 |
SUBROUTINE phystokenc(pdtphys, pt, pmfu, pmfd, pen_u, pde_u, pen_d, pde_d, & |
SUBROUTINE phystokenc(pdtphys, pt, pmfu, pmfd, pen_u, pde_u, pen_d, pde_d, & |
8 |
pfm_therm, pentr_therm, pcoefh, yu1, yv1, ftsol, pctsrf, frac_impa, & |
pfm_therm, pentr_therm, pcoefh, yu1, yv1, ftsol, pctsrf, frac_impa, & |
9 |
frac_nucl, pphis, paire, dtime, itap) |
frac_nucl, pphis, paire, dtime) |
10 |
|
|
11 |
! From phylmd/phystokenc.F, version 1.2 2004/06/22 11:45:35 |
! From phylmd/phystokenc.F, version 1.2 2004/06/22 11:45:35 |
12 |
! Author: Fr\'ed\'eric Hourdin |
! Author: Fr\'ed\'eric Hourdin |
19 |
USE indicesol, ONLY: nbsrf |
USE indicesol, ONLY: nbsrf |
20 |
use initphysto_m, only: initphysto |
use initphysto_m, only: initphysto |
21 |
USE dimphy, ONLY: klev, klon |
USE dimphy, ONLY: klev, klon |
22 |
|
use time_phylmdz, only: itap |
23 |
USE tracstoke, ONLY: istphy |
USE tracstoke, ONLY: istphy |
24 |
|
|
25 |
REAL, INTENT (IN):: pdtphys ! pas d'integration pour la physique (seconde) |
REAL, INTENT (IN):: pdtphys ! pas d'integration pour la physique (seconde) |
64 |
REAL, INTENT(IN):: pphis(klon) |
REAL, INTENT(IN):: pphis(klon) |
65 |
real paire(klon) |
real paire(klon) |
66 |
REAL, INTENT (IN):: dtime |
REAL, INTENT (IN):: dtime |
|
INTEGER, INTENT (IN):: itap |
|
67 |
|
|
68 |
! Variables local to the procedure: |
! Variables local to the procedure: |
69 |
|
|
109 |
|
|
110 |
IF (iadvtr==0) CALL initphysto('phystoke', dtime, dtime*istphy, dtime*istphy, physid) |
IF (iadvtr==0) CALL initphysto('phystoke', dtime, dtime*istphy, dtime*istphy, physid) |
111 |
|
|
112 |
i = itap |
CALL histwrite(physid, 'phis', itap, gr_phy_write(pphis)) |
113 |
CALL histwrite(physid, 'phis', i, gr_phy_write(pphis)) |
CALL histwrite(physid, 'aire', itap, gr_phy_write(paire)) |
|
i = itap |
|
|
CALL histwrite(physid, 'aire', i, gr_phy_write(paire)) |
|
114 |
iadvtr = iadvtr + 1 |
iadvtr = iadvtr + 1 |
115 |
|
|
116 |
IF (mod(iadvtr, istphy) == 1 .OR. istphy == 1) THEN |
IF (mod(iadvtr, istphy) == 1 .OR. istphy == 1) THEN |