[plplot]интерактивное обновление графиков

call plinit()
call plcol0(3)
call plenv(xmin,xmax,ymin,ymax,0,0)
do ...
    <самрасчёт>
    plclear()
    call plline(x_,y_)
end do
call plend()

use plplot
implicit none
integer :: j
call plsdev ('xwin')
call plinit()
do j=1,100
	call demo1(j)
end do
call plend
end

      subroutine demo1(j)
      use plplot
      implicit none
      real(plflt) vmin, vmax
      real(plflt), dimension(0:100):: x,y
      integer i,j
      call plbop
	do i=0,100
	x(i)=i
	end do
      call plssub(1,2)

        call plcol0(15)
        call pladv(0)
        vmin = 0.1_plflt
        vmax = 0.9_plflt
          call plvpor(vmin,vmax,vmin,vmax)
          call plwind(0.0_plflt, 100.0_plflt, -1.1_plflt, +1.1_plflt)
          call plbox('bctn', 0.0_plflt, 0, 'bctn', 0.0_plflt, 0)
         call plcol(1)
         y=sin((x+j)*0.25)
        call plline( x, y )
        
        call pladv(0)
        call plcol0(15)
        vmin = 0.1_plflt
        vmax = 0.9_plflt
          call plvpor(vmin,vmax,vmin,vmax)
          call plwind(0.0_plflt, 100.0_plflt, -1.1_plflt, +1.1_plflt)
          call plbox('bctn', 0.0_plflt, 0, 'bctn', 0.0_plflt, 0)
         call plcol(5)
         y=sin((x-j)*0.25)
		call plline( x, y )	
	call sleep(1)
	end