36 |
REAL qsol(klon) ! column-density of water in soil, in kg m-2 |
REAL qsol(klon) ! column-density of water in soil, in kg m-2 |
37 |
REAL snow(klon, nbsrf) |
REAL snow(klon, nbsrf) |
38 |
REAL albe(klon, nbsrf) |
REAL albe(klon, nbsrf) |
|
REAL evap(klon, nbsrf) |
|
39 |
REAL rain_fall(klon) |
REAL rain_fall(klon) |
40 |
REAL snow_fall(klon) |
REAL snow_fall(klon) |
41 |
real solsw(klon) |
real solsw(klon) |
76 |
pa = 5e4 |
pa = 5e4 |
77 |
call disvert |
call disvert |
78 |
p = ap + bp * preff |
p = ap + bp * preff |
79 |
call phyetat0(pctsrf, ftsol, ftsoil, qsurf, qsol, snow, albe, evap, & |
call phyetat0(pctsrf, ftsol, ftsoil, qsurf, qsol, snow, albe, rain_fall, & |
80 |
rain_fall, snow_fall, solsw, sollw, fder, radsol, frugs, agesno, zmea, & |
snow_fall, solsw, sollw, fder, radsol, frugs, agesno, zmea, zstd, zsig, & |
81 |
zstd, zsig, zgam, zthe, zpic, zval, t_ancien, q_ancien, ancien_ok, & |
zgam, zthe, zpic, zval, t_ancien, q_ancien, ancien_ok, rnebcon, ratqs, & |
82 |
rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01, ncid_startphy) |
clwcon, run_off_lic_0, sig1, w01, ncid_startphy) |
83 |
|
|
84 |
do julien = 1, 360 |
do julien = 1, 360 |
85 |
wo(:, :, julien) = ozonecm(REAL(julien), spread(p, dim=1, ncopies=jjm+1)) |
wo(:, :, julien) = ozonecm(REAL(julien), spread(p, dim=1, ncopies=jjm+1)) |