4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE alboc(rjour, rlat, albedo) |
SUBROUTINE alboc(jour, rlat, albedo) |
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 |
|
|
18 |
USE orbite_m, ONLY: orbite |
USE orbite_m, ONLY: orbite |
19 |
|
|
20 |
! Arguments; |
! Arguments; |
21 |
! rjour (in, R) : jour dans l'annee (a compter du 1 janvier) |
! jour (in) : jour dans l'annee (a compter du 1 janvier) |
22 |
! rlat (in, R) : latitude en degre |
! rlat (in, R) : latitude en degre |
23 |
! albedo (out, R): albedo obtenu (de 0 a 1) |
! albedo (out, R): albedo obtenu (de 0 a 1) |
24 |
|
|
31 |
INTEGER npts ! il controle la precision de l'integration |
INTEGER npts ! il controle la precision de l'integration |
32 |
PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes |
PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes |
33 |
|
|
34 |
REAL rlat(klon), rjour, albedo(klon) |
integer jour |
35 |
|
REAL rlat(klon), albedo(klon) |
36 |
REAL zdist, zlonsun, zpi, zdeclin |
REAL zdist, zlonsun, zpi, zdeclin |
37 |
REAL rmu, alb, srmu, salb, fauxo, aa, bb |
REAL rmu, alb, srmu, salb, fauxo, aa, bb |
38 |
INTEGER i, k |
INTEGER i, k |
46 |
zpi = 4.*atan(1.) |
zpi = 4.*atan(1.) |
47 |
|
|
48 |
! Calculer la longitude vraie de l'orbite terrestre: |
! Calculer la longitude vraie de l'orbite terrestre: |
49 |
CALL orbite(rjour, zlonsun, zdist) |
CALL orbite(real(jour), zlonsun, zdist) |
50 |
|
|
51 |
! Calculer la declinaison du soleil (qui varie entre + et - R_incl): |
! Calculer la declinaison du soleil (qui varie entre + et - R_incl): |
52 |
zdeclin = asin(sin(zlonsun*zpi/180.0)*sin(r_incl*zpi/180.0)) |
zdeclin = asin(sin(zlonsun*zpi/180.0)*sin(r_incl*zpi/180.0)) |
87 |
zpi = 4.*atan(1.) |
zpi = 4.*atan(1.) |
88 |
|
|
89 |
! Calculer la longitude vraie de l'orbite terrestre: |
! Calculer la longitude vraie de l'orbite terrestre: |
90 |
CALL orbite(rjour, zlonsun, zdist) |
CALL orbite(real(jour), zlonsun, zdist) |
91 |
|
|
92 |
! Calculer la declinaison du soleil (qui varie entre + et - R_incl): |
! Calculer la declinaison du soleil (qui varie entre + et - R_incl): |
93 |
zdeclin = asin(sin(zlonsun*zpi/180.0)*sin(r_incl*zpi/180.0)) |
zdeclin = asin(sin(zlonsun*zpi/180.0)*sin(r_incl*zpi/180.0)) |