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

Diff of /trunk/phylmd/orbite.f

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

trunk/libf/phylmd/orbite.f90 revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC trunk/phylmd/orbite.f revision 82 by guez, Wed Mar 5 14:57:53 2014 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 63  CONTAINS Line 63  CONTAINS
63      USE dimphy, ONLY : klon      USE dimphy, ONLY : klon
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 nr_util, only: assert, pi
     use numer_rec, only: assert  
67    
68      ! Author : O. Boucher (LMD/CNRS), d'après les routines "zenith" et      ! Author : O. Boucher (LMD/CNRS), d'après les routines "zenith" et
69      ! "angle" de Z.X. Li      ! "angle" de Z.X. Li
# Line 110  CONTAINS Line 109  CONTAINS
109    
110      !----------------------------------------------------------------------      !----------------------------------------------------------------------
111    
112      call assert((/size(pmu0), size(frac)/) == klon, "zenang")      if (present(frac)) call assert((/size(pmu0), size(frac)/) == klon, &
113             "zenang")
114        
115      deux_pi = 2*pi      deux_pi = 2*pi
116    
117      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.82
  ViewVC Help
Powered by ViewVC 1.1.21