SPINS Science Files

SPINS comes with a variety of Science files which help with simulation data analysis. Below is a list of currently available functions with a brief description of their purpose and the assumptions behind their implementation. These Science files are in spins/src/Science.

bottom_slope.cpp

Returns the vector of the bottom slope at each x-grid point. This is used for calculating the stresses along the bottom boundary (cf. bottom_stress_x.cpp and bottom_stress_y.cpp).

Assumes the bottom is at the 0-index for $z$.

bottom_stress_x.cpp

Returns the x-component of stress along the bottom boundary.

If the case is unmapped: $t_{x}=\mu u_{z}$ where $\mu$ is the dynamic viscosity.

If the case is mapped: $t_{x}={\frac {\mu }{\alpha ^{2}(x)}}\left(2h'(x)(w_{z}-u_{x})+(u_{z}+w_{x})(1-h'^{2}(x))\right)$ where $\alpha (x)=sqrt{1+h'^{2}(x)}$ and where $h'^{2}(x)$ is the bottom slope.