- #1

MelihAltunan

- 8

- 0

subroutine collesion(u,v,f,feq,rho,omega,w,cx,cy,n,m,tm,tminv,stmiv)

real f(0:8,0:n,0:m)

real feq(0:8,0:n,0:m),rho(0:n,0:m)

real w(0:8), cx(0:8),cy(0:8)

real u(0:n,0:m), v(0:n,0:m)

real tm(0:8,0:8),tminv(0:8,0:8),stmiv(0:8,0:8)

real fmom(0:8,0:n,0:m),fmeq(0:8,0:n,0:m)

!calculate equilibrium moments

do i=0,n

do j=0,m

fmeq(0,i,j)=rho(i,j)

fmeq(1,i,j)=rho(i,j)*(-2.0+3.0*rho(i,j)*(u(i,j)*u(i,j)+v(i,j)*v(i,j)))

fmeq(2,i,j)=rho(i,j)*(1.0-3.0*rho(i,j)*(u(i,j)*u(i,j)+v(i,j)*v(i,j)))

fmeq(3,i,j)=rho(i,j)*u(i,j)

fmeq(4,i,j)=-rho(i,j)*u(i,j)

fmeq(5,i,j)=rho(i,j)*v(i,j)

fmeq(6,i,j)=-rho(i,j)*v(i,j)

fmeq(7,i,j)=rho(i,j)*(u(i,j)*u(i,j)-v(i,j)*v(i,j))

fmeq(8,i,j)=rho(i,j)*u(i,j)*v(i,j)

end do

end do

the error is Warning: Rank mismatch in argument 'rho' at (1) (rank-2 and scalar)

i guess it's because fmeq and rho have different ranks but how i can write in a proper way?

Thank you for now.