Changeset 4194

Timestamp:

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

Author:

dubos

Message:

simple_physics : created plugin for logging

Location:

dynamico_lmdz/simple_physics/phyparam/physics

Files:

• astronomy.F90 (modified) (5 diffs)
• logging.F90 (added)
• use_logging.h (added)
• vdif_mod.F90 (modified) (4 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

dynamico_lmdz/simple_physics/phyparam/physics/vdif_mod.F90

@@ -1 +1 @@
 MODULE vdif_mod
+  
+#include "use_logging.h"
+
   IMPLICIT NONE
   SAVE
   PRIVATE
-  
+
   REAL, PARAMETER :: karman=0.4
   REAL :: lmixmin=100., emin_turb=1e-8  
@@ -249 +252 @@
     IF(lwrite) THEN
        ig=ngrid/2+1
-       PRINT*,'Pression (mbar) ,altitude (km),u,v,theta, rho dz'
+       WRITELOG(*,*) 'Pression (mbar) ,altitude (km),u,v,theta, rho dz'
        DO ilay=1,nlay
-          WRITE(*,*) .01*pplay(ig,ilay),.001*pzlay(ig,ilay), &
+          WRITELOG(*,*) .01*pplay(ig,ilay),.001*pzlay(ig,ilay), &
                pu(ig,ilay),pv(ig,ilay),ph(ig,ilay),za(ig,ilay)
        ENDDO
-       PRINT*,'Pression (mbar) ,altitude (km),zb'
+       WRITELOG(*,*) 'Pression (mbar) ,altitude (km),zb'
        DO ilev=1,nlay
-          WRITE(*,*) .01*pplev(ig,ilev),.001*pzlev(ig,ilev), &
+          WRITELOG(*,*) .01*pplev(ig,ilev),.001*pzlev(ig,ilev), &
                zb0(ig,ilev)
        ENDDO
+       CALL flush_log
     ENDIF
     
@@ -287 +291 @@
     
     IF(lwrite) THEN
-       PRINT*
-       PRINT*,'Diagnostique diffusion verticale'
-       print*,'LMIXMIN',lmixmin
+       
+       WRITELOG(*,*) 'Diagnostique diffusion verticale'
+       WRITELOG(*,*) 'LMIXMIN',lmixmin
        PRINT*,'coefficients Cd pour v et h'
        PRINT*,zcdv(ngrid/2+1),zcdh(ngrid/2+1)
@@ -296 +300 @@
           PRINT*,zkv(ngrid/2+1,ilev),zkh(ngrid/2+1,ilev)
        ENDDO
+       CALL flush_log
     ENDIF