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

Teach load_state_point about internal tallies

Browse source 
When reading in tally data from a statepoint file, a check
is performed on the "internal" attribute of the tally group.
If the attribute doesn't exist (previous statepoint files), then
it is assumed that the tally is writeable / not internal. Otherwise,
the internal attribute is read from the dataset.

If the tally is determined to be internal, the value of
tally.writeable_ is set to false and no tally data
is read
This commit is contained in:
Andrew Johnson 2019-09-20 08:48:27 -05:00
parent 2ead24e172
commit 1300ca2678
No known key found for this signature in database
GPG key ID: 253418E91B7F6FEB
1 changed files with 18 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

23
src/state_point.cpp
View file

@ -433,11 +433,23 @@ void load_state_point()
// Read sum, sum_sq, and N for each bin
std::string name = "tally " + std::to_string(tally->id_);
hid_t tally_group = open_group(tallies_group, name.c_str());
auto& results = tally->results_;
read_tally_results(tally_group, results.shape()[0],
results.shape()[1], results.data());
read_dataset(tally_group, "n_realizations", tally->n_realizations_);
close_group(tally_group);
int internal;
if (attribute_exists(tally_group, "internal")) {
read_attribute(tally_group, "internal", internal);
} else {
internal = 0;
}
if (internal) {
tally->writeable_ = false;
} else {
auto& results = tally->results_;
read_tally_results(tally_group, results.shape()[0],
results.shape()[1], results.data());
read_dataset(tally_group, "n_realizations", tally->n_realizations_);
close_group(tally_group);
}
}
close_group(tallies_group);
@ -705,6 +717,7 @@ void write_tally_results_nr(hid_t file_id)
for (const auto& t : model::tallies) {
// Skip any tallies that are not active
if (!t->active_) continue;
if (!t->writeable_) continue;
if (mpi::master && !object_exists(file_id, "tallies_present")) {
write_attribute(file_id, "tallies_present", 1);