Flip ID
This commit is contained in:
@@ -486,8 +486,8 @@ int main(int argc, char **argv)
|
||||
FlipID(id,Nx*Ny*Nz);
|
||||
}
|
||||
#else
|
||||
// If negative phi_s is chosen, flip the ID for the wetting and non-wetting phase
|
||||
if (phi_s < 0.0 && !pBC){
|
||||
// If positive phi_s is chosen, flip the ID for the wetting and non-wetting phase
|
||||
if (phi_s > 0.0 && !pBC){
|
||||
phi_s = -phi_s;
|
||||
das = (phi_s+1.0)*0.5;
|
||||
dbs = 1.0 - das;
|
||||
|
||||
@@ -507,8 +507,8 @@ int main(int argc, char **argv)
|
||||
FlipID(id,Nx*Ny*Nz);
|
||||
}
|
||||
#else
|
||||
// If negative phi_s is chosen, flip the ID for the wetting and non-wetting phase
|
||||
if (phi_s < 0.0 && !pBC){
|
||||
// If positive phi_s is chosen, flip the ID for the wetting and non-wetting phase
|
||||
if (phi_s > 0.0 && !pBC){
|
||||
phi_s = -phi_s;
|
||||
das = (phi_s+1.0)*0.5;
|
||||
dbs = 1.0 - das;
|
||||
|
||||
Reference in New Issue
Block a user