devs should still be set, add function to set pressure at top wall
add pressure set in jacobi loop
set wall devs value in consolidate function
remove comment
added Dirichlet and Neumann BC at the dynamic top wall
remove p debug output
remove dp_dz output
corrected BCs at top wall
old formulation of wall stress