	function resolv(e)
c     return resolution at energy e
c
	include '/u1/gcf/e350/crystal/cb/clouds'
c
	resolv=sqrt(sigma*e**2+sigmb*e)
	return
c
	end
