crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Add derivative info to tallies.out

Browse source
This commit is contained in:
Sterling Harper 2016-10-21 02:56:15 -04:00
parent 9313d8720b
commit f0467e35f4
1 changed files with 20 additions and 15 deletions
Download patch file Download diff file Expand all files Collapse all files

35
src/output.F90
View file

@ -778,21 +778,26 @@ contains
endif
! Write derivative information.
!if (allocated(t % deriv)) then
! select case (t % deriv % variable)
! case (DIFF_DENSITY)
! write(unit=unit_tally, fmt="(' Density derivative Material ',A)") &
! to_str(t % deriv % diff_material)
! case (DIFF_NUCLIDE_DENSITY)
! i_listing = nuclides(t % deriv % diff_nuclide) % listing
! write(unit=unit_tally, fmt="(' Nuclide density derivative Material '&
! &,A,' Nuclide ',A)") trim(to_str(t % deriv % diff_material)), &
! trim(xs_listings(i_listing) % alias)
! case default
! call fatal_error("Differential tally dependent variable for tally " &
! // trim(to_str(t % id)) // " not defined in output.F90.")
! end select
!end if
if (t % deriv /= NONE) then
associate(deriv => tally_derivs(t % deriv))
select case (deriv % variable)
case (DIFF_DENSITY)
write(unit=unit_tally, fmt="(' Density derivative Material ',A)") &
to_str(deriv % diff_material)
case (DIFF_NUCLIDE_DENSITY)
write(unit=unit_tally, fmt="(' Nuclide density derivative &
&Material ',A,' Nuclide ',A)") &
trim(to_str(deriv % diff_material)), &
trim(nuclides(deriv % diff_nuclide) % name)
case (DIFF_TEMPERATURE)
write(unit=unit_tally, fmt="(' Temperature derivative Material ',&
&A)") to_str(deriv % diff_material)
case default
call fatal_error("Differential tally dependent variable for tally "&
// trim(to_str(t % id)) // " not defined in output.F90.")
end select
end associate
end if
! Handle surface current tallies separately
if (t % type == TALLY_SURFACE_CURRENT) then