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