/[lmdze]/trunk/phylmd/Interface_surf/alboc.f
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/alboc.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 125 by guez, Fri Feb 6 15:00:28 2015 UTC revision 130 by guez, Tue Feb 24 15:43:51 2015 UTC
# Line 4  module alboc_m Line 4  module alboc_m
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    
# Line 18  contains Line 18  contains
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    
# Line 31  contains Line 31  contains
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
# Line 45  contains Line 46  contains
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))
# Line 86  contains Line 87  contains
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))
Legend:
Removed from v.125  
changed lines
  Added in v.130
  ViewVC Help
Powered by ViewVC 1.1.21