--- trunk/libf/phylmd/orbite.f90 2009/07/31 15:18:47 22 +++ trunk/libf/phylmd/orbite.f90 2010/12/02 17:11:04 36 @@ -9,7 +9,7 @@ SUBROUTINE orbite(xjour, longi, dist) USE yomcst, ONLY : r_ecc, r_peri - use numer_rec, only: pi + use nr_util, only: pi ! Auteur(s): Z.X. Li (LMD/CNRS) ! Date: 1993/08/18 @@ -64,7 +64,7 @@ USE yomcst, ONLY : r_incl USE phyetat0_m, ONLY : rlat, rlon USE comconst, ONLY : pi - use numer_rec, only: assert + use nr_util, only: assert ! Author : O. Boucher (LMD/CNRS), d'après les routines "zenith" et ! "angle" de Z.X. Li @@ -110,8 +110,9 @@ !---------------------------------------------------------------------- - call assert((/size(pmu0), size(frac)/) == klon, "zenang") - + if (present(frac)) call assert((/size(pmu0), size(frac)/) == klon, & + "zenang") + deux_pi = 2*pi lat_sun = asin(sin(longi * pi / 180.) * sin(r_incl * pi / 180.))