-
Notifications
You must be signed in to change notification settings - Fork 98
Open
Description
Describe the bug
See title.
To Reproduce
Steps to reproduce the behavior:
- Run
conj_htexample. - Restart
conj_htexample, for only one timestep.
Expected behavior
The solution looks roughly like from the previous run.
Actual behavior
The solid domain is very close to zero.
Fix:
diff --git a/src/app/nrs/nrs.cpp b/src/app/nrs/nrs.cpp
index cb6a14b8..3231024a 100644
--- a/src/app/nrs/nrs.cpp
+++ b/src/app/nrs/nrs.cpp
@@ -641,7 +641,7 @@ void nrs_t::restartFromFiles(const std::vector<std::string> &fileList)
for (int i = scalarStart; i < Nscalar; i++) {
const auto sid = scalarDigitStr(i - scalarStart);
if (checkOption("s" + sid) && isAvailable("scalar" + sid)) {
- auto o_Si = scalar->o_S.slice(scalar->fieldOffsetScan[i], meshV->Nlocal);
+ auto o_Si = scalar->o_S.slice(scalar->fieldOffsetScan[i], scalar->mesh(i)->Nlocal);
std::vector<occa::memory> o_iofldSi = {o_Si};
iofld->addVariable("scalar" + sid, o_iofldSi);
}Commit: b5b876d
Metadata
Metadata
Assignees
Labels
No labels