Sediment Transport Tools

MS698_03

Sediment Transport Processes in Coastal Environments

Instructor: Courtney K. Harris
Spring Semester, 2003

Tools:

This is intended to be a storage spot for useful matlab files (m-files, *.mat files, figures), fortran files, web-links, or other data files. Anyone associated with the class can feel free to provide things for this page- just email them to ckharris@vims.edu, and we'll try to put them up here. Requirement for contributed files is that the author is easily accessible in case someone discovers a bug.

email from february 27:

HI
I've set up a new page on the class web-site; http://www.vims.edu/~ckharris/MS698_03/tools.html

This can serve as a repository of helpful tools (matlab, fortran, data, etc.) for doing sediment transport calculations. Nothing will be put up on this page unless I have some confidence that it is working correctly, but there are no guarantees in the world of community software.

Whenever you use shared code you should verify that it is giving reasonable answers in time to have any problems corrected (say well before handing in a homework assignment).

Your contributions are welcomed to this page; just email me the matlab or fortran code that you want to contribute and I'll put it on the page if I agree that it is a good contribution. (Art- do you still have that settling velocity code?).

You can also make suggestions as to other tools that would be handy.

You should be able to use the sed_normplot.m routine to do problem 1 of the first homework set. If you prefer to not do the work in matlab; let me know and I'll scrounge up some probability paper for you to hand plot the sediment distribution. I tried to improve it this morning... and think that it is plotting correctly.

Thanks,
Courtney

M-files

Matlab m-files that have been useful for doing sediment-transport calculations. Maybe some of them will help you with your problem set or research. Never accept any software on blind faith, they could contain bugs. Let me (or the author of the script) know if you find any problems with them.

phitod.m and dtophi.m: simple matlab functions for going back and forth from phi sizes to d(mm) sizes.
clc_z0_nik.m: curve fit to the relationship z0/ks = f(R*) used to relate hydraulic roughness to physical roughness based on dependency on roughness Reynolds number.
sed_normplot.m: matlab plotting function that takes a sediment size distribution (phi_coarse, frac) and plots a normal-probability plot. Vectors phi_coarse and frac should contain the coarse-end of size bins, and the fraction within each bin (0-1). What you get is h_plot=sed_normplot(phi,frac) and a plot. The handle, h_plot lets you modify the colors of the lines.

Data Sets

FORTRAN code

Wiberg, P.L. and C.K. Harris, 1994. Ripple geometry in wave-dominated environments. Journal of Geophysical Research, 99 (C1), 775-789. FORTRAN subroutine hlpred.f available for download.