Changeset 4077 for trunk/LMDZ.GENERIC/libf/phygeneric/ephemeris_orbit_init.F

Timestamp:

Feb 18, 2026, 10:28:23 AM (6 weeks ago)

Author:

mturbet

Message:

renaming and merging of radiative routines in phygeneric

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/ephemeris_orbit_init.F (moved) (moved from trunk/LMDZ.GENERIC/libf/phygeneric/iniorbit.F) (5 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/ephemeris_orbit_init.F

@@ -1 @@
-      SUBROUTINE iniorbit
+      SUBROUTINE ephemeris_orbit_init
      $     (papoastr,pperiastr,pyear_day,pperi_day,pobliq)
       
@@ -32 +32 @@
       peri_day=pperi_day
 
-      PRINT*,'iniorbit: Periastron in AU  ',periastr
-      PRINT*,'iniorbit: Apoastron in AU  ',apoastr 
-      PRINT*,'iniorbit: Obliquity in degrees  :',obliquit
+      PRINT*,'ephemeris_orbit_init: Periastron in AU  ',periastr
+      PRINT*,'ephemeris_orbit_init: Apoastron in AU  ',apoastr 
+      PRINT*,'ephemeris_orbit_init: Obliquity in degrees  :',obliquit
 
 
@@ -40 +40 @@
       p_elips=0.5*(periastr+apoastr)*(1-e_elips*e_elips)
 
-      print*,'iniorbit: e_elips',e_elips
-      print*,'iniorbit: p_elips',p_elips
+      print*,'ephemeris_orbit_init: e_elips',e_elips
+      print*,'ephemeris_orbit_init: p_elips',p_elips
 
 !-----------------------------------------------------------------------
@@ -52 +52 @@
       zanom=2.*pi*(zz-nint(zz))
       zxref=abs(zanom)
-      PRINT*,'iniorbit: zanom  ',zanom
+      PRINT*,'ephemeris_orbit_init: zanom  ',zanom
 
 !  solve equation  zx0 - e * sin (zx0) = zxref for eccentric anomaly zx0
@@ -66 +66 @@
       zx0=zx0+zdx
       if(zanom.lt.0.) zx0=-zx0
-      PRINT*,'iniorbit: zx0   ',zx0
+      PRINT*,'ephemeris_orbit_init: zx0   ',zx0
 
       timeperi=2.*atan(sqrt((1.+e_elips)/(1.-e_elips))*tan(zx0/2.))
-      PRINT*,'iniorbit: Perihelion solar long. Ls (deg)=',
+      PRINT*,'ephemeris_orbit_init: Perihelion solar long. Ls (deg)=',
      &       360.-timeperi*180./pi