function[vmax,dmax]=ws_max_param(iunit,th,nt,dt)
%
    ymax=0.;
    vmax=0.;
    dmax=0.;
%
    area=0.; 
    dth=dt/2;
%
    vth=zeros(nt,1);
%
    for j=1:nt
%
        if( (abs(th(j))>ymax))
             ymax=(abs(th(j)));
        end
%
    end
%
    for i=1:nt
%
      if(i==0 || i==nt)
          area=area+th(i)*dth;
 
      else   
          area=area+th(i)*dt;
      end
 
      if( abs(area)>vmax)
          vmax=(abs(area));
      end    
%
      vth(i)=area;
%
    end
%
    area=0.;
%
    for i=1:nt
%        
      if(i==0 || i==nt)
          area=area+vth(i)*dth;
      else
          area=area+vth(i)*dt;
      end
%
      if( (abs(area))>dmax)
          dmax=(abs(area));
      end    
    end
%
    if(iunit==1)
        vmax=vmax*(386.);
        dmax=dmax*(386.);
    end
    if(iunit==2)
        vmax=vmax*(9.81);
        dmax=dmax*(9810.);
    end
    if(iunit==3)
        vmax=vmax*(1.);
        dmax=dmax*(1000.);
    end