Index: /issm/trunk-jpl/src/c/Makefile.am
===================================================================
--- /issm/trunk-jpl/src/c/Makefile.am	(revision 25631)
+++ /issm/trunk-jpl/src/c/Makefile.am	(revision 25632)
@@ -803,4 +803,5 @@
 issm_dakota_CXXFLAGS= $(CXXFLAGS) $(CXXOPTFLAGS) $(COPTFLAGS)
 issm_dakota_LDADD= $(LDADD)
+bin_PROGRAMS += issm_post
 issm_post_SOURCES = main/issm_post.cpp
 issm_post_CXXFLAGS= $(CXXFLAGS) $(CXXOPTFLAGS)
Index: /issm/trunk-jpl/src/m/classes/qmu.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/qmu.py	(revision 25631)
+++ /issm/trunk-jpl/src/m/classes/qmu.py	(revision 25632)
@@ -25,4 +25,5 @@
         self.method = OrderedDict()
         self.params = OrderedStruct()
+        self.statistics = OrderedStruct() # TODO: Replace this with 'qmustatistics()' once class is translated
         self.results = OrderedDict()
         self.numberofresponses = 0
@@ -217,3 +218,6 @@
             flag = False
         WriteData(fid, prefix, 'data', flag, 'name', 'md.qmu.mass_flux_segments_present', 'format', 'Boolean')
+        # TODO: Uncomment this and remove the following line once qmustatistics class is translated
+        #self.statistics.marshall
+        WriteData(fid,prefix,'name','md.qmu.statistics','data',0,'format','Boolean');
     # }}}
Index: /issm/trunk-jpl/src/m/classes/qmustatistics.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/qmustatistics.m	(revision 25631)
+++ /issm/trunk-jpl/src/m/classes/qmustatistics.m	(revision 25632)
@@ -9,5 +9,5 @@
 		%generic information: 
 		nfiles_per_directory = 5;  %number of files per output directory.
-		ndirectories         = 50; %number of output directories.  %ndirectories should be < numcpus
+		ndirectories         = 50; %number of output directories; should be < numcpus
 	
 		method = struct('name',{'None'},'fields',{[]},'steps',{[]},'nbins',{NaN},'indices',{[]});
@@ -55,5 +55,5 @@
 			md = checkfield(md,'fieldname','qmu.statistics.nfiles_per_directory','>=',1);
 			if self.ndirectories>md.cluster.np, 
-				error('qmustatistics consistency check: number of cluster cpus should be > to number of output directories');
+				error('qmustatistics consistency check: number of cluster CPUs should be > number of output directories');
 			end
 			if (self.ndirectories*self.nfiles_per_directory)~=md.qmu.method.params.samples,
@@ -86,5 +86,5 @@
 		function disp(self) % {{{
 
-			disp(sprintf('qmustatistics:  post-dakota run processing of QMU statistics:'));
+			disp(sprintf('qmustatistics: post-Dakota run processing of QMU statistics:'));
 				
 			if strcmpi(self.method(1).name,'None'),  return; end;
@@ -92,5 +92,5 @@
 			%generic information:
 			fielddisplay(self,'nfiles_per_directory','number of files per output directory');
-			fielddisplay(self,'ndirectories','number of output directories.  %ndirectories should be < numcpus');
+			fielddisplay(self,'ndirectories','number of output directories; should be < numcpus');
 
 			for i=1:length(self.method),
Index: /issm/trunk-jpl/src/m/classes/sealevelmodel.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/sealevelmodel.py	(revision 25631)
+++ /issm/trunk-jpl/src/m/classes/sealevelmodel.py	(revision 25632)
@@ -70,4 +70,6 @@
         s += '{}\n'.format(fielddisplay(self, 'cluster', 'cluster parameters (number of cpus...'))
         s += '{}\n'.format(fielddisplay(self, 'miscellaneous', 'miscellaneous fields'))
+
+        return s
     #}}}
 
Index: /issm/trunk-jpl/src/m/qmu/postqmu.py
===================================================================
--- /issm/trunk-jpl/src/m/qmu/postqmu.py	(revision 25631)
+++ /issm/trunk-jpl/src/m/qmu/postqmu.py	(revision 25632)
@@ -9,14 +9,13 @@
 
 def postqmu(md):
-    '''
-    POSTQMU - Deal with Dakota output results in files
+    """POSTQMU - Deal with Dakota output results in files
     
-        Usage:
-            md = postqmu(md)
+    Usage:
+        md = postqmu(md)
 
-        TODO:
-        - Run Dakota test to check that updates from 6/26 are working
-        - Add checks to Popen
-    '''
+    TODO:
+    - Run Dakota test to check that updates from 6/26 are working
+    - Add checks to Popen
+    """
 
     # check to see if dakota returned errors in the err file
@@ -56,8 +55,9 @@
                 dakotaresults.modelresults.append(md2.results)
 
-    if md.qmu.statistics.method[0].name != 'None':
-        md.qmu.isdakota = 0
-        md = loadresultsfromdisk(md, [md.miscellaneous.name,'.stats'])
-        md.qmu.isdakota = 1
+    # TODO: Uncomment the following once qmustatistics class is translated
+    # if md.qmu.statistics.method[0].name != 'None':
+    #     md.qmu.isdakota = 0
+    #     md = loadresultsfromdisk(md, [md.miscellaneous.name,'.stats'])
+    #     md.qmu.isdakota = 1
 
     # put dakotaresults in their right location.
Index: /issm/trunk-jpl/src/m/qmu/preqmu.py
===================================================================
--- /issm/trunk-jpl/src/m/qmu/preqmu.py	(revision 25631)
+++ /issm/trunk-jpl/src/m/qmu/preqmu.py	(revision 25632)
@@ -11,19 +11,20 @@
 
 def preqmu(md, options):
-    '''
-    QMU - apply Quantification of Margins and Uncertainties techniques
+    """QMU - apply Quantification of Margins and Uncertainties techniques
     to a solution sequence (like stressbalance.py, progonstic.py, etc ...),
     using the Dakota software from Sandia.
 
-    options come from the solve.py routine. They can include Dakota options:
+    Options come from the solve.py routine. They can include Dakota options:
 
-    qmufile: input file for Dakota
-
-    (ivap, iresp, imethod, and iparams are currently unimplemented)
-    ivar: selection number for variables input (if several are specified in variables)
-    iresp: same thing for response functions
-    imethod: same thing for methods
-    iparams: same thing for params
-    '''
+        qmufile : input file for Dakota
+        ivap    : <currently unimplemented>
+        iresp   : <currently unimplemented>
+        imethod : <currently unimplemented>
+        iparams : <currently unimplemented>
+        ivar    : selection number for variables input (if several are specified in variables)
+        iresp   : same thing for response functions
+        imethod : same thing for methods
+        iparams : same thing for params
+    """
 
     print('preprocessing dakota inputs')