Changeset 1368

Timestamp:

Nov 17, 2014, 10:14:16 AM (10 years ago)

Author:

emillour

Message:

Mars GCM:
Some cleanup in orbit.F
EM

File:

• trunk/LMDZ.MARS/libf/phymars/orbite.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/orbite.F

@@ -1 +1 @@
       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 +14 @@
 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 +37 @@
 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