Open ChristinaB opened 6 years ago
Notes on fixing recharge (by Bert) I added a units conversion for recharge and upwelling that may be a fix for the problem. A conversion constant named 'units' is read in from line 7 in 'modelspc.dat' that converts from millimeters to meters. That constant is one of the variables passed to topmod(). There is a comment note on line 324 in topmod() about unit conversions: // Convert inputs to m per time step.
The constant in question is used to convert r1 net precipitation (throughfall + stemflow) on line 338 just before soil() is called. r1 in turn is one of the factors in r2 which is subsequently passed into soil(). Computations in soil() are thus in meters as are all of the computations in that section of topmod() code. I used the 'units' conversion constant to convert upwelling and recharge to millimeters, so now the task is to figure out why recharge for drainage 125 is so large, not why it is so small. But it may not be unreasonably large; mean annual precipitation for drainage 125 is 1432 mm, while mean recharge is 470.547 mm, which is 38.286%.
Weighting for irrigation should come from fraction irrigated in basinpars. Weighting for drainage type should come from 1=Tile Fraction, 2 = ditch fraction, Type 0 weight = 1-tilefraction-ditchfraction
Result = Rows Timestep; Columns Drainage total R2 and Rd