search

edoddridge / aronnax

An idealised isopycnal model that can be run either with n+1/2 layers, or with n layers and variable bathymetry.
http://aronnax.readthedocs.io/en/latest/
MIT License
23 stars 5 forks source link

Bottom drag only applied if there is more than one layer #193

Closed edoddridge closed 5 years ago

edoddridge commented 6 years ago

Currently the bottom drag is applied in the same block of code as the linear drag between layers: https://github.com/edoddridge/aronnax/blob/61a1115383a9a69e55c6c98c7d956d87dfaa296b/src/momentum.f90#L67-L75

This code is within an if statement and is only evaluated if there is more than one layer. Therefore, in a 1 layer simulation the bottom drag will not be applied, even though it should be.