dvm_module.f90

圆柱绕流问题的涡方法数值模拟（fortran源码）

module const
integer,parameter:: N=10,Nx=10,Ny=10,Mpv=2,Ns=500
real,parameter   :: pi=3.141592653589793,U=10.0,fmiu=0.0002
real,parameter   :: xl=0.1,yl=0.1
integer,parameter:: knp=5,pnum=50
integer,parameter:: Nall=2*Nx+2*Ny,NMAX=2*Ny+2*Nx+1,nlimit=Ns+1
end module

module varible
use const
integer,save::KK,nstep
integer,save::ndp(nlimit,knp)
integer,save::nt(knp)
real,save   ::fx(ns),fy(ns),bnp(nall)
real,save   ::btao(nall),Matrx1(nmax,nall),bmatrx(nall,nall)
real,save   ::Gx(nlimit,knp),Gy(nlimit,knp),GVx(nlimit,knp),GVy(nlimit,knp)
real,save   ::bv(nmax,2),bc(nall,2),thta(nall)
real,save   ::alpha,re,dt,ddt,PN,dL,dnt
real,save   ::tao(nlimit,knp),sigma(nlimit),sigma0(nlimit,Mpv)
!----------------------------------------------------
real,save   ::px(nlimit,pnum),py(nlimit,pnum),pz(nlimit,pnum),pvx(nlimit,pnum),pvy(nlimit,pnum),pvz(nlimit,pnum)
real,save   ::ti
integer,save::pt
real,save   ::fvx(2,pnum),fvy(2,pnum),dvdx(pnum),dvdy(pnum),dudx(pnum),dudy(pnum)
end module