49 |
REAL pplay(klon, klev) ! pression au milieu de couche (Pa) |
REAL pplay(klon, klev) ! pression au milieu de couche (Pa) |
50 |
REAL delp(klon, klev) ! epaisseur de couche en pression (Pa) |
REAL delp(klon, klev) ! epaisseur de couche en pression (Pa) |
51 |
|
|
52 |
REAL, intent(inout):: radsol(:) ! (knon) |
REAL, intent(in):: radsol(:) ! (knon) |
53 |
! rayonnement net au sol (Solaire + IR) W / m2 |
! rayonnement net au sol (Solaire + IR) W / m2 |
54 |
|
|
55 |
REAL, intent(inout):: albedo(:) ! (knon) albedo de la surface |
REAL, intent(inout):: albedo(:) ! (knon) albedo de la surface |
243 |
p1lay(1:knon) = pplay(1:knon, 1) |
p1lay(1:knon) = pplay(1:knon, 1) |
244 |
|
|
245 |
CALL interfsurf_hq(dtime, julien, rmu0, nisurf, knindex, debut, tsoil, & |
CALL interfsurf_hq(dtime, julien, rmu0, nisurf, knindex, debut, tsoil, & |
246 |
qsol, u1lay, v1lay, temp_air, spechum, tq_cdrag(:knon), petAcoef, peqAcoef, & |
qsol, u1lay, v1lay, temp_air, spechum, tq_cdrag(:knon), petAcoef, & |
247 |
petBcoef, peqBcoef, precip_rain, precip_snow, rugos, rugoro, snow, & |
peqAcoef, petBcoef, peqBcoef, precip_rain, precip_snow, rugos, & |
248 |
qsurf, ts, p1lay, psref, radsol, evap, flux_t, fluxlat, dflux_l, & |
rugoro, snow, qsurf, ts, p1lay, psref, radsol, evap, flux_t, fluxlat, & |
249 |
dflux_s, tsurf_new, albedo, z0_new, pctsrf_new_sic, agesno, & |
dflux_l, dflux_s, tsurf_new, albedo, z0_new, pctsrf_new_sic, agesno, & |
250 |
fqcalving, ffonte, run_off_lic_0) |
fqcalving, ffonte, run_off_lic_0) |
251 |
|
|
252 |
flux_q = - evap |
flux_q = - evap |