%-----------------------------------------------------------
function model2d
global im jm km lm
im=64
jm=64
km=21
lm=21

nstep=10
h=1
m=1
nv=5
dnu=1e-13
deltat=h/m

sbinit(nv,dnu,deltat);
t=0
z=sbiniz(0);
sbgrph(t,z);

for istep=1:nstep
  [t,z]=tdrknu(m,h,t,z,'sbgdzl','sbgdzn');
  sbgrph(t,z);
end



