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	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	299	SUBROUTINE albsno(agesno, alb_neig, precip_snow)
10	guez	3
11	guez	299	use comconst, only: dtphys
12
13	guez	174	REAL, intent(inout):: agesno(:) ! (knon)
14			real, intent(out):: alb_neig(:) ! (knon)
15			real, intent(in):: precip_snow(:) !(knon)
16	guez	116
17			!------------------------------------------------------------------------
18
19	guez	174	! D\'esert partout:
20			alb_neig = 0.55 + 0.3 * EXP(- agesno / 5.)
21	guez	116
22	guez	174	! Modulation en fonction de l'\^age de la neige :
23	guez	299	agesno = max((agesno + (1. - agesno / 50.) * dtphys / 86400.) &
24			* EXP(- MAX(0., precip_snow) * dtphys / 0.3), 0.)
25	guez	116
26	guez	3	END SUBROUTINE albsno
27
28			end module albsno_m