real*8 a,b,c,reso
	integer*2 h,k,l,hmax,kmax,lmax
	integer*4 count
	write (6,*)'Enter cell edges for ortho cell:'
	read (5,*)a,b,c
	write (6,*)' Enter Dmin and Dmax (eof to quit):'
	write(6,*)' order of Dmin, Dmax irrelevant'
	write(6,*)' Caution- if D<1.0 it will be teated as S'
10	read(5,*,end=100) dmin,dmax
	if (dmin.lt.1.0) dmin=sqrt(1/dmin)
	if (dmax.lt.1.0) dmax=sqrt(1/dmax)
	if (dmax.lt.dmin) then
	x=dmax
	dmax=dmin
	dmin=x
	endif

	count=0
C*****This depends on hmax > a/dmin etc
	hmax=int(1.1*a/dmin)
	kmax=int(1.1*b/dmin)
	lmax=int(1.1*c/dmin)
c	write(6,*)'hkl limits: ',hmax,kmax,lmax
	
	Do 50 h=0,hmax
	Do 50 k=0,kmax
	Do 50 l=0,lmax
	Reso= 1./sqrt((h/a)**2+(k/b)**2 + (l/c)**2)  
C          (orthorhombic cells only)
c	write(6,*)h,k,l,dmin,reso,dmax,count
50 	if ((reso.lt.dmax).and.(reso.gt.dmin)) count=count+1
	write (6,*)'dmax, dmin, #refl:',dmax,dmin, count
	goto 10
100	continue
	end

c	169.59 182.52 240.57
c	3.0 12.0