I want to applly RLBM in lid driven cavity.
Is this format is correct?
DO i=0,n
DO j=0,m
sxx=0.0
sxy=0.0
syy=0.0
DO k=0,8
sxx=sxx+cx(k)cx(k)(f(k,i,j)-feq(k,i,j))
sxy=sxy+cx(k)cy(k)(f(k,i,j)-feq(k,i,j))
syy=syy+cy(k)cy(k)(f(k,i,j)-feq(k,i,j))
f1(k,i,j)=4.5w(k)((cx(k)*cx(k)-1./3.)*sxx+(cy(k)*cy(k)-1./3.)syy+2.0(cx(k)*cy(k))*sxy)
END DOEND DOEND DODO i=0,n
DO j=0,m
t1=u(i,j)*u(i,j)+v(i,j)*v(i,j)
DO k=0,8
t2=u(i,j)*cx(k)+v(i,j)*cy(k)
feq(k,i,j)=w(k)(rho(i,j)(1./3.)+( rhoo)(1./3.)(1.0+3.0t2+4.50t2t2-1.50t1))
f(k,i,j)=1*feq(k,i,j)+(1-omega)*f1(k,i,j)
END DOEND DOEND DOIs rhoo is constat?
Thanks for your attention in advance.