#!/bin/bash

# for hkl files, $2 is (0) or (1), $7 is 2theta, and #9 is intensity

awk '{print $8}' $1 > .hkl1a

# strip out parenthesis around phase identifier
sed 's/,//' .hkl1a > .hkl1b

awk '{print $1, "1"}' .hkl1b > .hklbar1

sed '1,5d' .hklbar1 > .hklbar2
sed '/^ 1$/d' .hklbar2 > .hklbar

# intensity of hkl isnt needed for line plot, but we will hide that dataset
rm .hkl1a
rm .hkl1b

cat > .bfile << !

with g0
    default font 4
    world xmin rint(min(g0.s0.x))
   world xmax rint(max(g0.s0.x))
    world ymax rint(1.1*maxof((max(g0.s0.y)),(max(g0.s1.y))))
    xaxis label "CuK\xa\f{} 2\xq\f{} (\So\N)"
    xaxis tick major 20
    yaxis label "counts"
    view .15, .15, .65, .62 
    xaxis  label font 4
    yaxis  label font 4
    xaxis  ticklabel font 4
    yaxis  ticklabel font 4
    s0 symbol 1
    s0 symbol size 0.2500000
    s0 line type 0
    s1 line linewidth 1  
    s2 line linewidth 1  
    s1 line color "red"
    s2 line color "navy"

    s1.y = s1.y + 0.1*(max(g0.s0.y))
    s0.y = s0.y + 0.1*(max(g0.s0.y))

with g1
    default font 4
    world xmin rint(min(g0.s0.x))
    world xmax rint(max(g0.s0.x))
    xaxis tick off
    xaxis ticklabel off
    yaxis ticklabel off
    yaxis tick major 200
    xaxis tick major 20
    view .15, .63, .65, .65
    world ymin 0.0
    world ymax 1.0
    s0 symbol size 0.000000
    xaxis  label font 4
    yaxis  label font 4
    xaxis  ticklabel font 4
    yaxis  ticklabel font 4 
!

awk '{print $2, $3, $4, $6}' $1 > .tempplg1

sed 's/,//' .tempplg1 > .tempplg2
sed '1,7d' .tempplg2 > .tempplg

xmgrace -pexec "arrange (2,1,.1,.1,.1,ON,ON,ON)" -nxy .tempplg -graph 1 -settype bar .hklbar -batch .bfile 

rm .grepout
rm .bfile
rm .tempplg
rm .tempplg1
rm .tempplg2
rm .hklbar
rm .hklbar1
rm .hklbar2