# first, save what we've got already
#cp -F saveStats.tar save
#tar -cvf saveStats.tar derived/* *.png
## scp files from cori

# first build utilities if needed
cd util
make all
cd ..

#scp dmartin@keep.lbl.gov:people/newRheology/newerMISMIP+/experimentStats/experimentStats.tar .
tar -xf experimentStats.tar
mv -f *.allOut data
ls 
set WD = $PWD


#set PLOTVAR="plot"
#cp -f VoF-template.gp VoFPlot.gp

# do stats for each experiment
foreach EXPT ( n3Control n3Control-noMuCoef  n4Control n3Experiment n3Experiment-ABUMIP  n4Experiment n4Experiment-ABUMIP )
set FILE = data/$EXPT.allOut
echo $FILE
awk/sumOver.awk < $FILE > derived/sumOver-$EXPT
awk/totalVol.awk < $FILE > derived/totalVol-$EXPT
awk/groundedArea.awk < $FILE > derived/groundedArea-$EXPT

util/computeDiffs  derived/sumOver-$EXPT derived/sumOverDiff-$EXPT
util/computeDiffs  derived/totalVol-$EXPT  derived/totalVolDiff-$EXPT
util/computeDiffs  derived/groundedArea-$EXPT derived/groundedAreaDiff-$EXPT

end

# process relative differences
util/subtractTwoRelative derived/sumOverDiff-n3Experiment derived/sumOverDiff-n4Experiment derived/sumOverDiff-relDiff

util/subtractTwoRelative derived/groundedAreaDiff-n3Experiment-ABUMIP derived/groundedAreaDiff-n4Experiment-ABUMIP derived/groundedAreaDiff-ABUMIP-relDiff

util/subtractTwoRelative derived/sumOverDiff-n3Experiment-ABUMIP derived/sumOverDiff-n4Experiment-ABUMIP derived/sumOverDiff-ABUMIP-relDiff

util/subtractTwoRelative derived/groundedAreaDiff-n3Experiment derived/groundedAreaDiff-n4Experiment derived/groundedAreaDiff-relDiff


cd plots
foreach PLOT (*.gp)
gnuplot $PLOT
end
cd ..
display *.png &



#gnuplot VoFPlot.gp
#gnuplot groundedAreaPlot.gp
#gnuplot totalVolPlot.gp

#gnuplot VoFPlot-small.gp
#gnuplot groundedAreaPlot-small.gp
#gnuplot totalVolPlot-small.gp

#gnuplot VoFvsBedmapAllPlot.gp
#gnuplot groundedAreaVsBedmapPlot.gp
#gnuplot totalVolVsBedmapPlot.gp

#display *.png &
 
#cd ..