function write_m3d(flnm,s)
fp=fopen(flnm,'w');
fprintf(fp,'%s\n',s.l1);
fprintf(fp,'%s\n',s.l2);
fprintf(fp,'%s\n',s.l3);
fprintf(fp,'%s\n',s.l4);
fprintf(fp,'%s\n',s.l5);
for k=1:length(s.yrday)
  fprintf(fp,'%f %f %f %f %f %f %f %f\n',s.yrday(k),s.x(k),s.y(k),s.u(k),s.v(k),s.zmin(k),s.zmax(k),s.zbottom(k));
  
end
