diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp index 8b29785a..4cbc3e2a 100644 --- a/models/ColorModel.cpp +++ b/models/ColorModel.cpp @@ -869,7 +869,7 @@ double ScaLBL_ColorModel::MorphOpenConnected(double target_volume_change){ signed char water=2; signed char notwater=1; double SW=-(target_volume_change)/count_connected; - MorphOpen(distance, id_connected, Dm, SW, notwater, water); + MorphOpen(distance, id_connected, Dm, SW, water, notwater); int count_morphopen=0.0; for (int k=0; kComm, count_morphopen); volume_change = double(count_morphopen - count_connected); - if (rank==0) printf(" opening of connected oil %f \n",volume_change); + if (rank==0) printf(" opening of connected oil %f \n",volume_change/count_connected); ScaLBL_CopyToDevice(Phi,phase.data(),N*sizeof(double)); ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm->LastExterior(), Np);