/[lmdze]/trunk/Sources/phylmd/Interface_surf/albsno.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/Interface_surf/albsno.f

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

revision 173 by guez, Wed Apr 29 15:47:56 2015 UTC revision 174 by guez, Wed Nov 25 20:14:19 2015 UTC
# Line 6  module albsno_m Line 6  module albsno_m
6    
7  contains  contains
8    
9    SUBROUTINE albsno(klon, knon, dtime, agesno, alb_neig_grid, precip_snow)    SUBROUTINE albsno(dtime, agesno, alb_neig, precip_snow)
10    
11      INTEGER klon, knon      REAL, intent(in):: dtime
12      REAL dtime      REAL, intent(inout):: agesno(:) ! (knon)
13      REAL, DIMENSION(klon):: alb_neig_grid, agesno, precip_snow      real, intent(out):: alb_neig(:) ! (knon)
14        real, intent(in):: precip_snow(:) !(knon)
     ! Local:  
   
     INTEGER, PARAMETER :: nvm = 8  
     REAL veget(klon, nvm)  
   
     INTEGER i, nv  
   
     REAL init(nvm), decay(nvm)  
     REAL as  
     DATA init /0.55, 0.14, 0.18, 0.29, 0.15, 0.15, 0.14, 0./  
     DATA decay /0.30, 0.67, 0.63, 0.45, 0.40, 0.14, 0.06, 1./  
15    
16      !------------------------------------------------------------------------      !------------------------------------------------------------------------
17    
18      veget = 0.      ! D\'esert partout:
19      veget(:, 1) = 1. ! desert partout      alb_neig = 0.55 + 0.3 * EXP(- agesno / 5.)
20      DO i = 1, knon  
21         alb_neig_grid(i) = 0.0      ! Modulation en fonction de l'\^age de la neige :
22      ENDDO      agesno = max((agesno + (1. - agesno / 50.) * dtime / 86400.) &
23      DO nv = 1, nvm           * EXP(- MAX(0., precip_snow) * dtime / 0.3), 0.)
        DO i = 1, knon  
           as = init(nv)+decay(nv)*EXP(-agesno(i)/5.)  
           alb_neig_grid(i) = alb_neig_grid(i) + veget(i, nv)*as  
        ENDDO  
     ENDDO  
   
     !! modilation en fonction de l'age de la neige  
   
     DO i = 1, knon  
        agesno(i) = (agesno(i) + (1.-agesno(i)/50.)*dtime/86400.)&  
             & * EXP(-1.*MAX(0.0, precip_snow(i))*dtime/0.3)  
        agesno(i) = MAX(agesno(i), 0.0)  
     ENDDO  
24    
25    END SUBROUTINE albsno    END SUBROUTINE albsno
26    
Legend:
Removed from v.173  
changed lines
  Added in v.174
  ViewVC Help
Powered by ViewVC 1.1.21