Changeset 2662 for trunk/LMDZ.GENERIC/libf/phystd/interpolateN2N2.F90

Timestamp:

Apr 19, 2022, 12:03:01 PM (3 years ago)

Author:

dbardet

Message:

Using of is_master keyword to allow nly the master processor to write the messages and avoid too heavy output files

File:

• trunk/LMDZ.GENERIC/libf/phystd/interpolateN2N2.F90 (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phystd/interpolateN2N2.F90

@@ -15 +15 @@
 
   use datafile_mod, only: datadir
+  use mod_phys_lmdz_para, only : is_master
+
   implicit none
 
@@ -56 +58 @@
 
   if(temp.gt.400)then
-     print*,'Your temperatures are too high for this N2-N2 CIA dataset.'
-     print*,'Currently, HITRAN provides data for this pair in the range 40-400 K.'
+     if (is_master) then
+       print*,'Your temperatures are too high for this N2-N2 CIA dataset.'
+       print*,'Currently, HITRAN provides data for this pair in the range 40-400 K.'
+     endif
      stop
   endif
@@ -64 +68 @@
 
   if(firstcall)then ! called by sugas_corrk only
-     print*,'----------------------------------------------------'
-     print*,'Initialising N2-N2 continuum from HITRAN database...'
+     if (is_master) print*,'----------------------------------------------------'
+     if (is_master) print*,'Initialising N2-N2 continuum from HITRAN database...'
 
      !     1.1 Open the ASCII files
@@ -73 +77 @@
      open(33,file=dt_file,form='formatted',status='old',iostat=ios)
      if (ios.ne.0) then        ! file not found
-        write(*,*) 'Error from interpolateN2N2'
-        write(*,*) 'Data file ',trim(dt_file),' not found.'
-        write(*,*) 'Check that your path to datagcm:',trim(datadir)
-        write(*,*) 'is correct. You can change it in callphys.def with:'
-        write(*,*) 'datadir = /absolute/path/to/datagcm'
-        write(*,*) 'Also check that the continuum data continuum_data/N2-N2_norm_2011.cia is there.'
+        if (is_master) then
+          write(*,*) 'Error from interpolateN2N2'
+          write(*,*) 'Data file ',trim(dt_file),' not found.'
+          write(*,*) 'Check that your path to datagcm:',trim(datadir)
+          write(*,*) 'is correct. You can change it in callphys.def with:'
+          write(*,*) 'datadir = /absolute/path/to/datagcm'
+          write(*,*) 'Also check that the continuum data continuum_data/N2-N2_norm_2011.cia is there.'
+        endif
         call abort
      else
@@ -86 +92 @@
            read(33,fmat1) bleh,blah,blah,nres,Ttemp
            if(nS.ne.nres)then
-              print*,'Resolution given in file: ',trim(dt_file)
-              print*,'is ',nres,', which does not match nS.'
-              print*,'Please adjust nS value in interpolateN2N2.F90'
+              if (is_master) then
+                print*,'Resolution given in file: ',trim(dt_file)
+                print*,'is ',nres,', which does not match nS.'
+                print*,'Please adjust nS value in interpolateN2N2.F90'
+              endif
               stop
            endif
@@ -103 +111 @@
 !$OMP END MASTER
 !$OMP BARRIER
-
-     print*,'interpolateN2N2: At wavenumber ',wn,' cm^-1'
-     print*,'   temperature ',temp,' K'
-     print*,'   pressure ',pres,' Pa'
-
+     if (is_master) then
+       print*,'interpolateN2N2: At wavenumber ',wn,' cm^-1'
+       print*,'   temperature ',temp,' K'
+       print*,'   pressure ',pres,' Pa'
+     endif
   endif
      call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,temp,abcoef,ind)