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(agesno, alb_neig, precip_snow)

    use comconst, only: dtphys

    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.) * dtphys / 86400.) &
         * EXP(- MAX(0., precip_snow) * dtphys / 0.3), 0.)

  END SUBROUTINE albsno

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