!!!!!!!!!!!变量定义!!!!!!!!!!!!!
module param
implicit none
double precision,parameter::pi=3.14159265358979
double precision,parameter::dt=0.0001
integer,parameter::dim=2
integer,parameter::max=1600
integer:: i,j,d,k,i$mestep,n%le, niac, scale_k, maxn,nvirt,ntotal
integer::pair_i(max),pair_j(max),scal_k
double precision::vx(dim,4000),x(dim,4000),v_min(dim,4000),mass(4000)
double precision::dvxdt(dim,4000),indvxdt(dim,400),exdvxdt(dim,4000),av(dim,4000)
double precision::hsml, eta,rho,factor
double precision::dx(dim),dvx(dim),dr,q,tw,tdwdr
double precision::dwdx(dim,max),w(max),r(max),dwdr(max)
end
program poiseuille
use param
implicit none
character(len=100)::%lename
rho=1000.0
eta=0.001
mass=2.5e-5*2.5e-5*rho
hsml=2.5e-5*1.1000
scal_k=2
factor=15.0/(7.0*PI*hsml*hsml)
n%le=0
!!!!!!!!!生成真实粒子!!!!!!!!!!!!
k=0
do i=1,40
do j=1,3
k=k+1
x(1,k)=2.5e-5*(j-1)+1.25e-5
x(2,k)=2.5e-5*(i-1)+1.25e-5
vx(1,k)=0
vx(2,k)=0
enddo
enddo
ntotal = 120
!!!!!!!!!!!生成真实粒子完毕!!!!!!