

NUMERICA,
A Library of Source Codes for 
Teaching, Research and Applications
by E. F. Toro
NUMERITEK Limited, 1999, UK


HYPER-LIN Sub-library


INTRODUCTION

Welcome to the HYPER-LIN sub-library. All files in these directories 
are DOS formatted, i.e. they conform to the extended DOS character 
set and can thus be used directly on any IBM compatible PC equipped 
with a FORTRAN compiler. 

This sub-library contains 20 source FORTRAN 77 codes for solving the 
linear advection equation, the inviscid Burgers equation, the linearised 
equations of gas dynamics and the linearised shallow water equations. 
The programming language is STANDARD FORTRAN 77, the only exception 
being the statement IMPLICIT NONE. These source codes originate from early 
research and teaching codes developed by the author over a period of 
many years of involvement in research, applications and teaching 
activities. The codes have been used as demonstration codes in regular 
courses on numerical methods given to MSc students in two British 
Universities and in short courses on Computational Fluid Dynamics given 
in several countries.

THE PROGRAMS

The 20 programs of HYPER-LIN are found in the 4 directories:

burger:  the inviscid Burgers equation
linadv:  the linear advection equation
lingas:  the linearised equations of gas dynamics 
linwat:  the linearised shallow water equations

Each FORTRAN 77 program  "progra.f" is self-contained and only requires 
an input data file "progra.ini". Numerical results are printed to an 
output file "progra.out". The programs can be copied directly to any 
computer with a FORTRAN 77 compiler and used immediately. The input file 
"progra.ini" can be edited and changed as desired. 

VISUALISATION

There is a selection of commercial plotting packages that can be used to 
visualise the numerical results. For simplicity we assume the user has 
access to the public domain plotting package "GNUPLOT". This can be 
downloaded from the web, visit http://www.gnu.org for further details. 

The one-dimensional codes produce output of the solution at a specified 
time TIMEOU. The output is printed in  the form of columns:

x1      y1    z1   ........

where x1 is a position in the spatial domain, y is the corresponding value 
of a quantity y (eg density) and z1 is that of a quantity z (eg pressure). 
The profile (x, y) can be plotted by first typing <gnuplot>. Then, within 
the GNUPLOT environment type:

<plot 'progra.out' u 1:2> This will plot the profile using symbols.

Similarly,  to plot the profile (x, z) one types:

<plot 'progra.out' u 1:3>

To produce both symbols and a line one types:

<plot 'progra.out' u 1:2,'progra.out' u 1:2 w l>

 
The user is encouraged to utilise more sophisticated visualisation 
software  available in their systems.


SUPPORT

Users on support agreement may contact NUMERITEK LIMITED regarding any 
matter connected to the NUMERICA programs by emailing: 

support@numeritek.com.  

All users are encouraged to report bugs or any problems encountered by 
emailing: 

bugs@numeritek.com.


IMPORTANT NOTICE:

NUMERICA is the copyright of NUMERITEK LIMITED and its use is governed 
by the terms of the License Agreement accompanying this package.  
 
If you can not locate the License Agreement please write to NUMERITEK 
LIMITED at 

PO Box 68 
Bramhall 
Cheshire 
SK7 3FT 
United Kingdom 

or email inquiries@numeritek.com for a copy.  


 NUMERITEK LTD., 1999
