CFD Utilities Files

program finterp1d

Purpose:

   This adaptation of earlier basic utilities was prompted by a set of sparse
   radiative heat flux calculations along a capsule forebody for multiple
   flight conditions along a trajectory.  Using SMOOTH to interpolate each
   coarse dataset to the full generatrix for one case was cumbersome enough
   that automating the interpolations for a list of datasets was called for.
   Application to other comparable scenarios should be possible.

Sample List of Dataset Files to Process:

   /nobackup/smuppidi/DAVINCI/SIMS/NOM/t100/FOREBODY-TANGENT-SLAB/qrad.dat
   /nobackup/smuppidi/DAVINCI/SIMS/NOM/t105/FOREBODY-TANGENT-SLAB/qrad.dat
   /nobackup/smuppidi/DAVINCI/SIMS/NOM/t110/FOREBODY-TANGENT-SLAB/qrad.dat
   :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::

Sample Dataset, where (y, qrad) data are to be interpolated at dense ys::

   VARIABLES = "x" "y" "z" "qrad W/cm2"
   0 0 0 42.82057
   0.025 0.1773849 0 42.30668
   0.05 0.24842058 0 41.31608
   0.1 0.344178092 0 37.03088
   0.2 0.466090926 0 21.86707
   0.3 0.566090908 0 15.90821
   0.4 0.666090930 0 14.52199
   0.5 0.766090930 0 12.39096
   0.6 0.866090909 0 10.56066
   0.7 0.966090916 0 9.39016
   0.8 1.066090910 0 8.57758
   0.9 1.166090870 0 8.03349
   1.0 1.282534878 0 7.64047

Sample target abscissa file (x, y):

   0.000000000E+00 0.000000000E+00
   4.764895912E-05 7.824138284E-03
   1.895706537E-04 1.560527311E-02
   4.246242209E-04 2.335331111E-02
   :::::::::::::::::::::::::::::::
   1.074885464E+00 1.234345831E+00
   1.078287486E+00 1.230943809E+00
   1.081779012E+00 1.227452283E+00
   1.085003328E+00 1.224227968E+00

Sample Output (interpolations.dat, Tecplottable):

   TITLE = "FINTERP1D interpolations"
   VARIABLES = "x" "y" "z" "qrad W/cm2"
   0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 4.2820571234E+01

Programming Considerations:

   o  Namelist control (standard input) is most flexible/extensible.
   o  One or more dataset file names are read from a dataset list file
      indicated on standard input.
   o  Datasets are in Tecplot form, multi-column, with (abscissa, ordinate)
      columns selectable at run time.
   o  More than one function/ordinate may be interpolated in a single run.
      All of the functions present are interpolated in order.
   o  A single target file (one column or more, no header records) applies to
      all datasets.  Its abscissa column number is specified at run time.
   o  Minor extrapolation is tolerable during the interpolations.
   o  Outputs are in Tecplot point-order form, with as many zones as input
      datasets.
   o  (Initially,) a single interpolation method is provided: the local cubic
      spline method of LCSFIT.  Some day, a smoothing option may be called
      for to deal with noisy data.

Namelist Control ('finterp1d.inp' on standard input; defaults shown):

   $FINTERP1D_INPUTS
   verbose = F
   dataset_list = 'finterp1d.list'
   target_data = 'target.dat'
   interpolated_data = 'finterp1d.dat'
   data_abscissa_col = 1
   target_abscissa_col = 1
   interpolation_method = 'lcsfit' ! Only choice initially
   lcs_method = 'M'                ! L|M|B = linear|monotonic|"Bessel"=loose
   $END

History:

   11/22/2024  D.A.Saunders  Initial design and implementation, for padding
                             forebody tangent-slab radiation results at
                             multipole trajectory points. (SURFACE_INTERP_2D
                             has somewhat related functionality.)
   11/23/2024    "     "     Interpolating all functions present is easier
                             than one or more in specified order. Interpolat-
                             ing one or two other coordinates is messy!
   11/24/2024    "     "     Initial testing is done.
   11/29/2024    "     "     Allowed for verbose mode while handling datasets
                             with missing file header and/or zone header.

Author:  David Saunders, AMA, Inc. at NASA Ames Research Center, CA.