phylmd/Interface_surf/alboc_cd.f

module alboc_cd_m

  IMPLICIT NONE

contains

  pure function alboc_cd(rmu0)

    ! 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):: rmu0(:) ! cosinus de l'angle solaire z\'enithal
    real alboc_cd(size(rmu0)) ! alb\'edo de surface de l'oc\'ean

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

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

  END function alboc_cd

end module alboc_cd_m
1	module alboc_cd_m
2
3	IMPLICIT NONE
4
5	contains
6
7	pure function alboc_cd(rmu0)
8
9	! From LMDZ4/libf/phylmd/albedo.F, version 1.2, 2005/02/07 15:00:52
10
11	! Author: Z. X. Li (LMD/CNRS)
12	! Date: 1994/06/24
13
14	! 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
20	REAL, intent(in):: rmu0(:) ! cosinus de l'angle solaire z\'enithal
21	real alboc_cd(size(rmu0)) ! alb\'edo de surface de l'oc\'ean
22
23	!----------------------------------------------------------
24
25	alboc_cd = max(min(0.058 / (max(rmu0, 0.) + 0.3), 0.6), 0.04)
26
27	END function alboc_cd
28
29	end module alboc_cd_m