5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE clqh(dtime, itime, jour, debut, rlat, knon, nisurf, knindex, & |
SUBROUTINE clqh(dtime, itime, jour, debut, rlat, knon, nisurf, knindex, & |
8 |
pctsrf, soil_model, tsoil, qsol, ok_veget, ocean, rmu0, co2_ppm, & |
pctsrf, soil_model, tsoil, qsol, rmu0, co2_ppm, rugos, rugoro, u1lay, & |
9 |
rugos, rugoro, u1lay, v1lay, coef, t, q, ts, paprs, pplay, delp, & |
v1lay, coef, t, q, ts, paprs, pplay, delp, radsol, albedo, alblw, & |
10 |
radsol, albedo, alblw, snow, qsurf, precip_rain, precip_snow, fder, & |
snow, qsurf, precip_rain, precip_snow, fder, swnet, fluxlat, & |
11 |
swnet, fluxlat, pctsrf_new, agesno, d_t, d_q, d_ts, z0_new, flux_t, & |
pctsrf_new, agesno, d_t, d_q, d_ts, z0_new, flux_t, flux_q, dflux_s, & |
12 |
flux_q, dflux_s, dflux_l, fqcalving, ffonte, run_off_lic_0, flux_o, & |
dflux_l, fqcalving, ffonte, run_off_lic_0, flux_o, flux_g, tslab, seaice) |
|
flux_g, tslab, seaice) |
|
13 |
|
|
14 |
! Author: Z. X. Li (LMD/CNRS) |
! Author: Z. X. Li (LMD/CNRS) |
15 |
! Date: 1993/08/18 |
! Date: 1993/08/18 |
53 |
integer, intent(in):: jour ! jour de l'annee en cours |
integer, intent(in):: jour ! jour de l'annee en cours |
54 |
real, intent(in):: rmu0(klon) ! cosinus de l'angle solaire zenithal |
real, intent(in):: rmu0(klon) ! cosinus de l'angle solaire zenithal |
55 |
real rugos(klon) ! rugosite |
real rugos(klon) ! rugosite |
56 |
integer knindex(klon) |
integer, intent(in):: knindex(klon) |
57 |
real, intent(in):: pctsrf(klon, nbsrf) |
real, intent(in):: pctsrf(klon, nbsrf) |
58 |
real, intent(in):: rlat(klon) |
real, intent(in):: rlat(klon) |
59 |
logical ok_veget |
REAL, intent(in):: co2_ppm ! taux CO2 atmosphere |
|
REAL co2_ppm ! taux CO2 atmosphere |
|
|
character(len=*), intent(in):: ocean |
|
60 |
|
|
61 |
REAL d_t(klon, klev) ! incrementation de "t" |
REAL d_t(klon, klev) ! incrementation de "t" |
62 |
REAL d_q(klon, klev) ! incrementation de "q" |
REAL d_q(klon, klev) ! incrementation de "q" |
270 |
endif |
endif |
271 |
ccanopy = co2_ppm |
ccanopy = co2_ppm |
272 |
|
|
273 |
CALL interfsurf_hq(itime, dtime, jour, rmu0, iim, jjm, & |
CALL interfsurf_hq(itime, dtime, jour, rmu0, iim, jjm, nisurf, knon, & |
274 |
nisurf, knon, knindex, pctsrf, rlat, debut, & |
knindex, pctsrf, rlat, debut, soil_model, nsoilmx, tsoil, qsol, & |
275 |
ok_veget, soil_model, nsoilmx, tsoil, qsol, u1lay, v1lay, & |
u1lay, v1lay, temp_air, spechum, tq_cdrag, petAcoef, peqAcoef, & |
276 |
temp_air, spechum, tq_cdrag, petAcoef, peqAcoef, & |
petBcoef, peqBcoef, precip_rain, precip_snow, fder, rugos, rugoro, & |
277 |
petBcoef, peqBcoef, precip_rain, precip_snow, & |
snow, qsurf, ts, p1lay, psref, radsol, evap, fluxsens, & |
278 |
fder, rugos, rugoro, & |
fluxlat, dflux_l, dflux_s, tsurf_new, alb_new, alblw, z0_new, & |
279 |
snow, qsurf, ts, p1lay, psref, radsol, ocean, & |
pctsrf_new, agesno, fqcalving, ffonte, run_off_lic_0, flux_o, & |
280 |
evap, fluxsens, fluxlat, dflux_l, dflux_s, tsurf_new, & |
flux_g, tslab, seaice) |
|
alb_new, alblw, z0_new, pctsrf_new, agesno, fqcalving, & |
|
|
ffonte, run_off_lic_0, flux_o, flux_g, tslab, seaice) |
|
281 |
|
|
282 |
do i = 1, knon |
do i = 1, knon |
283 |
flux_t(i, 1) = fluxsens(i) |
flux_t(i, 1) = fluxsens(i) |