tweak morphlogy routine

This commit is contained in:
James E McClure 2018-12-17 12:33:02 -05:00
parent 3e1ed52929
commit 00dec1714a

View File

@ -359,6 +359,7 @@ double MorphGrow(DoubleArray &BoundaryDist, DoubleArray &Dist, Array<char> &id,
if (rank == 0) printf(" delta=%f, growth=%f, max. displacement = %f \n",morph_delta, GrowthEstimate, MAX_DISPLACEMENT);
// Now adjust morph_delta
morph_delta *= TargetGrowth/GrowthEstimate;
MAX_DISPLACEMENT *= TargetGrowth/GrowthEstimate
if (MAX_DISPLACEMENT > 2.0 ){
morph_delta /= 0.5*MAX_DISPLACEMENT;
COUNT_FOR_LOOP = 100; // exit loop if displacement is too large