vplanet
vplanet copied to clipboard
Initial ice sheets heights shouldn't be negative!
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.
This issue is in the file src/poise.c.