Changeset 4192 for dynamico_lmdz/simple_physics/phyparam/physics/astronomy.F90

Timestamp:

Dec 20, 2019, 12:56:43 AM (5 years ago)

Author:

dubos

Message:

simple_physics : cleanup astronomy

File:

• dynamico_lmdz/simple_physics/phyparam/physics/astronomy.F90 (modified) (2 diffs)

dynamico_lmdz/simple_physics/phyparam/physics/astronomy.F90

@@ -7 +7 @@
   REAL, PARAMETER :: unitastr=149.597927, & ! millions of km
        pi=2.*ASIN(1.)
+
 CONTAINS
  
@@ -166 +167 @@
   END SUBROUTINE iniorbit
 
+  PURE SUBROUTINE orbite(pls,pdist_sol,pdecli)
+    !=======================================================================
+    !
+    !   Objet:
+    !   ------
+    !
+    !   Distance from sun and declimation as a function of the solar
+    !   longitude Ls
+    !
+    !   Arguments:
+    !   ----------
+    !
+    !   Input:
+    !   ------
+    !   pls          Ls
+    !
+    !   Output:
+    !   -------
+    !   pdist_sol     Distance Sun-Planet in UA
+    !   pdecli        declinaison ( en radians )
+    !
+    !=======================================================================
+    !-----------------------------------------------------------------------
+    !   Declarations:
+    !   -------------
+    
+    ! arguments:
+    ! ----------
+
+    REAL, INTENT(IN) :: pls
+    REAL, INTENT(OUT) :: pdist_sol,pdecli
+    
+    !-----------------------------------------------------------------------
+    
+    ! Distance Sun-Planet
+    
+    pdist_sol=p_elips/(1.+e_elips*cos(pls+timeperi))
+    
+    ! Solar declination
+    
+    pdecli= asin (sin(pls)*sin(obliquit*pi/180.))
+    
+  END SUBROUTINE orbite
+
 END MODULE astronomy