2D lid- driven cavity flow divergence problem
 

I implemented 2D lid driven flow.

simple algorithm,
10x10 grid , collocated grid
reynolds number 100

wall condition,
left , down , right u, v = 0 ,
The upper wall is u = 1.e-4 ,v =0.

Density = 1000 , Viscosity Modulus = 0.001

convection term is
A hybrid scheme was used,

When calculating the correction pressure, rhie-chow interpolation was performed.

But I keep diverge.

I really really don't know.

Please take a look at my code and let me know what the problem is.