Author: William Chan
ifar = 0 do not prompt for a variable far field distance
= 1 constant far field distance
= 2 piece-wise linear 1-D function in J and K respectively,
distance specified at nodes.
= 3 compute distance to centroid in J or K stations,
good for internal grids.
= 4 transfinite interpolation in J and K, distance specified
at nodes on the four boundaries
= 10 read in bilinear surface and do shooting (to be implemented
in future)
idzs = 0 do not prompt for a variable initial grid spacing
= 1 constant initial grid spacing
= 2 piece-wise linear 1-D function in J and K respectively,
spacing specified at nodes.
= 4 transfinite interpolation in J and K, distance specified
at nodes on the four boundaries.
= 10 read in volume grid and use final spacing in L as initial
spacing (jmax and kmax must match of course).
idze = 0 do not prompt for a variable end grid spacing
= 1 constant end grid spacing
= 2 piece-wise linear 1-D function in J and K respectively,
spacing specified at nodes.
= 4 transfinite interpolation in J and K, distance specified
at nodes on the four boundaries
zeta(j,k,1) - variable far field distance
zeta(j,k,2) - variable initial spacing
zeta(j,k,3) - variable end spacing
Ascii/Unformatted.
Little/Big Endian for unformatted files.
Single/Double precision for unformatted files.
Single/Multiple grid file.
2D/3D.
Only files with 1 grid are processed.
Output file format follows the endian and precision of the input file.