diff --git a/CEP/BB/BBS/test/makeDemoParmDB b/CEP/BB/BBS/test/makeDemoParmDB
new file mode 100755
index 0000000000000000000000000000000000000000..ceeda37443614e355b8a3cd71036168d8076ebd5
--- /dev/null
+++ b/CEP/BB/BBS/test/makeDemoParmDB
@@ -0,0 +1,57 @@
+#!/bin/sh
+
+pert="1e-6"
+
+echo "creating instrument_model (aips++)..."
+parmdb > $1_tmp.outpdb1 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_instr'
+adddef gain:11:real values=1, perturbation=${pert}
+adddef gain:12:real values=0, perturbation=${pert}
+adddef gain:21:real values=0, perturbation=${pert}
+adddef gain:22:real values=1, perturbation=${pert}
+adddef gain:11:imag values=1, perturbation=${pert}
+adddef gain:12:imag values=0, perturbation=${pert}
+adddef gain:21:imag values=0, perturbation=${pert}
+adddef gain:22:imag values=1, perturbation=${pert}
+adddef gain:11:ampl values=1, perturbation=${pert}
+adddef gain:12:ampl values=0, perturbation=${pert}
+adddef gain:21:ampl values=0, perturbation=${pert}
+adddef gain:22:ampl values=0, perturbation=${pert}
+adddef gain:11:phase values=0, perturbation=${pert}
+adddef gain:12:phase values=0, perturbation=${pert}
+adddef gain:21:phase values=0, perturbation=${pert}
+adddef gain:22:phase values=0, perturbation=${pert}
+adddef frot values=0, perturbation=${pert}
+adddef drot values=0, perturbation=${pert}
+adddef dell values=0, perturbation=${pert}
+adddef dgain:X values=1, perturbation=${pert}
+adddef dgain:Y values=1, perturbation=${pert}
+quit
+EOF
+if [ $? != 0 ]; then
+  exit 1
+fi
+
+echo "creating sky_model (aips++)..."
+parmdb > $1_tmp.outpdb2 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_sky'
+adddef RA:SRC1  values=2.734003,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC1 values=0.4537895, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC1 values=1, perturbation=${pert}
+adddef StokesQ:SRC1 values=0, perturbation=${pert}
+adddef StokesU:SRC1 values=0, perturbation=${pert}
+adddef StokesV:SRC1 values=0, perturbation=${pert}
+adddef RA:SRC2  values=2.7340205,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC2 values=0.45368975, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC2 values=1, perturbation=${pert} 
+adddef StokesQ:SRC2 values=0, perturbation=${pert} 
+adddef StokesU:SRC2 values=0, perturbation=${pert}
+adddef StokesV:SRC2 values=0, perturbation=${pert}
+adddef RA:SRC3  values=2.733981,   perturbation=${pert}, pert_rel=F
+adddef DEC:SRC3 values=0.45375025, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC3 values=1, perturbation=${pert} 
+adddef StokesQ:SRC3 values=0, perturbation=${pert} 
+adddef StokesU:SRC3 values=0, perturbation=${pert}
+adddef StokesV:SRC3 values=0, perturbation=${pert}
+quit
+EOF
diff --git a/CEP/BB/BBS/test/makeDemoParmDBCorr b/CEP/BB/BBS/test/makeDemoParmDBCorr
new file mode 100755
index 0000000000000000000000000000000000000000..78a352d4f917fcb5568e5839ae56835a6a625122
--- /dev/null
+++ b/CEP/BB/BBS/test/makeDemoParmDBCorr
@@ -0,0 +1,57 @@
+#!/bin/sh
+
+pert="1e-6"
+
+echo "creating instrument_model (aips++)..."
+parmdb > $1_tmp.outpdb1 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_instr'
+adddef gain:11:real values=1, perturbation=${pert}
+adddef gain:12:real values=0, perturbation=${pert}
+adddef gain:21:real values=0, perturbation=${pert}
+adddef gain:22:real values=1, perturbation=${pert}
+adddef gain:11:imag values=0, perturbation=${pert}
+adddef gain:12:imag values=0, perturbation=${pert}
+adddef gain:21:imag values=0, perturbation=${pert}
+adddef gain:22:imag values=0, perturbation=${pert}
+adddef gain:11:ampl values=1, perturbation=${pert}
+adddef gain:12:ampl values=0, perturbation=${pert}
+adddef gain:21:ampl values=0, perturbation=${pert}
+adddef gain:22:ampl values=1, perturbation=${pert}
+adddef gain:11:phase values=0, perturbation=${pert}
+adddef gain:12:phase values=0, perturbation=${pert}
+adddef gain:21:phase values=0, perturbation=${pert}
+adddef gain:22:phase values=0, perturbation=${pert}
+adddef frot values=0, perturbation=${pert}
+adddef drot values=0, perturbation=${pert}
+adddef dell values=0, perturbation=${pert}
+adddef dgain:X values=1, perturbation=${pert}
+adddef dgain:Y values=1, perturbation=${pert}
+quit
+EOF
+if [ $? != 0 ]; then
+  exit 1
+fi
+
+echo "creating sky_model (aips++)..."
+parmdb > $1_tmp.outpdb2 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_sky'
+adddef RA:SRC1  values=2.734003,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC1 values=0.4537895, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC1 values=1, perturbation=${pert}
+adddef StokesQ:SRC1 values=0, perturbation=${pert}
+adddef StokesU:SRC1 values=0, perturbation=${pert}
+adddef StokesV:SRC1 values=0, perturbation=${pert}
+adddef RA:SRC2  values=2.7340205,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC2 values=0.45368975, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC2 values=1, perturbation=${pert} 
+adddef StokesQ:SRC2 values=0, perturbation=${pert} 
+adddef StokesU:SRC2 values=0, perturbation=${pert}
+adddef StokesV:SRC2 values=0, perturbation=${pert}
+adddef RA:SRC3  values=2.733981,   perturbation=${pert}, pert_rel=F
+adddef DEC:SRC3 values=0.45375025, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC3 values=1, perturbation=${pert} 
+adddef StokesQ:SRC3 values=0, perturbation=${pert} 
+adddef StokesU:SRC3 values=0, perturbation=${pert}
+adddef StokesV:SRC3 values=0, perturbation=${pert}
+quit
+EOF
diff --git a/CEP/BB/BBS/test/makeDemoParmDBPred b/CEP/BB/BBS/test/makeDemoParmDBPred
new file mode 100755
index 0000000000000000000000000000000000000000..0352547183b5109b0b18e46fbd733f6a4311d0fe
--- /dev/null
+++ b/CEP/BB/BBS/test/makeDemoParmDBPred
@@ -0,0 +1,57 @@
+#!/bin/sh
+
+pert="1e-6"
+
+echo "creating instrument_model (aips++)..."
+parmdb > $1_tmp.outpdb1 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_instr'
+adddef gain:11:real values=1, perturbation=${pert}
+adddef gain:12:real values=0, perturbation=${pert}
+adddef gain:21:real values=0, perturbation=${pert}
+adddef gain:22:real values=1, perturbation=${pert}
+adddef gain:11:imag values=1, perturbation=${pert}
+adddef gain:12:imag values=0, perturbation=${pert}
+adddef gain:21:imag values=0, perturbation=${pert}
+adddef gain:22:imag values=1, perturbation=${pert}
+adddef gain:11:ampl values=1, perturbation=${pert}
+adddef gain:12:ampl values=0, perturbation=${pert}
+adddef gain:21:ampl values=0, perturbation=${pert}
+adddef gain:22:ampl values=0, perturbation=${pert}
+adddef gain:11:phase values=0, perturbation=${pert}
+adddef gain:12:phase values=0, perturbation=${pert}
+adddef gain:21:phase values=0, perturbation=${pert}
+adddef gain:22:phase values=0, perturbation=${pert}
+adddef frot values=0, perturbation=${pert}
+adddef drot values=0, perturbation=${pert}
+adddef dell values=0, perturbation=${pert}
+adddef dgain:X values=1, perturbation=${pert}
+adddef dgain:Y values=1, perturbation=${pert}
+quit
+EOF
+if [ $? != 0 ]; then
+  exit 1
+fi
+
+echo "creating sky_model (aips++)..."
+parmdb > $1_tmp.outpdb2 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_sky'
+adddef RA:SRC1  values=2.734,   perturbation=${pert}, pert_rel=F
+adddef DEC:SRC1 values=0.45379, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC1 values=1, perturbation=${pert}
+adddef StokesQ:SRC1 values=0, perturbation=${pert}
+adddef StokesU:SRC1 values=0, perturbation=${pert}
+adddef StokesV:SRC1 values=0, perturbation=${pert}
+adddef RA:SRC2  values=2.73402, perturbation=${pert}, pert_rel=F
+adddef DEC:SRC2 values=0.45369, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC2 values=0.5, perturbation=${pert} 
+adddef StokesQ:SRC2 values=0, perturbation=${pert} 
+adddef StokesU:SRC2 values=0, perturbation=${pert}
+adddef StokesV:SRC2 values=0, perturbation=${pert}
+adddef RA:SRC3  values=2.73398, perturbation=${pert}, pert_rel=F
+adddef DEC:SRC3 values=0.45375, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC3 values=0.3, perturbation=${pert} 
+adddef StokesQ:SRC3 values=0, perturbation=${pert} 
+adddef StokesU:SRC3 values=0, perturbation=${pert}
+adddef StokesV:SRC3 values=0, perturbation=${pert}
+quit
+EOF
diff --git a/CEP/BB/BBSKernel/test/makeDemoParmDB b/CEP/BB/BBSKernel/test/makeDemoParmDB
new file mode 100755
index 0000000000000000000000000000000000000000..ceeda37443614e355b8a3cd71036168d8076ebd5
--- /dev/null
+++ b/CEP/BB/BBSKernel/test/makeDemoParmDB
@@ -0,0 +1,57 @@
+#!/bin/sh
+
+pert="1e-6"
+
+echo "creating instrument_model (aips++)..."
+parmdb > $1_tmp.outpdb1 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_instr'
+adddef gain:11:real values=1, perturbation=${pert}
+adddef gain:12:real values=0, perturbation=${pert}
+adddef gain:21:real values=0, perturbation=${pert}
+adddef gain:22:real values=1, perturbation=${pert}
+adddef gain:11:imag values=1, perturbation=${pert}
+adddef gain:12:imag values=0, perturbation=${pert}
+adddef gain:21:imag values=0, perturbation=${pert}
+adddef gain:22:imag values=1, perturbation=${pert}
+adddef gain:11:ampl values=1, perturbation=${pert}
+adddef gain:12:ampl values=0, perturbation=${pert}
+adddef gain:21:ampl values=0, perturbation=${pert}
+adddef gain:22:ampl values=0, perturbation=${pert}
+adddef gain:11:phase values=0, perturbation=${pert}
+adddef gain:12:phase values=0, perturbation=${pert}
+adddef gain:21:phase values=0, perturbation=${pert}
+adddef gain:22:phase values=0, perturbation=${pert}
+adddef frot values=0, perturbation=${pert}
+adddef drot values=0, perturbation=${pert}
+adddef dell values=0, perturbation=${pert}
+adddef dgain:X values=1, perturbation=${pert}
+adddef dgain:Y values=1, perturbation=${pert}
+quit
+EOF
+if [ $? != 0 ]; then
+  exit 1
+fi
+
+echo "creating sky_model (aips++)..."
+parmdb > $1_tmp.outpdb2 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_sky'
+adddef RA:SRC1  values=2.734003,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC1 values=0.4537895, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC1 values=1, perturbation=${pert}
+adddef StokesQ:SRC1 values=0, perturbation=${pert}
+adddef StokesU:SRC1 values=0, perturbation=${pert}
+adddef StokesV:SRC1 values=0, perturbation=${pert}
+adddef RA:SRC2  values=2.7340205,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC2 values=0.45368975, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC2 values=1, perturbation=${pert} 
+adddef StokesQ:SRC2 values=0, perturbation=${pert} 
+adddef StokesU:SRC2 values=0, perturbation=${pert}
+adddef StokesV:SRC2 values=0, perturbation=${pert}
+adddef RA:SRC3  values=2.733981,   perturbation=${pert}, pert_rel=F
+adddef DEC:SRC3 values=0.45375025, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC3 values=1, perturbation=${pert} 
+adddef StokesQ:SRC3 values=0, perturbation=${pert} 
+adddef StokesU:SRC3 values=0, perturbation=${pert}
+adddef StokesV:SRC3 values=0, perturbation=${pert}
+quit
+EOF
diff --git a/CEP/BB/BBSKernel/test/makeDemoParmDBCorr b/CEP/BB/BBSKernel/test/makeDemoParmDBCorr
new file mode 100755
index 0000000000000000000000000000000000000000..78a352d4f917fcb5568e5839ae56835a6a625122
--- /dev/null
+++ b/CEP/BB/BBSKernel/test/makeDemoParmDBCorr
@@ -0,0 +1,57 @@
+#!/bin/sh
+
+pert="1e-6"
+
+echo "creating instrument_model (aips++)..."
+parmdb > $1_tmp.outpdb1 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_instr'
+adddef gain:11:real values=1, perturbation=${pert}
+adddef gain:12:real values=0, perturbation=${pert}
+adddef gain:21:real values=0, perturbation=${pert}
+adddef gain:22:real values=1, perturbation=${pert}
+adddef gain:11:imag values=0, perturbation=${pert}
+adddef gain:12:imag values=0, perturbation=${pert}
+adddef gain:21:imag values=0, perturbation=${pert}
+adddef gain:22:imag values=0, perturbation=${pert}
+adddef gain:11:ampl values=1, perturbation=${pert}
+adddef gain:12:ampl values=0, perturbation=${pert}
+adddef gain:21:ampl values=0, perturbation=${pert}
+adddef gain:22:ampl values=1, perturbation=${pert}
+adddef gain:11:phase values=0, perturbation=${pert}
+adddef gain:12:phase values=0, perturbation=${pert}
+adddef gain:21:phase values=0, perturbation=${pert}
+adddef gain:22:phase values=0, perturbation=${pert}
+adddef frot values=0, perturbation=${pert}
+adddef drot values=0, perturbation=${pert}
+adddef dell values=0, perturbation=${pert}
+adddef dgain:X values=1, perturbation=${pert}
+adddef dgain:Y values=1, perturbation=${pert}
+quit
+EOF
+if [ $? != 0 ]; then
+  exit 1
+fi
+
+echo "creating sky_model (aips++)..."
+parmdb > $1_tmp.outpdb2 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_sky'
+adddef RA:SRC1  values=2.734003,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC1 values=0.4537895, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC1 values=1, perturbation=${pert}
+adddef StokesQ:SRC1 values=0, perturbation=${pert}
+adddef StokesU:SRC1 values=0, perturbation=${pert}
+adddef StokesV:SRC1 values=0, perturbation=${pert}
+adddef RA:SRC2  values=2.7340205,  perturbation=${pert}, pert_rel=F
+adddef DEC:SRC2 values=0.45368975, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC2 values=1, perturbation=${pert} 
+adddef StokesQ:SRC2 values=0, perturbation=${pert} 
+adddef StokesU:SRC2 values=0, perturbation=${pert}
+adddef StokesV:SRC2 values=0, perturbation=${pert}
+adddef RA:SRC3  values=2.733981,   perturbation=${pert}, pert_rel=F
+adddef DEC:SRC3 values=0.45375025, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC3 values=1, perturbation=${pert} 
+adddef StokesQ:SRC3 values=0, perturbation=${pert} 
+adddef StokesU:SRC3 values=0, perturbation=${pert}
+adddef StokesV:SRC3 values=0, perturbation=${pert}
+quit
+EOF
diff --git a/CEP/BB/BBSKernel/test/makeDemoParmDBPred b/CEP/BB/BBSKernel/test/makeDemoParmDBPred
new file mode 100755
index 0000000000000000000000000000000000000000..0352547183b5109b0b18e46fbd733f6a4311d0fe
--- /dev/null
+++ b/CEP/BB/BBSKernel/test/makeDemoParmDBPred
@@ -0,0 +1,57 @@
+#!/bin/sh
+
+pert="1e-6"
+
+echo "creating instrument_model (aips++)..."
+parmdb > $1_tmp.outpdb1 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_instr'
+adddef gain:11:real values=1, perturbation=${pert}
+adddef gain:12:real values=0, perturbation=${pert}
+adddef gain:21:real values=0, perturbation=${pert}
+adddef gain:22:real values=1, perturbation=${pert}
+adddef gain:11:imag values=1, perturbation=${pert}
+adddef gain:12:imag values=0, perturbation=${pert}
+adddef gain:21:imag values=0, perturbation=${pert}
+adddef gain:22:imag values=1, perturbation=${pert}
+adddef gain:11:ampl values=1, perturbation=${pert}
+adddef gain:12:ampl values=0, perturbation=${pert}
+adddef gain:21:ampl values=0, perturbation=${pert}
+adddef gain:22:ampl values=0, perturbation=${pert}
+adddef gain:11:phase values=0, perturbation=${pert}
+adddef gain:12:phase values=0, perturbation=${pert}
+adddef gain:21:phase values=0, perturbation=${pert}
+adddef gain:22:phase values=0, perturbation=${pert}
+adddef frot values=0, perturbation=${pert}
+adddef drot values=0, perturbation=${pert}
+adddef dell values=0, perturbation=${pert}
+adddef dgain:X values=1, perturbation=${pert}
+adddef dgain:Y values=1, perturbation=${pert}
+quit
+EOF
+if [ $? != 0 ]; then
+  exit 1
+fi
+
+echo "creating sky_model (aips++)..."
+parmdb > $1_tmp.outpdb2 2>&1 << EOF
+create db='$USER', dbtype='aips', tablename='$1_tmp.demo_sky'
+adddef RA:SRC1  values=2.734,   perturbation=${pert}, pert_rel=F
+adddef DEC:SRC1 values=0.45379, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC1 values=1, perturbation=${pert}
+adddef StokesQ:SRC1 values=0, perturbation=${pert}
+adddef StokesU:SRC1 values=0, perturbation=${pert}
+adddef StokesV:SRC1 values=0, perturbation=${pert}
+adddef RA:SRC2  values=2.73402, perturbation=${pert}, pert_rel=F
+adddef DEC:SRC2 values=0.45369, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC2 values=0.5, perturbation=${pert} 
+adddef StokesQ:SRC2 values=0, perturbation=${pert} 
+adddef StokesU:SRC2 values=0, perturbation=${pert}
+adddef StokesV:SRC2 values=0, perturbation=${pert}
+adddef RA:SRC3  values=2.73398, perturbation=${pert}, pert_rel=F
+adddef DEC:SRC3 values=0.45375, perturbation=${pert}, pert_rel=F
+adddef StokesI:SRC3 values=0.3, perturbation=${pert} 
+adddef StokesQ:SRC3 values=0, perturbation=${pert} 
+adddef StokesU:SRC3 values=0, perturbation=${pert}
+adddef StokesV:SRC3 values=0, perturbation=${pert}
+quit
+EOF