fini
/cle
/batch,list
*CFOPEN,'efi-rpa6-b0','dat'
!*dim,frac,array,5,5
NMAX=3000
N0=25
/input,hydra6b-bm,inp

!--          Voltages
v_g = 0.0   ! Ground Potential
v_c = -12.0   ! Collector Potential

pi=acos(-1.)
qq = 1.602e-19
mm = 15.9994*1.672e-027


*do,i1,50,50
	v_r = 0.1*i1
	/input,hydra6b-av,inp
	*do,k1,0,5,5
 		alpha = k1*pi/180.
		*do,l1,0,1
			*if,k1,eq,0,and,l1,eq,1,cycle
		beta = l1*pi*45/180.
 		xl0 = xl - g_loc1*tan(alpha)*cos(beta)
	 	yl0 = yl - g_loc1*tan(alpha)*sin(beta)
	 	xu0 = xu - g_loc1*tan(alpha)*cos(beta)
	 	yu0 = yu - g_loc1*tan(alpha)*sin(beta)
		reg=1
		*do,j1,0,5*(i1+1)
			*dim,x,array,NMAX
			*dim,y,array,NMAX
			*vfill,x,rand,xl0,xu0
			*vfill,y,rand,yl0,yu0
			energy = 0.02*(j1+0.5)
			/input,hydra6b-ic,inp
			*set,x
			*set,y
		*enddo
		reg=2
		*do,j1,1,5
			*dim,x,array,NMAX
			*dim,y,array,NMAX
			*vfill,x,rand,xl0,xu0
			*vfill,y,rand,yl0,yu0
			energy = 0.1*(i1+1+j1)
			/input,hydra6b-ic,inp
			*set,x
			*set,y
		*enddo
		*do,j1,1,20
			*dim,x,array,NMAX
			*dim,y,array,NMAX
			*vfill,x,rand,xl0,xu0
			*vfill,y,rand,yl0,yu0
			energy = 0.1*(i1+6)+0.5*(j1)
			/input,hydra6b-ic,inp
			*set,x
			*set,y
			*if,energy,gt,10.,exit
		*enddo
		*enddo !-l1
	*enddo !-k1
*enddo !-i1
*CFCLOS