New version of the code where the projection merhod is implemented correctly. In addition to that, the predictor corrector method is firstly applied. The intermediate velocities are computed implicitly along the y