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

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

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

revision 299 by guez, Thu Aug 2 14:27:11 2018 UTC revision 301 by guez, Thu Aug 2 17:23:07 2018 UTC
# Line 5  module fonte_neige_m Line 5  module fonte_neige_m
5  contains  contains
6    
7    SUBROUTINE fonte_neige(nisurf, precip_rain, precip_snow, snow, qsol, &    SUBROUTINE fonte_neige(nisurf, precip_rain, precip_snow, snow, qsol, &
8         tsurf_new, evap, fqcalving, ffonte, run_off_lic_0)         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
# Line 45  contains Line 45  contains
45      real, intent(INOUT):: run_off_lic_0(:) ! (knon)      real, intent(INOUT):: run_off_lic_0(:) ! (knon)
46      ! run off glacier du pas de temps pr\'ecedent      ! run off glacier du pas de temps pr\'ecedent
47    
48        REAL, intent(OUT):: run_off_lic(:) ! (knon) ruissellement total
49    
50      ! Local:      ! Local:
51    
52      integer knon ! nombre de points \`a traiter      integer knon ! nombre de points \`a traiter
# Line 60  contains Line 62  contains
62      REAL, parameter:: chaice = 3.334E5 / (2.3867E6 * 0.15)      REAL, parameter:: chaice = 3.334E5 / (2.3867E6 * 0.15)
63      real, parameter:: max_eau_sol = 150. ! in kg m-2      real, parameter:: max_eau_sol = 150. ! in kg m-2
64      real coeff_rel      real coeff_rel
     REAL, ALLOCATABLE, SAVE:: run_off_lic(:) ! ruissellement total  
65    
66      !--------------------------------------------------------------------      !--------------------------------------------------------------------
67    
# Line 111  contains Line 112  contains
112      IF (nisurf == is_ter) then      IF (nisurf == is_ter) then
113         qsol = MIN(qsol + bil_eau_s, max_eau_sol)         qsol = MIN(qsol + bil_eau_s, max_eau_sol)
114      else if (nisurf == is_lic) then      else if (nisurf == is_lic) then
        if (.not. allocated(run_off_lic)) allocate(run_off_lic(knon))  
        ! assumes that the fraction of land-ice does not change during the run  
   
115         do i = 1, knon         do i = 1, knon
116            run_off_lic(i) = (coeff_rel * fqcalving(i)) + &            run_off_lic_0(i) = (coeff_rel * fqcalving(i)) + &
117                 (1. - coeff_rel) * run_off_lic_0(i)                 (1. - coeff_rel) * run_off_lic_0(i)
118            run_off_lic_0(i) = run_off_lic(i)            run_off_lic(i) = run_off_lic_0(i) + bil_eau_s(i) / dtphys
           run_off_lic(i) = run_off_lic(i) + bil_eau_s(i) / dtphys  
119         enddo         enddo
120      endif      endif
121    
Legend:
Removed from v.299  
changed lines
  Added in v.301
  ViewVC Help
Powered by ViewVC 1.1.21