Flex2D.m
Solves the 2D flexural equation (e.g. Turcotte and Schubert) using a centered-finite difference technique that allows for spatially variable loads and flexural rigidity (or elastic thickness). Nestor Cardozo has a cool Mac program that does something similar with a nicer interface, but it's useful to dig into code occasionally.
Solves the 2D flexural equation (e.g. Turcotte and Schubert) using a centered-finite difference technique that allows for spatially variable loads and flexural rigidity (or elastic thickness). Nestor Cardozo has a cool Mac program that does something similar with a nicer interface, but it's useful to dig into code occasionally.
subsidence.m
Calculates and plots sediment accumulation (compacted/uncompacted) and subsidence (total/tectonic) curves from stratigraphic data. with uncertainties based on paleobathymetry. Included is a sample dataset.
Calculates and plots sediment accumulation (compacted/uncompacted) and subsidence (total/tectonic) curves from stratigraphic data. with uncertainties based on paleobathymetry. Included is a sample dataset.
strain_lattice_D.m
Finds the best-fit lattice strain model (Blundy and Wood, 1994; 2003) to a given set of mineral/melt or mineral/bulk rock (e.g. zircon/whole rock) REE partition coefficients. Returns the theoretical strain compensated (i.e. "strain free") partition coefficient (D) and Young's Modulus (E) for the mineral. Included is a sample dataset to experiment with. Useful for examining the internal consistency among a set of partition coefficients.
Finds the best-fit lattice strain model (Blundy and Wood, 1994; 2003) to a given set of mineral/melt or mineral/bulk rock (e.g. zircon/whole rock) REE partition coefficients. Returns the theoretical strain compensated (i.e. "strain free") partition coefficient (D) and Young's Modulus (E) for the mineral. Included is a sample dataset to experiment with. Useful for examining the internal consistency among a set of partition coefficients.
dem_contact.m
Finds the intersection line between a topographic surface and a planar feature using a DEM and user defined strike/dip. Especially useful for mapping out low-angle contacts or structures. I commonly use this with the companion script load_DEM.m to format a .tiff raster DEM exported from ArcMap.
Finds the intersection line between a topographic surface and a planar feature using a DEM and user defined strike/dip. Especially useful for mapping out low-angle contacts or structures. I commonly use this with the companion script load_DEM.m to format a .tiff raster DEM exported from ArcMap.
threept.m
Simple three-point problem solver. Returns the strike/dip of a plane with user entered x,y,z points (UTMs). I use this all the time in conjunction with Google Earth where you can read coordinates right off the screen.
Simple three-point problem solver. Returns the strike/dip of a plane with user entered x,y,z points (UTMs). I use this all the time in conjunction with Google Earth where you can read coordinates right off the screen.
points2plane.m
A 3-point problem on steroids . Loads a text file with unlimited x,y,z data points and inverts to find the best-fit plane. Good for geomorphic surfaces as well as structural features.
A 3-point problem on steroids . Loads a text file with unlimited x,y,z data points and inverts to find the best-fit plane. Good for geomorphic surfaces as well as structural features.
disloc2D.m
A basic elastic edge dislocation model in 2D that returns surface displacement with user defined fault geometry and magnitude of slip. There are much better and more rigorous codes and software out there, but this is fun to play around with.
A basic elastic edge dislocation model in 2D that returns surface displacement with user defined fault geometry and magnitude of slip. There are much better and more rigorous codes and software out there, but this is fun to play around with.
iso_mix.m
A simple GUI for playing around with isotopic mixing, this one is set up for 87Sr/86Sr vs. εNd. Again, there are lots of programs that do this already, but I like having a basic stand-alone app to play with. Included is a sample data file if you want to try fitting mixing curves to data.
A simple GUI for playing around with isotopic mixing, this one is set up for 87Sr/86Sr vs. εNd. Again, there are lots of programs that do this already, but I like having a basic stand-alone app to play with. Included is a sample data file if you want to try fitting mixing curves to data.