mrayson / iwaves

Internal wave equation solvers (KdV, normal modes, TG)
Other
6 stars 7 forks source link

Bottom drag formulation #3

Open iosonobert opened 4 years ago

iosonobert commented 4 years ago

Issue in the bottom drag formulation of the kdvimex class. Current code has self.H^2 in the denominator, where self.H is the greatest depth input in the background density profile i.e. spatially constant and potentially arbitrary (kdv classes will extrapolate past this max depth id need be). In Holloway 97 the vertical scale (h) was a spatially variable function of Beta which is much smaller than the depth.