%-----------------------------------------------------------
function z=sbgdzl(t,dt,z)
  global dl
  global im jm km lm
  
  %for k=-km:km
  %  for l=-lm:lm
  %    z(1+lm+l,1+km+k)=dl(1+lm+l,1+km+k)*z(1+lm+l,1+km+k);
  %  end
  %end
  z = dl.*z;
