ST_example

 Example usage of Sediment Toolbox

 function example

 DESCRIPTION:
    Demonstrate ST Toolbox functionality 

 INPUT:

 OUTPUT:
    
 Author(s):  
    Geoff Cowles (University of Massachusetts Dartmouth)

 Revision history
   
==============================================================================

0001 function ST_example
0002 % Example usage of Sediment Toolbox
0003 %
0004 % function example
0005 %
0006 % DESCRIPTION:
0007 %    Demonstrate ST Toolbox functionality
0008 %
0009 % INPUT:
0010 %
0011 % OUTPUT:
0012 %
0013 % Author(s):
0014 %    Geoff Cowles (University of Massachusetts Dartmouth)
0015 %
0016 % Revision history
0017 %
0018 %==============================================================================
0019 
0020 global ftbverbose
0021 ftbverbose = false;
0022 
0023 close all
0024 fprintf('       phi          class       d(mm)    Dstar    wset(mm/s)  taucr (Pa)  erate x1e-3(kg/(m^2-s))\n')
0025 i = 0;
0026 for phi=-8:11
0027     i = i + 1;
0028     phiclass    = ST_wentworth(phi);
0029     d(i)        = ST_phi2d(phi);
0030     Dstar(i)    = ST_Dstar(d(i));
0031     Wset(i)     = ST_wset(d(i));
0032     Taucr(i)    = ST_taucr(d(i));
0033     erate(i)    = ST_erate(d(i));
0034     fprintf('%10d %20s %8.4f %8.2f %9.4f %8.3f %8.3f\n',phi,phiclass,d(i)*1000,Dstar(i),Wset(i)*1000,Taucr(i),1000*erate(i))
0035 end;
0036 
0037 loglog(d*1000,Taucr)
0038 title('critical shear stress')
0039 xlabel('Grain diameter (mm)')
0040 ylabel('critical shear Pa')
0041 axis([.01,10,.01,10])