Changeset 3504 for trunk/LMDZ.PLUTO/libf/phypluto/setspi.F90

Timestamp:

Nov 8, 2024, 10:57:14 AM (13 days ago)

Author:

afalco

Message:

Pluto: print only on master process.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/setspi.F90 (modified) (14 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/setspi.F90

@@ -2 +2 @@
 
 !==================================================================
-!     
+!
 !     Purpose
 !     -------
 !     Set up spectral intervals and Planck function in the longwave.
-!     
+!
 !     Authors
-!     ------- 
+!     -------
 !     Adapted from setspi in the NASA Ames radiative code by
 !     Robin Wordsworth (2009).
-!     
+!
 !     Called by
 !     ---------
 !     callcorrk.F
-!     
+!
 !     Calls
 !     -----
 !     none
-!     
+!
 !==================================================================
 
@@ -26 +26 @@
       use datafile_mod, only: datadir
       use comcstfi_mod, only: pi
+      use mod_phys_lmdz_para, only : is_master
 
       implicit none
@@ -42 +43 @@
       real*8 :: c1 = 3.741832D-16 ! W m^-2
       real*8 :: c2 = 1.438786D-2  ! m K
-      
+
       real*8 :: lastband(2), plancksum
 
@@ -75 +76 @@
       write(temp1,'(i2.2)') L_NSPECTI
       !file_id='/corrk_data/' // corrkdir(1:LEN_TRIM(corrkdir)) // '/narrowbands_IR.in'
-      file_id='/corrk_data/'//trim(adjustl(banddir))//'/narrowbands_IR.in' 
+      file_id='/corrk_data/'//trim(adjustl(banddir))//'/narrowbands_IR.in'
       file_path=TRIM(datadir)//TRIM(file_id)
 
@@ -81 +82 @@
       inquire(FILE=file_path,EXIST=file_ok)
       if(.not.file_ok) then
-         write(*,*)'The file ',TRIM(file_path)
-         write(*,*)'was not found by setspi.F90, exiting.'
-         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 corrkdir you chose in callphys.def exists.'
-         call abort
-      endif
-    
-!$OMP MASTER    
+         if (is_master) then
+            write(*,*)'The file ',TRIM(file_path)
+            write(*,*)'was not found by setspi.F90, exiting.'
+            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 corrkdir you chose in callphys.def exists.'
+         endif
+         call abort_physic("setspi", "File not found by setspi", 1)
+      endif
+
+!$OMP MASTER
       nb=0
       ierr=0
-      ! check that the file contains the right number of bands 
+      ! check that the file contains the right number of bands
       open(131,file=file_path,form='formatted')
       read(131,*,iostat=ierr) file_entries
@@ -103 +106 @@
       close(131)
 
-      write(*,*) 'setspi: L_NSPECTI = ',L_NSPECTI, 'in the model '
-      write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)
+      if (is_master) then
+         write(*,*) 'setspi: L_NSPECTI = ',L_NSPECTI, 'in the model '
+         write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)
+      endif
       if(nb.ne.L_NSPECTI) then
-         write(*,*) 'MISMATCH !! I stop here'
-         call abort
+         call abort_physic("setspi",'MISMATCH !! I stop here',1)
       endif
 
       ! load and display the data
       open(111,file=file_path,form='formatted')
-      read(111,*) 
+      read(111,*)
       do M=1,L_NSPECTI-1
          read(111,*) BWNI(M)
@@ -123 +127 @@
 !$OMP BARRIER
 
-      print*,''
-      print*,'setspi: IR band limits:'
-      do M=1,L_NSPECTI+1
-         print*,m,'-->',BWNI(M),' cm^-1'
-      end do
-
-!     Set up mean wavenumbers and wavenumber deltas.  Units of 
+      if (is_master)then
+         print*,''
+         print*,'setspi: IR band limits:'
+         do M=1,L_NSPECTI+1
+            print*,m,'-->',BWNI(M),' cm^-1'
+         end do
+      endif
+
+!     Set up mean wavenumbers and wavenumber deltas.  Units of
 !     wavenumbers is cm^(-1); units of wavelengths is microns.
 
@@ -136 +142 @@
          DWNI(M)  = BWNI(M+1)-BWNI(M)
          WAVEI(M) = 1.0D+4/WNOI(M)
-         BLAMI(M) = 0.01D0/BWNI(M)         
+         BLAMI(M) = 0.01D0/BWNI(M)
       end do
       BLAMI(M) = 0.01D0/BWNI(M)
@@ -147 +153 @@
 !     original planck.f; W m^-2 wavenumber^-1, where wavenumber is in CM^-1.
 
-      print*,''
-      print*,'setspi: Current Planck integration range:'
-      print*,'T = ',dble(NTstart)/NTfac, ' to ',dble(NTstop)/NTfac,' K.'
+      if (is_master)then
+         print*,''
+         print*,'setspi: Current Planck integration range:'
+         print*,'T = ',dble(NTstart)/NTfac, ' to ',dble(NTstop)/NTfac,' K.'
+      endif
 
       IF(.NOT.ALLOCATED(planckir)) ALLOCATE(planckir(L_NSPECTI,NTstop-NTstart+1))
@@ -159 +167 @@
          bma = (b-a)/2.0D0
          ! if (nw .eq. 25) then !LT debug
-         !    print*, "a = ",a 
-         !    print*, "b= ",b 
-         !    print*,"bpa = ",bpa 
+         !    print*, "a = ",a
+         !    print*, "b= ",b
+         !    print*,"bpa = ",bpa
          !    print*, "bma = ",bma
          ! endif
@@ -170 +178 @@
                y    = bma*x(mm)+bpa
                !to avoid floating overflow when T is low and optical wavelength
-               if ((c2/(y*T)) .lt. 700.0D0) then 
+               if ((c2/(y*T)) .lt. 700.0D0) then
                   ans  = ans + w(mm)*c1/(y**5*(exp(c2/(y*T))-1.0D0))
-               else 
+               else
                   ans = ans +0.0D0
                endif
@@ -179 +187 @@
          end do
       end do
-         
+
       ! 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=NTstart,NTstop
             plancksum=0.0D0
             T=dble(NT)/NTfac
-       
+
             do NW=1,L_NSPECTI
                plancksum=plancksum+  &
@@ -207 +215 @@
             plancksum=plancksum+planckir(NW,nt-NTstart+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) then
+            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'
+         endif
          ! check energy conservation at upper temperature boundary
          plancksum=0.0D0
@@ -217 +226 @@
             plancksum=plancksum+planckir(NW,nt-NTstart+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) then
+            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*,''
+         endif
       endif