program demo_spitpgm
  use spitpgm

  integer      :: fancydatas(128, 128)
  integer      :: i, j

  fancydatas = 0

  do i=1, 128
    do j=1, 128
    fancydatas(i,j) = mod(i*j, 65000)
    enddo
  enddo

  call spit_as_pgm_16(fancydatas, 'fancydatas.pgm')

end program