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

Timestamp:

Dec 20, 2019, 10:02:32 AM (5 years ago)

Author:

dubos

Message:

simple_physics : created plugin for logging

File:

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

dynamico_lmdz/simple_physics/phyparam/physics/astronomy.F90

@@ -1 +1 @@
 MODULE astronomy
+  
+#include "use_logging.h"
+
   IMPLICIT NONE
   SAVE
@@ -79 +82 @@
     
     IF (lwrite) THEN
-       PRINT*,'jour de l"annee   :',pday
-       PRINT*,'distance au soleil (en unite astronomique) :',pdist_sol
-       PRINT*,'declinaison (en degres) :',pdecli*180./pi
+       WRITELOG(*,*) 'jour de l"annee   :',pday
+       WRITELOG(*,*) 'distance au soleil (en unite astronomique) :',pdist_sol
+       WRITELOG(*,*) 'declinaison (en degres) :',pdecli*180./pi
+       CALL flush_log
     ENDIF
     
@@ -124 +128 @@
     !-----------------------------------------------------------------------
     
-    PRINT*,'Perihelie en Mkm  ',periheli
-    PRINT*,'Aphelise  en Mkm  ',aphelie 
-    PRINT*,'obliquite en degres  :',obliquit
+    WRITELOG(*,*) 'Perihelie en Mkm  ',periheli
+    WRITELOG(*,*) 'Aphelise  en Mkm  ',aphelie 
+    WRITELOG(*,*) 'obliquite en degres  :',obliquit
     
     e_elips=(aphelie-periheli)/(periheli+aphelie)
     p_elips=0.5*(periheli+aphelie)*(1-e_elips*e_elips)/unitastr
     
-    print*,'e_elips',e_elips
-    print*,'p_elips',p_elips
-    
+    WRITELOG(*,*) 'e_elips',e_elips
+    WRITELOG(*,*) 'p_elips',p_elips
+
     !-----------------------------------------------------------------------
     ! calcul de l'angle polaire et de la distance au soleil :
@@ -143 +147 @@
     zanom=2.*pi*(zz-nint(zz))
     zxref=abs(zanom)
-    PRINT*,'zanom  ',zanom
+    WRITELOG(*,*) 'zanom  ',zanom
     
     !  resolution de l'equation horaire  zx0 - e * sin (zx0) = zxref
@@ -157 +161 @@
     zx0=zx0+zdx
     if(zanom.lt.0.) zx0=-zx0
-    PRINT*,'zx0   ',zx0
+    WRITELOG(*,*) 'zx0   ',zx0
     
     ! zteta est la longitude solaire
     
     timeperi=2.*atan(sqrt((1.+e_elips)/(1.-e_elips))*tan(zx0/2.))
-    PRINT*,'longitude solaire du perihelie timeperi = ',timeperi
-    
+    WRITELOG(*,*) 'longitude solaire du perihelie timeperi = ',timeperi
+    
+    CALL flush_log
     
   END SUBROUTINE iniorbit