
subroutine init
sumv=0
sumv2=0
do i=1,npart
x(i)=lattice-pos(i)
v(i)=(ranf( )-05)
sumv=sumv+v(i)
sumv2=sumv2+v(i)2
end do
sum = sumv/npart
fs = sqrt(3temp/sumv2) ! 速度矫正因子
do i=1,npart
v(i)=(v(i)-sumv)fs ! 速度矫正
x(i)=x(i)-v(i)dt
end do
return
end
从上面的代码看,公式还是很简单的,你可以对照书上的公式一一对照。
使用编译器将源代码编译成可执行程序。首先,你得确认你已经正确安装了任意一个FORTRAN编译器;然后,要保证源代码正确有效;第三,建立项目,开始编译,如果一切顺利就可以得到可执行文件。
以上就是关于用Fortran编程实现分子动力学模拟方法,研究微流体流动特性全部的内容,包括:用Fortran编程实现分子动力学模拟方法,研究微流体流动特性、如何编译fortran语言、等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)