Changeset 3974 for trunk/LMDZ.MARS/libf/phymars/tabfi.F

Timestamp:

Nov 21, 2025, 4:28:06 PM (5 weeks ago)

Author:

jbclement

Message:

Mars PCM:
Simplification of the orbital parameters initialization (resolve ticket #109): removing redundant actions in subroutine 'iniorbit' in regard of what is done in "tabfi.F" and moving the subroutine 'lsp2solp' from "tabfi.F" to "planete_h.F90" + some cleanings.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/tabfi.F (modified) (4 diffs)

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

@@ -61 +61 @@
       use time_phylmdz_mod, only: daysec, dtphys
       use planete_h, only: aphelie, emin_turb, lmixmin, obliquit,
-     &                     peri_day, periheli, year_day
+     &                     peri_day, periheli, year_day,
+     &                     lsp2solp, iniorbit
+
       implicit none
  
@@ -460 +462 @@
           write(*,*) 'perihelion =',periheli
           write(*,*) 'and',year_day,'sols/year'
-          call lsp2solp(peri_ls,peri_day,aphelie,periheli,year_day)
+          call lsp2solp(peri_ls,peri_day)
           write(*,*) 'peri_day (new value):',peri_day
 
@@ -543 +545 @@
 
 c-----------------------------------------------------------------------
+c       Initialization of orbital parameters
+c-----------------------------------------------------------------------
+      call iniorbit()
+
+c-----------------------------------------------------------------------
 c       Case when using a start file from before March 1996 (without iceradius...
 c-----------------------------------------------------------------------
@@ -555 +562 @@
        end if
 
-c-----------------------------------------------------------------------
       END SUBROUTINE tabfi
 
-
-
-      
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-! gives sol at perihelion for ls at perihelion (for precession cycles)
-      subroutine lsp2solp(lsp,solp,aphelie,periheli,year_day)
-
-      implicit none
-!  Arguments:
-      real lsp     ! Input: ls at perihelion
-      real solp    ! Output: sol at perihelion
-      real aphelie,periheli,year_day ! Input: parameters
- 
-!  Local:
-      double precision zx0 ! eccentric anomaly at Ls=0
-      double precision e_elips
-      double precision pi,degrad 
-      
-      parameter (pi=4.d0*atan(1.d0))
-      parameter (degrad=180.d0/pi)
-
-      e_elips=(aphelie-periheli)/(aphelie+periheli)      
-      zx0 = -2.0*datan(dtan(0.5*lsp/degrad)
-     .          *dsqrt((1.-e_elips)/(1.+e_elips)))
-      if (zx0 .le. 0.) zx0 = zx0 + 2.*pi
-      
-      solp  = year_day*(1.-(zx0-e_elips*dsin(zx0))/(2.*pi))
-
-
-      end subroutine lsp2solp
-
-
       END MODULE tabfi_mod
-