PC-Progress Discussion Forums

It's irrelevant in such case. J.

here is the entire routine. Figure it out. I do not have time for that any more. J. ************************************************************************* subroutine Weight(HO,FO,WT,iType,iPos,NOBb,NOB,NOB1,NOB2,iWeight, ! Type,iConcType,ierr) double precision FO(NOBb),WT(NOBb) dimension HO(NOBb)...

The easiest way to explain is to give you the code (when the standard deviation is used): else if(iWeight.eq.2) then if(nQ.gt.0) QAve=QAve/nQ if(nH.gt.0) hAve=hAve/nh if(nTh.gt.0) ThAve=ThAve/nTh if(nV.gt.0) vAve=vAve/nV if(nC.gt.0) cAve=cAve/nC if(nK.gt.0) rKAve=rKAve/nK QSig2=0. hSig2=0. ThSig2=0....

cBot should be directly the concentration of groundwater, without any adjustments and conversions. You need to have consistent units for all input concentrations, i.e., initial and boundary concentrations. In US, we usually express EC in dS/m. J.

No problem. Good luck with the simulations. Jirka

You are using a dual-porosity model and specify the saturated water content of the immobile (Ths_im) and mobile (Ths_m) domains equal to 0.806 and 0.428, respectively. That is obviously possible since then the porosity would be larger than 1 (Ths_im + Ths_m = 1.334). That's the first problem. Second...

That cannot be true. If you have at the bottom Variable Pressure Head and Concentration BC, the program will honor the bottom pressure head as well as the bottom concentration. J.

That depends what bottom boundary condition you use. If you use the pressure head BC (i.e., you specify the position of groundwater table, then you should use the Dirichlet BC and groundwater EC for cBot. If you use a free drainage BC, then cBot is not used.

No problem. Also, I usually simply zip the project and post it; and I have not encountered any problems with that. J.

It seems that you do all the right things. Obviously, since HYDRUS-1D uses a local gradient optimization method (and not a global optimization method), you could do some more comprehensive analysis. See, for example, some recent papers by Brunetti: Brunetti, G., J. Šimůnek, and P. Piro, A comprehens...

This is not particularly useful. Can you just post the entire project (with the inverse file)? J.

HYDRUS-1D prints solute fluxes for up to 4 observation nodes at the end of the Solute.out file. See Table 13.3. SOLUTE.OUT in the Manual (the last item in the table). J.

It would be good to see the project when this happens. For a uniform flow model, this is pretty straightforward, but for a dual-porosity model it is more complicated. J.

It is a standard version. You can select this option ("Print Fluxes (instead of Temp) for Observation Nodes") on the Print Information dialog window. J/

Use the Cauchy (third-type) BC at both ends. At the outflow BC, the Cauchy BC will automatically be transformed into the Neumann BC (as is explain in the manual). J.