5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE clqh(julien, nisurf, knindex, tsoil, qsol, mu0, rugos, rugoro, & |
SUBROUTINE clqh(julien, nisurf, knindex, tsoil, qsol, mu0, rugos, rugoro, & |
8 |
u1lay, v1lay, coef, tq_cdrag, t, q, ts, paprs, pplay, delp, radsol, & |
u1lay, v1lay, coef, cdragh, t, q, ts, paprs, pplay, delp, radsol, & |
9 |
albedo, snow, qsurf, rain_fall, snow_fall, fluxlat, pctsrf_new_sic, & |
albedo, snow, qsurf, rain_fall, snow_fall, fluxlat, pctsrf_new_sic, & |
10 |
agesno, d_t, d_q, d_ts, z0_new, flux_t, flux_q, dflux_s, dflux_l, & |
agesno, d_t, d_q, d_ts, z0_new, flux_t, flux_q, dflux_s, dflux_l, & |
11 |
fqcalving, ffonte, run_off_lic_0, run_off_lic) |
fqcalving, ffonte, run_off_lic_0, run_off_lic) |
39 |
! Le coefficient d'echange (m**2 / s) multiplie par le cisaillement |
! Le coefficient d'echange (m**2 / s) multiplie par le cisaillement |
40 |
! du vent (dV / dz) |
! du vent (dV / dz) |
41 |
|
|
42 |
REAL, intent(in):: tq_cdrag(:) ! (knon) sans unite |
REAL, intent(in):: cdragh(:) ! (knon) sans unite |
43 |
|
|
44 |
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
45 |
REAL, intent(in):: q(:, :) ! (knon, klev) humidite specifique (kg / kg) |
REAL, intent(in):: q(:, :) ! (knon, klev) humidite specifique (kg / kg) |
55 |
! epaisseur de couche en pression (Pa) |
! epaisseur de couche en pression (Pa) |
56 |
|
|
57 |
REAL, intent(in):: radsol(:) ! (knon) |
REAL, intent(in):: radsol(:) ! (knon) |
58 |
! rayonnement net au sol (Solaire + IR) W / m2 |
! surface net downward radiative flux, in W / m2 |
59 |
|
|
60 |
REAL, intent(inout):: albedo(:) ! (knon) albedo de la surface |
REAL, intent(inout):: albedo(:) ! (knon) albedo de la surface |
61 |
REAL, intent(inout):: snow(:) ! (knon) ! hauteur de neige |
REAL, intent(inout):: snow(:) ! (knon) ! hauteur de neige |
114 |
|
|
115 |
call climb_hq_down(pkf, cq, dq, ch, dh, paprs, pplay, t, coef, delp, q) |
call climb_hq_down(pkf, cq, dq, ch, dh, paprs, pplay, t, coef, delp, q) |
116 |
CALL interfsurf_hq(julien, mu0, nisurf, knindex, tsoil, qsol, u1lay, & |
CALL interfsurf_hq(julien, mu0, nisurf, knindex, tsoil, qsol, u1lay, & |
117 |
v1lay, t(:, 1), q(:, 1), tq_cdrag, ch(:, 1), cq(:, 1), dh(:, 1), & |
v1lay, t(:, 1), q(:, 1), cdragh, ch(:, 1), cq(:, 1), dh(:, 1), & |
118 |
dq(:, 1), rain_fall, snow_fall, rugos, rugoro, snow, qsurf, ts, & |
dq(:, 1), rain_fall, snow_fall, rugos, rugoro, snow, qsurf, ts, & |
119 |
pplay(:, 1), paprs(:, 1), radsol, evap, flux_t, fluxlat, dflux_l, & |
pplay(:, 1), paprs(:, 1), radsol, evap, flux_t, fluxlat, dflux_l, & |
120 |
dflux_s, tsurf_new, albedo, z0_new, pctsrf_new_sic, agesno, & |
dflux_s, tsurf_new, albedo, z0_new, pctsrf_new_sic, agesno, & |