********************************************************************

	real function dhitrost(e)
**-	real function dhitrost(e)

**- izracun potovalne hitrosti za elektrone v metanu
**- vhodni podatek je reducirano polje v v/cm/torr
**- izhodni podatek je v cm/s	
**- podatki so vzeti iz ?
	real e,el,xf,yf,k1,k2,c1,c2,dv,vdat(29)
	integer i

	data vdat /5.6,8.3,11.1,13.8,16.4,19,21.3,23.4,24.8,26,26.8
     &		  ,27,27,26.8,26,25.2,24.1,22.9,21.8,20.3,18.9,17.5
     &		  ,16.3,15.8,16.9,18.9,21.5,24,26.5/
	data xf,yf,k1,c1,k2,c2 /26.67,26.9,1.125,5.6,0.97,-41.4/

	el=log10(abs(e))
	el=xf*(el+1)
	if (el.le.0) then
	  dv=k1*(el)+c1
	 elseif (el.ge.70) then
	  dv=k2*(el)+c2
	 else
	  i=int(el/2.5)+1
	  el=el-2.5*(i-1)
	  dv=(vdat(i+1)-vdat(i))*el/2.5+vdat(i)
	endif
	dhitrost=10.**(6.+dv/yf)
	return
	end

