I hope you all are doing well. I am simulating a 100cmx190cm a sloping lysimeter in Hydrus 2D. The domain is subject to a long time series of rainfall and for that I am using its upper boundary as Atmospheric boundary condition. The right side is a seepage face. The figure

*fig_domain*(attached) depicts my domain.

I have a long sequence of rainfall and discharge for which I have observed values and will try to calibrate the internal parameters of flow. I am trying to obtain the discharge for each time interval (L/T) by getting the difference on the cumulative seepage and dividing by the difference in time for which each cumulative discharge was calculated. I also want to plot the actual flow into the domain as a way of verifying if the rainfall rates I provided throught the ATMOSPH.IN file are correct.

My problem arises when I plot the results for a 10 day experiment with varying precip rates: After a certain point my precipitation rate becomes coarser, so does the discharge rates. The attached figure

*input_output*(attached) depicts what I am saying.

I realized the Cum_Q.OUT file writes the cumulative flows with 5 digits only and at a certain point my cumulative fluxes become too big to be able to be written with enough precision. The same thing happens when I use mm, cm or m as my unit.

Is there any way I can increase the number of digits of these output files? Any other suggestion would be highly appreciated.