8 |
tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic) |
tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic) |
9 |
|
|
10 |
! Routine de traitement de la fonte de la neige dans le cas du traitement |
! Routine de traitement de la fonte de la neige dans le cas du traitement |
11 |
! de sol simplifi\'e |
! de sol simplifi\'e. |
12 |
|
|
13 |
! Laurent Fairhead, March, 2001 |
! Laurent Fairhead, March, 2001 |
14 |
|
|
88 |
! Y a-t-il fonte de neige ? |
! Y a-t-il fonte de neige ? |
89 |
|
|
90 |
do i = 1, knon |
do i = 1, knon |
91 |
if ((snow(i) > epsfra .OR. nisurf == is_sic & |
if ((snow(i) > epsfra .OR. nisurf == is_sic .OR. nisurf == is_lic) & |
92 |
.OR. nisurf == is_lic) .AND. tsurf_new(i) >= RTT) then |
.AND. tsurf_new(i) >= RTT) then |
93 |
fq_fonte = MIN(MAX((tsurf_new(i) - RTT) / chasno, 0.), snow(i)) |
fq_fonte = MIN(MAX((tsurf_new(i) - RTT) / chasno, 0.), snow(i)) |
94 |
ffonte(i) = fq_fonte * RLMLT / dtphys |
ffonte(i) = fq_fonte * RLMLT / dtphys |
95 |
snow(i) = max(0., snow(i) - fq_fonte) |
snow(i) = max(0., snow(i) - fq_fonte) |