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

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

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

trunk/libf/phylmd/albsno.f90 revision 68 by guez, Wed Nov 14 16:59:30 2012 UTC trunk/phylmd/Interface_surf/albsno.f revision 305 by guez, Tue Sep 11 11:08:38 2018 UTC
# Line 1  Line 1 
1  module albsno_m  module albsno_m
2    
3    ! From phylmd/interface_surf.F90,v 1.8 2005/05/25 13:10:09    ! From phylmd/interface_surf.F90, version 1.8 2005/05/25 13:10:09
4    
5    IMPLICIT none    IMPLICIT none
6    
7  contains  contains
8    
9    SUBROUTINE albsno(klon, knon,dtime,agesno,alb_neig_grid, precip_snow)    SUBROUTINE albsno(agesno, alb_neig, snow_fall)
10    
11      INTEGER :: klon, knon      use comconst, only: dtphys
12      INTEGER, PARAMETER :: nvm = 8  
13      REAL   :: dtime      REAL, intent(inout):: agesno(:) ! (knon)
14      REAL, dimension(klon,nvm) :: veget      real, intent(out):: alb_neig(:) ! (knon)
15      REAL, DIMENSION(klon) :: alb_neig_grid, agesno, precip_snow      real, intent(in):: snow_fall(:) !(knon)
16    
17      INTEGER :: i, nv      !------------------------------------------------------------------------
18    
19      REAL, DIMENSION(nvm),SAVE :: init, decay      ! D\'esert partout:
20      REAL :: as      alb_neig = 0.55 + 0.3 * EXP(- agesno / 5.)
21      DATA init /0.55, 0.14, 0.18, 0.29, 0.15, 0.15, 0.14, 0./  
22      DATA decay/0.30, 0.67, 0.63, 0.45, 0.40, 0.14, 0.06, 1./      ! Modulation en fonction de l'\^age de la neige :
23        agesno = max((agesno + (1. - agesno / 50.) * dtphys / 86400.) &
24      veget = 0.           * EXP(- MAX(0., snow_fall) * dtphys / 0.3), 0.)
     veget(:,1) = 1.     ! desert partout  
     DO i = 1, knon  
        alb_neig_grid(i) = 0.0  
     ENDDO  
     DO nv = 1, nvm  
        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  
25    
26    END SUBROUTINE albsno    END SUBROUTINE albsno
27    
Legend:
Removed from v.68  
changed lines
  Added in v.305
  ViewVC Help
Powered by ViewVC 1.1.21