Changeset 298 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/calc_rayleigh.F90

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/calc_rayleigh.F90 (modified) (5 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/calc_rayleigh.F90

@@ -28 +28 @@
       use radcommon_h, only: WAVEV, BWNV, DWNV, tstellar, tauray, taurayvar, scalep
       use gases_h
+      use mod_phys_lmdz_para, only : is_master
 
       implicit none
@@ -55 +56 @@
       do igas=1,ngasmx
          if(igas.eq.vgas)then
-            print*,'variable gas is ',trim(gnom(igas)),' in Rayleigh scattering '
+            if (is_master) print*,'variable gas is ',trim(gnom(igas)),' in Rayleigh scattering '
          endif
          if((igas/=vgas).and.(gfrac(igas).lt.5.e-2))then
-            print*,'Ignoring ',trim(gnom(igas)),' in Rayleigh scattering '// &
+            if (is_master) print*,'Ignoring ',trim(gnom(igas)),' in Rayleigh scattering '// &
             'as its mixing ratio is less than 0.05.' 
             ! ignore variable gas in Rayleigh calculation
@@ -84 +85 @@
 
             if((gfrac(igas).eq.1.0).and.(vgas.eq.0))then
-               print*,'Rayleigh scattering is for a pure ',trim(gnom(igas)),' atmosphere.'
+               if (is_master) print*,'Rayleigh scattering is for a pure ',trim(gnom(igas)),' atmosphere.'
                typeknown=.true.
             endif
@@ -92 +93 @@
 
       if(.not.typeknown)then
-         print*,'Rayleigh scattering is for a mixed gas atmosphere.'
+         if (is_master) print*,'Rayleigh scattering is for a mixed gas atmosphere.'
          typeknown=.true.
       endif
@@ -159 +160 @@
       IF (L_NSPECTV > 6) THEN
         icantbewrong = L_NSPECTV-6
-        print*,'At 1 atm and lambda = ',WAVEV(icantbewrong),' um'
-        print*,'tau_R = ',TAURAY(icantbewrong)*1013.25
-        print*,'sig_R = ',TAURAY(icantbewrong)*g*mugaz*1.67e-27*100, &
+        if (is_master) print*,'At 1 atm and lambda = ',WAVEV(icantbewrong),' um'
+        if (is_master) print*,'tau_R = ',TAURAY(icantbewrong)*1013.25
+        if (is_master) print*,'sig_R = ',TAURAY(icantbewrong)*g*mugaz*1.67e-27*100, &
                'cm^2 molecule^-1'
       ENDIF