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

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

@@ -30 +30 @@
 !      use ioipsl_getincom
       use ioipsl_getincom_p 
+      use mod_phys_lmdz_para, only : is_master
       implicit none
 
@@ -105 +106 @@
       do igas=1,ngas
          read(111,*) gastype(igas)
-         print*,'Gas ',igas,' is ',gastype(igas)
+         if (is_master) print*,'Gas ',igas,' is ',gastype(igas)
       enddo
 
@@ -135 +136 @@
 
       ! display the values
+      if (is_master) then 
       print*,'Variable gas volume mixing ratios:'
       do n=1,L_REFVAR
@@ -141 +143 @@
       end do
       print*,''
+      end if
 
 !=======================================================================
@@ -172 +175 @@
  
       ! display the values
+      if (is_master) then
       print*,'Correlated-k g-space grid:'
       do n=1,L_NGAUSS
@@ -177 +181 @@
       end do
       print*,''
+      end if
 
 !=======================================================================
@@ -210 +215 @@
 
       ! display the values
+      if (is_master) then
       print*,'Correlated-k pressure grid (mBar):'
       do n=1,L_NPREF
@@ -215 +221 @@
       end do
       print*,''
+      end if
 
       ! save the min / max matrix values
@@ -255 +262 @@
 
       ! display the values
+      if (is_master) then
       print*,'Correlated-k temperature grid:'
       do n=1,L_NTREF
          print*,n,'.',tgasref(n),' K'
       end do
+      end if
 
       ! save the min / max matrix values
@@ -275 +284 @@
 
       ! display the values
-      print*,''
-      print*,'Correlated-k matrix size:' 
-      print*,'[',L_NTREF,',',L_NPREF,',',L_REFVAR,',',L_NGAUSS,']' 
+      if (is_master) print*,''
+      if (is_master) print*,'Correlated-k matrix size:' 
+      if (is_master) print*,'[',L_NTREF,',',L_NPREF,',',L_REFVAR,',',L_NGAUSS,']' 
 
 !=======================================================================
@@ -285 +294 @@
 !        wavelength used to separate IR from VI in graybody. We will need that anyway
          IR_VI_wnlimit=3000.
-         write(*,*)"graybody: Visible / Infrared separation set at",10000./IR_VI_wnlimit,"um"
+         if (is_master) write(*,*)"graybody: Visible / Infrared separation set at",10000./IR_VI_wnlimit,"um"
          
          nVI_limit=0
@@ -304 +313 @@
       if (graybody) then
 !        constant absorption coefficient in visible
-         write(*,*)"graybody: constant absorption coefficient in visible:"
+         if (is_master) write(*,*)"graybody: constant absorption coefficient in visible:"
          kappa_VI=-100000.
          call getin_p("kappa_VI",kappa_VI)
-         write(*,*)" kappa_VI = ",kappa_VI
+         if (is_master) write(*,*)" kappa_VI = ",kappa_VI
          kappa_VI=kappa_VI*1.e4* mugaz * 1.672621e-27    ! conversion from m^2/kg to cm^2/molecule         
       
 !        constant absorption coefficient in IR
-         write(*,*)"graybody: constant absorption coefficient in InfraRed:"
+         if (is_master) write(*,*)"graybody: constant absorption coefficient in InfraRed:"
          kappa_IR=-100000.
          call getin_p("kappa_IR",kappa_IR)
-         write(*,*)" kappa_IR = ",kappa_IR       
+         if (is_master) write(*,*)" kappa_IR = ",kappa_IR        
          kappa_IR=kappa_IR*1.e4* mugaz * 1.672621e-27    ! conversion from m^2/kg to cm^2/molecule 
 
-         write(*,*)"graybody: Visible / Infrared separation set at band: IR=",nIR_limit,", VI=",nVI_limit
+         if (is_master) write(*,*)"graybody: Visible / Infrared separation set at band: IR=",nIR_limit,", VI=",nVI_limit
                
       Else
@@ -330 +339 @@
          if ((corrkdir(1:4).eq.'null'))then   !(TRIM(corrkdir).eq.'null_LowTeffStar')) then
             gasv8(1:L_NTREF,1:L_NPREF,1:L_REFVAR,1:L_NSPECTV,1:L_NGAUSS)=0.0
-            print*,'using no corrk data'
-            print*,'Visible corrk gaseous absorption is set to zero if graybody=F'
+            if (is_master) print*,'using no corrk data'
+            if (is_master) print*,'Visible corrk gaseous absorption is set to zero if graybody=F'
          else
             file_id='/corrk_data/'//trim(adjustl(banddir))//'/corrk_gcm_VI.dat' 
@@ -363 +372 @@
          end if
       else
-         print*,'Visible corrk gaseous absorption is set to zero.'
+         if (is_master) print*,'Visible corrk gaseous absorption is set to zero.'
          gasv8(1:L_NTREF,1:L_NPREF,1:L_REFVAR,1:L_NSPECTV,1:L_NGAUSS)=0.0
       endif
@@ -371 +380 @@
       ! INFRA-RED
       if ((corrkdir(1:4).eq.'null'))then       !.or.(TRIM(corrkdir).eq.'null_LowTeffStar')) then
-         print*,'Infrared corrk gaseous absorption is set to zero if graybody=F'
+         if (is_master) print*,'Infrared corrk gaseous absorption is set to zero if graybody=F'
 !$OMP MASTER         
          gasi8(1:L_NTREF,1:L_NPREF,1:L_REFVAR,1:L_NSPECTI,1:L_NGAUSS)=0.0
@@ -662 +671 @@
       endif
 
+      if (is_master) then
       print*,'----------------------------------------------------'
       print*,'And that`s all we have. It`s possible that other'
@@ -667 +677 @@
       print*,'don`t yet have data for it...'
       print*,''
+      end if
 
 !     Deallocate local arrays