39 |
integer, intent(IN):: nisurf ! index de la surface a traiter |
integer, intent(IN):: nisurf ! index de la surface a traiter |
40 |
integer, intent(IN):: knon ! nombre de points de la surface a traiter |
integer, intent(IN):: knon ! nombre de points de la surface a traiter |
41 |
|
|
42 |
integer, intent(in):: knindex(klon) |
integer, intent(in):: knindex(:) ! (knon) |
43 |
! index des points de la surface a traiter |
! index des points de la surface a traiter |
44 |
|
|
45 |
real, intent(IN):: pctsrf(klon, nbsrf) |
real, intent(IN):: pctsrf(klon, nbsrf) |
90 |
! ps pression au sol |
! ps pression au sol |
91 |
REAL, DIMENSION(klon), INTENT(INOUT):: radsol |
REAL, DIMENSION(klon), INTENT(INOUT):: radsol |
92 |
! radsol rayonnement net aus sol (LW + SW) |
! radsol rayonnement net aus sol (LW + SW) |
93 |
real, dimension(klon), intent(INOUT):: evap |
real, intent(INOUT):: evap(klon) ! evaporation totale |
|
! evap evaporation totale |
|
94 |
real, dimension(klon), intent(OUT):: fluxsens, fluxlat |
real, dimension(klon), intent(OUT):: fluxsens, fluxlat |
95 |
! fluxsens flux de chaleur sensible |
! fluxsens flux de chaleur sensible |
96 |
! fluxlat flux de chaleur latente |
! fluxlat flux de chaleur latente |
206 |
|
|
207 |
! calcul albedo: lecture albedo fichier boundary conditions |
! calcul albedo: lecture albedo fichier boundary conditions |
208 |
! puis ajout albedo neige |
! puis ajout albedo neige |
209 |
call interfsur_lim(itime, dtime, jour, nisurf, knon, knindex, & |
call interfsur_lim(itime, dtime, jour, nisurf, knindex, debut, & |
210 |
debut, alb_new, z0_new) |
alb_new, z0_new) |
211 |
|
|
212 |
! calcul snow et qsurf, hydrol adapté |
! calcul snow et qsurf, hydrol adapté |
213 |
CALL calbeta(nisurf, snow(:knon), qsol(:knon), beta(:knon), & |
CALL calbeta(nisurf, snow(:knon), qsol(:knon), beta(:knon), & |
244 |
case (is_oce) |
case (is_oce) |
245 |
! Surface "ocean" appel à l'interface avec l'océan |
! Surface "ocean" appel à l'interface avec l'océan |
246 |
! lecture conditions limites |
! lecture conditions limites |
247 |
call interfoce_lim(itime, dtime, jour, klon, knon, knindex, debut, & |
call interfoce_lim(itime, dtime, jour, knindex, debut, tsurf_temp, & |
248 |
tsurf_temp, pctsrf_new) |
pctsrf_new) |
249 |
|
|
250 |
cal = 0. |
cal = 0. |
251 |
beta = 1. |
beta = 1. |
286 |
! Surface "glace de mer" appel a l'interface avec l'ocean |
! Surface "glace de mer" appel a l'interface avec l'ocean |
287 |
|
|
288 |
! ! lecture conditions limites |
! ! lecture conditions limites |
289 |
CALL interfoce_lim(itime, dtime, jour, klon, knon, knindex, & |
CALL interfoce_lim(itime, dtime, jour, knindex, debut, tsurf_new, & |
290 |
debut, tsurf_new, pctsrf_new) |
pctsrf_new) |
291 |
|
|
292 |
DO ii = 1, knon |
DO ii = 1, knon |
293 |
tsurf_new(ii) = tsurf(ii) |
tsurf_new(ii) = tsurf(ii) |