function trm=vels(ffreqs,dirs,wns,z,depth,xpos,ypos)

trm=(-i*ffreqs)*ones(size(dirs));