Initial ice sheets heights shouldn't be negative!

VirtualPlanetaryLaboratory / vplanet

The Virtual Planet Simulator

MIT License

136 stars 55 forks source link

Initial ice sheets heights shouldn't be negative! #150

Open RoryBarnes opened 3 years ago

RoryBarnes commented 3 years ago

The option dinitIceSheetHeight in POISE allows negative arguments, which shouldn't be allowed. It'd be better to change it so that if the user inputs a negative value, it defaults to meters since the ice sheet height is positive definite.

RoryBarnes commented 3 years ago

This issue is in the file src/poise.c.