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

Allow momentum forcing to be given by quadratic drag from wind field #163

Closed edoddridge closed 6 years ago

edoddridge commented 6 years ago

Currently the momentum forcing is in N/m^2, and the momentum input does not depend on the flow field. It would be helpful in some situations to be able to specify a wind field in m/s and use a quadratic drag law to force the momentum equations.

I expect it should be possible to use the same wind input files, and just repurpose them internally. This means that the RelativeWind logical will change the units and meaning of the input wind files. It also means that the wind magnitude time series machinery can be used with this functionality.

The new parameters required are: