diff --git a/DFTFringe_Dale.pro b/DFTFringe_Dale.pro
index b0208d0f..c4307ac2 100644
--- a/DFTFringe_Dale.pro
+++ b/DFTFringe_Dale.pro
@@ -423,7 +423,7 @@ RC_FILE = DFTFringe.rc
 QMAKE_CXXFLAGS += -std=c++11
 
 # The application version
-VERSION = Dale7.3.0
+VERSION = Dale7.3.2
 
 # Define the preprocessor macro to get the application version in our application.
 DEFINES += APP_VERSION=\\\"$$VERSION\\\"
diff --git a/RevisionHistory.html b/RevisionHistory.html
index f6b877b6..2d6af3c9 100644
--- a/RevisionHistory.html
+++ b/RevisionHistory.html
@@ -425,9 +425,18 @@
 <li>Added polar plot of the astig of selected wave fronts under View menu.</li>
 <li>Changed create movie feature to add user provided prefix to each from created.</li>
 <li>Added hot keys to import interferogram</li>
-<li>Added hot key to help</li>
+<li>Added hot key help</li>
 <li>User can move mouse cursor over any profile shown,  click and drag it up or down.  Useful for comparing two work sessions results so that they match at zero height.</li>
 <li>If mouse is over the profile plot the scroll wheel can increase or decrease the y axis range.</li>
 <li>Added auto collimation setting to Ronchi and Foucault feature</li>
 <li>Remembered last ROC offset value in Ronchi and Foucault feature to remember last setting when wave front is changed to a different value</li>
 </ul>
+</ul>
+d
+<ul><li>Version 7.3.2</li>
+<ul>
+<li>Added hover and click to astig polar plot to select the wave front and show it as the currently selected wave front.  Added click to astig polar table to highlight the line on the plot and display as current wave front.</li>
+<li>Corrected test stand astig removal feature to display the correct image sizes on the report no matter what the screen resolution is.</li>
+<li>Corrected bug when subtracted wave front result is saved and then loaded so that it does not lose the fact to use or not use the null.</LI>
+<li>Updated the subtract dialog to display a reasonable size based on screen resolution.  Improved help information.</li>
+</ul>
diff --git a/astigpolargraph.cpp b/astigpolargraph.cpp
index ef88e2c8..aa880750 100644
--- a/astigpolargraph.cpp
+++ b/astigpolargraph.cpp
@@ -3,17 +3,20 @@
 #include "surfacemanager.h"
 
 
-astigPolargraph::astigPolargraph(    QList<int>list, QWidget *parent) :
-    QDialog(parent),
-    ui(new Ui::astigPolargraph), m_list(list)
+astigPolargraph::astigPolargraph(    QList<astigSample>list, QWidget *parent) :
+
+        QDialog(parent),ui(new Ui::astigPolargraph)
 {
     ui->setupUi(this);
-    
-    QPolarChart *chart = new QPolarChart();
+
+    chart = new QPolarChart();
 
     // process each wave front and place astig on the chart
-    ui->waveFrontTable->setRowCount(m_list.size());
+    ui->waveFrontTable->setRowCount(list.size());
+    ui->waveFrontTable->setSelectionBehavior(QAbstractItemView::SelectRows);
 
+    // Configure selection mode (e.g., single selection)
+    ui->waveFrontTable->setSelectionMode(QAbstractItemView::SingleSelection);
     QValueAxis *angularAxis = new QValueAxis();
     angularAxis->setTickCount(9); // First and last ticks are co-located on 0/360 angle.
     angularAxis->setLabelFormat("%.0f");
@@ -28,21 +31,30 @@ astigPolargraph::astigPolargraph(    QList<int>list, QWidget *parent) :
     double maxAstig = 1.;
 
     QVector<wavefront *>  wavefronts =SurfaceManager::get_instance()->m_wavefronts;
-
-    for(int ndx = 0; ndx < m_list.length(); ++ndx){
-        wavefront *wf = wavefronts[m_list[ndx]];
-        QString name = wf->name;
+    QScreen *screen = QGuiApplication::primaryScreen();
+    qreal screenDPI = screen->physicalDotsPerInchX();
+    int pensize = 5 * screenDPI/256.;  // adjust pen size to screen resolution.  256 is DPI of my 4K 17 inch laptop
+    if (pensize < 1)
+        pensize = 1;
+    for(int ndx = 0; ndx < list.length(); ++ndx){
+
+        QString name = list[ndx].m_name;
         int slashndx = name.lastIndexOf('/');
         QString shortName = name.mid(name.lastIndexOf('/',slashndx-1));
         QTableWidgetItem *item = new QTableWidgetItem(shortName, 0);
         ui->waveFrontTable->setItem(ndx,0,item);
-        double xastig = wf->InputZerns[4];
-        double yastig = wf->InputZerns[5];
+        double xastig = list[ndx].m_xastig;
+        double yastig = list[ndx].m_yastig;
         double mag = sqrt(xastig * xastig + yastig * yastig);
         if (mag > maxAstig) maxAstig = mag;
 
         double angle = (atan2(yastig,xastig)/2.) * 180./M_PI;
+
         angle = 90 - angle;
+
+        astigSample sample(name, angle, mag);
+
+        m_list << sample;
         QScatterSeries *series = new QScatterSeries();
 
         int lastndx = name.lastIndexOf('/');
@@ -61,9 +73,12 @@ astigPolargraph::astigPolargraph(    QList<int>list, QWidget *parent) :
         chart->addSeries(line);
         line->attachAxis(radialAxis);
         line->attachAxis(angularAxis);
+        line->setName(name);
+            connect(line, &QLineSeries::hovered, this, &astigPolargraph::tooltip);
+            connect(line, &QLineSeries::clicked, this, &astigPolargraph::clicked);
         chart->legend()->markers(line)[0]->setVisible(false);
 
-        line->setPen(QPen(series->brush(),5));
+        line->setPen(QPen(series->brush(),pensize));
 
         QTableWidgetItem *pv = new QTableWidgetItem(QString().number(mag), 0);
         item->setForeground(series->brush());
@@ -73,18 +88,110 @@ astigPolargraph::astigPolargraph(    QList<int>list, QWidget *parent) :
     }
 
     chart->setTitle("Magnitude and axis of high edge");
-    if (m_list.length() > 4)
+    if (list.length() > 4)
         chart->legend()->setAlignment(Qt::AlignRight);
     else chart->legend()->setAlignment(Qt::AlignBottom);
 
     maxAstig = ceil(maxAstig);
     radialAxis->setRange(0, maxAstig);
     angularAxis->setRange(0, 360);
-
+    ui->waveFrontTable->resizeColumnsToContents();
     ui->polarChart->setChart(chart);
+
 }
+void astigPolargraph::tooltip(QPointF point, bool state)
+{
+    if (state) {
+        findClosestPoint(point);
+    }
 
+}
+void astigPolargraph::clicked(QPointF point)
+{
+
+   findClosestPoint(point);
+
+}
 astigPolargraph::~astigPolargraph()
 {
     delete ui;
+    delete chart;
 }
+
+void astigPolargraph::hideHoverHelp(){
+    ui->hoverText->hide();
+}
+int astigPolargraph::findClosestPoint(const QPointF clickedPoint){
+
+    QPointF closest(INT_MAX, INT_MAX);
+    qreal distance(INT_MAX);
+    int closeNdx = -1;
+    int ndx = 0;
+    for (auto sample : m_list) {
+        QPointF currentPoint(sample.m_xastig, sample.m_yastig);
+        qreal currentDistance = qSqrt((currentPoint.x() - clickedPoint.x())
+                                      * (currentPoint.x() - clickedPoint.x())
+                                      + (currentPoint.y() - clickedPoint.y())
+                                      * (currentPoint.y() - clickedPoint.y()));
+
+        if (currentDistance < distance) {
+            distance = currentDistance;
+            closest = currentPoint;
+            closeNdx = ndx;
+
+        }
+        ++ndx;
+    }
+    QString name = m_list[closeNdx].m_name;
+    emit waveSeleted(name);
+    int slashndx = name.lastIndexOf('/');
+    QString shortName = name.mid(name.lastIndexOf('/',slashndx-1));
+
+    QList<QTableWidgetItem*> items = ui->waveFrontTable->findItems(shortName, Qt::MatchEndsWith);
+
+
+    if (items.length() > 0){
+
+        ui->waveFrontTable->selectRow(items[0]->row());
+    }
+    return closeNdx;
+}
+
+void astigPolargraph::on_waveFrontTable_itemClicked(QTableWidgetItem *item)
+{
+    QString name = ui->waveFrontTable->item(item->row(),0)->text();
+    emit waveSeleted(name);
+
+    int lastndx = name.lastIndexOf('/');
+    if (lastndx != -1)
+        name = name.mid(lastndx);
+
+
+    int seriesCount = chart->series().count();
+
+    for (int i = 0; i < seriesCount; ++i) {
+      QAbstractSeries* series = chart->series().at(i);
+      if (series) {
+
+           if (series->type()== QAbstractSeries::SeriesTypeLine){
+              if (series->name() == name){
+                  QLineSeries *line = static_cast<QLineSeries*>(series);
+                  // Create a pen object to get the current pen attributes.
+                  QPen pen = line->pen();
+                  if (pen.style() == Qt::DotLine){
+                      pen.setStyle(Qt::SolidLine);
+                      pen.setWidth(pen.width()/4);
+                  }
+                  else {
+                    pen.setStyle(Qt::DotLine);
+                    pen.setWidth(pen.width() * 4);
+                  }
+
+                  line->setPen(pen);
+              }
+           }
+
+      }
+    }
+}
+
diff --git a/astigpolargraph.h b/astigpolargraph.h
index ea8033aa..9a1e7013 100644
--- a/astigpolargraph.h
+++ b/astigpolargraph.h
@@ -9,22 +9,43 @@
 #include <QtCharts/QValueAxis>
 #include <QtCharts/QPolarChart>
 #include "wavefront.h"
+#include <QTableWidgetItem>
 
+QT_CHARTS_USE_NAMESPACE
 namespace Ui {
 class astigPolargraph;
 }
+class astigSample{
+
+public:
+    QString m_name;
+    double m_xastig;
+    double m_yastig;
+    astigSample(QString name, double xastig, double yastig): m_name(name), m_xastig(xastig), m_yastig(yastig){};
+};
 
 class astigPolargraph : public QDialog
 {
     Q_OBJECT
 
 public:
-    explicit astigPolargraph(QList<int> list,QWidget *parent = nullptr);
+    explicit astigPolargraph(QList<astigSample> list,QWidget *parent = nullptr);
     ~astigPolargraph();
 
+signals:
+    void waveSeleted(QString);
+private slots:
+    void tooltip(QPointF point, bool state);
+    void clicked(QPointF point);
+    void on_waveFrontTable_itemClicked(QTableWidgetItem *item);
+
 private:
+    int findClosestPoint(const QPointF clickedPoint );
     Ui::astigPolargraph *ui;
-    QList<int> m_list;  // list index of selected wave fronts in surface manager's list.
+    QPolarChart *chart;
+    QList<astigSample> m_list;
+public:
+        void hideHoverHelp();
 };
 
 #endif // ASTIGPOLARGRAPH_H
diff --git a/astigpolargraph.ui b/astigpolargraph.ui
index 566f4f73..7cddb452 100644
--- a/astigpolargraph.ui
+++ b/astigpolargraph.ui
@@ -16,6 +16,30 @@
   <layout class="QVBoxLayout" name="verticalLayout">
    <item>
     <layout class="QHBoxLayout" name="horizontalLayout">
+     <item>
+      <spacer name="horizontalSpacer">
+       <property name="orientation">
+        <enum>Qt::Horizontal</enum>
+       </property>
+       <property name="sizeHint" stdset="0">
+        <size>
+         <width>40</width>
+         <height>20</height>
+        </size>
+       </property>
+      </spacer>
+     </item>
+     <item>
+      <widget class="QLabel" name="hoverText">
+       <property name="text">
+        <string>Hover or left click over right side end point to highlight associated wave front.</string>
+       </property>
+      </widget>
+     </item>
+    </layout>
+   </item>
+   <item>
+    <layout class="QHBoxLayout" name="horizontalLayout_2">
      <item>
       <widget class="QTableWidget" name="waveFrontTable">
        <property name="editTriggers">
@@ -55,7 +79,7 @@
       <enum>Qt::Horizontal</enum>
      </property>
      <property name="standardButtons">
-      <set>QDialogButtonBox::Cancel|QDialogButtonBox::Ok</set>
+      <set>QDialogButtonBox::NoButton</set>
      </property>
     </widget>
    </item>
diff --git a/mainwindow.cpp b/mainwindow.cpp
index 7efc47eb..2d5d60af 100644
--- a/mainwindow.cpp
+++ b/mainwindow.cpp
@@ -2074,8 +2074,21 @@ void MainWindow::on_actionastig_in_polar_triggered()
 {
     surfaceAnalysisTools *saTools = surfaceAnalysisTools::get_Instance();
     QList<int> list = saTools->SelectedWaveFronts();
-    astigPolargraph * graph = new astigPolargraph( list, this);
-    graph->resize(2000,1000);
+
+    QList<astigSample> samples;
+    foreach(int i, list){
+        wavefront *wf = m_surfaceManager->m_wavefronts[i];
+        astigSample sample(wf->name, wf->InputZerns[4], wf->InputZerns[5]);
+        samples << sample;
+    }
+    astigPolargraph * graph = new astigPolargraph(samples);
+    connect(graph, SIGNAL(waveSeleted(QString)), m_surfaceManager, SLOT(wavefrontDClicked(QString)));
+    QScreen *screen = QGuiApplication::primaryScreen();
+    QSizeF physicalSize = screen->availableSize();
+    graph->resize(physicalSize.width()/2,physicalSize.height()/2);
+
     graph->exec();
+    graph->disconnect();
+    delete graph;
 }
 
diff --git a/standastigwizard.cpp b/standastigwizard.cpp
index 6c617378..adecf2dc 100644
--- a/standastigwizard.cpp
+++ b/standastigwizard.cpp
@@ -42,8 +42,12 @@ standAstigWizard::standAstigWizard(SurfaceManager *sm, QWidget *parent, Qt::Wind
     setWindowTitle(tr("Stand Astig analysis"));
     setPixmap(QWizard::WatermarkPixmap, QPixmap(":/res/wats2.png"));
     setPixmap(QWizard::LogoPixmap, QPixmap(":/res/mirror_stand.png"));
-    //AstigReportTitle = mirrorDlg::get_Instance()->m_name;
-    //resize( QSize(600, 489));
+
+
+    QScreen *screen = QGuiApplication::primaryScreen();
+
+    resize(screen->availableSize().width() * .5 ,screen->availableSize().height() * .5);
+
 }
 
 standAstigWizard::~standAstigWizard()
@@ -122,6 +126,7 @@ makeAverages::makeAverages(QWidget *parent)
                                    );
     txt->setReadOnly(true);
     QVBoxLayout *layout = new QVBoxLayout;
+
     layout->addWidget(txt);
     setLayout(layout);
 }
@@ -227,9 +232,15 @@ define_input::define_input(QWidget *parent)
     QPushButton *browsePath = new QPushButton("...");
     CWRb = new QRadioButton("Rotate CW");
     CCWRb = new QRadioButton("Rotate CCW");
-    CCWRb->setChecked(true);
+    if (set.value("stand astig ccw", true).toBool())
+        CCWRb->setChecked(true);
+    else
+        CWRb->setChecked(true);
+
     connect(browsePath, SIGNAL(pressed()), this, SLOT(setBasePath()));
 
+
+
     browsePath->setStyleSheet("");
     QGridLayout *l = new QGridLayout();
 
@@ -246,6 +257,12 @@ define_input::define_input(QWidget *parent)
     listDisplay = new QListWidget();
     listDisplay->setSelectionMode( QAbstractItemView::MultiSelection);
     listDisplay->setContextMenuPolicy(Qt::CustomContextMenu);
+//    int size = set.beginReadArray("stand astig removal files");
+//    for (int i = 0; i < size; ++i) {
+//        set.setArrayIndex(i);
+//        listDisplay->addItem(set.value("item").toString());
+//    }
+    set.endArray();
     connect(listDisplay, SIGNAL(customContextMenuRequested(QPoint)), this,
             SLOT(showContextMenu(QPoint)));
 
@@ -260,6 +277,7 @@ define_input::define_input(QWidget *parent)
     m_log = new QTextEdit();
     m_log->append("Ready to add wave fronts to process.");
     l->addWidget(m_log,21,0,1,0);
+
     setLayout(l);
 
 }
@@ -270,7 +288,15 @@ void define_input::compute(){
         QMessageBox::warning(0, "", "You must first add at least 2 wavefront files to the list.");
         return;
     }
-
+    QSettings set;
+    set.setValue("stand astig ccw",CCWRb->isChecked() );
+    set.beginWriteArray("stand astig removal files");
+    // save the listDisplay of files to process
+    for (int i = 0; i < listDisplay->count(); ++i) {
+        set.setArrayIndex(i);
+        set.setValue("item", listDisplay->item(i)->text());
+    }
+    set.endArray();
     AstigReportTitle = title->text();
     runpb->setText("Working");
     runpb->setEnabled(false);
diff --git a/subtractwavefronatsdlg.ui b/subtractwavefronatsdlg.ui
index 171959ac..86274b58 100644
--- a/subtractwavefronatsdlg.ui
+++ b/subtractwavefronatsdlg.ui
@@ -17,7 +17,7 @@
    <item>
     <widget class="QLabel" name="label">
      <property name="text">
-      <string>&lt;html&gt;&lt;head/&gt;&lt;body&gt;&lt;p&gt;Select wavefront to be subtracted from the currently viewed wavefront. Remember to disable the software null in the mirror configuration in order to view the result properly.&lt;/p&gt;&lt;p&gt;When you subtract two wavefronts that both needed to have SA software null subtracted from them because the interferogram that made them was taken at ROC of the mirror under test, the result will have not have the original SA in it.&lt;/p&gt;&lt;p&gt;As a result subtracting the artificial null will artificially change the SA of the result. Usually you don't want that to happen.&lt;/p&gt;&lt;p&gt;So in this case do not check this check box.&lt;/p&gt;&lt;p&gt;If you are using the subtraction operation to remove system errors for example then you want to have this check box checked.&lt;/p&gt;&lt;p&gt;&lt;br/&gt;&lt;/p&gt;&lt;/body&gt;&lt;/html&gt;</string>
+      <string>&lt;html&gt;&lt;head/&gt;&lt;body&gt;&lt;p&gt;Select wavefront to be subtracted from the currently viewed wavefront. &lt;/p&gt;&lt;p&gt;&lt;span style=&quot; font-weight:700;&quot;&gt;When to check the &amp;quot;Use softwar Null after subtraction&amp;quot;.&lt;/span&gt;&lt;/p&gt;&lt;p&gt;When you subtract two wavefronts that both needed to have SA software null subtracted from them because for example the interferogram that made them was taken at ROC of the mirror under test, the result will have not have the original SA in it. Usually that is what you want. You just want to see the difference. In that case do not check the box.&lt;/p&gt;&lt;p&gt;If you are using the subtraction operation to remove system errors, for example if you have a test stand only induced wave front you want to subtract from the selected wave front  then you want to &lt;span style=&quot; font-weight:700;&quot;&gt;check that box&lt;/span&gt;.&lt;/p&gt;&lt;p&gt;&lt;br/&gt;&lt;/p&gt;&lt;/body&gt;&lt;/html&gt;</string>
      </property>
      <property name="wordWrap">
       <bool>true</bool>
diff --git a/surfacemanager.cpp b/surfacemanager.cpp
index bfbf9eb6..b186dccc 100644
--- a/surfacemanager.cpp
+++ b/surfacemanager.cpp
@@ -72,6 +72,7 @@
 #include "oglrendered.h"
 #include "ui_oglrendered.h"
 #include "spdlog/spdlog.h"
+#include "astigpolargraph.h"
 
 cv::Mat theMask;
 cv::Mat deb;
@@ -948,6 +949,9 @@ void SurfaceManager::writeWavefront(QString fname, wavefront *wf, bool saveNulle
         file << "DIAM " << wf->diameter << std::endl;
         file << "ROC " << wf->roc << std::endl;
         file << "Lambda " << wf->lambda << std::endl;
+        if (!wf->useSANull){
+            file << "Do Not use null" <<  1 << std::endl;
+        }
         mirrorDlg &md = *mirrorDlg::get_Instance();
         if (md.isEllipse()){
             file << "ellipse_vertical_axis " << md.m_verticalAxis;
@@ -1145,7 +1149,7 @@ wavefront * SurfaceManager::readWaveFront(QString fileName){
             md->m_outlineShape = ELLIPSE;
             iss >> dummy >> md->m_verticalAxis;
         }
-        if (l.startsWith("nulled")){
+        if (l.startsWith("Do Not use null")){
             wf->useSANull = false;
         }
     }
@@ -1559,7 +1563,7 @@ void SurfaceManager::average(QList<wavefront *> wfList){
 
 
     QApplication::setOverrideCursor(Qt::WaitCursor);
-qDebug() << "maxkey" << maxkey << rrows << rcols << sizes[maxkey];
+
     cv::Mat mask = wfList[sizes[maxkey][0]]->workMask.clone();
     if (mask.cols != rcols || mask.rows != rrows){
         cv::resize(mask,mask,Size(rrows,rcols));
@@ -1780,10 +1784,15 @@ void SurfaceManager::subtract(wavefront *wf1, wavefront *wf2, bool use_null){
 
 void SurfaceManager::subtractWavefronts(){
     QList<QString> list;
+    QList<int> doThese =  m_surfaceTools->SelectedWaveFronts();
     for (int i = 0; i < m_wavefronts.size(); ++i){
-        list.append(m_wavefronts[i]->name);
+        if (!m_wavefronts[i]->name.contains(doThese[0]))
+             list.append(m_wavefronts[i]->name);
     }
     subtractWavefronatsDlg dlg(list);
+    QScreen *screen = QGuiApplication::primaryScreen();
+    QSizeF physicalSize = screen->availableSize();
+    dlg.resize(physicalSize.width()/2,physicalSize.height()/2);
     if (dlg.exec()){
         int ndx2 = dlg.getSelected();
         if (ndx2 == -1){
@@ -1909,31 +1918,41 @@ void plotlineThruAstigs(QwtPlot *plt, cv::Mat xastig, cv::Mat yastig, QList<rota
 
 }
 
+int getImageSize(QPrinter &printer){
+    QScreen *screen = QGuiApplication::primaryScreen();
+
+    qreal screenDPI = screen->physicalDotsPerInchX();
+
+    int Width;
+
+    printer.setPageSize(QPageSize(printer.pageLayout().pageSize()));
+    QPainter painterPDF( &printer );
+
+    // I want to create images that are about 1/2 as wide as the printer paper.
+    // So get the printer paper size
+    // THen get the screen resolution
+    // then calculate how many pixels create that size image on the screen
+
+    // get the size of the printed page in inches.
+    QSize printerSizePixels = printer.pageLayout().paintRectPixels(printer.resolution()).size();
+    qreal pageWidthInches = printerSizePixels.width()/printer.resolution();
+    qreal imageWidthInches= .8 * pageWidthInches/2;
+
+    Width = imageWidthInches * screenDPI;
+    return Width;
+}
 // calculate stand astig for each input
 // for each input rotate the average by the input angle and subtract it from the input
 // plot the astig of each of the inputs which will be the stand only astig.
 // on input list is the list of wavefront files and thier rotation angle.
 //     inputs are the actual wavefronts at original rotations.
 //     avgNdx is the index in the m_wavefronts list of the average with stand removed.
-textres SurfaceManager::Phase2(QList<rotationDef *> list, QList<wavefront *> inputs, int avgNdx ){
+textres SurfaceManager::Phase2(QList<rotationDef *> list, QList<wavefront *> inputs, int avgNdx , int Width, QPrinter &printer){
     QTextEdit *editor = new QTextEdit;
 
     QTextDocument *doc = editor->document();
     textres results;
-    results.Edit = editor;
-    const int Width = 800 * .6;
-    const int Height = 600 * .6;
-    QImage contour(Width ,Height, QImage::Format_ARGB32 );
 
-    QPrinter printer(QPrinter::HighResolution);
-    printer.setColorMode( QPrinter::Color );
-    printer.setFullPage( true );
-    printer.setOutputFileName( "stand.pdf" );
-    printer.setOutputFormat( QPrinter::PdfFormat );
-    //printer.setResolution(85);
-    printer.setPageSize(QPageSize(QPageSize::A4));
-
-    QPainter painterPDF( &printer );
     cv::Mat xastig = cv::Mat::zeros(list.size(), 1, numType);
     cv::Mat yastig = cv::Mat::zeros(list.size(), 1, numType);
     cv::Mat standxastig = cv::Mat::zeros(list.size(), 1, numType);
@@ -1941,8 +1960,15 @@ textres SurfaceManager::Phase2(QList<rotationDef *> list, QList<wavefront *> inp
     cv::Mat standastig = cv::Mat::zeros(list.size(), 1, numType);
     QVector<cv::Mat> standwfs;
     QVector<double> astigMag;
-    editor->resize(printer.pageLayout().paintRectPixels(printer.resolution()).size());
+
+    int Height = Width;
+
+    QImage contour(Width , Height, QImage::Format_ARGB32 );
+
+    results.Edit = editor;
+    editor->resize(Width * 1.2, Height);
     doc->setPageSize(printer.pageLayout().paintRectPixels(printer.resolution()).size());
+
     cv::Mat standavg = cv::Mat::zeros(inputs[0]->workData.size(), numType);
     cv::Mat standavgZernMat = cv::Mat::zeros(inputs[0]->workData.size(), numType);
     double mirrorXaverage = 0;
@@ -2021,7 +2047,7 @@ textres SurfaceManager::Phase2(QList<rotationDef *> list, QList<wavefront *> inp
         }
         fittedcircle2 = Circle(xmean/xastig.rows, ymean/xastig.rows, avgRadius/xastig.rows );
         fittedcircle = fittedcircle2;
-qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
+
     fittedcircle = fittedcircle1;
     if (xastig.rows == 2 || fittedcircle1.r > 1.5 * fittedcircle2.r ){
         fittedcircle = fittedcircle2;
@@ -2118,6 +2144,7 @@ qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
     QPolygonF mirrorAstigAtEachRotation;
     double mirrorAstigRadius =0;
     plotlineThruAstigs(pl1, xastig, yastig, list);
+    QList<astigSample> samples;
     for (int i = 0; i < list.size(); ++i){
         if (cnt++ == 0)
             imagesHtml.append("<tr>");
@@ -2166,15 +2193,21 @@ qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
         wf->name = QString("%1").arg(list[i]->angle, 6, 'f', 2, QLatin1Char('0'));
 
         cp->setSurface(wf);
-        cp->resize(Width,Height);
+        cp->resize(Width,.8 * Height);
         cp->replot();
+        astigSample sample(wf->name, wf->InputZerns[4], wf->InputZerns[5]);
+        samples << sample;
+        QSize s = cp->size();
         contour.fill( QColor( Qt::white ).rgb() );
-        renderer.render( cp, &painter, QRect(0,0,Width,Height) );
+        renderer.render( cp, &painter, QRect(0,0,s.width(),s.height() ));
+
         QString imageName = QString("mydata://zern%1.png").arg(wf->name);
         imageName.replace("-","CCW");
         doc->addResource(QTextDocument::ImageResource,  QUrl(imageName), QVariant(contour));
         results.res.append (imageName);
-        imagesHtml.append("<td><p> <img src='" +imageName + "' /></p></td>");
+
+        imagesHtml.append(QString("<td style=\"text-align: center;\"> ") + wf->name +  "<br><img src='" +imageName + "'><br></td>");
+
         if (cnt == 2){
             cnt = 0;
             imagesHtml.append("</tr>");
@@ -2200,20 +2233,20 @@ qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
     minMaxIdx(wf2->data, &mmin,&mmax);
     wf2->min = mmin;
     wf2->max = mmax;
-    //wf2->workMask = m_wavefronts[0]->workMask.clone();
+
     wf2->name = QString("Average Stand zernike based");
     ContourPlot *cp1 = new ContourPlot();
-    //cp1->m_zRangeMode = "Min/Max";
+
     cp1->setSurface(wf2);
-    cp1->resize(Width, Height);
+    cp1->resize(Width, .8 * Height);
     cp1->replot();
-    QImage contour2(500, 500, QImage::Format_ARGB32 );
+    QImage contour2(Width, Width, QImage::Format_ARGB32 );
     contour2.fill( QColor( Qt::white ).rgb() );
     QPainter painter2( &contour2 );
 
     renderer.setDiscardFlag(QwtPlotRenderer::DiscardLegend, false);
     renderer.render( cp1, &painter2, QRect(0,0,Width,Height) );
-    QString imageName = "mydata://StandCotourZerns.png";
+    QString imageName = "mydata://StandContourZerns.png";
     doc->addResource(QTextDocument::ImageResource,  QUrl(imageName), QVariant(contour2));
 
     wf2->data = wf2->workData = standavg;
@@ -2232,7 +2265,7 @@ qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
     contour2.fill( QColor( Qt::white ).rgb() );
     renderer.render( cp1, &painter2, QRect(0,0,Width,Height) );
     cp1->m_zRangeMode = "Auto"; // restore contour plot to auto scaling.
-    imageName = QString("mydata://StandCotourMat.png");
+    imageName = QString("mydata://StandContourMat.png");
     doc->addResource(QTextDocument::ImageResource,  QUrl(imageName), QVariant(contour2));
 
 
@@ -2271,11 +2304,6 @@ qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
     mirrorXastig /= count;
     mirrorYastig /= count;
 
-    //QwtPlotMarker *mirrorMeanMark = new QwtPlotMarker();
-    //mirrorMeanMark->setValue(QPointF(bestFitAstig.a, bestFitAstig.b));
-    //mirrorMeanMark->setSymbol(new QwtSymbol(QwtSymbol::Triangle,QColor(0,255,0,40), QColor(0,0,0), QSize(30,30)));
-    //mirrorMeanMark->attach(pl1);
-
     QPolygonF stdCircle;
     double SE = standStd[0]/sqrt(standastig.rows);
     for (double rho = 0; rho <= 2 * M_PI; rho += M_PI/32.){
@@ -2309,15 +2337,16 @@ qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
     curveAvgMirror->setPen(pen);
     curveAvgMirror->setSamples(stdCircle);
     curveAvgMirror->attach(pl1);
-    pl1->resize(500,500);
+    pl1->resize(1.5 * Width,1.5 * Width);
     pl1->replot();
-    QwtPlotRenderer renderer3;
-    QImage contour3(500, 500, QImage::Format_ARGB32 );
+
+    QImage contour3(1.5 * Width, 1.5 * Width , QImage::Format_ARGB32 );
     contour3.fill( QColor( Qt::white ).rgb() );
     QPainter painter3( &contour3 );
 
     renderer.setDiscardFlag(QwtPlotRenderer::DiscardLegend, false);
-    renderer.render( pl1, &painter3, QRect(0,0,500,500) );
+    renderer.render( pl1, &painter3, QRect(0,0,1.5  * Width,1.5 * Width) );
+
     imageName = QString("mydata://plot.png");
     doc->addResource(QTextDocument::ImageResource,  QUrl(imageName), QVariant(contour3));
     results.res.append (imageName);
@@ -2336,18 +2365,38 @@ qDebug() << "circle fit"<< avgRadius << fittedcircle1.r << fittedcircle2.r;
                 "The stand removal was good.  Idealy the STD (standard deviation) should be"
                 " less than .1 which means less than .1 wave pv on the surface of the mirror"
                 "<br>The colored plus signs are what is calculated for test stand induced astig."
-                "at each rotation angle.</p><br>"
+                "at each rotation angle.</p><br>The plot below is a polar plot of each rotations test stand astig."
+                "The ideal would be they are all the same magnitude and angle. If any one is much different than the rest then the astig removal may not be good.");
+            astigPolargraph *polar = new astigPolargraph(samples);
+            polar->hideHoverHelp();
+            polar->setWindowFlags(Qt::Dialog | Qt::CustomizeWindowHint | Qt::WindowTitleHint);
+            polar->resize(Width * 2, 2 *  Width);
+            QImage polarImage(2 * Width, 2 * Width , QImage::Format_ARGB32 );
+            polarImage.fill( QColor( Qt::white ).rgb() );
+            QPainter painterPolar( &polarImage );
+
+
+            polar->render(  &painterPolar);
+            delete polar;
+            imageName = QString("mydata://polarplot.png");
+            doc->addResource(QTextDocument::ImageResource,  QUrl(imageName), QVariant(polarImage));
+            results.res.append (imageName);
+            html.append("<p> <img src='" + imageName + "' /></p>");
+
+                html.append("<table><tr><td><img src='mydata://StandContourZerns.png' /></td><td><img src='mydata://StandContourMat.png' /></td></tr></table>"
+                "<p style=\"page-break-before:always\"; font-size:12pt>The contours above are the average system induced forces derived from the average of all rotations.<br>"
+                "The left contour is based on the zernike values and the contour on the right is based on the wavefront.</p>");
 
 
-                "<table><tr><td><img src='mydata://StandCotourZerns.png' /></td><td><img src='mydata://StandCotourMat.png' /></td></tr></table>"
-                "<p font-size:12pt>The contours above are the average system induced forces derived from the average of all rotations.<br>"
-                "The left contour is based on the zernike values and the contour on the right is based on the wavefront.</p>"
-                "<p font-size:12pt>The contour plots below are of what is beleived to be test stand only induced errors at each rotation. "
+
+    html.append( "<p style=\"page-break-before:always\"; font-size:12pt>The contour plots below are of what is beleived to be test stand only induced errors at each rotation. "
                 "Check that they are similar at each rotation.  If not then maybe "
                 "stand (system) induced error is not same at each rotation then the stand removal is not reliable. "
                 "However it is unlikely that they will all look exactly the same.</p>");
 
 
+
+
     html.append(imagesHtml);
     html.append("</body></html>");
     editor->setHtml(html);
@@ -2392,21 +2441,19 @@ void SurfaceManager::computeStandAstig(define_input *wizPage, QList<rotationDef
             break;
         }
     }
-        QApplication::setOverrideCursor(Qt::WaitCursor);
+    QApplication::setOverrideCursor(Qt::WaitCursor);
     QPrinter printer(QPrinter::HighResolution);
     printer.setColorMode( QPrinter::Color );
     printer.setFullPage( true );
     printer.setOutputFileName( "stand.pdf" );
     printer.setOutputFormat( QPrinter::PdfFormat );
-    printer.setResolution(85);
-    printer.setPageSize(QPageSize(QPageSize::A4));
 
+    int Width = getImageSize(printer);
 
     QTextEdit *editor = new QTextEdit;
     editor->resize(printer.pageLayout().paintRectPixels(printer.resolution()).size());
-    const int contourWidth = 2 * 340/3;
-    const int contourHeight = 2 * 360/3;
-    QImage contour = QImage(contourWidth ,contourHeight, QImage::Format_ARGB32 );
+
+    QImage contour = QImage(Width ,Width, QImage::Format_ARGB32 );
 
 
     QwtPlotRenderer renderer;
@@ -2455,8 +2502,9 @@ void SurfaceManager::computeStandAstig(define_input *wizPage, QList<rotationDef
         inputs.append(wf);
         unrotatedNdxs.append(m_currentNdx);
         plot->setSurface(wf);
+        plot->resize(Width, .8 * Width);
         plot->replot();
-        renderer.render( plot, &painter, QRect(0,0,contourWidth,contourHeight) );
+        renderer.render( plot, &painter, QRect(0,0,Width,.8 * Width) );
 
         QString imageName = QString("mydata://%1.png").arg(list[i]->fname);
         QString angle = QString("%1 Deg").arg(-list[i]->angle, 6, 'f', 2);
@@ -2472,18 +2520,17 @@ void SurfaceManager::computeStandAstig(define_input *wizPage, QList<rotationDef
         l.append(ndx);
         ndx = m_wavefronts.size();
 
-
-        qDebug() << "Counter rotating " << list[i]->fname << " angle " << list[i]->angle;
         rotateThese(wrapAngle(list[i]->angle),l);
         rotated.append(m_wavefronts[ndx]);
         wf = m_wavefronts[ndx];
 
         loadComplete();
         plot->setSurface(wf);
+        plot->resize(Width, .8 * Width);
         plot->replot();
 
         contour.fill( QColor( Qt::white ).rgb() );
-        renderer.render( plot, &painter, QRect(0,0,contourWidth,contourHeight) );
+        renderer.render( plot, &painter, QRect(0,0,Width, .8 * Width) );
         angle = QString("%1 Deg").arg(-list[i]->angle, 6, 'f', 2);
         imageName = QString("mydata://CR%1%2.png").arg(list[i]->fname).arg(angle); // clazy:exclude=qstring-arg
 
@@ -2525,11 +2572,11 @@ void SurfaceManager::computeStandAstig(define_input *wizPage, QList<rotationDef
     ContourPlot *plotAvg =new ContourPlot(0,0,false);//m_contourPlot;
 
     plotAvg->setSurface(m_wavefronts[m_currentNdx]);
-    contour = QImage(550,450, QImage::Format_ARGB32 );
+    contour = QImage(1.5 * Width,1.5 * Width, QImage::Format_ARGB32 );
     contour.fill( QColor( Qt::white ).rgb() );
     QPainter painter( &contour );
-
-    renderer.render( plotAvg, &painter, QRect(0,0,550,450) );
+    plotAvg->resize(1.5 * Width, 1.5 * .8 * Width);
+    renderer.render( plotAvg, &painter, QRect(0,0,1.5 * Width, 1.5 *  Width) );
 
     QString imageName = "mydata://AvgAstigremoved.png";
     doc2->addResource(QTextDocument::ImageResource,  QUrl(imageName), QVariant(contour));
@@ -2540,6 +2587,7 @@ void SurfaceManager::computeStandAstig(define_input *wizPage, QList<rotationDef
     page2->setHtml(html);
     //page2->show();
 
+
     /**************************************************************************************************/
     // PHASE 2
     // calculate stand astig for each input
@@ -2548,13 +2596,18 @@ void SurfaceManager::computeStandAstig(define_input *wizPage, QList<rotationDef
 
     wizPage->m_log->setText("computing stand astigs");
     QApplication::processEvents();
-    textres page3res = Phase2(list, inputs, avgNdx);
+    textres page3res = Phase2(list, inputs, avgNdx, Width, printer);
     QTabWidget *tabw = new QTabWidget();
     tabw->setTabShape(QTabWidget::Triangular);
     tabw->addTab(editor, "Page 1 input analysis");
     tabw->addTab(page2, "Page 2 Stand removed.");
     tabw->addTab(page3res.Edit, "Page 3 stand analysis");
-    tabw->resize(800,600);
+    QScreen *screen = QGuiApplication::primaryScreen();
+    QSizeF physicalSize = screen->availableSize();
+
+
+    //tabw->resize(Width,2 * Width);
+    tabw->setGeometry(Width, 200, 2 * Width, physicalSize.height()* .7);
     tabw->show();
 
     // build pdf doc from the three textedits
@@ -2562,10 +2615,10 @@ void SurfaceManager::computeStandAstig(define_input *wizPage, QList<rotationDef
     foreach( QString res, page3res.res){
         pdfDoc.addResource(QTextDocument::ImageResource, res, page3res.Edit->document()->resource(QTextDocument::ImageResource,res));
     }
-    pdfDoc.addResource(QTextDocument::ImageResource, QString("mydata://StandCotourZerns.png"),
-                       page3res.Edit->document()->resource(QTextDocument::ImageResource,QString("mydata://StandCotourZerns.png") ));
-    pdfDoc.addResource(QTextDocument::ImageResource, QString("mydata://StandCotourMat.png"),
-                       page3res.Edit->document()->resource(QTextDocument::ImageResource,QString("mydata://StandCotourMat.png") ));
+    pdfDoc.addResource(QTextDocument::ImageResource, QString("mydata://StandContourZerns.png"),
+                       page3res.Edit->document()->resource(QTextDocument::ImageResource,QString("mydata://StandContourZerns.png") ));
+    pdfDoc.addResource(QTextDocument::ImageResource, QString("mydata://StandContourMat.png"),
+                       page3res.Edit->document()->resource(QTextDocument::ImageResource,QString("mydata://StandContourMat.png") ));
     foreach(QString res, doc1Res){
         pdfDoc.addResource(QTextDocument::ImageResource, res, editor->document()->resource(QTextDocument::ImageResource,res));
     }
@@ -2606,6 +2659,7 @@ void SurfaceManager::computeTestStandAstig(){
         spdlog::get("logger")->trace("new standAstigWizard");
         m_standAstigWizard = new standAstigWizard(this);
         m_standAstigWizard->setAttribute( Qt::WA_DeleteOnClose, true );
+
         m_standAstigWizard->show();
     }
     else{
diff --git a/surfacemanager.h b/surfacemanager.h
index 7e064f79..c1182160 100644
--- a/surfacemanager.h
+++ b/surfacemanager.h
@@ -135,7 +135,7 @@ class SurfaceManager : public QObject
 
     explicit SurfaceManager(QObject *parent=0, surfaceAnalysisTools *tools = 0, ProfilePlot *profilePlot =0,
                    contourView *contourView = 0, SurfaceGraph *glPlot = 0, metricsDisplay *mets = 0);
-    textres Phase2(QList<rotationDef *> list, QList<wavefront *> inputs, int avgNdx);
+    textres Phase2(QList<rotationDef *> list, QList<wavefront *> inputs, int avgNdx, int Width, QPrinter &printer);
 
 signals:
     void currentNdxChanged(int);