316 |
coefh(:knon, 2:), coefm(:knon, 1), coefh(:knon, 1)) |
coefh(:knon, 2:), coefm(:knon, 1), coefh(:knon, 1)) |
317 |
|
|
318 |
IF (iflag_pbl == 1) THEN |
IF (iflag_pbl == 1) THEN |
319 |
CALL coefkz2(nsrf, knon, ypaprs, ypplay, yt, ycoefm0(:knon, :), & |
CALL coefkz2(nsrf, knon, ypaprs, ypplay, yt, ycoefm0(:knon, 2:), & |
320 |
ycoefh0(:knon, :)) |
ycoefh0(:knon, 2:)) |
321 |
|
ycoefm0(:knon, 1) = 0. |
322 |
|
ycoefh0(:knon, 1) = 0. |
323 |
coefm(:knon, :) = max(coefm(:knon, :), ycoefm0(:knon, :)) |
coefm(:knon, :) = max(coefm(:knon, :), ycoefm0(:knon, :)) |
324 |
coefh(:knon, :) = max(coefh(:knon, :), ycoefh0(:knon, :)) |
coefh(:knon, :) = max(coefh(:knon, :), ycoefh0(:knon, :)) |
325 |
END IF |
END IF |
393 |
! calculer la diffusion de "q" et de "h" |
! calculer la diffusion de "q" et de "h" |
394 |
CALL clqh(dtime, julien, firstcal, nsrf, ni(:knon), & |
CALL clqh(dtime, julien, firstcal, nsrf, ni(:knon), & |
395 |
ytsoil(:knon, :), yqsol(:knon), mu0, yrugos, yrugoro, & |
ytsoil(:knon, :), yqsol(:knon), mu0, yrugos, yrugoro, & |
396 |
yu(:knon, 1), yv(:knon, 1), coefh(:knon, :), yt, yq, & |
yu(:knon, 1), yv(:knon, 1), coefh(:knon, 2:), coefh(:knon, 1), & |
397 |
yts(:knon), ypaprs, ypplay, ydelp, yrads(:knon), yalb(:knon), & |
yt, yq, yts(:knon), ypaprs, ypplay, ydelp, yrads(:knon), & |
398 |
snow(:knon), yqsurf, yrain_f, ysnow_f, yfluxlat(:knon), & |
yalb(:knon), snow(:knon), yqsurf, yrain_f, ysnow_f, & |
399 |
pctsrf_new_sic, yagesno(:knon), y_d_t, y_d_q, y_d_ts(:knon), & |
yfluxlat(:knon), pctsrf_new_sic, yagesno(:knon), y_d_t, y_d_q, & |
400 |
yz0_new, y_flux_t(:knon), y_flux_q(:knon), y_dflux_t(:knon), & |
y_d_ts(:knon), yz0_new, y_flux_t(:knon), y_flux_q(:knon), & |
401 |
y_dflux_q(:knon), y_fqcalving, y_ffonte, y_run_off_lic_0) |
y_dflux_t(:knon), y_dflux_q(:knon), y_fqcalving, y_ffonte, & |
402 |
|
y_run_off_lic_0) |
403 |
|
|
404 |
! calculer la longueur de rugosite sur ocean |
! calculer la longueur de rugosite sur ocean |
405 |
yrugm = 0. |
yrugm = 0. |