never rescale unsteady cases

This commit is contained in:
James McClure committed 2021-11-07 07:56:18 -05:00
1 parent ce57491bce
commit 1fe39b8223
1 file changed
+2 -4
+2 -4
View File
@@ -707,7 +707,7 @@ double ScaLBL_ColorModel::Run(int returntime){
analysis.basic(timestep, current_db, *Averages, Phi, Pressure, Velocity, fq, Den ); // allow initial ramp-up to get closer to steady state
CURRENT_TIMESTEP += 2;
if (CURRENT_TIMESTEP > MIN_STEADY_TIMESTEPS){
if (CURRENT_TIMESTEP > MIN_STEADY_TIMESTEPS && BC == 0){
analysis.finish();
double volB = Averages->gwb.V;
@@ -744,9 +744,7 @@ double ScaLBL_ColorModel::Run(int returntime){
isSteady = true;
if (CURRENT_TIMESTEP >= MAX_STEADY_TIMESTEPS)
isSteady = true;
if (BC != 0 ) isSteady = false; // Never rescale unsteady simulation cases
if (isSteady && (Ca > maxCa || Ca < minCa) && SET_CAPILLARY_NUMBER ){
/* re-run the point if the actual Ca is too far from the target Ca */
isSteady = false;