crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Moved default Watt spectrum.

Browse source
This commit is contained in:
Paul Romano 2012-07-28 22:16:02 -04:00
parent 1133b69011
commit 8564a203dc
2 changed files with 15 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

3
src/input_xml.F90
View file

@ -234,7 +234,10 @@ contains
call fatal_error()
end select
else
! Set default energy distribution to Watt fission spectrum
external_source % type_energy = SRC_ENERGY_WATT
allocate(external_source % params_energy(2))
external_source % params_energy = (/ 0.988_8, 2.249_8 /)
end if
! Survival biasing

19
src/source.F90
View file

@ -102,12 +102,14 @@ contains
type(Bank), intent(inout) :: site ! source site
integer :: i ! dummy loop index
real(8) :: r(3) ! sampled coordinates
real(8) :: phi ! azimuthal angle
real(8) :: mu ! cosine of polar angle
real(8) :: p_min(3) ! minimum coordinates of source
real(8) :: p_max(3) ! maximum coordinates of source
integer :: i ! dummy loop index
real(8) :: r(3) ! sampled coordinates
real(8) :: phi ! azimuthal angle
real(8) :: mu ! cosine of polar angle
real(8) :: p_min(3) ! minimum coordinates of source
real(8) :: p_max(3) ! maximum coordinates of source
real(8) :: Watt_a ! Watt spectrum parameter 'a'
real(8) :: Watt_b ! Watt spectrum parameter 'b'
! Sample position
select case (external_source % type_space)
@ -136,9 +138,12 @@ contains
! Sample energy
select case (external_source % type_energy)
case (SRC_ENERGY_WATT)
Watt_a = external_source % params_energy(1)
Watt_b = external_source % params_energy(2)
do
! Sample Watt fission spectrum
site % E = watt_spectrum(0.988_8, 2.249_8)
site % E = watt_spectrum(Watt_a, Watt_b)
! resample if energy is >= 20 MeV
if (site % E < 20) exit