1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
| subroutine col(u,v,g,geq,th,theta,rho,
&omegat,t_in,tw,w,cx,cy,n,m)
real g(0:8,0:n,0:m),geq(0:8,0:n,0:m),th(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),theta(0:n,0:m)
real rho(0:n,0:m)
real t_in(m/2),tw
tw=0.0
t_in(m/2)=1
do i=0,n
do j=0,m
t1=u(i,j)*u(i,j)+v(i,j)*v(i,j)
g(0,i,j)= -(2.d0/3.d0)*th(i,j)*t1(0,i,j)*rho(i,j)
do k=1,4
t2=u(i,j)*cx(k)+v(i,j)*cy(k)
g(k,i,j)= w(k)* th(i,j)*rho(i,j)*(3.d0/2.d0
&+(3.d0/2.d0)*t1+(9.d0/2.d0)*t2*t2-(3.d0/2.d0)*t2)
end do
do k=5,8
t2=u(i,j)*cx(k)+v(i,j)*cy(k)
g(k,i,j)= w(k)* th(i,j)*rho(i,j)*(3+(2.d0/3.d0)*t1+
&(9.d0/2.d0)*t2*t2-(3.d0/2.d0)*t2)
end do
end do
end do
return
end |
Partager