53 |
real, dimension(klon), intent(IN):: temp_air, spechum |
real, dimension(klon), intent(IN):: temp_air, spechum |
54 |
! temp_air temperature de l'air 1ere couche |
! temp_air temperature de l'air 1ere couche |
55 |
! spechum humidite specifique 1ere couche |
! spechum humidite specifique 1ere couche |
56 |
real, dimension(klon), intent(INOUT):: tq_cdrag ! coefficient d'echange |
real, intent(IN):: tq_cdrag(:) ! (knon) coefficient d'echange |
57 |
|
|
58 |
real, dimension(klon), intent(IN):: petAcoef, peqAcoef |
real, dimension(klon), intent(IN):: petAcoef, peqAcoef |
59 |
! coefficients A de la r\'esolution de la couche limite pour t et q |
! coefficients A de la r\'esolution de la couche limite pour t et q |
175 |
cal = RCPD * capsol(:knon) |
cal = RCPD * capsol(:knon) |
176 |
ENDIF |
ENDIF |
177 |
|
|
178 |
CALL calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag(:knon), & |
CALL calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag, & |
179 |
ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), & |
ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), & |
180 |
temp_air(:knon), spechum(:knon), u1_lay, v1_lay, & |
temp_air(:knon), spechum(:knon), u1_lay, v1_lay, & |
181 |
petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |
petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |
197 |
beta = 1. |
beta = 1. |
198 |
dif_grnd = 0. |
dif_grnd = 0. |
199 |
call calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, & |
call calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, & |
200 |
tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol, & |
tq_cdrag, ps(:knon), qsurf(:knon), radsol, & |
201 |
dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay, & |
dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay, & |
202 |
v1_lay, petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |
v1_lay, petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |
203 |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
234 |
beta = 1. |
beta = 1. |
235 |
|
|
236 |
CALL calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, & |
CALL calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, & |
237 |
tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol, & |
tq_cdrag, ps(:knon), qsurf(:knon), radsol, & |
238 |
dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay, & |
dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay, & |
239 |
v1_lay, petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |
v1_lay, petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |
240 |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
264 |
beta = 1. |
beta = 1. |
265 |
dif_grnd = 0. |
dif_grnd = 0. |
266 |
|
|
267 |
call calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag(:knon), & |
call calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag, & |
268 |
ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), & |
ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), & |
269 |
temp_air(:knon), spechum(:knon), u1_lay, v1_lay, & |
temp_air(:knon), spechum(:knon), u1_lay, v1_lay, & |
270 |
petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |
petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), & |