Changeset 2662 for trunk/LMDZ.GENERIC/libf/phystd/interpolateCH4CH4.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/interpolateCH4CH4.F90 (modified) (6 diffs)

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

@@ -16 +16 @@
   use callkeys_mod, only: strictboundcia
   use datafile_mod, only: datadir
+  use mod_phys_lmdz_para, only : is_master
+
   implicit none
 
@@ -57 +59 @@
   if(temp.gt.400)then
     if (strictboundcia) then
-      print*,'Your temperatures are too high for this CH4-CH4 CIA dataset.'
-      print*,'Please run mixed CH4-CH4 atmospheres below T = 400 K.'      
+      if (is_master) then
+        print*,'Your temperatures are too high for this CH4-CH4 CIA dataset.'
+        print*,'Please run mixed CH4-CH4 atmospheres below T = 400 K.'      
+      endif
       stop
     else
-      print*,'Your temperatures are too high for this CH4-CH4 CIA dataset'
-      print*,'you have chosen strictboundcia = ', strictboundcia
-      print*,'*********************************************************'
-      print*,' we allow model to continue but with temp = 400          '
-      print*,'  ...       for CH4-CH4 CIA dataset        ...           '
-      print*,'  ... we assume we know what you are doing ...           '
-      print*,'*********************************************************'
+      if (is_master) then
+        print*,'Your temperatures are too high for this CH4-CH4 CIA dataset'
+        print*,'you have chosen strictboundcia = ', strictboundcia
+        print*,'*********************************************************'
+        print*,' we allow model to continue but with temp = 400          '
+        print*,'  ...       for CH4-CH4 CIA dataset        ...           '
+        print*,'  ... we assume we know what you are doing ...           '
+        print*,'*********************************************************'
+      endif
       temp = 400
     endif
   elseif(temp.lt.40)then
     if (strictboundcia) then
-      print*,'Your temperatures are too low for this CH4-CH4 CIA dataset.'
-      print*,'Please run mixed CH4-CH4 atmospheres above T = 40 K.'      
+      if (is_master) then
+        print*,'Your temperatures are too low for this CH4-CH4 CIA dataset.'
+        print*,'Please run mixed CH4-CH4 atmospheres above T = 40 K.'      
+      endif
       stop
     else
-      print*,'Your temperatures are too low for this CH4-CH4 CIA dataset'
-      print*,'you have chosen strictboundcia = ', strictboundcia
-      print*,'*********************************************************'
-      print*,' we allow model to continue but with temp = 40           '
-      print*,'  ...       for CH4-CH4 CIA dataset        ...           '
-      print*,'  ... we assume we know what you are doing ...           '
-      print*,'*********************************************************'
+      if (is_master) then
+        print*,'Your temperatures are too low for this CH4-CH4 CIA dataset'
+        print*,'you have chosen strictboundcia = ', strictboundcia
+        print*,'*********************************************************'
+        print*,' we allow model to continue but with temp = 40           '
+        print*,'  ...       for CH4-CH4 CIA dataset        ...           '
+        print*,'  ... we assume we know what you are doing ...           '
+        print*,'*********************************************************'
+      endif
       temp = 40
     endif
@@ -90 +100 @@
 
   if(firstcall)then ! called by sugas_corrk only
-     print*,'----------------------------------------------------'
-     print*,'Initialising CH4-CH4 continuum from HITRAN database...'
+     if (is_master) print*,'----------------------------------------------------'
+     if (is_master) print*,'Initialising CH4-CH4 continuum from HITRAN database...'
 
      !     1.1 Open the ASCII files
@@ -99 +109 @@
      open(33,file=dt_file,form='formatted',status='old',iostat=ios)
      if (ios.ne.0) then        ! file not found
-        write(*,*) 'Error from interpolateCH4CH4'
-        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/CH4-CH4_2011.cia is there.'
+        if (is_master) then
+          write(*,*) 'Error from interpolateCH4CH4'
+          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/CH4-CH4_2011.cia is there.'
+        endif
         call abort
      else
@@ -112 +124 @@
            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 interpolateCH4CH4.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 interpolateCH4CH4.F90'
+              endif
               stop
            endif
@@ -129 +143 @@
 !$OMP END MASTER
 !$OMP BARRIER
-
-     print*,'interpolateCH4CH4: At wavenumber ',wn,' cm^-1'
-     print*,'   temperature                 ',temp,' K'
-     print*,'   pressure                    ',pres,' Pa'
-
+     if (is_master) then
+       print*,'interpolateCH4CH4: At wavenumber ',wn,' cm^-1'
+       print*,'   temperature                 ',temp,' K'
+       print*,'   pressure                    ',pres,' Pa'
+     endif
   endif