我正在尝试在Fortran OpenGl中编写一个简单的代码来动画线性聚合物的动力学。在正方形晶格上,最初直线型聚合物是连续的左上和右上键合的水平序列。同样,这是一系列的上下扭结。偏置动态是向下移动的序列 - >>扭结,向下扭结是随机选择的,具有统一的概率。为了制作动画,我需要擦除向下扭结的两个键,然后绘制向上扭结的两个键。 “擦除”键是通过使用背景颜色“绘制”其上的另一个键来执行的。问题是当我运行我的代码一段时间的蒙特卡洛时间时,它给出了模拟结束时聚合物最终配置的正确图片,但没有显示扭结翻转的中间微观移动序列。我想,我必须缺少一个或多个图形命令。如果有人试图提供帮助,我可以发送代码。最后,我在Linux(Ubuntu)环境中工作。非常感谢提前。
部分代码:
!----用于删除下行链接
call glColor3f(1.0,1.0,0.0) ! Back ground colour: Yellow
call drawline((i-1)*f,s(im1,j)*f,i*f,s(i,j)*f)
call drawline(i*f,s(i,j)*f,(i+1)*f,s(ip1,j)*f)
s(i,j)=s(i,j)+2
!----用于绘制向上链接
call glColor3f(1.0,0.0,0.0) ! Polymer colour: Red
if(i.ne.1)call drawline((i-1)*f,s(im1,j)*f,i*f,s(i,j)*f)
if(i.ne.Lx)call drawline(i*f,s(i,j)*f,(i+1)*f,s(ip1,j)*f)