Changeset 3504 for trunk/LMDZ.PLUTO/libf/phypluto/setspv.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/setspv.F90 (modified) (9 diffs)

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

@@ -2 +2 @@
 
 !==================================================================
-!     
+!
 !     Purpose
 !     -------
-!     Set up spectral intervals, stellar spectrum and Rayleigh 
-!     opacity in the shortwave. 
-!     
+!     Set up spectral intervals, stellar spectrum and Rayleigh
+!     opacity in the shortwave.
+!
 !     Authors
-!     ------- 
+!     -------
 !     Adapted from setspv in the NASA Ames radiative code by
 !     Robin Wordsworth (2009).
@@ -16 +16 @@
 !     ---------
 !     callcorrk.F
-!     
+!
 !     Calls
 !     -----
 !     ave_stelspec.F
-!     
+!
 !==================================================================
 
@@ -28 +28 @@
       use datafile_mod, only: datadir
       use callkeys_mod, only: Fat1AU,rayleigh
+      use mod_phys_lmdz_para, only : is_master
 
       implicit none
@@ -53 +54 @@
 
       write(temp1,'(i2.2)') L_NSPECTV
-      file_id='/corrk_data/'//trim(adjustl(banddir))//'/narrowbands_VI.in' 
+      file_id='/corrk_data/'//trim(adjustl(banddir))//'/narrowbands_VI.in'
       file_path=TRIM(datadir)//TRIM(file_id)
 
@@ -65 +66 @@
          write(*,*)' datadir = /absolute/path/to/datagcm'
          write(*,*)'Also check that the corrkdir you chose in callphys.def exists.'
-         call abort
+         call abort_physic("setspv", "File not found by setspi", 1)
       endif
-        
-!$OMP MASTER        
+
+!$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
@@ -80 +81 @@
       close(131)
 
-      write(*,*) 'setspv: L_NSPECTV = ',L_NSPECTV, 'in the model '
-      write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)
+      if (is_master) then
+         write(*,*) 'setspv: L_NSPECTV = ',L_NSPECTV, 'in the model '
+         write(*,*) '        there are   ',nb, 'entries in ',TRIM(file_path)
+      endif
       if(nb.ne.L_NSPECTV) then
-         write(*,*) 'MISMATCH !! I stop here'
-         call abort
+         call abort_physic("setspv",'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_NSPECTV-1
          read(111,*) BWNV(M)
@@ -100 +102 @@
 !$OMP BARRIER
 
-      print*,'setspv: VI band limits:'
-      do M=1,L_NSPECTV+1
-         print*,m,'-->',BWNV(M),' cm^-1'
-      end do
-      print*,' '
+      if (is_master) then
+         print*,'setspv: VI band limits:'
+         do M=1,L_NSPECTV+1
+            print*,m,'-->',BWNV(M),' cm^-1'
+         end do
+         print*,' '
+      endif
 
-!     Set up mean wavenumbers and wavenumber deltas.  Units of 
+!     Set up mean wavenumbers and wavenumber deltas.  Units of
 !     wavenumbers is cm^(-1); units of wavelengths is microns.
 
@@ -121 +125 @@
 !     Set up stellar spectrum
 
-      write(*,*)'setspv: Interpolating stellar spectrum from the hires data...'
+      if (is_master) write(*,*)'setspv: Interpolating stellar spectrum from the hires data...'
       call ave_stelspec(STELLAR)
 
-!     Sum the stellar flux, and write out the result.  
-      sum = 0.0  
+!     Sum the stellar flux, and write out the result.
+      sum = 0.0
       do N=1,L_NSPECTV
          STELLARF(N) = STELLAR(N) * Fat1AU
          sum         = sum+STELLARF(N)
       end do
-      write(6,'("setspv: Stellar flux at 1 AU = ",f9.2," W m-2")') sum
-      print*,' '
+      if (is_master)then
+         write(6,'("setspv: Stellar flux at 1 AU = ",f9.2," W m-2")') sum
+         print*,' '
+      endif
 
 
@@ -142 +148 @@
          call calc_rayleigh
       else
-         print*,'setspv: No Rayleigh scattering, check for NaN in output!'
+         if (is_master) print*,'setspv: No Rayleigh scattering, check for NaN in output!'
          do N=1,L_NSPECTV
             TAURAY(N) = 1E-16