Commit a62d979a authored by Philippe Delandmeter's avatar Philippe Delandmeter
Browse files

slim3d: bug fixed for bottom friction in slim3d

parent 0335a833
Pipeline #1763 failed with stage
in 43 minutes and 1 second
......@@ -274,7 +274,9 @@ void dgConservationLawSW3dMomentum::boundaryBottom(functorCache &m, fullMatrix<d
double uBot = zBotDist(i,0)*uvgrad(i,2) + uv(i,0);
double vBot = zBotDist(i,0)*uvgrad(i,5) + uv(i,1);
double Cd = slim3dParameters::vonKarman/log(zBotDist(i,0)/z0(i,0));
Cd = Cd * Cd * hypot(uBotE,vBotE)*normals(i,2);
Cd = - Cd * Cd * hypot(uBotE,vBotE)*normals(i,2);
if (normals(i,2) > 0)
Msg::Fatal("nz should be negative here!!");
val(i,0) = -Cd * uBot;
val(i,1) = -Cd * vBot;
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment