Changeset 298 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/setspi.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/setspi.F90 (modified) (8 diffs)

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

@@ -25 +25 @@
       use radcommon_h, only: BWNI,BLAMI,WNOI,DWNI,WAVEI,planckir,sigma
       use datafile_mod, only: datadir
+      use mod_phys_lmdz_para, only : is_master
 
       implicit none
@@ -105 +106 @@
       close(131)
 
-      write(*,*) 'setspi: L_NSPECTI = ',L_NSPECTI, 'in the model '
-      write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)
+      if (is_master) write(*,*) 'setspi: L_NSPECTI = ',L_NSPECTI, 'in the model '
+      if (is_master) write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)
       if(nb.ne.L_NSPECTI) then
          write(*,*) 'MISMATCH !! I stop here'
@@ -125 +126 @@
 !$OMP BARRIER
 
+      if (is_master) then
       print*,''
       print*,'setspi: IR band limits:'
@@ -130 +132 @@
          print*,m,'-->',BWNI(M),' cm^-1'
       end do
+      end if
 
 !     Set up mean wavenumbers and wavenumber deltas.  Units of 
@@ -149 +152 @@
 !     original planck.f; W m^-2 wavenumber^-1, where wavenumber is in CM^-1.
 
+      if (is_master) then
       print*,''
       print*,'setspi: Current Planck integration range:'
       print*,'T = ',dble(NTstar)/NTfac, ' to ',dble(NTstop)/NTfac,' K.'
+      end if
 
       do NW=1,L_NSPECTI
@@ -173 +178 @@
       ! force planck=sigma*eps*T^4 for each temperature in array
       if(forceEC)then
-         print*,'setspi: Force F=sigma*eps*T^4 for all values of T!'
+         if (is_master) print*,'setspi: Force F=sigma*eps*T^4 for all values of T!'
          do nt=NTstar,NTstop
             plancksum=0.0D0
@@ -198 +203 @@
             plancksum=plancksum+planckir(NW,nt-NTstar+1)*DWNI(NW)*pi
          end do
-         print*,'setspi: At lower limit:'
-         print*,'in model sig*T^4 = ',plancksum,' W m^-2'
-         print*,'actual sig*T^4   = ',sigma*(dble(nt)/NTfac)**4,' W m^-2'
+         if (is_master) print*,'setspi: At lower limit:'
+         if (is_master) print*,'in model sig*T^4 = ',plancksum,' W m^-2'
+         if (is_master) print*,'actual sig*T^4   = ',sigma*(dble(nt)/NTfac)**4,' W m^-2'
          
          ! check energy conservation at upper temperature boundary
@@ -208 +213 @@
             plancksum=plancksum+planckir(NW,nt-NTstar+1)*DWNI(NW)*pi
          end do
-         print*,'setspi: At upper limit:'
-         print*,'in model sig*T^4 = ',plancksum,' W m^-2'
-         print*,'actual sig*T^4   = ',sigma*(dble(nt)/NTfac)**4,' W m^-2'
-         print*,''
+         if (is_master) print*,'setspi: At upper limit:'
+         if (is_master) print*,'in model sig*T^4 = ',plancksum,' W m^-2'
+         if (is_master) print*,'actual sig*T^4   = ',sigma*(dble(nt)/NTfac)**4,' W m^-2'
+         if (is_master) print*,''
       endif