natural convection trhrough an 2D enclosure with partially active wall

Hello LBM Friends,
I need a code matlab software based on lattice boltzmann method(D2Q9 coupled with D2Q5) to simulate heat transfer (free convection) in an 2D enclosure with partially active wall
the code i need should compute : steamlines contour , isotherms and local Nusselt number and also average Nusselt number
please add a comments to the lines of the code ,i hope understand it .
i am waiting please thank you so much.

hi friends,
i am waiting plz give me the code or suggest me a website ^^