91 |
|
|
92 |
! Flux d'eau "perdue" par la surface et n\'ecessaire pour limiter la |
! Flux d'eau "perdue" par la surface et n\'ecessaire pour limiter la |
93 |
! hauteur de neige, en kg / m2 / s |
! hauteur de neige, en kg / m2 / s |
94 |
real, dimension(klon), intent(INOUT):: fqcalving |
real, intent(OUT):: fqcalving(:) ! (knon) |
95 |
|
|
96 |
! Flux thermique utiliser pour fondre la neige |
! Flux thermique utiliser pour fondre la neige |
97 |
real, dimension(klon), intent(INOUT):: ffonte |
real, dimension(klon), intent(INOUT):: ffonte |
144 |
! Initialisations diverses |
! Initialisations diverses |
145 |
|
|
146 |
ffonte(1:knon) = 0. |
ffonte(1:knon) = 0. |
|
fqcalving(1:knon) = 0. |
|
147 |
dif_grnd = 999999. |
dif_grnd = 999999. |
148 |
z0_new = 999999. |
z0_new = 999999. |
149 |
|
|
172 |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
173 |
CALL fonte_neige(is_ter, dtime, precip_rain(:knon), & |
CALL fonte_neige(is_ter, dtime, precip_rain(:knon), & |
174 |
precip_snow(:knon), snow, qsol, tsurf_new, evap, & |
precip_snow(:knon), snow, qsol, tsurf_new, evap, & |
175 |
fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon)) |
fqcalving, ffonte(:knon), run_off_lic_0(:knon)) |
176 |
|
|
177 |
call albsno(dtime, agesno, alb_neig, precip_snow(:knon)) |
call albsno(dtime, agesno, alb_neig, precip_snow(:knon)) |
178 |
where (snow < 0.0001) agesno = 0. |
where (snow < 0.0001) agesno = 0. |
194 |
agesno = 0. |
agesno = 0. |
195 |
albedo = alboc_cd(rmu0(knindex)) * fmagic |
albedo = alboc_cd(rmu0(knindex)) * fmagic |
196 |
z0_new = sqrt(rugos**2 + rugoro**2) |
z0_new = sqrt(rugos**2 + rugoro**2) |
197 |
|
fqcalving = 0. |
198 |
case (is_sic) |
case (is_sic) |
199 |
! Surface "glace de mer" appel a l'interface avec l'ocean |
! Surface "glace de mer" appel a l'interface avec l'ocean |
200 |
|
|
221 |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
222 |
CALL fonte_neige(is_sic, dtime, precip_rain(:knon), & |
CALL fonte_neige(is_sic, dtime, precip_rain(:knon), & |
223 |
precip_snow(:knon), snow, qsol, tsurf_new, evap, & |
precip_snow(:knon), snow, qsol, tsurf_new, evap, & |
224 |
fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon)) |
fqcalving, ffonte(:knon), run_off_lic_0(:knon)) |
225 |
|
|
226 |
! Compute the albedo: |
! Compute the albedo: |
227 |
|
|
246 |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l) |
247 |
call fonte_neige(is_lic, dtime, precip_rain(:knon), & |
call fonte_neige(is_lic, dtime, precip_rain(:knon), & |
248 |
precip_snow(:knon), snow, qsol, tsurf_new, evap, & |
precip_snow(:knon), snow, qsol, tsurf_new, evap, & |
249 |
fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon)) |
fqcalving, ffonte(:knon), run_off_lic_0(:knon)) |
250 |
|
|
251 |
! calcul albedo |
! calcul albedo |
252 |
CALL albsno(dtime, agesno, alb_neig, precip_snow(:knon)) |
CALL albsno(dtime, agesno, alb_neig, precip_snow(:knon)) |