dvm_parameters.f90

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

subroutine parameters()
use varible 
integer i,j,k
real tk,ts

entry para_sigma()
do k=1,Ns+1
 tk=float(k)*dt
 sigma(k)=2.24*sqrt(fmiu*tk)
 do m=1,Mpv
    ts=float((k-1)*Mpv+m)*ddt
    sigma0(k,m)=2.24*sqrt(fmiu*ts)
 end do
end do

return
end subroutine