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

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

@@ -2 +2 @@
 
 !==================================================================
-!     
+!
 !     Purpose
 !     -------
 !     Average the chosen high resolution stellar spectrum over the
 !     visible bands in the model.
-!     
+!
 !     Authors
-!     ------- 
+!     -------
 !     Robin Wordsworth (2010).
 !     Generalized to very late spectral types (and Brown dwarfs) Jeremy Leconte (2012)
@@ -16 +16 @@
 !     ---------
 !     setspv.F
-!     
+!
 !     Calls
 !     -----
 !     none
-!     
+!
 !==================================================================
 
@@ -28 +28 @@
       use callkeys_mod, only: stelbbody,stelTbb,startype
       use ioipsl_getin_p_mod, only: getin_p
+      use mod_phys_lmdz_para, only : is_master
 
       implicit none
 
       real*8 STELLAR(L_NSPECTV)
-!      integer startype 
+!      integer startype
 
       integer Nfine
@@ -48 +49 @@
       real lam_temp
       double precision stel_temp
-      
+
       integer :: ios ! file opening/reading status
 
       STELLAR(:)=0.0
 
-      print*,'enter ave_stellspec'
+      if (is_master) print*,'enter ave_stellspec'
       if(stelbbody)then
          tstellar=stelTbb
@@ -65 +66 @@
                call blackl(dble(lam_temp*1e-6),dble(tstellar),stel_temp)
                STELLAR(band)=STELLAR(band)+stel_temp*dl
-            enddo           
+            enddo
          end do
          STELLAR(1:L_NSPECTV)=STELLAR(1:L_NSPECTV)/sum(STELLAR(1:L_NSPECTV))
@@ -99 +100 @@
            Case(6)
             file_id='/stellar_spectra/BD_Teff-1600K.txt'
-            tstellar=1600. 
+            tstellar=1600.
             file_id_lam='/stellar_spectra/lamBD.txt'
             Nfine=5000
            Case(7)
             file_id='/stellar_spectra/BD_Teff-1000K.txt'
-            tstellar=1000. 
+            tstellar=1000.
             file_id_lam='/stellar_spectra/lamBD.txt'
             Nfine=5000
            Case(8)
             file_id='/stellar_spectra/Flux_K5_Teff4700_logg4.5_Met-0.5_BTsettle.dat'
-            tstellar=4700. 
+            tstellar=4700.
             file_id_lam='/stellar_spectra/lambda_K5_Teff4700_logg4.5_Met-0.5_BTsettle.dat'
             Nfine=3986
            Case(9)
             file_id='/stellar_spectra/Flux_TRAPPIST1.dat'
-            tstellar=2550. 
+            tstellar=2550.
             file_id_lam='/stellar_spectra/lambda_TRAPPIST1.dat'
             Nfine=5000
            Case(10)
             file_id='/stellar_spectra/Flux_Proxima.dat'
-            tstellar=3050. 
+            tstellar=3050.
             file_id_lam='/stellar_spectra/lambda_Proxima.dat'
             Nfine=5000
@@ -139 +140 @@
               stop
             end if
-   
+
             ! Opening file
             file_path = trim(datadir)//'/stellar_spectra/'//stelspec_file
             print*, 'stellar flux : ', file_path
             OPEN(UNIT=110,FILE=file_path,STATUS='old',iostat=ios)
-   
+
             if (ios /= 0) THEN
               write(*,*)'Error: cannot open stelspec_file file ', trim(stelspec_file)
@@ -219 +220 @@
 !$OMP END MASTER
 !$OMP BARRIER
-         
+
          ! sum data by band
          band=1
@@ -236 +237 @@
             STELLAR(band)=STELLAR(band)+stel_f(ifine)*dl
          end do
-               
-         
+
+
          STELLAR(1:L_NSPECTV)=STELLAR(1:L_NSPECTV)/sum(STELLAR(1:L_NSPECTV))
 !$OMP BARRIER
@@ -244 +245 @@
          if (allocated(stel_f)) deallocate(stel_f)
 !$OMP END MASTER
-!$OMP BARRIER         
+!$OMP BARRIER
       endif