84 |
! rayonnement net au sol (LW + SW) |
! rayonnement net au sol (LW + SW) |
85 |
|
|
86 |
real, intent(OUT):: evap(:) ! (knon) evaporation totale |
real, intent(OUT):: evap(:) ! (knon) evaporation totale |
87 |
|
|
88 |
real, intent(OUT):: flux_t(:) ! (knon) flux de chaleur sensible |
real, intent(OUT):: flux_t(:) ! (knon) flux de chaleur sensible |
89 |
real, dimension(klon), intent(OUT):: fluxlat ! flux de chaleur latente |
! (Cp T) à la surface, positif vers le bas, W / m2 |
90 |
|
|
91 |
|
real, intent(OUT):: fluxlat(:) ! (knon) flux de chaleur latente |
92 |
real, dimension(klon), intent(OUT):: dflux_l, dflux_s |
real, dimension(klon), intent(OUT):: dflux_l, dflux_s |
93 |
real, intent(OUT):: tsurf_new(:) ! (knon) temp\'erature au sol |
real, intent(OUT):: tsurf_new(:) ! (knon) temp\'erature au sol |
94 |
real, intent(OUT):: albedo(:) ! (knon) albedo |
real, intent(OUT):: albedo(:) ! (knon) albedo |
189 |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
190 |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
191 |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
192 |
fluxlat(:knon), flux_t, dflux_s(:knon), dflux_l(:knon)) |
fluxlat, flux_t, dflux_s(:knon), dflux_l(:knon)) |
193 |
|
|
194 |
CALL fonte_neige(is_ter, dtime, tsurf, p1lay(:knon), beta(:knon), & |
CALL fonte_neige(is_ter, dtime, tsurf, p1lay(:knon), beta(:knon), & |
195 |
tq_cdrag(:knon), ps(:knon), precip_rain(:knon), & |
tq_cdrag(:knon), ps(:knon), precip_rain(:knon), & |
216 |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
217 |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
218 |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
219 |
fluxlat(:knon), flux_t, dflux_s(:knon), dflux_l(:knon)) |
fluxlat, flux_t, dflux_s(:knon), dflux_l(:knon)) |
220 |
fder = fder + dflux_s + dflux_l |
fder = fder + dflux_s + dflux_l |
221 |
albedo = alboc_cd(rmu0(knindex)) * fmagic |
albedo = alboc_cd(rmu0(knindex)) * fmagic |
222 |
z0_new = sqrt(rugos**2 + rugoro**2) |
z0_new = sqrt(rugos**2 + rugoro**2) |
254 |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
255 |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
256 |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
257 |
fluxlat(:knon), flux_t, dflux_s(:knon), dflux_l(:knon)) |
fluxlat, flux_t, dflux_s(:knon), dflux_l(:knon)) |
258 |
|
|
259 |
CALL fonte_neige(is_sic, dtime, tsurf_temp, p1lay(:knon), beta(:knon), & |
CALL fonte_neige(is_sic, dtime, tsurf_temp, p1lay(:knon), beta(:knon), & |
260 |
tq_cdrag(:knon), ps(:knon), precip_rain(:knon), & |
tq_cdrag(:knon), ps(:knon), precip_rain(:knon), & |
296 |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), & |
297 |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), & |
298 |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap, & |
299 |
fluxlat(:knon), flux_t, dflux_s(:knon), dflux_l(:knon)) |
fluxlat, flux_t, dflux_s(:knon), dflux_l(:knon)) |
300 |
|
|
301 |
call fonte_neige(is_lic, dtime, tsurf, p1lay(:knon), beta(:knon), & |
call fonte_neige(is_lic, dtime, tsurf, p1lay(:knon), beta(:knon), & |
302 |
tq_cdrag(:knon), ps(:knon), precip_rain(:knon), & |
tq_cdrag(:knon), ps(:knon), precip_rain(:knon), & |