/[lmdze]/trunk/phylmd/Interface_surf/fonte_neige.f90
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/fonte_neige.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 305 by guez, Tue Sep 11 11:08:38 2018 UTC revision 311 by guez, Mon Dec 3 17:52:21 2018 UTC
# Line 8  contains Line 8  contains
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    
# Line 88  contains Line 88  contains
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)
Legend:
Removed from v.305  
changed lines
  Added in v.311
  ViewVC Help
Powered by ViewVC 1.1.21