/[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 300 by guez, Thu Aug 2 15:55:01 2018 UTC
# Line 60  contains Line 60  contains
60      REAL, parameter:: chaice = 3.334E5 / (2.3867E6 * 0.15)      REAL, parameter:: chaice = 3.334E5 / (2.3867E6 * 0.15)
61      real, parameter:: max_eau_sol = 150. ! in kg m-2      real, parameter:: max_eau_sol = 150. ! in kg m-2
62      real coeff_rel      real coeff_rel
63      REAL, ALLOCATABLE, SAVE:: run_off_lic(:) ! ruissellement total      REAL run_off_lic(size(precip_rain)) ! (knon) ruissellement total
64    
65      !--------------------------------------------------------------------      !--------------------------------------------------------------------
66    
# Line 111  contains Line 111  contains
111      IF (nisurf == is_ter) then      IF (nisurf == is_ter) then
112         qsol = MIN(qsol + bil_eau_s, max_eau_sol)         qsol = MIN(qsol + bil_eau_s, max_eau_sol)
113      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  
   
114         do i = 1, knon         do i = 1, knon
115            run_off_lic(i) = (coeff_rel * fqcalving(i)) + &            run_off_lic_0(i) = (coeff_rel * fqcalving(i)) + &
116                 (1. - coeff_rel) * run_off_lic_0(i)                 (1. - coeff_rel) * run_off_lic_0(i)
117            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  
118         enddo         enddo
119      endif      endif
120    
Legend:
Removed from v.299  
changed lines
  Added in v.300
  ViewVC Help
Powered by ViewVC 1.1.21