# Time-stamp: </home/wrf/git/all/c/overprop/fns, Tue, 20 Nov 2018, 16:34:33 EST, w540, wrf, W. Randolph Franklin, http://wrfranklin.org/>

# Run program and save output 

function doit() {
./overlay${V:-10} $* | tee log
}

# Simple example

function do1() {
doit db/sq/1 db/sq/2 1
}

function do2() {
doit db/sq/1 db/sq/2 5
}

# Bigger example.  Output has errors, perhaps due to bad input..

function do5() {
doit db/ch/county db/ch/hydrog ${NG:-1000}
}

# Bigger examples from Salles.

function do6() {
A=db/salles
doit $A/brc $A/brs ${NG:-1000}
}

function do7() {
A=db/salles
doit $A/usa $A/usc ${NG:-1000}
}




# Select a polygon from county into a separate map file, extract data, and plot.
# E.g., selplot1 12345

function selplot1() {
  selectpol $1 < db/ch/county >! db/ch/county-$1
  plotpol $1 < db/ch/county-$1 >! county-${1}.dat
  echo "plot \"county-${1}.dat\" with linespoints" | gnuplot -p
}

function selplot2() {
  selectpol $1 < db/ch/hydrog >! db/ch/hydrog-$1
  plotpol $1 < db/ch/hydrog-$1 >! hydrog-${1}.dat
  echo "plot \"hydrog-${1}.dat\" with linespoints" | gnuplot -p
}

# Do both county and hydrog polygons
function selplot() {
  selectpol $1 < db/ch/county >! db/ch/county-$1
  plotpol $1 < db/ch/county-$1 >! county-${1}.dat
  selectpol $2 < db/ch/hydrog >! db/ch/hydrog-$2
  plotpol $2 < db/ch/hydrog-$2 >! hydrog-${2}.dat
  echo "plot \"county-${1}.dat\" with linespoints, \"hydrog-${2}.dat\" with linespoints" | gnuplot -p
}


# The original floating version
function doitf() {
F=logf
echo '==========================='  >>! $F
echo Maps: $1 $2  >>! $F
echo Grid size: $3  >>! $F
echo >>! $F
echo db/$1 db/$2 $3 | ./overpropf >>! $F
tail -20 $F
}
