implicit none integer nx,ny,ic,jc parameter (nx=90,ny=90) ; integer m(nx,ny) character file*20,titolo*20 call cellularrule (-1,0) m=0 call mvsetflags('X size',600.0) do ic=18,72,9 ; do jc=18,72,9 call cellularinit(m,nx,ny,1,ic,jc) enddo ; enddo call cellularplot(m,nx,ny) call produciplot ('tmp000.gif','Cellular Automata package', 2 'Cycle n: 0') do ic=1,13 call cellularcycle(m,nx,ny) call cellularplot(m,nx,ny) call makefilename('tmp',ic,'.gif',file) write(titolo,'(a,i3)') 'Cycle n:',ic call produciplot (file,'Cellular Automata package',titolo) enddo call viacolvento('tmp*.gif','example381.gif') end