phylmd/Interface_surf/alboc_cd.f

module alboc_cd_m

  IMPLICIT NONE

contains

  pure function alboc_cd(mu0)

    ! From LMDZ4/libf/phylmd/albedo.F, version 1.2, 2005/02/07 15:00:52

    ! Author: Z. X. Li (LMD/CNRS)
    ! Date: 1994/06/24

    ! Calculer l'alb\'edo sur l'oc\'ean en fonction de l'angle
    ! z\'enithal moyen. Formule due \`a Larson and Barkstrom,
    ! Proceedings of the symposium on radiation in the atmosphere,
    ! 19-28 August 1976, science Press, 1977, pages 451-453, ou
    ! th\`ese de 3\`eme cycle de Sylvie Joussaume.

    REAL, intent(in):: mu0(:) ! cosinus de l'angle solaire z\'enithal
    real alboc_cd(size(mu0)) ! alb\'edo de surface de l'oc\'ean

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

    alboc_cd = max(min(0.058 / (max(mu0, 0.) + 0.3), 0.6), 0.04)

  END function alboc_cd

end module alboc_cd_m
1	guez	154	module alboc_cd_m
2
3	guez	117	IMPLICIT NONE
4
5	guez	154	contains
6	guez	117
7	guez	300	pure function alboc_cd(mu0)
8	guez	117
9	guez	208	! From LMDZ4/libf/phylmd/albedo.F, version 1.2, 2005/02/07 15:00:52
10	guez	117
11	guez	154	! Author: Z. X. Li (LMD/CNRS)
12			! Date: 1994/06/24
13	guez	117
14	guez	154	! Calculer l'alb\'edo sur l'oc\'ean en fonction de l'angle
15			! z\'enithal moyen. Formule due \`a Larson and Barkstrom,
16			! Proceedings of the symposium on radiation in the atmosphere,
17			! 19-28 August 1976, science Press, 1977, pages 451-453, ou
18			! th\`ese de 3\`eme cycle de Sylvie Joussaume.
19	guez	117
20	guez	300	REAL, intent(in):: mu0(:) ! cosinus de l'angle solaire z\'enithal
21			real alboc_cd(size(mu0)) ! alb\'edo de surface de l'oc\'ean
22	guez	117
23	guez	154	!----------------------------------------------------------
24	guez	117
25	guez	300	alboc_cd = max(min(0.058 / (max(mu0, 0.) + 0.3), 0.6), 0.04)
26	guez	117
27	guez	208	END function alboc_cd
28	guez	117
29	guez	154	end module alboc_cd_m