phylmd/Interface_surf/albsno.f90

module albsno_m

  ! From phylmd/interface_surf.F90, version 1.8 2005/05/25 13:10:09

  IMPLICIT none

contains

  SUBROUTINE albsno(dtime, agesno, alb_neig, precip_snow)

    REAL, intent(in):: dtime
    REAL, intent(inout):: agesno(:) ! (knon)
    real, intent(out):: alb_neig(:) ! (knon)
    real, intent(in):: precip_snow(:) !(knon)

    !------------------------------------------------------------------------

    ! D\'esert partout:
    alb_neig = 0.55 + 0.3 * EXP(- agesno / 5.)

    ! Modulation en fonction de l'\^age de la neige :
    agesno = max((agesno + (1. - agesno / 50.) * dtime / 86400.) &
         * EXP(- MAX(0., precip_snow) * dtime / 0.3), 0.)

  END SUBROUTINE albsno

end module albsno_m
1	guez	3	module albsno_m
2
3	guez	116	! From phylmd/interface_surf.F90, version 1.8 2005/05/25 13:10:09
4	guez	3
5			IMPLICIT none
6
7			contains
8
9	guez	174	SUBROUTINE albsno(dtime, agesno, alb_neig, precip_snow)
10	guez	3
11	guez	174	REAL, intent(in):: dtime
12			REAL, intent(inout):: agesno(:) ! (knon)
13			real, intent(out):: alb_neig(:) ! (knon)
14			real, intent(in):: precip_snow(:) !(knon)
15	guez	116
16			!------------------------------------------------------------------------
17
18	guez	174	! D\'esert partout:
19			alb_neig = 0.55 + 0.3 * EXP(- agesno / 5.)
20	guez	116
21	guez	174	! Modulation en fonction de l'\^age de la neige :
22			agesno = max((agesno + (1. - agesno / 50.) * dtime / 86400.) &
23			* EXP(- MAX(0., precip_snow) * dtime / 0.3), 0.)
24	guez	116
25	guez	3	END SUBROUTINE albsno
26
27			end module albsno_m