1 |
guez |
154 |
module alboc_cd_m |
2 |
|
|
|
3 |
guez |
117 |
IMPLICIT NONE |
4 |
|
|
|
5 |
guez |
154 |
contains |
6 |
guez |
117 |
|
7 |
guez |
154 |
SUBROUTINE alboc_cd(rmu0, albedo) |
8 |
guez |
117 |
|
9 |
guez |
154 |
! 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 |
154 |
REAL, intent(in):: rmu0(:) ! cosinus de l'angle solaire z\'enithal |
21 |
|
|
real, intent(out):: albedo(:) ! alb\'edo de surface de l'oc\'ean |
22 |
guez |
117 |
|
23 |
guez |
154 |
! Local: |
24 |
|
|
REAL, PARAMETER:: fmagic = 1. ! facteur magique pour r\'egler l'alb\'edo |
25 |
guez |
117 |
|
26 |
guez |
154 |
!---------------------------------------------------------- |
27 |
guez |
117 |
|
28 |
guez |
154 |
albedo = max(min(fmagic * 0.058 / (max(rmu0, 0.) + 0.3), 0.6), 0.04) |
29 |
guez |
117 |
|
30 |
guez |
154 |
END SUBROUTINE alboc_cd |
31 |
guez |
117 |
|
32 |
guez |
154 |
end module alboc_cd_m |