Hi all, I would like to add forces at the INS equation (a volume force as gravity, and also a force proportional to the velocity filed), and define new boundary conditions. In particular, I need to impose at the inflow and outflow boundaries (top and bottom respectively) the velocity profiles to be u(x,t)=0; v(x,t) = a*(cosh(b*x)/cosh(b*c) - 1); being a,b,c constants. I've found in ins/src the whole source code, but what and where I have to modify exactly? (lineSolveINS.h? Cgins.C? ins(bC).C? ...) Thanks, Regards Alessandro