fortran - 尝试使用 do 循环和打开文件在 Fortran 中制作 gif

Question

我目前正在尝试使用 Fortran 制作 gif，但我仍然是初学者，我不知道该怎么做。我的程序是关于一个受干扰的钟摆的相位图，我想用它制作一个 gif。这是我的代码：

integer :: i,j,ij
real(kind=8) :: x,xp,xpp,a,dt,eps,t,pii
character*12 :: filename
dt=0.01
pii=3.14159265359
eps=1
x=modulo((-pii+2*pii)/2+pii,2*pii)-pii
xp=0
t=0
do i=1,300
    x=modulo((-pii+2*i*pii)/(2*i)+pii,2*pii)-pii
    xp=0
    t=0
    write(filename,'("graphespoincare",i2.2,"txt")')i
    ij=30+i
    open(unit=ij,file=filename,form='unformatted',status='replace',action='readwrite')

    do j=1,300         !Euler
        a=x
        xpp=acc(a,t)   !acceleration function
        xp=xp+xpp*dt
        x=x+xp*dt
        t=t+dt
        write(ij,*) x,xp,t,en(x,xp)
    enddo
    close(ij)
enddo

现在的问题是当我使用它时，我没有得到任何文件。我究竟做错了什么？

标签： fortran