Index: wflow-py/wflow/wflow_hbv.py =================================================================== diff -u -r8e380ffbfaedb9df9fccee3c03ebb53e346d1a71 -rc7476868c506252d099ab57c9909e317b43aeeee --- wflow-py/wflow/wflow_hbv.py (.../wflow_hbv.py) (revision 8e380ffbfaedb9df9fccee3c03ebb53e346d1a71) +++ wflow-py/wflow/wflow_hbv.py (.../wflow_hbv.py) (revision c7476868c506252d099ab57c9909e317b43aeeee) @@ -480,6 +480,9 @@ self.Slope= slope(self.Altitude) self.Slope=ifthen(boolean(self.TopoId),max(0.001,self.Slope*celllength()/self.reallength)) Terrain_angle=scalar(atan(self.Slope)) + temp = catchmenttotal(cover(1.0), self.TopoLdd) * self.reallength * 0.001 * 0.001 * self.reallength + self.QMMConvUp = cover(self.timestepsecs * 0.001)/temp + # Multiply parameters with a factor (for calibration etc) -P option in command line for k, v in multpars.iteritems(): @@ -865,7 +868,9 @@ Runoff=self.SurfaceRunoff - + self.QCatchmentMM = self.SurfaceRunoff * self.QMMConvUp + self.RunoffCoeff = self.QCatchmentMM/catchmenttotal(self.Precipitation, self.TopoLdd)/catchmenttotal(cover(1.0), self.TopoLdd) + self.sumprecip=self.sumprecip + self.Precipitation #accumulated rainfall for water balance self.sumevap=self.sumevap + ActEvap #accumulated evaporation for water balance self.sumpotevap=self.sumpotevap + self.PotEvaporation