Changeset 3040 for trunk/LMDZ.MARS/libf/phymars/orbite.F

Timestamp:

Sep 12, 2023, 11:41:22 AM (15 months ago)

Author:

jbclement

Message:

Mars PCM:
The variable 'timeperi' (defined in "planete_h.F90" and computed in "iniorbit.F") is renamed into 'lsperi' and thus slightly changed to be coherent to the solar longitude of perihelion in radian. It can now be used out of the box by other subroutines/programs like the PEM.
JBC

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, ONLY: e_elips, p_elips, obliquit, timeperi
+      USE planete_h, ONLY: e_elips, p_elips, obliquit, lsperi
       USE comcstfi_h, ONLY: pi
       IMPLICIT NONE
@@ -19 +19 @@
 c   -------
 c   pdist_sol     Distance Sun-Planet in UA
-c   pdecli        declinaison ( in radians )
+c   pdecli        Declinaison ( in radians )
 c
 c=======================================================================
@@ -26 +26 @@
 c ----------
 
-      REAL,INTENT(IN) :: pls
-      REAL,INTENT(OUT) :: pdist_sol,pdecli
+      REAL,INTENT(IN)  :: pls
+      REAL,INTENT(OUT) :: pdist_sol, pdecli
 
 c-----------------------------------------------------------------------
 
 c Distance Sun-Planet
-
-      pdist_sol=p_elips/(1.+e_elips*cos(pls+timeperi))
+      pdist_sol = p_elips/(1.+e_elips*cos(pls-lsperi))
 
 c Solar declination
-
-      pdecli= asin (sin(pls)*sin(obliquit*pi/180.))
+      pdecli = asin(sin(pls)*sin(obliquit*pi/180.))
 
       END