Calculate settling velocity of particle diameter d (m) in m/s
function [wset] = ST_wset(d,varargin)
DESCRIPTION:
Calculate settling velocity of particle diameter d (m) in m/s
INPUT:
d: sediment grain size in m
[optional] 'temperature' = temperature of the seawater in C [default=10]
[optional] 'salinity' = salinity of seawater in PSU [default=35]
[optional] 'sdens' = sediment density in kg/m^3 [default=2650]
OUTPUT:
wset: settling velocity in m/s
EXAMPLE USAGE
wset = ST_wset(.0005,'temperature',10,'salinity',35,'sdens',2650)
Author(s):
Geoff Cowles (University of Massachusetts Dartmouth)
References
Soulsby DMS (102)
Revision history
==============================================================================0001 function [wset] = ST_wset(d,varargin) 0002 % Calculate settling velocity of particle diameter d (m) in m/s 0003 % 0004 % function [wset] = ST_wset(d,varargin) 0005 % 0006 % DESCRIPTION: 0007 % Calculate settling velocity of particle diameter d (m) in m/s 0008 % 0009 % INPUT: 0010 % d: sediment grain size in m 0011 % [optional] 'temperature' = temperature of the seawater in C [default=10] 0012 % [optional] 'salinity' = salinity of seawater in PSU [default=35] 0013 % [optional] 'sdens' = sediment density in kg/m^3 [default=2650] 0014 % 0015 % OUTPUT: 0016 % wset: settling velocity in m/s 0017 % 0018 % EXAMPLE USAGE 0019 % wset = ST_wset(.0005,'temperature',10,'salinity',35,'sdens',2650) 0020 % 0021 % Author(s): 0022 % Geoff Cowles (University of Massachusetts Dartmouth) 0023 % 0024 % References 0025 % Soulsby DMS (102) 0026 % 0027 % Revision history 0028 % 0029 %============================================================================== 0030 0031 subname = 'ST_wset'; 0032 %fprintf('\n') 0033 %fprintf(['begin : ' subname '\n']) 0034 0035 % constants 0036 grav = 9.8106; %g 0037 T = 10; %T (C) 0038 S = 35; %S (PSU) 0039 sdens = 2650; %sediment density in kg/m^3 0040 0041 % parse arguments 0042 for i=1:2:length(varargin)-1 0043 keyword = lower(varargin{i}); 0044 if( ~ischar(keyword) ) 0045 error('incorrect usage of ST_wset') 0046 end; 0047 0048 switch(keyword(1:3)) 0049 0050 case 'tem' 0051 T = varargin{i+1}; 0052 case 'sal' 0053 S = varargin{i+1}; 0054 case 'sde' 0055 sdens = varargin{i+1}; 0056 otherwise 0057 error(['Can''t understand value for:' keyword]); 0058 end; %switch keyword 0059 end; 0060 0061 0062 % calculate nu 0063 nu = SW_Kviscosity(T,S); 0064 0065 % calculate rho 0066 dens = SW_Density(T,S); 0067 0068 s = sdens; 0069 dstar = ([grav*(s-1)/(nu^2)])^(1/3)*d; 0070 0071 0072 % calculate wset 0073 wset = (nu/d)*( sqrt(10.36^2 + 1.049*(dstar^3)) - 10.36); 0074 0075 0076 %fprintf(['end : ' subname '\n'])