Vadose Zone Leaching Model, Version 2.2c

Version:  2.2c
 May 2006, Thomas Harter, University of California, Davis

Based on:
 Version (2.2b) updated by Thomas Harter in June 1998 (see Notes below)

 Version (2.2a) updated by Rashid Islam (CDSI) in June 1996
 (fixed the bugs in ground water impact calculations)

 Version (2.2) updated by Rashid Islam (CDSI) in October 1995
 Modified exponetial constant in the Millington equation for calculating 
 free-air diffusion coefficient from 10/3 to 7/3.  Also corrected a bug in 
 the Total ground water impact output.

 Version (1.0) of VLEACH developed by CH2M Hill for USEPA, Region IX,
 version updated by Varadhan Ravi (Dynamac Corp.) in March 1994

---------------------------------------------------------------------

Notes on version 2.2c:
 This version is identical to version 2.2b. But it has been recompiled for
 use in the Windows XP environment and the array dimensions have been adjusted.

The program allows for as much as 100,000 cells per polygon and as many as
500 printout times for the profile. Minimum time-step size is still 0.0001 years.

To run the program in Windows XP:

 * start Windows Explorer
 * make sure that "vl22c.exe" and "vleach-input.txt" are in the same working directory
 * find the file "vleach-input.txt" in your vleach working directory
 * double-click on "vleach-input.txt" to edit the input file in Notepad(tm)
 * click on FILE, then click on SAVE to save the input file.

 * click "Start"
 * click "Run"
 * browse to find "vl22c.exe"
 * double click on "vl22c.exe"
 * click "ok"
 A window will briefly appear and disappear.  Check the file "vleach.prm" to
 make sure that the input file was read correctly. If you have problems
 you may want to open a DOS Window ("Start" -> "Programs" -> "Accessories"
 -> "Command Prompt") and run vl22c.exe from there.

For further instructions, see Notes on Version 2.2b below.

FOR A SUMMARY OF THE OUTPUT FILES GENERATED, SEE BOTTOM OF THIS NOTE.

Signed, THomas Harter, May 9, 2006.

---------------------------------------------------------------------

Notes on Version 2.2b:

 This version (2.2b) updated by Thomas Harter (Dep. of Land, Air,
 and Water Resources, Univ. of California, Davis;  thharter@ucdavis.edu).

 This version performs the identical computations to VLEACH version 2.2a,
 which I downloaded in May 1998 from EPA's Center for Subsurface Modeling
 Support (CSMOS) at EPA's Research Laboratory in Ada, Oklahoma.

 The differences to version 2.2a are the following:

1.The input file is different from the originial input files to make it
 easier to read the input data.  The program cannot read the original
 input files (not upwardly compatible)!  The input file is divided into
 three sections: general project information, contaminant information, and
 polygon information.  The entire polygon section (including the initial
 concentration distribution subsection) must be repeated for each polygon.

 The format of each input line is:
 40 characters (including blank spaces), number

 The format of the initial concentration profile information is:
 40 characters,integer, comma, integer, comma, real number (no spaces)

 Be sure that the number starts on the 41st position of the line!
 An example input file is included.  The easiest way to create an input
 file is to modify the example input file.  Note that the input file MUST
 be named "vleach-input.txt" and must reside in the same directory as the
 program "vl22b.exe".  I recommend that you always make a copy of your
 input file to another file.

2.The program has been modified to accomodate much larger columns.  The
 program now accepts up to 9999 model cells.  It is recommended to use
 a vertical cell size of 0.1 feet or less to reduce the amount of numerical
 dispersion to a minimum.  The smaller cell size will generally require a
 time-step less than 1 year or even less than 1 day.  The program has been
 altered to be able to handle time-steps as small as 0.0001 years.  In
 earlier versions, time-steps smaller than one year caused output problems.
 The total time, gw impact output and profile output times must be exact
 multiples of the time-steps!  The time-step should be small enough that
 the advective transport distance within one time-step is less than the
 cell thickness (Courant criterion).

3.A new transient profile output file is created to allow for easy plotting
 in spreadsheet programs such as Microsoft Excel(tm), or Corel Quattro Pro (tm).
 The new output file is called "trnsprof.dat".  For each polygon, it collects
 the profile information at each requested profile output time interval.  The
 output is the concentration in the liquid phase (concentrations in the gas
 and sorbed phase are linearly proportional to the liquid phase, i.e. the 
 shape of the profiles are identical in all three phases). The output under
 each polygon contains:

 1st row:               number_of_cells time0, time1, time2 .... endtime
 other rows:    depth_of_center_of_cell conc(time0) conc(time1)... conc(endtime)
 (1 row per cell)

 This file can be imported into a spreadsheet as a space delimited file (with
 repeated spaces counted as one delimiter).  By highlighting ALL rows below
 the first row (up to the next polygon), starting with the column containing
 the depth information, a x-y graph containing all the profiles can easily be
 created in the spreadsheet.  Some programs will give a choice of "data in rows"
 or "data in columns".  Since each column presents one profile (at a given time)
 the user should select "data in columns".

 Alternatively, if ALL rows below the first row are highlighted WITHOUT the
 first data column (which contains the depth information), and the user
 selects "data in rows", the result will be the breakthrough curves for each
 cell, beginning with the top-most cell.

4.The program was compiled for Windows95 and should run in Windows95 (tm),
 Windows NT (tm), and Windows98 (tm).  The recommended minimum configuration
 is an Intel 486 compatible processor and at least 24 MB of memory (the
 program uses probably around 1 - 4 MB).

5.To run the program in Win95:

 start Windows Explorer
 find the file "vleach-input.txt" in your vleach working directory
 double-click on "vleach-input.txt" to edit the input file in Notepad(tm)
 click on FILE, then click on SAVE to save the input file.

 click "Start"
 click "Run"
 browse to find "vl22b.exe"
 double click on "vl22b.exe"
 click "ok"
 A window will briefly appear and disappear.  Check the file "vleach.prm" to
 make sure that the input file was read correctly.

5.No warranties attached.  You are using the program at your own risk.
 Please refer to the VLEACH manual published by EPA for details on VLEACH.

Thomas Harter, June 29, 1998.

===========================================================================

Output Files (all are ASCII files that can be opened with NOTEPAD):

VLEACH.PRM  - This file echoes the input. Useful to ensure that input file is
              read correctly.

VLEACH.PRF  - This file is always the largest file. It contains the solid,
              liquid, and gas phase concentrations for each polygon at all
              depth, and for each time interval in the input line "gw impact
              output interval"

VLEACH.OUT  - File that contains the mass balance report for each time interval
              in the input line "gw impact output interval"

TRNSPROF.DAT - File that contains the soil aqueous phase profile for each time
              interval in the input line "transient profile plot interval".
              This file intended to be used for plotting multiple profiles
              with a spreadsheet program.

SOILIMP.DAT - File with the soil profile at the time "plot time" in the input file.
              The file contains two columns: The first is the concentration, the
              second is the depth.

GWIMP.DAT   - File with the groundwater impact concentrations for each polygon at
              each time interval in the input line "gw impact output interval".

Note that the last two files are only printed if "y" is chosen in the input line
"plot SOILIMP, GWIMP [y,n]".


