New Features and Developments in
HYDRUS Software Packages
Jirka Šimůnek
Department of Environmental Sciences
University of California Riverside, CA, USA

Miroslav Šejna
PC-Progress, s.r.o., Anglicka 28, Prague 120 00, Czech Republic

Rien van Genuchten
U.S. Salinity Laboratory, USDA, ARS, Riverside, CA 92507, USA

OUTLINEHYDRUS (2D/3D)- new features- new GUI
- ParSWMS
 HYDRUS-1D- new features
 HP1
 HYDRUS Package for MODFLOW

Hydrus workshop, Prague, March 28 2008

HYDRUS - History
Israel: Neuman [1972] - UNSAT
U. of Arizona: Davis and Neuman [1983]

2006 – Beta Version
2007 – Version 1.0
2008 – Version 1.05
Šimůimůnek,nek, J., M. Th. van Genuchten, and M. Šejna,
ejna, TheHYDRUS Software Package for Simulating TwoTwo- and ThreeThreeDimensional Movement of Water, Heat, and Multiple Solutes in
VariablyVariably-Saturated Media, Technical Manual, Version 1.0, PC
Progress, Prague, Czech Republic, pp. 241, 2006.
http://www.pchttp://www.pc-progress.com
http://www.hydrus2d.com

Agr. U. in Wageningen:
Feddes et al. [1978]
Vogel [1987] - SWMII

Princeton U.:
van Genuchten [1978]
MIT:
Celia et al. [1990]

USSL - SWMS-2D
Šimůnek et al. [1992]

USSL - CHAIN-2D
Šimůnek and van Genuchten [1994]

IGWMC - HYDRUS-2D (1.0)Šimůnek et al. [1996]

IGWMC - HYDRUS-2D (2.0)Šimůnek et al. [1999]

USSL – SWMS_3D
Šimůnek et al. [1995]

PC-Progress – HYDRUS (2/3D)Šimůnek et al. [2007]1

Spatial Root Distribution Function(Vrugt et al., 2001)

HYDRUS (2D/3D) – New features
TwoTwo- and threethree-dimensional water flow and solute transport
Water flow and solute transport in a dualdual-porosity system Root water uptake with compensation
 Spatial root distribution functions of Vrugt et al. (2001)+
 Soil hydraulic property models of Kosugi (1996) and Durner(1994)+
 Hysteresis model of Lenhard et al. (1991) (no pumping effect)+ Transport of viruses, colloids, and/or bacteria (attachment/
detachment model, filtration theory, and blocking functions) Dynamic, systemsystem-dependent boundary conditions+ Flowing particles in twotwo-dimensional applications+ Calculations of actual and cumulative fluxes across internal
meshlines+ Constructed wetland module (only in 2D)+
 New print management options



⎛ p

k
θ ( h) - θ r
1= ∑ win m
θ s -θ ri =1(1+ α i h )eθr , θs - residual and saturated water contents, respectivelyk- number of overlapping subregionswi- weighting factors for the sub-curves
αi, ni, mi (= 1 - 1/ni), and l - empirical parameters of the sub-curves
The hydraulic characteristics contain 4+2k unknown parameters: θr , θs , αi , ni , l, and Ks.
Of these, θr, θs, and Ks have a clear physical meaning, whereas αi, ni and l are essentiallyempirical parameters determining the shape of the retention and hydraulic conductivity
functions [van Genuchten, 1980].

2

Log(Conductiv ity [cm/days])

Se

⎞

0
Total

0.5

m
⎛
⎞
wiα i ⎡1- (1- S 1/ei mi ) i ⎤ ⎟
∑
⎜
k
⎣
⎦
⎛
⎞
⎠
K (θ ) = K s ⎜ ∑ wi S eli ⎟ ⎝ i =12k
⎝ i =1⎠⎛⎞⎜ ∑ wiα i ⎟⎝ i =1⎠- effective water content

*

⎛ p
⎞
p
− ⎜ z z * − z + r x* − x ⎟
Xm
⎝ Zm
⎠

0.6

i

k

p

*

Multi-Porosity Hydraulic Property Model
Durner (1994)

Wate r Conte nt [-]

i

*

⎛
z ⎞⎛
x ⎞
b ( x, z ) = ⎜ 1 −⎟ ⎜1 −⎟Xm ⎠⎝ Zm ⎠ ⎝Multi-Porosity Hydraulic Property Model
Durner (1994)
S e (h) =

py

x
z
⎛
x ⎞⎛
y ⎞⎛
z ⎞ −⎜⎜⎝ X m x − x + Ym y − y + Z m z − z ⎟⎟⎠−−b ( x, y , z ) = ⎜ 1 −
e
1
1
⎟⎜
⎟⎜
⎟
⎝ X m ⎠⎝ Ym ⎠⎝ Z m ⎠MatrixFracture

0.4
0.3
0.2
0.1
0

-2
Total
Matrix-4

Fracture
-6

-8

-10
-1

0

1

2

3

Log(|Pressure Head [cm]|)

4

5

-1

0

1

2

3

4

5

Log(|Pressure Head [cm]|)

Example of composite retention (left) and hydraulic conductivity(right) functions (θr=0.00, θs=0.50, α1=0.01 cm-1, n1=1.50, l=0.5,
Ks=1 cm d-1, w1=0.975, w2=0.025, α2=1.00 cm-1, n2=5.00).

2

Dynamic, System-Dependent
Boundary Conditions

Hysteresis Model (Lenhard et al., 1991)(pumping effect)(no pumping effect)

Old: only Atmospheric and Seepage Face boundary conditions
Interpolate variable pressure head and flux boundary conditions
smoothly with time Switch the boundary condition from variable pressure head to zero
zero
flux (e.g., disc permeameter)
permeameter) Switch the boundary conditions from timetime-variable pressure head tozero flux when the specified nodal pressure head is negative (e.g.,(e.g.,
above the water table) As above, except that an atmospheric boundary condition is assigned
assigned
to nodes with negative calculated pressure heads
 As above, except that a seepage face boundary condition is assigned
assigned
to nodes with negative calculated pressure heads
 Treat the timetime-variable flux boundary conditions similarly as
atmospheric, i.e., with limiting pressure heads (hCritS(hCritS and hCritA)
hCritA) Apply atmospheric boundary conditions on nonnon-active seepage face Snow accumulation on top of the soil surface when temperatures arearenegative

Kool and Parker (1987)

Se ( h) = [1 + |α w h| ]− mnLenhard et al. (1991)n −m

Se ( h) = [1 + |α d h| ]α w ≈ 2α d

Flowing Particles in TwoDimensional Applications

Actual and Cumulative Fluxes
Across Internal Meshlines

3

Constructed Wetlands Module(CW2D: Langergraber and Simunek, 2005)
Components

Processes

dissolved oxygen [mg O2/l], readily and slowly biodegradable,
and inert OM [mg COD/l],
 ammonia, nitrite, nitrate, and N2[mg N/l] inorganic phosphorus [mg P/l], heterotrophic MO [mg COD/l],
 autotrophic MO (Nitrosomonas
and Nitrobacter) [mg COD/l], Organic N and P are modelled as
N and P content of the COD.












Constructed Wetlands Module
1. Hydrolysis - conversion of slowly biodegradable organic matter CS into
readily biodegradable organic matter CR, with a small fraction being converted into
inert organic matter CI. Ammonium (NH4+) and inorganic phosphorus (IP) arereleased.

hydrolysis (slowly Î readily
biodegradable)
aerobic growth of heterotrophic
MO (degradation of OM)
anoxic growth of heterotrophic
MO using nitrate and nitrite(denitrification)
aerobic growth of Nitrosomonas
and Nitrobacter (nitrification)lysis (for heterotrophic andautotrophic MO)

rc1 = Kh

2. Aerobic growth of heterotrophic bacteria - consumes oxygen (O2)and readily biodegradable organic matter (CR), while ammonium (NH4+) andinorganic phosphorus (IP) are incorporated in the biomass

rc2 = μH

cO2
cCR
f N,Het cXH
KHet,O2 + cO2 KHet,CR + cCR

3. NO3-based growth of heterotrophs on readily biodegradable COD(denitrification) - consumes nitrate (NO3-) and readily biodegradable organic

matter (CR). Nitrate is reduced to dinitrogen (N2). Again, ammonium (NH4+) andinorganic phosphorus (IP) are incorporated in the biomass

rc3 = μDN

Two-stage subsurface verticalflow constructed wetlandcCS cXH
cXH
KX + cCS cXH

KDN ,O2
KDN ,NO2
cNO3
cCR
f N ,DN cXH
KDN ,O2 + cO2 KDN ,NO3 + cNO3 KDN ,NO2 + cNO2 KDN ,CR + cCR

HYDRUS (2D/3D) – New GUI
Completely new GUI based on Hi-End 3D graphics libraries
MDI architecture – multiple projects and multiple views
 New organization of geometric objects
 Navigator window with an object explorer
 Many new functions improving the user-friendliness, such as
drag-and-drop and context sensitive pop-up menus Improved interactive tools for graphical input
 Saving Cross-Sections and Mesh-Lines for charts within agiven project Display Options – all colors, line styles, fonts and otherparameters of graphical objects can be customized
 Extended print options
 Extended information in the Project Manager (including
project preview)



0

50

100

150

200

250

Simulated steady-state distribution of heterotrophic organisms XHMany additional improvements

4

ProjectManager

Navigator Bar
The Navigator Bar is by default located on the left side of the
HYDRUS main window. The Navigator Bar has three Tabs:• A Data Tab to allow quick access to all input and output
data.
Input data include:

Navigator
Bars

Domain Geometry
Flow Parameters
FE-MeshDomain Properties
Initial and Boundary Conditions
Auxiliary Objects

Output data include various Results.
Data are organized in a tree-like structure.
• A View Tab to specify what and how information will be
displayed in the View window, and• A Sections Tab to show various Sections

5

Edit
Bars

Edit BarThe Edit Bar is by default located on the rightside of the HYDRUS main window.
The Edit Bar is very dynamic since it changes
depending upon the process being carried out.

Toolbars

Edit
Bars


Standard Toolbar

Tools Toolbar

View Toolbar



GUI Toolbar

Time Layer Toolbar

6

Menus
File

Edit

View

InsertMenus
Calculation Results Tools Options Windows

View Window
Capillary Barrier
Material Distributions

Tabs: Geometry, FE-Mesh, Domain Properties, Initial
Conditions, Boundary Conditions, Results

7Capillary Barrier
Velocity Vectors

Cut-off Wall
Finite Element MeshCut-off Wall
Solute Plume

Plume Movement in
a Transect with Stream

8

HYDRUS (2D/3D)
129,461 FE
I415
411 m

61 m

1

2

3

4

5

Finite element mesh and material distribution in the transport domain

9

HYDRUS (2D/3D)

HYDRUS (2D/3D)

-4
0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

5 years

10 years

15 years

25 years

-3

-2

-1

0

1

2

3

4

5

6

log (Conc [mg/L])

Velocity vectors

Contours of MTBE concentrations at 5, 10, 15, and 25 years

Flow and Transport Under the
Banana Tree (Sansoulet et al., 2008)

ParSWMSParSWMS (Hardelauf et al., 2007) - Parallelizedversion of SWMS_3D (Simunek et al., 1995).
1995).
 Developed by Forschungszentrum Jülich,lich, Germany.
 SWMS_3D – earlier and simpler version of HydrusHydrus-3D
 MPI (Message(Message-Passing Interface) - a libraryspecification for messagemessage-passing between the different
processors. MPI is free software for LINUX or UNIX
operating systems.
 Test- Water flow and solute transport problem - 492264 FE nodes- Supercomputer with 41 SMP nodes with 32 processors each (total
1312 processors - Power4+ 1.7 GHz)

10ParSWMS

1998 – Version 2.0 (USSL)
2005 – Version 3.0 (UCR)
2008 – Version 4.0 (UCR+PC Progress)

Time gain as compared to the one processor run (in log2) as a function of thenumber of processors np (in log2) for solute transport scenario with 492264 nodes(open circles) and water flow with atmospheric upper boundary conditions(diamonds) (Hardelauf et al., 2007).
2007)

Šimů
imůnek,nek, J., M. Šejna,
ejna, H. Saito, M. Sakai, and M. Th. van
Genuchten, The HYDRUSHYDRUS-1D Software Package for Simulating
the Movement of Water, Heat, and Multiple Solutes in Variably
Saturated Media, Version 4.0, HYDRUS Software Series 3,3,
Department of Environmental Sciences, University of California
Riverside, Riverside, California, USA, pp. 315, 2008.
http://www.pchttp://www.pc-progress.com
http://www.hydrus2d.com

Coupled Water, Vapor, andEnergy Transport

HYDRUS -1D – New features
Coupled water, vapor, and energy transport DualDual-porosity water flow and solute transport(two(two-site sorption in the mobile zone)
 DualDual-permeability type water flow and solutetransport (two(two-site sorption in both zones) Potential evapotranspiration calculated using
PenmanPenman-Monteith combination or Hargreaves
equations Daily variations in the evaporation, transpiration,and precipitation rates
 Support for the HP1 codeModified Richards’
Richards’ equation for water flow:flow:∂θ ∂ ⎡
∂h
∂T
∂h
∂T ⎤
= ⎢ KLh (h) + KLh (h) + KLT (h) + Kvh + KvT ⎥ − S∂t ∂z ⎣
∂z
∂z
∂z
∂z ⎦
∂θ ∂ ⎡∂h∂T ⎤∂ ⎡ ∂h∂T ⎤=( KLh + Kvh ) + KLh + ( KLT + KvT ) ⎥ − S = ⎢ KTh + KLh + KTT ⎥ − S∂t ∂z ⎢⎣∂z∂z ⎦∂z ⎣ ∂z
∂z ⎦

KLh
KLT
Kvh
KvT

- hydraulic conductivity for liquid phase fluxes due to gradient of h [LT-1]- hydraulic conductivity for liquid phase fluxes due to gradient of T [L2T-1K-1]- isothermal vapor hydraulic conductivity [LT-1]- thermal vapor hydraulic conductivity [L2K-1T-1]

11

Coupled Water, Vapor, andEnergy Transport

Coupled Water, Vapor, andEnergy Transport

Energy Transport:
10

∂CpT

∂θ
∂⎡
∂T ⎤
∂qT
∂q
∂q T
+ L0 v = ⎢λ (θ ) ⎥ − Cw l − CwST − L0 v − Cv v∂t∂t ∂z ⎣
∂z ⎦
∂z
∂z
∂z
(1)
(2)
(3)
(4)
(5)

10

10

8

T=0

T=0

t=0.25 dt=0.25 dt=0.25 d

8

t=1

t=1

D e p th [cm ]t=5(1) Soil heat flow by conduction(2) Convection of sensible heat by water flow(3) Heat removed by root water uptake(4) Transfer of latent heat by diffusion of watervapor(5) Transfer of sensible heat by diffusion of water
vapor

t=25

6

10

T=0

6

8

T=0
t=0.25 d
8

t=1

t=5

t=5

t=25

t=25

6

t=5

4

4

4

2

2

2

2

0

0
0.05

0.1

0.15

Water Content [-]

0.2

0

0

0.02

0.04

Total Flux [cm/d]

t=25

6

4

0

t=1

0

0.06 01020Temperature [C]

30

0

2

4

6

8

10

Concentration [-]Total flux=water flux+vapor flux

Physical Nonequilibrium Water Flowand Solute Transport Models
a)

b)
Water
Soluteθc)
Water

Immob.

SoluteImmob.

Mobileθ = θim + θ mo

d)

Water

MobileSoluteImmob.

Mobileθ = θim + θ mo

Water

Slow

e)
Slow

Fast

Im. Slow

Solute
Slow

Water

Fast

Physical Nonequilibrium Water Flowand Solute Transport Modelsa)b) + c)

d)

e)

Fast

Soluteθ = θm + θ fFastθ = θ m ,im + θ m ,mo + θ fa) Uniform Flowb)Mobile-Immobile Waterc) Dual-Porosityd)Dual-Permeabilitye) Dual-Permeability with MIM in the Matrix Domaina) Uniform Flowd)Dual-Permeabilityb)Mobile-Immobile Water e) Dual-Permeability withc) Dual-Porosity
MIM

12Chemical Nonequilibrium Solute
Transport Models
a)

b)

c)

se

s1k

d)

θ

θ

sk

The PenmanPenman-Monteith combination method for calculating ofpotential evapotranspiration [FAO, 1990]:e)
Slow Fast

Immob. Mob.

θ

Penman-Monteith Equation

θim

θmo

cim

cmo

smoe

sme

smok

smk

θm

θf

cm

cf

sfe

ET0 = ETrad + ETaero =

sime
c

sk

c

s2k

c

s fk

a) One-Site Kinetic Modelb) Two-Site Model (kinetic and instantaneous sorption)c) Two-Site Kinetic Modeld) Dual-Porosity with One Kinetic Site Modele) Dual-permeability with Two-Site Model

Penman-Monteith Equation

ρ c p (ea - ed ) / ra ⎤1 ⎡ Δ( Rn - G )+⎢⎥λ ⎣ Δ + γ (1 + rc / ra ) Δ + γ (1 + rc / ra ) ⎦- reference crop evapotranspiration[mm d-1]
ETrad - radiation term [mm d-1]
ETaero - aerodynamic term [mm d-1]
Rn
- net radiation at crop surface[MJ m-2d-1]
G- soil heat flux [MJ m-2d-1](ea-ed) - vapour pressure deficit [kPa]Δ- slope vapour pressure curve[kPa oC-1]y- psychrometric constant [kPa oC-1]
λ- latent heat of vaporization [MJ kg-1]
ETo

ρ
cp
ea
ed
rc
ra

- atmospheric density [kg m-3]
- specific heat of moist air[i.e., 1.013 kJ kg-1 oC-1]- saturation vapor pressure attemperature T [kPa],- actual vapor pressure [kPa]- crop canopy resistance [s m-1]
- aerodynamic resistance [s m-1]

Penman-Monteith Equation

13

Daily Variations in Evaporation, andTranspiration Rates

Hargreaves Equation
ETp = 0.0023Ra (Tm + 17.8 ) TR

Tp (t ) = 0.24Tp

Ra - extraterrestrial radiation in the same units as ETp [e.g., mm d-1 or J m-2s-1]
Tm - daily mean air temperature [oC]
TR - temperature range between the mean daily maximum and minimum [oC]t < 0.264d, t > 0.736d⎛ 2π t π ⎞− ⎟ t ∈ (0.264d, 0.736d)
Tp (t ) = 2.75Tp sin ⎜⎝ 1day 2 ⎠1.2

Hourly values between 0-6
a.m. and 18-24 p.m.
represent 1% of the total
daily value and a sinusoidal
shape is followed during
the rest of the day(Fayer, 2000)

0.8

0.6

0.4

0.2

0
1/1/1999

3/2/1999

5/1/1999

6/30/1999

8/29/1999

10/28/1999

2.5

E, T [cm/d]

Evapotranspiration [cm/d]

3

Hargreaves
Penman-Montheith
Measured Bare Lysimeter
Measured Vegetated Lysimeter

1

2
1.5
1
0.5
0
0

0.2

0.4

0.6

0.8

1

Time [d]12/27/1999

Time [d]HP1 (Jacques and Šimůnek, 2005)

HYDRUS-1D

Coupled HYDRUS-1D and PHREEQC

Cumulative Fluxes

0
-500

60

-1000

40

-1500

20

-2000

0

potTopHYDRUS-1D [ŠŠimů
imůnek et al., 2008]:

potRoot
actTopactRoot

actBot




-2500

-20

-3000
0.0

0.1

0.2

0.3

0.4

0.5

50

Water Contents [-]All Fluxes



-40
100

150

200

250

300

Time [days]

Observation Nodes

PHREEQC [Parkhurst and Appelo, 1999]:
Available chemical reactions:

1.0
200



0.5
100



-0.5

0



-1.0

-100



-1.5

-200

0.0

-2.0



-300

50

100

150

200

Time [days]

250

300

50

100

Variably Saturated Water Flow
Solute Transport
Heat transport
Root water uptake

150

200

Time [days]

250

300




Aqueous complexation
Redox reactions
Ion exchange (Gains-Thomas)
Surface complexation – diffuse double-layer model and nonelectrostatic surface complexation model
Precipitation/dissolution
Chemical kinetics
Biological reactions

14

HP1 Examples

HYDRUS-1D GUI for HP1

Transport of heavy metals (Zn2+, Pb2+, and Cd2+) subject tomultiple cation exchange Transport with mineral dissolution of amorphous SiO2 andgibbsite (Al(OH)3)
 Heavy metal transport in a medium with a pHpH-dependent
cation exchange complex
 Infiltration of a hyperalkaline solution in a clay sample (thisexample considers kinetic precipitationprecipitation-dissolution of
kaolinite,kaolinite, illite,
illite, quartz, calcite, dolomite, gypsum, hydrotalcite,hydrotalcite,and sepiolite)
sepiolite) LongLong-term transient flow and transport of major cations (Na+,K+, Ca2+, and Mg2+) and heavy metals (Cd2+, Zn2+, and Pb2+) ina soil profile.
 Kinetic biodegradation of NTA (biomass, cobalt)

HP1 Publications
Jacques, D., and J. Šimůimůnek,nek, User Manual of the Multicomponent VariablyVariablySaturated Flow and Transport Model HP1, Description, Verification
Verification and
Examples, Version 1.0, SCK•SCK•CENCEN-BLGBLG-998, Waste and Disposal,
SCK•
SCK•CEN, Mol, Belgium, 79 pp., 2005.
 Jacques, D., J. Šimůimůnek,nek, D. Mallants,
Mallants, and M. Th. van Genuchten, OperatorOperatorsplitting errors in coupled reactive transport codes for transienttransient variably
saturated flow and contaminant transport in layered soil profiles,
profiles, J.
Contam.
Contam. Hydrology,
Hydrology, 88,
88, 197197-218, 2006.
 Šimů
imůnek,nek, J., D. Jacques, M. Th. van Genuchten, and D. Mallants,
Mallants,
Multicomponent geochemical transport modeling using the HYDRUS
computer software packages, J. Am. Water Resour.
Resour. Assoc.,
Assoc., 42(6),42(6), 15371537-1547,
2006.
 Jacques, D., J. Šimůimůnek,nek, D. Mallants,
Mallants, and M. Th. van Genuchten, Modelingcoupled hydrological and chemical processes in the vadose zone: A casestudy on long term uranium migration following mineral phosphorus
phosphorus
fertilization, Vadose Zone Journal,
Journal, Special Issue “Vadose Zone Modeling”
Modeling”,
2008a (in press).
 Jacques, D., J. Šimůimůnek,nek, D. Mallants and M. Th. van Genuchten, Modelling
coupled water flow, solute transport and geochemical reactions affection
affection
heavy metal migration in a Podzol soil, Geoderma,
Geoderma, 2008b (in press).


The HYDRUS (Unsaturated Flow)
Package for MODFLOW-2000
Hyeyoung Sophia Seo,
Seo, Jirka Šimůnek,nek, Eileen P. Poeter
Colorado School of Mines, 1500 Illinois St. Golden CO 80401
University of California, Riverside
Seo, H. S., J. Šimůnek, and E. P. Poeter, Documentation of the HYDRUS Package
for MODFLOW-2000, the U.S. Geological Survey Modular Ground-Water Model,
GWMI 2007-01, International Ground Water Modeling Center, Colorado School ofMines, Golden, Colorado, 96 pp., 2007.
Twarakavi, N. K. C., J. Šimůnek, and H. S. Seo, Evaluating interactions between
groundwater and vadose zone using HYDRUS-based flow package for MODFLOW,
Vadose Zone Journal, Special Issue “Vadose Zone Modeling”, (submitted April 25
2007, revised September 16 2007, accepted October 5 2007).

15

HYDRUS Package for MODFLOWHYDRUS Package for MODFLOW(Seo et al., 2006)

ZSURF

UNSF Soil Profile1

a: Ground Surface
b: Bottom of Soil Column

Layer 1
Columns
2

3

4

5

6

7

8

9

10 11 12 13

Explanation

Layer 2

K: Hydraulic
Conductivity

Flux (q)Explanation

Depth toGroundWater

1

MODFLOW Grid2

4

Zone 2

5
6

Z1

Zone 1

3

Rows

K1

7

b

Layer 3Layer 1
Layer 2

Z1

a

Depth

K1
K2
K3
K4
K5
K6
K7
K8
K9
K10
K11
K12

b

Layer 3

ZSURF

Explanation

Depth

a
Depth toGroundWater

Flux (q)a: Ground Surface
b: Bottom of Soil Column

K3
2

K: Hydraulic Conductivity

K: Hydraulic
Conductivity

UNSF Soil Profile

K2
1

Layers

K1
K2
K3
K4
K5
K6
K7
K8
K9
K10
K11
K12

MODFLOW Grid

3

HYDRUS - Web Site - Discussion ForumHYDRUS - Web Site – Tutorials

16

HYDRUS – Useful Links
http://www.pchttp://www.pc-progress.cz/Fr_Hydrus3D.htm
 http://www.pchttp://www.pc-progress.cz/Fr_Hydrus3D_Ordering.htm
 http://www.pchttp://www.pc-progress.cz/Fr_Hydrus3D_Downloads.htm
 http://www.pchttp://www.pc-progress.cz/Fr_Hydrus3D_Tutorials.htm
 http://www.pchttp://www.pc-progress.cz/Fr_Hydrus3D_Resellers.htm
 http://www.pchttp://www.pcprogress.cz/Fr_Hydrus3D_Troubleshooting.htm
 http://www.pchttp://www.pc-progress.cz/Pg_Hydrus_References.htm
 http://www.pchttp://www.pc-progress.cz/Pg_Hydrus1D_References.htm
Introduction to Soil Physics
with HYDRUS:
Modeling and Applications
David Radcliffe and Jirka Šimůnek
CRC Press, Taylor & Francis Group
ISBN-10: 142007380X, ISBN-13: 9781420073805
due 4/15/200917