Source:
http://people.bu.edu/robinc/kzk

This set of files constitutes the time-domain computer code developed at
The University of Texas at Austin to model axisymmetric sound beams in
fluids.  The code is based on an augmented KZK equation that accounts for
nonlinearity, diffraction, thermoviscous absorption, and absorption and
dispersion due to an arbitrary number of relaxation phenomena.  The
principal references for this code are the following:
 
Y.-S. Lee, "Numerical solution of the KZK equation for pulsed finite
amplitude sound beams in thermoviscous fluids," Ph.D. dissertation, The
University of Texas at Austin (1993).

Y.-S. Lee and M. F. Hamilton, "Time-domain modeling of pulsed
finite-amplitude sound beams," J. Acoust. Soc. Am., Vol. 97, pp. 906-917
(1995).

R. O. Cleveland, M. F. Hamilton, and D. T. Blackstock,  "Time-domain
modeling of finite-amplitude sound in relaxing fluids," J. Acoust. Soc.
Am., Vol. 99, 3312-3318 (1996).

The dissertation by Lee (1993) contains the original listing of the code,
description of algorithm, and test cases.  Lee's dissertation is summarized
in the article by Lee and Hamilton (1995).  The article by Cleveland et al.
(1996) describes the inclusion of relaxation and other modifications of the
code.

  The files that should be in this bundle are:

kzktexas1.f    the source code

GetCols.f     code that pulls the output file "pulse.out" apart so that individual
                  waveforms can be viewed

yslmfhFig4    input file to generate Fig 4 from the Lee and Hamilton JASA paper

relaxFig3     input file to generate Fig 3 from Cleveland et al JASA paper. Note to get
              the TV case only set the relaxation flag to 0.

relaxFig3input  This does the same as the above file except it reads the
                inital waveform in from a file rather than generating it
                itself.

C4M6P240      The file with the initial pressure waveform for the above 
              input file. The cryptic name stands for cycles=4, menvelope=6,
              number of tau points = 240.