HP1 History

History of the HP1 Code Development - The List of Changes and Fixed Errors

HP1 Version 2.3.001 (to be released)

To get a pre-release of Version 2.3.001, please contact djacques@sckcen.be.

New Processes, Models and Options

1. Diffusion and advection in the gas phase
2. Update physical transport properties depending on geochemical state variables by user-defined relations: scaling factors of water content, pressure head and hydraulic conductivity; tortuosity in the water and gas phase, dispersivity for solute transport, heat capacity, thermal conductivity and thermal dispersivity
3. Include model of Moldrup for calculating tortuosity coefficients

Input Features and Changes:

1. New data block keyword DATABLOCK to define vectors and matrices of user-defined variables and parameters. Information can be used in BASIC-statements.
2. New data block keyword REACTIVE_TRANSPORT to define different options for update of physical properties

Output Features and Changes:

1. Print overview of aqueous and gaseous components and of the chosen update options in phreeqc.out

New Identifiers and BASIC Statements:

1. Identifier –bound_layer in data block GAS_PHASE to indicate the thickness of the boundary layer for upper boundary condition
2. Identifier –state in data block CALCULATE_VALUES to indicate the specific use of the BASIC function
3. BASIC statement change_sfcond(number) to change the scaling factor of the hydraulic conductivity for a given node number.
4. BASIC statement change_sfpressure(number) to change the scaling factor of the pressure head for a given node number.
5. BASIC statement change_sfwater(number) to change the scaling factor of the water content for a given node number
6. BASIC statement conv_t("unit") to convert time units to the unit defined in the GUI.
7. BASIC statement conv_x("unit") to convert length units to the unit defined in the GUI.
8. BASIC statement get_def(…) to get the value of a defined value or a value from the data block DATABLOCK
9. BASIC statement get_sfcond(number) to get the scaling factor of the hydraulic conductivity for a given node number
10. BASIC statement get_sfpressure(number) to get the scaling factor of the pressure head for a given node number
11. BASIC statement get_sfwater(number) to get the scaling factor of the water content for a given node nmber
12. BASIC statement hyd_cond(number) to get the hydraulic conductivity for a given node number
13. BASIC statement pressure_head(number) to get the pressure head for a given node number
14. BASIC statement put_def(…) to put the value of a defined variable
15. BASIC statement sum_solid(…) to obtain the sum of all solid phases multiplied by a phase-specific value
16. BASIC statement water_flux(number) to get the water flux for a given node number
17. BASIC statement water_sink(number) to get the root water uptake for a given node number

HP1 Version 2.2.002 - September 2009

Input Features and Changes:

1. Possibility to define initial and boundary solutions using solution types or compositions.
2. Full integration in HYDRUS-1D Graphical User Interface; a complete project can be defined and executed from the HYDRUS-1D Graphical User Interface.
3. No need anymore to give components, which are not present in the initial solution, a small concentration of, e.g., 1E-24 mol/kgw.

Output Features and Changes:

1. A warning when a component is present in PHREEQC, but not transported by HYDRUS was added.
2. Selected-output: changes in the amount of minerals: recalculate to units in mol/1000cm3 soil/time.
3. Possibility to create dump files.
4. phreeqc.out: Indication of time step levels when printing during transport calculations.
5. Printing of templates for plotting profile data and time series of variables defined in the SELECTED_OUTPUT. These plots can be generated using the free software GNUPLOT.
6. An output file createdfiles.out, listing all output files created by HP1.

New Identifiers and BASIC Statements:

1. Identifier –gnuplot for use in TRANSPORT data block. If –gnuplot is followed by true, GNUPLOT templates are generated. Default –gnuplot false
2. BASIC statement bulkdensity(number): returns the bulk density, which is defined in the Solute Transport Parameters window for node number.

Fixed Errors:

1. Transferring the components Total_O, Total_H and Charge in immobile water cells used Concentration instead of Im_Cocentration.
2. No memory was allocated for cell_data in case of MIM.

HP1 Version 2.1.001 - May 2008

Source Codes:

1. Includes HYDRUS-1D 4.0
2. Includes PHREEQC-2-15.0-2697

HP1 Version 2.0 - January 2008

Source Codes:

1. Includes HYDRUS-1D 3.003 (in double precision)
2. Includes PHREEQC-2-14.3-2411

New Processes:

1. Surface complexation and oxidation/reduction processes were added by defining the components Total_O, Total_H, and Charge

Output Features and Changes:

1. Add identifiers to TRANSPORT data block to control saving of output in additional files: -punch_hydrus (selected output variables are saved for observation nodes and print times defined by the variables ), -separate_tsfiles (save output for each observation node in a separate output file), -separate_mim (save output for mobile and immobile phase in a separate output file)

HP1 Version 1.0 - November 2004

Source Codes:

1. Includes HYDRUS-1D version 2.0
2. Includes PHREEQC-2.4

Note on Version Numbering: HP1 Versions A.B.C