--- trunk/libf/phylmd/orbite.f90	2009/07/31 15:18:47	22
+++ trunk/libf/phylmd/orbite.f90	2009/12/14 15:25:16	23
@@ -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.))