- Timestamp:
- 12/02/14 19:21:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/setspv.F90
r227 r298 27 27 STELLARF,TAURAY 28 28 use datafile_mod, only: datadir 29 use mod_phys_lmdz_para, only : is_master 29 30 30 31 implicit none … … 82 83 close(131) 83 84 84 write(*,*) 'setspv: L_NSPECTV = ',L_NSPECTV, 'in the model '85 write(*,*) ' there are ',nb, 'entries in ',TRIM(file_path)85 if (is_master) write(*,*) 'setspv: L_NSPECTV = ',L_NSPECTV, 'in the model ' 86 if (is_master) write(*,*) ' there are ',nb, 'entries in ',TRIM(file_path) 86 87 if(nb.ne.L_NSPECTV) then 87 88 write(*,*) 'MISMATCH !! I stop here' … … 102 103 !$OMP BARRIER 103 104 105 if (is_master) then 104 106 print*,'setspv: VI band limits:' 105 107 do M=1,L_NSPECTV+1 … … 107 109 end do 108 110 print*,' ' 111 end if 109 112 110 113 ! Set up mean wavenumbers and wavenumber deltas. Units of … … 123 126 ! Set up stellar spectrum 124 127 125 write(*,*)'setspv: Interpolating stellar spectrum from the hires data...'128 if (is_master) write(*,*)'setspv: Interpolating stellar spectrum from the hires data...' 126 129 call ave_stelspec(STELLAR) 127 130 … … 132 135 sum = sum+STELLARF(N) 133 136 end do 134 write(6,'("setspv: Stellar flux at 1 AU = ",f7.2," W m-2")') sum135 print*,' '137 if (is_master) write(6,'("setspv: Stellar flux at 1 AU = ",f7.2," W m-2")') sum 138 if (is_master) print*,' ' 136 139 137 140 … … 144 147 call calc_rayleigh 145 148 else 146 print*,'setspv: No Rayleigh scattering, check for NaN in output!'149 if (is_master) print*,'setspv: No Rayleigh scattering, check for NaN in output!' 147 150 do N=1,L_NSPECTV 148 151 TAURAY(N) = 1E-16
Note: See TracChangeset
for help on using the changeset viewer.