search

hamishHay / GeodesicODIS

Finite volume fluid dynamics code to simulate ocean flow in icy satellites and other bodies.
1 stars 1 forks source link

Interpolation of ocean height needs updating #4

Open hamishHay opened 1 month ago

hamishHay commented 1 month ago

Currently, ocean height is interpolated to a control volume edge using 0.5*(h1+h2), but this leads to poor performance if control volume edge center does not lie on the line directly between two nodes.

hamishHay commented 1 month ago

Will modify so that ocean height will be interpolated using a barycentric weighting following Eq. 6 in Piexoto (2016). This should be updated in momentumAdvection.cpp as well.