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