subroutine me_slit_init

	implicit none

	include 'sieve.inc'

	integer i, j
	real*8 dxsum, dx, xcount
	real*8 dysum, dy, ycount,inch_cm

        inch_cm=2.54
	slit_info.x = slit_info.holes(1,1).x
	slit_info.y = slit_info.holes(1,1).y
	slit_info.num_horz = slit_info.num_holes / slit_info.num_vert
	dxsum  = 0
	xcount = 0

	do i = 1, slit_info.num_vert
	  if (slit_info.num_horz .gt. 1) then
	    do j = 2, slit_info.num_horz
	      dx = slit_info.holes(j,i).x - slit_info.holes(j-1,i).x
	      dxsum = dxsum + dx
	      xcount = xcount + 1
	    enddo
	  endif
	enddo

	if (xcount .ne. 0) then
	  slit_info.dx = dxsum / xcount
	else
	  slit_info.dx = 1
	endif

	dysum  = 0
	ycount = 0

	do i = 1, slit_info.num_horz
	  if (slit_info.num_vert .gt. 1) then
	    do j = 2, slit_info.num_vert
	      dy = slit_info.holes(i,j).y - slit_info.holes(i,j-1).y
	      dysum = dysum + dy
	      ycount = ycount + 1
	    enddo
	  endif
	enddo

	if (ycount .ne. 0) then
	  slit_info.dy = dysum / ycount
	else
	  slit_info.dy = 1
	endif

	write (*,*) ' Total number of holes ', slit_info.num_holes
	write (*,*) ' Number of vertical holes ', slit_info.num_vert
	write (*,*) ' Number of horizontal holes ', slit_info.num_horz
	write (*,*) ' average spacing between horizontal holes ',
     >              slit_info.dx
	write (*,*) ' average spacing between vertical holes ',
     >              slit_info.dy
	write (*,*) ' hole 1,1 location ', slit_info.x, slit_info.y

	end