fix growth estimation

This commit is contained in:
James E McClure 2018-12-17 14:05:53 -05:00
parent 0553309fc5
commit ad553399d7

View File

@ -335,7 +335,7 @@ double MorphGrow(DoubleArray &BoundaryDist, DoubleArray &Dist, Array<char> &id,
double MAX_DISPLACEMENT = 0.0;
int COUNT_FOR_LOOP = 0;
if (rank == 0) printf("Estimate delta for growth=%f \n",TargetGrowth);
while ( fabs(GrowthEstimate - TargetGrowth)/TargetGrowth < 0.01) && COUNT_FOR_LOOP < 10 ){
while ( fabs(GrowthEstimate - TargetGrowth)/TargetGrowth < 0.01 && COUNT_FOR_LOOP < 10 ){
COUNT_FOR_LOOP++;
count = 0.0;
for (int k=1; k<Nz-1; k++){