/[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

trunk/phylmd/Interface_surf/fonte_neige.f revision 305 by guez, Tue Sep 11 11:08:38 2018 UTC trunk/phylmd/Interface_surf/fonte_neige.f90 revision 332 by guez, Tue Aug 13 09:19:22 2019 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 29  contains Line 29  contains
29      ! precipitation, solid water mass flux (kg / m2 / s), positive down      ! precipitation, solid water mass flux (kg / m2 / s), positive down
30    
31      real, intent(INOUT):: snow(:) ! (knon)      real, intent(INOUT):: snow(:) ! (knon)
32      ! column-density of mass of snow, in kg m-2      ! column-density of mass of snow at the surface, in kg m-2
33    
34      real, intent(INOUT):: qsol(:) ! (knon)      real, intent(INOUT):: qsol(:) ! (knon)
35      ! column-density of water in soil, in kg m-2      ! column-density of water in soil, in kg m-2
# 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.332
  ViewVC Help
Powered by ViewVC 1.1.21