phylmd/Interface_surf/albsno.f90

module albsno_m

  IMPLICIT none

contains

  SUBROUTINE albsno(agesno, alb_neig, snow_fall)

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

    use comconst, only: dtphys

    REAL, intent(inout):: agesno(:) ! (knon) age of snow, in days
    real, intent(out):: alb_neig(:) ! (knon)

    real, intent(in):: snow_fall(:) !(knon)
    ! precipitation, solid water mass flux (kg / m2 / s), positive down

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

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

    ! Modulation en fonction de l'\^age de la neige :
    agesno = (agesno + (1. - agesno / 50.) * dtphys / 86400.) &
         * EXP(- snow_fall * dtphys / 0.3)

  END SUBROUTINE albsno

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