HISTORY OF HST PROGRAMS (Status as of 20 Oct 2010)

Version 2.2.16 is released:

Patch levels 14 and 15 have never been released. Starting at patch level 16, a separate distribution package for Unix platforms is no longer provided. Corrected the software to handle leading zeros in the 9-digit boundary condition codes. Enabled skipping of columns of dry cells for all the flux types of boundary conditions. Improved some output messages and unit labels. Fixed a deallocation bug for boundary conditions when the values change for only some of the simulation periods. Added checking for presence of new boundary-condition data for subsequent simulation periods. Improved checking of well-screen location data. Corrected some Fortran 90 usage errors including argument intent declarations, argument target declarations, and explicit interfaces. Nested some concatenated if tests to avoid subscript range errors. Modified the precipitation flux and river boundary conditions to work with a variable land-surface configuration. Added Liang-Barsky line clipping algorithm to identify river segments in cells. Changed some output labels from nodes to cells for clarity. Changed the flow rate calculation at each time step to use the same temporal weighting as chosen for the finite-difference system equations. Removed a redundant saturation fraction calculation for specified pressure cells. Added warnings of excessive movement of the water table (more than one cell layer per time step). Increased the number of digits in some output formats. The version of BCFLOW that accompanies HST3D Version 2.2.16 is BCFLOW Version 2.5.3

Version 2.2.13 is released (May 2006):

Made the dump/restart option functional under dynamic array allocation of Fortran 90. Improved some output message formats. Repaired some uninitialized variables. Fixed a bug in handling of printout frequency. Fixed an error in calculation of the balance tables under conditions of unconfined flow. Corrected some deallocation errors. Free surface location is only printed now for unconfined flow.

Version 2.2.12 is released (June 2005):

Corrected the conversion factors for output of changes in temperature. Fixed a bug in cumulative well amounts for case of both injection and production periods for a given well. Changed the step total and cumulative amounts to be calculated using the same time weighting as for the finite difference equations. Fixed print control to avoid printing out unused variables. Increased the number of digits in output of time values for plotting. Improved the time unit labels and balance table labels. Fixed an uninitialized printout variable. Included source code with the Linux distribution file.

Version 2.2.11 is released (March 2004):

Corrected the potentiometric head output for the case of tilted coordinates. Improved some error messages. Added a digit to the output of time values. Initialized some undefined, yet unused, variables in argument lists.

Version 2.0.12 is released (February 2004):

Fixed initialization of interpolation routine. Improved two error messages. Modified makefiles to work on Sun UltraSparc instead of Data General platform. The documentation for HST3D Version 1 is now available electronically as a pdf file. There will be no further development of the FORTRAN 77 source.

Version 2.2.8 is released (April 2003):

Source code is now converted to Fortran 90 enabling dynamic array allocation. Converted Bcnfr output file to binary. Corrected cross dispersion for variable mesh and boundaries and cylindrical coordinates. Improved well bore and solver error messages. Fixed round-off error for time level. Removed extra spaces in printouts of units of measure. Modified to allow free surface to be located below excluded cells and allow for keywords in the input file to be an arbitrary mix of upper and lower case letters. Reduced the potential for round-off error in computed nodal coordinates. Reformatted overflow fractional errors in balance tables to print as asterisks not zeros. Enabled current version and patch number to be printed in the problem definition file, Out.*. Modified format of output to screen. Fixed print control to handle a print interval that is not a multiple of a fixed-length time step.

Version 2.0.10 is released (August 2001):

Fixed bugs affecting input of river head data, printout when there is an entire column of excluded cells with unconfined flow, and the direct linear equation solver when there are excluded cells. Corrected the lengths of two common blocks for dump and restart. A header was added indicating when the printout of the identical second slice of a 2-slice region is being skipped. Restricted the redirection of flux to the cells containing the water table to unconfined flow only. Initialized some variables even when they are unused because heat and/or solute are not being simulated. Fixed an undefined variable for the well riser calculation.

Version 2.0.9 is released (May 2001):

Improved the print control and fixed a bug so that print intervals that are not multiples of the smallest print interval now work properly. A print event time no longer causes the time step to revert to the initial, smallest value under automatic time stepping. Solute concentrations in empty cells are now set to zero. Convergence criteria for the semi-implicit well flow allocation have been improved. An unconfined simulation now handles excluded cells above the free surface properly. Improved some number formats.

Version 2.0.8 is released (December 2000):

Print control frequency is now a real number. Revised initial condition pressure calculation. Now uses trivial equations for excluded cells in direct solver to fix a bug. New convergence criterion for iterative solver. Bug fix for flux b.c. into water-table cells. Bug fix for handling of dry cells in unconfined systems.

Version 2.0.7 is released (November 1999):

Restore routine ITER that was left out of Version 2.0.6. Documentation now available only as a separate package (pdf or postscript).

Version 2.0.6 is released (October 1999):

Patches to fix well iteration in HST routine INIT, zero well flow condition in HST routine WELLSR. Documentation now available only as a separate package (pdf or postscript).

Documentation of Version 2 is released as a pdf file (July 1999):

A pdf file of the documentation report, USGS/WRIR-97-4157, is available. Difficulties obtaining the Postscript version have been resolved.

Version 2.0.5 is released for Windows (June 1999):

A self-extracting file for Windows 95, 98, and NT4.0 is available. It includes executable files dimensioned sufficiently large to execute all the examples.

Version 2.0.5 is released (March 1999):

Patches to fix ascending node range in HST routine HUNT, cell saturation fraction in HST routine INIT2, and printout cosmetic changes. Patch to fix bug in DIMEN routine READ1.

Version 2.0.4 is released (February 1999):

Patch to fix a conflict in temporary storage for rivers and head printouts in HST routines ERROR4, INIT2, READ2, WRITE2. Reduce stripped file record lengths in DIMEN routine READ1.

Version 2.0.3 is released (January 1999):

Patch to add some undeclared variables, output equation name to screen, remove some obsolete solution method indexing in HST routine APLBCI.

Version 2.0.2 is released (December 1998):

Patch to fix divide by zero bug in DIMEN, restart bug in HST, include static information in restart output files, increase input data file name length to 20 characters in HST, DIMEN, BCFLOW.

Version 2.0.1 is released (October 1998):

Patch to fix some format errors in routine REWI3.

Version 2.0 is released (August 1998):

Documentation: Guide to the revised heat and solute transport simulator: HST3D, by Kenneth L. Kipp, Jr., U.S. Geological Survey Water-Resources Investigations Report 97-4157, 1997, 149 p.

Version 1.4 is released (1989):

Documentation: HST3D: A computer code for simulation of heat and solute transport in three-dimensional ground-water flow systems, by Kenneth L. Kipp, Jr., U.S. Geological Survey Water-Resources Investigations Report 86-4095, 1987, 597 p.