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 global ftbverbose 0032 [~, subname] = fileparts(mfilename('fullpath')); 0033 if ftbverbose 0034 fprintf('\nbegin : %s\n', subname) 0035 end 0036 0037 % constants 0038 grav = 9.8106; %g 0039 T = 10; %T (C) 0040 S = 35; %S (PSU) 0041 sdens = 2650; %sediment density in kg/m^3 0042 0043 % parse arguments 0044 for i=1:2:length(varargin)-1 0045 keyword = lower(varargin{i}); 0046 if( ~ischar(keyword) ) 0047 error('incorrect usage of ST_wset') 0048 end; 0049 0050 switch(keyword(1:3)) 0051 0052 case 'tem' 0053 T = varargin{i+1}; 0054 case 'sal' 0055 S = varargin{i+1}; 0056 case 'sde' 0057 sdens = varargin{i+1}; 0058 otherwise 0059 error(['Can''t understand value for:' keyword]); 0060 end; %switch keyword 0061 end; 0062 0063 0064 % calculate nu 0065 nu = SW_Kviscosity(T,S); 0066 0067 % calculate rho 0068 dens = SW_Density(T,S); 0069 %dens/1000. 0070 0071 % calculate dstar 0072 dstar = ST_Dstar(d,'temp',T,'sal',S,'sdens',sdens); 0073 0074 % calculate wset 0075 if ismatrix(d) 0076 wset = (nu./d).*( sqrt(10.36^2 + 1.049*(dstar.^3)) - 10.36); 0077 else 0078 wset = (nu/d)*( sqrt(10.36^2 + 1.049*(dstar^3)) - 10.36); 0079 end 0080 0081 if ftbverbose 0082 fprintf('end : %s\n', subname) 0083 end