4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE alboc_cd(rmu0, albedo) |
pure function alboc_cd(mu0) |
8 |
|
|
9 |
! From LMDZ4/libf/phylmd/albedo.F, version 1.2 2005/02/07 15:00:52 |
! From LMDZ4/libf/phylmd/albedo.F, version 1.2, 2005/02/07 15:00:52 |
10 |
|
|
11 |
! Author: Z. X. Li (LMD/CNRS) |
! Author: Z. X. Li (LMD/CNRS) |
12 |
! Date: 1994/06/24 |
! Date: 1994/06/24 |
17 |
! 19-28 August 1976, science Press, 1977, pages 451-453, ou |
! 19-28 August 1976, science Press, 1977, pages 451-453, ou |
18 |
! th\`ese de 3\`eme cycle de Sylvie Joussaume. |
! th\`ese de 3\`eme cycle de Sylvie Joussaume. |
19 |
|
|
20 |
REAL, intent(in):: rmu0(:) ! cosinus de l'angle solaire z\'enithal |
REAL, intent(in):: mu0(:) ! cosinus de l'angle solaire z\'enithal |
21 |
real, intent(out):: albedo(:) ! alb\'edo de surface de l'oc\'ean |
real alboc_cd(size(mu0)) ! alb\'edo de surface de l'oc\'ean |
|
|
|
|
! Local: |
|
|
REAL, PARAMETER:: fmagic = 1. ! facteur magique pour r\'egler l'alb\'edo |
|
22 |
|
|
23 |
!---------------------------------------------------------- |
!---------------------------------------------------------- |
24 |
|
|
25 |
albedo = max(min(fmagic * 0.058 / (max(rmu0, 0.) + 0.3), 0.6), 0.04) |
alboc_cd = max(min(0.058 / (max(mu0, 0.) + 0.3), 0.6), 0.04) |
26 |
|
|
27 |
END SUBROUTINE alboc_cd |
END function alboc_cd |
28 |
|
|
29 |
end module alboc_cd_m |
end module alboc_cd_m |