#! /bin/bash

# Define variables for grid, solution, and chemistry file
grid=run.cgd
sol=run.cfl
chm=h2-air-7sp-mbody-3287.chm

# Create a CFPOST input script
echo "
  clear all

  grid $grid
  sol  $sol
  chem file $chm

  units english

  subset i all j 8 k 1
  zone 1 to 10

  variable x in; y in; C-H2; C-O2; C-H2O 
  list output massfrac.lis raw lines 55555
  exit
" > massfrac.jou

# Run CFPOST and supply the input script
cfpost_prod.exe < massfrac.jou

# Strip header from CFPOST output file
grep -v -e 'Subset' -e '-x-' -e 'in' massfrac.lis > massfrac.dat

# Add new header for use with Tecplot
sed -e '1i\ TITLE     = "Str. Mass Frac."
1i\ VARIABLES = "x [in]"
1i\ "y [in]"
1i\ "C-H2"
1i\ "C-O2"
1i\ "C-H2O"
1i\ ZONE T="Str. ", F=POINT' massfrac.dat > mf.dat

# Cleanup
mv mf.dat ./results/mf.dat
rm mass*.lis mass*.jou mass*.dat cfpost.jou