%-----------------------------------------------------------
function sbinit(nv,dnu,deltat)
  global dl
  global im jm km lm  
  for k=-km:km
    for l=-lm:lm
      dl(1+lm+l,1+km+k)=exp(-dnu*deltat/2*(k*k+l*l)^nv);
    end
  end
