Ocean Modeling with Fortran

Ocean Modeling with the Navier-Stokes Equations: or as I think of it, 'my missed oppurtunity.'

Tyler Tucker

Project outline

This work borrows heavily from Jochen Kaempf's book Advanced Ocean Modelling: Using Open-Source Software. Equations, terminology, and code are used to model physical ocean processes. Kaempf's numerical method of choice is the Successive Over-Relaxation (SOR) method. Although there are faster more efficient methods using Krylov Subspace methods, the SOR can be as an exclusive solver for non-symmetric matrices, where the latter methods may need some tweaking. Nevertheless, I am set out to see how much more efficient these methods can be.

Bayesian Predictive Posterior for Normal Distribution using SIR prior

Bayesian Predictive Posterior for a Normal Distribution using SIR Prior

Hello all,

My Bayesian statistics class had me compute the predictive posterior distribution of a normal distribution. You can do this by hand if you use the standard improper reference (SIR) prior. This was painfull to do on my own, and I thought it would be helpful to include this on my blog, in case anyone wants to compute this, and wants to see how I did it.

Taylor Series

Recursive Taylor Series Using Python's Sympy Library

Greetings everyone, I ran into some homework problems involving taking a multivariate Taylor expansion and noticed that Sympy's doesn't have this feature. They do have a series expansion at one point. However, things go sour for additional terms. I managed to write a multivariate Taylor series expanded about the origin up to the second order and provided the code below for those who are interested. If you see something that can be improved, let me know at Thanks!

Tibet Snow-Man


Snow and Ice cover has been collected in .asc files here by the National Snow and Ice Data Center (NSIDC). One product of interest is the Ice Mapping System (IMS) which takes satellite and field measurements and puts them on a grid. Currently, they have resolutions of about 24x24km 4x4km and 1x1 km. Their website is located here

IMS Polar Stereographic Map Projections

An explanation of IMS's map projection:

I feel the need to clarify how the National Ice Center's Interactive Multisensor Snow and Ice Mapping System (IMS) generates their maps.

But first a little background. The IMS project has been uploading zipped .asc files daily since 1997. From the start, they were using a 24x24 km grid. Since 2004, they have improved their measurement methods and process and implemented a 4x4 km resolution grid. Lastly, 2014 saw further improvements with a 1x1 km resolution map. Below is a picture was taken from their website.