73 |
REAL sn(klon) |
REAL sn(klon) |
74 |
REAL qsolsrf(klon, nbsrf), snsrf(klon, nbsrf) |
REAL qsolsrf(klon, nbsrf), snsrf(klon, nbsrf) |
75 |
REAL albe(klon, nbsrf), evap(klon, nbsrf) |
REAL albe(klon, nbsrf), evap(klon, nbsrf) |
|
REAL alblw(klon, nbsrf) |
|
76 |
REAL tsoil(klon, nsoilmx, nbsrf) |
REAL tsoil(klon, nsoilmx, nbsrf) |
77 |
REAL radsol(klon), rain_fall(klon), snow_fall(klon) |
REAL radsol(klon), rain_fall(klon), snow_fall(klon) |
78 |
REAL solsw(klon), sollw(klon), fder(klon) |
REAL solsw(klon), sollw(klon), fder(klon) |
331 |
albe(:, is_lic) = 0.6 |
albe(:, is_lic) = 0.6 |
332 |
albe(:, is_oce) = 0.5 |
albe(:, is_oce) = 0.5 |
333 |
albe(:, is_sic) = 0.6 |
albe(:, is_sic) = 0.6 |
|
alblw = albe |
|
334 |
evap = 0. |
evap = 0. |
335 |
qsolsrf = 150. |
qsolsrf = 150. |
336 |
tsoil = spread(spread(pack(tsol_2d, dyn_phy), 2, nsoilmx), 3, nbsrf) |
tsoil = spread(spread(pack(tsol_2d, dyn_phy), 2, nsoilmx), 3, nbsrf) |
357 |
|
|
358 |
call phyredem("startphy.nc", pctsrf, tsoil(:, 1, :), tsoil, & |
call phyredem("startphy.nc", pctsrf, tsoil(:, 1, :), tsoil, & |
359 |
tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, & |
tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, & |
360 |
albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, & |
albe, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, frugs, & |
361 |
frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, & |
agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, & |
362 |
q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01) |
q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01) |
363 |
|
|
364 |
END SUBROUTINE etat0 |
END SUBROUTINE etat0 |