Add boundary conditions velocity, pressure and symmetry. Add loops to read...

Add boundary conditions velocity, pressure and symmetry. Add loops to read implicit cells and nodes to fill the vector grad_J and the matrix Hess_J