Search found 21 matches

by mvcallaghan
Tue Jan 09, 2007 10:18 pm
Forum: STANMOD
Topic: Use of 3DADE for Simple GW Transport
Replies: 1
Views: 5296

I have since found many examples of the use of 3DADE to horizontal GW transport. Including its implementation in the SOLUTRANS program. So I guess that I've answered my own question...yes, 3DADE can be used to simulate 1D horizontal groundwater transport with 3D dispersion.

Mike
by mvcallaghan
Tue Jan 09, 2007 9:13 pm
Forum: STANMOD
Topic: Use of 3DADE for Simple GW Transport
Replies: 1
Views: 5296

Use of 3DADE for Simple GW Transport

Hi All, I have been testing 3DADE for simple 1D groundwater transport of chloride. The way that the source dimensions are set out in the program are for vertical transport in a one dimensional flow field. Is there any conceptual problem with rotating the frame of reference so that the flow is modele...
by mvcallaghan
Wed Nov 29, 2006 7:44 pm
Forum: Hydrus 1D
Topic: high intensity rainfall events
Replies: 3
Views: 4515

I'd like to reopen the discussion on this topic. Is it possible that the Hydrus code could be modified to apportion excess water at the atmospheric boundary to runoff when the soil profile becomes saturated? I have many cases where I am examining shallow water tables 0.5 to 1 m depth (fixed pressure...
by mvcallaghan
Mon Sep 18, 2006 7:42 pm
Forum: Hydrus 1D
Topic: multiple runs
Replies: 5
Views: 6903

Hi Bruno, I do alot of batch runs with Hydrus. I needed to create three files. The first was a batch file, which I simply called run.bat. The run.bat file contains the following text: copy Path1.txt level_01.dir H1d_calc<return.txt copy Path2.txt level_01.dir H1d_calc<return.txt ...etc. ,where Path1...
by mvcallaghan
Thu Aug 17, 2006 8:58 pm
Forum: Hydrus 1D
Topic: How to eliminate crashes in batch runs
Replies: 3
Views: 3803

Thanks Jirka, Enjoy your vacation! It would be be a great aid to my work if you were able to reprogram Hydrus to terminate on non-convergence. I'll contact you by e-mail at the beginning of September for an update. I was also wondering if it would be an easy task to increase the allowable number of ...
by mvcallaghan
Wed Aug 09, 2006 4:41 am
Forum: Hydrus 1D
Topic: more than 1 variable boundry in the simulation
Replies: 1
Views: 2330

A1. What kind of variable boundary condition do you want on the bottom boundary? I have sucessfully used a variable pressure head boundary condition to represent a fluctuating water table. It that case it is important to make sure that the function used to define the boundary pressure head is smooth...
by mvcallaghan
Mon Jul 31, 2006 9:31 pm
Forum: Hydrus 1D
Topic: How to eliminate crashes in batch runs
Replies: 3
Views: 3803

How to eliminate crashes in batch runs

Hi Everyone, Does anyone have experience programming batch runs of Hydrus 1D so that the simulation continues past a model locked in a non-convergence loop? My problem is that everytime the model enters non-convergence in a run, the batch run gets stuck in a loop that requires manual stoppage. This ...
by mvcallaghan
Wed May 17, 2006 5:24 pm
Forum: Hydrus 1D
Topic: convergence problem
Replies: 1
Views: 2459

There are so many variations in model setup, that really all I can do is offer broad suggestions. The top two problems for convergence that I have encountered are the 1) the runoff BC, and 2) high intensity rainfall. Both seem to be related to saturating the upper BC. The model does not seem to like...
by mvcallaghan
Wed May 17, 2006 5:13 pm
Forum: Hydrus 1D
Topic: UNSATCHEM Mass Units
Replies: 4
Views: 4646

It's been awhile since your response, but I wanted to add that your advice was very useful. Many thanks. FYI for others: The use of activities in the "Gapon" type equation (which normally uses concentrations), has been referred to as being used to calculate the Vanselow coefficient. The paper "Sodiu...
by mvcallaghan
Thu Apr 20, 2006 12:09 am
Forum: Hydrus 1D
Topic: rainfall
Replies: 17
Views: 17232

It could be one of many problems. So I can only give some suggestions at this point. 1) make sure your minimum time step is small enough e.g. 1E-5 days 2) try a node spacing as small as < 1cm 3) turn off interpolation tables by setting upper and lower limit of tension interval to zero (in the Iterat...
by mvcallaghan
Wed Apr 05, 2006 4:54 pm
Forum: Hydrus 1D
Topic: UNSATCHEM Mass Units
Replies: 4
Views: 4646

Many thanks for the quick reply, it was very helpful. My checks on the mass balance for solute concentrations check out fine with the addition of the factor of 100. For adsorbed species, the mass units in CHEMBAL.OUT don't match up with my hand calculations. It is best illustrated with the following...
by mvcallaghan
Wed Apr 05, 2006 3:17 pm
Forum: Hydrus 1D
Topic: Vadose Zone modelling
Replies: 2
Views: 3278

When you say "the model continues to run", what do you mean by that? Are you getting output to the screen in the DOS window? If you are getting output, is the model converging, or are you getting a convergence error message? Or is the model stopping at 1 sec. of runtime? More details would be helpful.
by mvcallaghan
Tue Apr 04, 2006 5:33 pm
Forum: Hydrus 1D
Topic: UNSATCHEM Mass Units
Replies: 4
Views: 4646

UNSATCHEM Mass Units

Jirka, I have a few questions about the mass units used in the UNSATCHEM module of HYDRUS-1D. The mass and length units I have specified in my simulation are mmol and cm. In the UNSATCHEM manual the solute concentration units are listed as mmol/l and for adsorbed species as mmol/kg (in Table 12.8) e...
by mvcallaghan
Mon Jan 30, 2006 5:11 pm
Forum: Hydrus 1D
Topic: rainfall
Replies: 17
Views: 17232

Likely, your input files are buried too deep in the file structure. There is a limit to the length of the file pathway. Beyond that limit, Hydrus can't read the filename.

I would suggest putting your input files into a directory right off your C: root. For example C:\Hydrus\....

Mike
by mvcallaghan
Fri Jan 27, 2006 7:51 pm
Forum: Hydrus 1D
Topic: rainfall
Replies: 17
Views: 17232

I've seen the non-convergence problem before. Are you getting through a partial run before a crash or is it crashing at the start. 1. Problem could be with your time step, node spacing or iteration criteria. Try a time step of 1E-5 days, node spacing of 10 mm, and max iterations of say 200. Set the ...