Program to derive harmonic constants from water level observations