#!/bin/csh

set Cols  = (Black=1860-1899 Red=1900-1939 Green=1940-1979 Blue=1980-2019)
set MIL   = (Millions of Barrels a Day)
set BIL   = (Billions of Barrels So Far)
set GRAPH = (xgraph style=n width=800 height=400 linewidth=2 grid1=dot grid2=dot label1=Year \
             linecolor=0,2,3,4,0,2,3,4 mark=8 marksize=6)
set pause = (sleep 1)
set pause = pause

@ n = `cat US_crude_oil_production_per_day.txt | wc -l`
@ no4 = $n / 4

set T = /tmp/j.txt

wc US_crude_oil_production_per_day.txt $T

<US_crude_oil_production_per_day.txt awk '{sum += $2; print $1,365.25*sum/1000000,0.001*$2}' >$T

<$T awk '{print $1,$2}' | a2b | $GRAPH n=$n label2="$BIL" title="Production So Far" &; $pause
<$T awk '{print $1,$3}' | a2b | $GRAPH n=$n label2="$MIL" title="Production a Day"  &; $pause

<$T awk '{print $1,$2}' | a2b | $GRAPH n=$no4 nplot=4 label2="$BIL" title="Production So Far $Cols"&; $pause
<$T awk '{print $1,$3}' | a2b | $GRAPH n=$no4 nplot=4 label2="$MIL" title="Production a Day $Cols"&; $pause

<$T awk '{print $1,log($2)/log(10)}' | a2b \
| $GRAPH n=$no4 nplot=4 label2="log $BIL" title="ProductionSoFar -- logarithmic" &; $pause
<$T awk '{print $1,log($3)/log(10)}' | a2b \
| $GRAPH n=$no4 nplot=4 label2="log $MIL" title="Production -- logarithmic" &; $pause

<$T awk '{print log($2)/log(10),log($3)/log(10)}' | a2b \
| $GRAPH n=$no4 nplot=4 label2="log $MIL" label1="log $BIL" title="log-log X-plot $Cols" &; $pause

<$T awk '{print $2,$3}' | a2b \
| $GRAPH n=$no4 nplot=4 label2="$MIL" label1="$BIL" title="X-plot $Cols" &; $pause