Index: /trunk/LMDZ.MARS/libf/phymars/orbite.F
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/orbite.F	(revision 1367)
+++ /trunk/LMDZ.MARS/libf/phymars/orbite.F	(revision 1368)
@@ -1,19 +1,11 @@
       SUBROUTINE orbite(pls,pdist_sol,pdecli)
-      use planete_h
-      USE comcstfi_h
+      USE planete_h, ONLY: e_elips, p_elips, obliquit, timeperi
+      USE comcstfi_h, ONLY: pi
       IMPLICIT NONE
 
 c=======================================================================
 c
-c   Objet:
-c   ------
-c
-c   Distance from sun and declimation as a function of the solar
+c   Compute distance to Sun and declination as a function of the solar
 c   longitude Ls
-c
-c   Interface:
-c   ----------
-c
-c
 c
 c   Arguments:
@@ -22,20 +14,18 @@
 c   Input:
 c   ------
-c   pls          Ls
+c   pls          Ls (radians)
 c
 c   Output:
 c   -------
 c   pdist_sol     Distance Sun-Planet in UA
-c   pdecli        declinaison ( en radians )
+c   pdecli        declinaison ( in radians )
 c
 c=======================================================================
-c-----------------------------------------------------------------------
-c   Declarations:
-c   -------------
 
 c arguments:
 c ----------
 
-      REAL pdist_sol,pdecli,pls,i
+      REAL,INTENT(IN) :: pls
+      REAL,INTENT(OUT) :: pdist_sol,pdecli
 
 c-----------------------------------------------------------------------
@@ -47,16 +37,5 @@
 c Solar declination
 
-c ********************* version before 01/01/2000 *******
-
       pdecli= asin (sin(pls)*sin(obliquit*pi/180.))
 
-c ********************* version after 01/01/2000 *******
-c     i=obliquit*pi/180.
-c     pdecli=asin(sin(pls)*sin(i)/sqrt(sin(pls)**2+
-c    & cos(pls)**2*cos(i)**2))
-c ******************************************************
-
-
-
-      RETURN
       END