Fix zs computation while using DEM with nodata

Having nodata into DEM leads to introduce false information when compute zs.

1. define valid range for zs value #65 (closed)
2. fix return value, switching "-1" to "-1000" #66 (closed)
3. add information into documentation