4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
pure function alboc_cd(rmu0) |
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 |
|
|
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 alboc_cd(size(rmu0)) ! alb\'edo de surface de l'oc\'ean |
real alboc_cd(size(mu0)) ! 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) |
alboc_cd = max(min(0.058 / (max(mu0, 0.) + 0.3), 0.6), 0.04) |
26 |
|
|
27 |
END function alboc_cd |
END function alboc_cd |
28 |
|
|