PC-Progress Discussion Forums

I suggest that you read the following two papers, where we discuss these issues. Note that some simulations from these papers are posted on the HYDRUS website. J. Ramos, T. B., J. Šimůnek, M. C. Gonçalves, J. C. Martins, A. Prazeres, N. L. Castanheira, and L. S. Pereira, Field evaluation of a multic...

No. HYDRUS does not keep track of that. Such option would need to be newly coded. You could get actual uptake from different subregions for a particular print time if you integrate the sink term in the Nod_inf.out file. J.

The parameter α (denoted as h0 in Kosugi [1996]) has units of [L] and should be entered as such in GUI. J.

The results look reasonable to me for the conditions that you described. I assume that the saturated water content is about 0.27 (the water content of the bottom node). Further I would assume that precipitation is smaller than Ks, and as a result, the profile will only reach the water content corres...

I do not think that it is possible to comment on your project without knowing more details, at the minimum the boundary conditions. J.

I do not see how it could be done in one run with the existing model. J.

The Beta function describing the root distribution with depth (in the GUI denoted as "Root [1/cm]") is a function that describes the proportional distribution of potential root water uptake (root distribution) over depth. It does not have any actual values, i.e., it is a relative function. The progr...

You need to use the modified van Genuchten model. This model then has additional parameters Qa (<Qr) and Qm (>Qs). Different values of the Qm parameter will give you different values of the air-entry value. The larger the difference between Qm and Qs, the larger the air-entry value. J.

I do not see how you could do all this in one run. You will likely have to split the simulation into multiple simulations for different time periods, while importing the final results of one simulation as an initial condition for the second simulations, etc, and changing the type of BCs. J.

Yes, it is the same thing. The l parameter is sometimes referred to as the "tortuosity and pore connectivity" parameter. Mualem recommended a value of 0.5, although other people recommended other values (e.g., Schaap et al. in his papers on Rosetta recommended a value of -1).

J.

PC-Progress could implement this function(s) into the computational module as a paid service. However, please note that having the analytical function for the retention curve is not sufficient to use HYDRUS to simulate water flow. You also need to have definitions of the capacity function (a derivat...

HYDRUS does not recognize any upper limit on the pressure head. It has upper and lower limits only for the Atmospheric BC. J.

I cannot see how this could happen. Can you post your project? J.

The code will always assigned as default cValue1 to all boundaries. You need to overwrite it with cValue2 and cValue3. Simply pick the Cauchy BC, overwrite the boundary, and you will be required to enter a code (pointer). That should do it. J.

Is you can choose any of the standard BCs for the bottom of your profile, why don't you extend the profile all the way to the groundwater level and then fix the pressure head there to keep GWT constant or use no flow to allow water to accumulate there? J.