	subroutine init_array(u,nx,ny,nv,xfill)
	dimension u(nx,ny,nv)
	do 10 ix=1,nx
	do 10 iy=1,ny
	do 10 iv=1,nv
10	u(ix,iy,iv)=xfill
	return
	end
