Reverting trimdist

This commit is contained in:
James E McClure
2015-11-16 09:36:55 -05:00
parent e1fced3818
commit 2e874fc41a
2 changed files with 4 additions and 4 deletions

View File

@@ -248,7 +248,7 @@ void TwoPhase::ComputeDelPhi()
void TwoPhase::Initialize()
{
trimdist=-1.0;
trimdist=1.0;
fluid_isovalue=solid_isovalue=0.0;
// Initialize the averaged quantities
awn = aws = ans = lwns = 0.0;
@@ -273,7 +273,6 @@ void TwoPhase::Initialize()
wwndnw = 0.0;
}
/*
void TwoPhase::SetupCubes(Domain &Dm)
{

View File

@@ -572,8 +572,9 @@ int main(int argc, char **argv)
MPI_Barrier(comm);
//.......................................................................
// Once phase has been initialized, map solid to account for 'smeared' interface
for (i=0; i<N; i++) Averages->SDs(i) -= (1.0); //
for (i=0; i<N; i++) Dm.id[i] = Mask.id[i]; //
for (i=0; i<N; i++) Averages->SDs(i) -= (1.0);
// Make sure the id match for the two domains
for (i=0; i<N; i++) Dm.id[i] = Mask.id[i];
//.......................................................................
// Finalize setup for averaging domain
//Averages->SetupCubes(Mask);