Solving 2D Navier Stokes using Finite Difference for staggered mesh
 

Hello Everyone. Before asking this, I do searched the code for the 2D navier stokes equation using finite difference for staggered mesh but I couldn't find any code.
I have the background theoretical knowledge to some extent but my coding knowledge is not really good. I coded 2D heat problem for both the steady and unsteady conditions in FORTRAN.Now, I am trying to solve the 2D NS problem using FORTRAN but I can't even start and can't figure out where to start. It seems really tough. I have written the discretized equation for the problem.I can share it too.
If someone has the code, I request to share it . As a beginner ,it will help me to learn things quickly.

Thank You.

In the book of Peric & Ferziger you find the infos to download some Fortran code

Thank you very much