#! /bin/bash -f
# Last edited on 2009-12-11 23:35:32 by stolfi

echo "PLOTTING MAJOR GLITCHES AND ADJUSTMENTS" 1>&2 

D=610.0
R=0.0006550

A1=40818
K1=5170.0
B1=1.0000
echo "B1 = ${B1}" 1>&2 

AA=`echo "${A1} + 36700" | bc -lq`
KA=`echo "${K1}" | bc -lq`

plot-counter-glitch.sh \
    2.00 ${D} ${R} \
    ${A1} ${K1} \
    ${AA} ${KA} \
    data/counter-glitch-2003-xx-irr.txt \
  > rslt/counter-glitch-2003-xx-irr.eps

A2=`echo "${A1}+ 4650" | bc -lq`
K2="${K1}"
B2=`echo "${K2}/${K1}" | bc -lq`
echo "B2 = ${B2}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A1} ${K1} \
    ${A2} ${K2} \
    data/counter-glitch-2002-09-irr.txt \
  > rslt/counter-glitch-2002-09-irr.eps

A3=`echo "${A2}+31850" | bc -lq`
K3="${K2}"
B3=`echo "${K3}/${K1}" | bc -lq`
echo "B3 = ${B3}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A2} ${K2} \
    ${A3} ${K3} \
    data/counter-glitch-2002-10-irr.txt \
  > rslt/counter-glitch-2002-10-irr.eps

A4=`echo "${A3}*0.994 +   0" | bc -lq`
K4=`echo "${K3}*0.994" | bc -lq`
B4=`echo "${K4}/${K1}" | bc -lq`
echo "B4 = ${B4}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A3} ${K3} \
    ${A4} ${K4} \
    data/counter-glitch-2002-11-irr.txt \
  > rslt/counter-glitch-2002-11-irr.eps

A5=`echo "${A4}*0.978 +    0" | bc -lq`
K5=`echo "${K4}*0.978" | bc -lq`
B5=`echo "${K5}/${K1}" | bc -lq`
echo "B5 = ${B5}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A4} ${K4} \
    ${A5} ${K5} \
    data/counter-glitch-2002-12-irr.txt \
  > rslt/counter-glitch-2002-12-irr.eps

A6=`echo "${A5}*0.986 +   0" | bc -lq`
K6=`echo "${K5}*0.986" | bc -lq`
B6=`echo "${K6}/${K1}" | bc -lq`
echo "B6 = ${B6}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A5} ${K5} \
    ${A6} ${K6} \
    data/counter-glitch-2003-01-irr.txt \
  > rslt/counter-glitch-2003-01-irr.eps

A7=`echo "${A6}*0.983 +    0" | bc -lq`
K7=`echo "${K6}*0.983" | bc -lq`
B7=`echo "${K7}/${K1}" | bc -lq`
echo "B7 = ${B7}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A6} ${K6} \
    ${A7} ${K7} \
    data/counter-glitch-2003-04-irr.txt \
  > rslt/counter-glitch-2003-04-irr.eps

A8=`echo "${A7}*1.046 + 0" | bc -lq`
K8=`echo "${K7}*1.046" | bc -lq`
B8=`echo "${K8}/${K1}" | bc -lq`
echo "B8 = ${B8}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A7} ${K7} \
    ${A8} ${K8} \
    data/counter-glitch-2003-05-irr.txt \
  > rslt/counter-glitch-2003-05-irr.eps

A9=`echo "${A8}*1.014 + 0" | bc -lq`
K9=`echo "${K8}*1.014" | bc -lq`
B9=`echo "${K9}/${K1}" | bc -lq`
echo "B9 = ${B9}" 1>&2 

plot-counter-glitch.sh \
    1.00 ${D} ${R} \
    ${A8} ${K8} \
    ${A9} ${K9} \
    data/counter-glitch-2003-07-irr.txt \
  > rslt/counter-glitch-2003-07-irr.eps