crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 13:45:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

read in photon heating number

Browse source
This commit is contained in:
liangjg 2019-03-18 13:56:12 -04:00
parent ce54673f1d
commit 0a18d85a13
2 changed files with 7 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

1
include/openmc/photon.h
View file

@ -67,6 +67,7 @@ public:
xt::xtensor<double, 1> pair_production_total_;
xt::xtensor<double, 1> pair_production_electron_;
xt::xtensor<double, 1> pair_production_nuclear_;
xt::xtensor<double, 1> heating_;
// Form factors
Tabulated1D incoherent_form_factor_;

6
src/photon.cpp
View file

@ -96,6 +96,11 @@ PhotonInteraction::PhotonInteraction(hid_t group, int i_element)
read_dataset(rgroup, "xs", photoelectric_total_);
close_group(rgroup);
// Read heating
rgroup = open_group(group, "heating");
read_dataset(rgroup, "xs", heating_);
close_group(rgroup);
// Read subshell photoionization cross section and atomic relaxation data
rgroup = open_group(group, "subshells");
std::vector<std::string> designators;
@ -280,6 +285,7 @@ PhotonInteraction::PhotonInteraction(hid_t group, int i_element)
xt::log(photoelectric_total_), -500.0);
pair_production_total_ = xt::where(pair_production_total_ > 0.0,
xt::log(pair_production_total_), -500.0);
heating_ = xt::where(heating_ > 0.0, xt::log(heating_), -500.0);
}
void PhotonInteraction::compton_scatter(double alpha, bool doppler,