function s=read_m3d(flnm)
fp=fopen(flnm,'r');
s.l1=fgetl(fp);
s.l2=fgetl(fp);
s.l3=fgetl(fp);
s.l4=fgetl(fp);
s.l5=fgetl(fp);
k=0;
a=fscanf(fp,'%f',[8,inf])';
whos
    s.yrday=a(:,1);
    s.x=a(:,2);
    s.y=a(:,3);
    s.u=a(:,4);
    s.v=a(:,5);
    s.zmin=a(:,6);
    s.zmax=a(:,7);
    s.zbottom=a(:,8);
fclose(fp);