Changeset 3040 for trunk/LMDZ.MARS/libf

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

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• call_dayperi.F (modified) (3 diffs)
• iniorbit.F (modified) (2 diffs)
• orbite.F (modified) (3 diffs)
• planete_h.F90 (modified) (1 diff)
• solarlong.F (modified) (1 diff)

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

@@ -15 +15 @@
 c   Input:
 c   ------
-c   Lsperi      solar longitude (Ls) of perohelion (rad)
-c   e_elips       Excentricity
-c   real year_day      ! number of sols per Mars yar
+c   Lsperi        ! Solar longitude (Ls) of perihelion (rad)
+c   e_elips       ! Excentricity
+c   real year_day ! Number of sols per Mars year
 c
 c   output
 c   ------
-c   dayperi       Martian date at perihelion (sol)
+c   dayperi       ! Martian date at perihelion (sol)
 
 c-----------------------------------------------------------------------
@@ -27 +27 @@
 c arguments:
 c ----------
-       real Lsperi,e_elips,dayperi
-       real year_day      ! number of sols per Mars yar
-
+       real Lsperi, e_elips, dayperi
+       real year_day ! Number of sols per Mars year
 
 c      Local
@@ -43 +42 @@
      &   ( 2*atan(x1*tan(0.5*Lsperi))
      &     -x2*sin(Lsperi)/(1+e_elips*cos(Lsperi)) )
-       if(dayperi.lt.0) dayperi=dayperi+year_day
+       if (dayperi < 0) dayperi=dayperi+year_day
        return
+
        end

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

@@ -2 +2 @@
      $     (paphelie,pperiheli,pyear_day,pperi_day,pobliq)
       use planete_h, only: aphelie, periheli, year_day, peri_day,
-     &                     obliquit, unitastr, e_elips, p_elips,
-     &                     timeperi
+     &                     obliquit, unitastr, e_elips, p_elips, lsperi
       use comcstfi_h, only: pi
       IMPLICIT NONE
@@ -67 +66 @@
       PRINT*,'iniorbit: zx0   ',zx0
 
-      timeperi=2.*atan(sqrt((1.+e_elips)/(1.-e_elips))*tan(zx0/2.))
+      lsperi=-2.*atan(sqrt((1.+e_elips)/(1.-e_elips))*tan(zx0/2.))
+      lsperi = modulo(lsperi,2.*pi)
       PRINT*,'iniorbit: Perihelion solar long. Ls (deg)=',
-     &       modulo(-timeperi*180./pi,360.)
+     &       lsperi*180./pi
 
       END

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

trunk/LMDZ.MARS/libf/phymars/planete_h.F90

@@ -2 +2 @@
       IMPLICIT NONE
 
-      REAL aphelie  ! Aphelion, in Mkm
-      REAL periheli ! Perihelion, in Mkm
-      REAL year_day ! number of days in the year
-      REAL peri_day ! date of perihelion, in days
-      REAL obliquit ! obliquity of the planet, in degrees
+      REAL aphelie   ! Aphelion, in Mkm
+      REAL periheli  ! Perihelion, in Mkm
+      REAL year_day  ! Number of days in the year
+      REAL peri_day  ! Date of perihelion, in days
+      REAL obliquit  ! Obliquity of the planet, in degrees
       REAL lmixmin
       REAL emin_turb
       REAL coefvis
       REAL coefir
-      REAL timeperi ! angle of the perihelion, in rad
-      REAL e_elips ! orbit eccentricity
-      REAL p_elips
-      REAL unitastr ! Astronomical unit AU, in Mkm
+      REAL lsperi    ! Solar longitude of the perihelion, angle in rad
+      REAL e_elips   ! Orbit eccentricity
+      REAL p_elips   ! Ellipse semi-latus rectum
+      REAL unitastr  ! Astronomical unit AU, in Mkm
 
       END MODULE planete_h

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

@@ -84 +84 @@
       zteta=2.*atan(sqrt((1.+e_elips)/(1.-e_elips))*tan(zx0/2.))
 
-      psollong=zteta-timeperi
+      psollong=zteta+lsperi
 
       IF(psollong.LT.0.) psollong=psollong+2.*pi