Changeset 298 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/iniorbit.F

Timestamp:

12/02/14 19:21:00 (10 years ago)

Author:

milmd

Message:

Less output messages are written. On 20000 cores it is better. In LMDZ, only master of MPI and OpenMP can write.

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/iniorbit.F (modified) (7 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/iniorbit.F

@@ -4 +4 @@
       USE planete_mod, only: apoastr, periastr, year_day, obliquit,
      &                       peri_day, e_elips, p_elips, timeperi
+      USE mod_phys_lmdz_para, only : is_master
       IMPLICIT NONE
 
@@ -67 +68 @@
       !!!! We hope that all cases are above 25 Mkm [OK with Gliese 581d]
       IF ( apoastr .gt. 25.) THEN
+        IF (is_master) THEN
         PRINT*,'!!!!! WARNING !!!!!'
         PRINT*,'!!!!! YOU ARE ABOUT TO WITNESS A DIRT HACK !!!!!'
@@ -73 +75 @@
         PRINT*,'So I am assuming units are in Mkm here'
         PRINT*,'and I am performing a conversion towards AU.'
+        END IF
         periastr = periastr / 149.598 ! Mkm to AU
         apoastr = apoastr / 149.598 ! Mkm to AU
@@ -80 +83 @@
 
  
-      PRINT*,'Periastron in AU  ',periastr
-      PRINT*,'Apoastron in AU  ',apoastr 
-      PRINT*,'Obliquity in degrees  :',obliquit
+      IF (is_master) PRINT*,'Periastron in AU  ',periastr
+      IF (is_master) PRINT*,'Apoastron in AU  ',apoastr 
+      IF (is_master) PRINT*,'Obliquity in degrees  :',obliquit
 
 
@@ -88 +91 @@
       p_elips=0.5*(periastr+apoastr)*(1-e_elips*e_elips)
 
-      print*,'e_elips',e_elips
-      print*,'p_elips',p_elips
+      IF (is_master) print*,'e_elips',e_elips
+      IF (is_master) print*,'p_elips',p_elips
 
 c-----------------------------------------------------------------------
@@ -100 +103 @@
       zanom=2.*pi*(zz-nint(zz))
       zxref=abs(zanom)
-      PRINT*,'zanom  ',zanom
+      IF (is_master) PRINT*,'zanom  ',zanom
 
 c  resolution de l'equation horaire  zx0 - e * sin (zx0) = zxref
@@ -114 +117 @@
       zx0=zx0+zdx
       if(zanom.lt.0.) zx0=-zx0
-      PRINT*,'zx0   ',zx0
+      IF (is_master) PRINT*,'zx0   ',zx0
 
 c zteta est la longitude solaire
 
       timeperi=2.*atan(sqrt((1.+e_elips)/(1.-e_elips))*tan(zx0/2.))
-      PRINT*,'Solar longitude of periastron timeperi = ',timeperi
+      IF (is_master) PRINT*,'Solar longitude of periastron',
+     &          ' timeperi = ',timeperi
 
       RETURN