17 |
USE conf_phys_m, ONLY: iflag_pbl |
USE conf_phys_m, ONLY: iflag_pbl |
18 |
USE dimphy, ONLY: klev, klon |
USE dimphy, ONLY: klev, klon |
19 |
USE dimsoil, ONLY: nsoilmx |
USE dimsoil, ONLY: nsoilmx |
|
USE indicesol, ONLY: nbsrf |
|
20 |
USE interfsurf_hq_m, ONLY: interfsurf_hq |
USE interfsurf_hq_m, ONLY: interfsurf_hq |
21 |
USE suphec_m, ONLY: rcpd, rd, rg, rkappa |
USE suphec_m, ONLY: rcpd, rd, rg, rkappa |
22 |
|
|
45 |
|
|
46 |
REAL t(klon, klev) ! temperature (K) |
REAL t(klon, klev) ! temperature (K) |
47 |
REAL q(klon, klev) ! humidite specifique (kg / kg) |
REAL q(klon, klev) ! humidite specifique (kg / kg) |
48 |
REAL, intent(in):: ts(klon) ! temperature du sol (K) |
REAL, intent(in):: ts(:) ! (knon) temperature du sol (K) |
49 |
REAL paprs(klon, klev+1) ! pression a inter-couche (Pa) |
REAL paprs(klon, klev+1) ! pression a inter-couche (Pa) |
50 |
REAL pplay(klon, klev) ! pression au milieu de couche (Pa) |
REAL pplay(klon, klev) ! pression au milieu de couche (Pa) |
51 |
REAL delp(klon, klev) ! epaisseur de couche en pression (Pa) |
REAL delp(klon, klev) ! epaisseur de couche en pression (Pa) |
66 |
REAL, intent(inout):: agesno(:) ! (knon) |
REAL, intent(inout):: agesno(:) ! (knon) |
67 |
REAL d_t(klon, klev) ! incrementation de "t" |
REAL d_t(klon, klev) ! incrementation de "t" |
68 |
REAL d_q(klon, klev) ! incrementation de "q" |
REAL d_q(klon, klev) ! incrementation de "q" |
69 |
REAL, intent(out):: d_ts(:) ! (knon) incrementation de "ts" |
REAL, intent(out):: d_ts(:) ! (knon) incr\'ementation de "ts" |
70 |
real z0_new(klon) |
real z0_new(klon) |
71 |
|
|
72 |
REAL, intent(out):: flux_t(:) ! (knon) |
REAL, intent(out):: flux_t(:) ! (knon) |
246 |
CALL interfsurf_hq(dtime, jour, rmu0, nisurf, knon, knindex, rlat, debut, & |
CALL interfsurf_hq(dtime, jour, rmu0, nisurf, knon, knindex, rlat, debut, & |
247 |
nsoilmx, tsoil, qsol, u1lay, v1lay, temp_air, spechum, tq_cdrag, & |
nsoilmx, tsoil, qsol, u1lay, v1lay, temp_air, spechum, tq_cdrag, & |
248 |
petAcoef, peqAcoef, petBcoef, peqBcoef, precip_rain, precip_snow, & |
petAcoef, peqAcoef, petBcoef, peqBcoef, precip_rain, precip_snow, & |
249 |
fder, rugos, rugoro, snow, qsurf, ts(:knon), p1lay, psref, radsol, & |
fder, rugos, rugoro, snow, qsurf, ts, p1lay, psref, radsol, & |
250 |
evap, flux_t, fluxlat, dflux_l, dflux_s, tsurf_new, albedo, & |
evap, flux_t, fluxlat, dflux_l, dflux_s, tsurf_new, albedo, & |
251 |
z0_new, pctsrf_new_sic, agesno, fqcalving, ffonte, run_off_lic_0) |
z0_new, pctsrf_new_sic, agesno, fqcalving, ffonte, run_off_lic_0) |
252 |
|
|
253 |
flux_q = - evap |
flux_q = - evap |
254 |
d_ts = tsurf_new - ts(:knon) |
d_ts = tsurf_new - ts |
255 |
|
|
256 |
!==== une fois on a zx_h_ts, on peut faire l'iteration ======== |
!==== une fois on a zx_h_ts, on peut faire l'iteration ======== |
257 |
DO i = 1, knon |
DO i = 1, knon |