/[lmdze]/trunk/phylmd/orbite.f
ViewVC logotype

Diff of /trunk/phylmd/orbite.f

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

revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC revision 36 by guez, Thu Dec 2 17:11:04 2010 UTC
# Line 9  CONTAINS Line 9  CONTAINS
9    SUBROUTINE orbite(xjour, longi, dist)    SUBROUTINE orbite(xjour, longi, dist)
10    
11      USE yomcst, ONLY : r_ecc, r_peri      USE yomcst, ONLY : r_ecc, r_peri
12      use numer_rec, only: pi      use nr_util, only: pi
13    
14      ! Auteur(s): Z.X. Li (LMD/CNRS)      ! Auteur(s): Z.X. Li (LMD/CNRS)
15      ! Date: 1993/08/18      ! Date: 1993/08/18
# Line 64  CONTAINS Line 64  CONTAINS
64      USE yomcst, ONLY : r_incl      USE yomcst, ONLY : r_incl
65      USE phyetat0_m, ONLY : rlat, rlon      USE phyetat0_m, ONLY : rlat, rlon
66      USE comconst, ONLY : pi      USE comconst, ONLY : pi
67      use numer_rec, only: assert      use nr_util, only: assert
68    
69      ! Author : O. Boucher (LMD/CNRS), d'après les routines "zenith" et      ! Author : O. Boucher (LMD/CNRS), d'après les routines "zenith" et
70      ! "angle" de Z.X. Li      ! "angle" de Z.X. Li
# Line 110  CONTAINS Line 110  CONTAINS
110    
111      !----------------------------------------------------------------------      !----------------------------------------------------------------------
112    
113      call assert((/size(pmu0), size(frac)/) == klon, "zenang")      if (present(frac)) call assert((/size(pmu0), size(frac)/) == klon, &
114             "zenang")
115        
116      deux_pi = 2*pi      deux_pi = 2*pi
117    
118      lat_sun = asin(sin(longi * pi / 180.) * sin(r_incl * pi / 180.))      lat_sun = asin(sin(longi * pi / 180.) * sin(r_incl * pi / 180.))
Legend:
Removed from v.22  
changed lines
  Added in v.36
  ViewVC Help
Powered by ViewVC 1.1.21