Index: /issm/trunk-jpl/src/m/classes/qmu.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/qmu.m	(revision 25040)
+++ /issm/trunk-jpl/src/m/classes/qmu.m	(revision 25041)
@@ -9,4 +9,5 @@
 		output                      = 0;
 		variables                   = struct();
+		correlation_matrix          = [];
 		responses                   = struct();
 		method                      = struct();
Index: /issm/trunk-jpl/src/m/qmu/dakota_in_data.m
===================================================================
--- /issm/trunk-jpl/src/m/qmu/dakota_in_data.m	(revision 25040)
+++ /issm/trunk-jpl/src/m/qmu/dakota_in_data.m	(revision 25041)
@@ -2,5 +2,5 @@
 %  define the data to write the dakota .in and .m files.
 %
-%  []=dakota_in_data(dmeth,variables,responses,dparams,filei,varargin)
+%  []=dakota_in_data(dmeth,variables,responses,dparams,filei,corrin,varargin)
 %
 %  where the required input is:
@@ -10,4 +10,5 @@
 %    dparams       (structure array, method-independent parameters)
 %    filei         (character, name of .in and .m files)
+%    corrin        (correlation matrix)
 %
 %  params may be empty, in which case defaults will be used.
@@ -46,5 +47,5 @@
 %  to foreign countries or providing access to foreign persons."
 %
-function []=dakota_in_data(dmeth,variables,responses,dparams,filei,varargin)
+function []=dakota_in_data(dmeth,variables,responses,dparams,filei,corrin,varargin)
 
 if ~nargin
@@ -118,5 +119,5 @@
 
 %Write in file
-dakota_in_write(dmeth,dvar,dresp,params,filei,varargin{:});
+dakota_in_write(dmeth,dvar,dresp,params,filei,corrin,varargin{:});
 
 end
Index: /issm/trunk-jpl/src/m/qmu/dakota_in_write.m
===================================================================
--- /issm/trunk-jpl/src/m/qmu/dakota_in_write.m	(revision 25040)
+++ /issm/trunk-jpl/src/m/qmu/dakota_in_write.m	(revision 25041)
@@ -40,5 +40,5 @@
 %  to foreign countries or providing access to foreign persons."
 %
-function []=dakota_in_write(method,dvar,dresp,params,filei,varargin)
+function []=dakota_in_write(method,dvar,dresp,params,filei,corrin,varargin)
 
 if ~nargin
@@ -99,5 +99,5 @@
 %  write the variables section
 
-variables_write(fidi,dmeth,dvar);
+variables_write(fidi,dmeth,dvar,corrin);
 
 %  write the interface section
@@ -181,5 +181,5 @@
 %%  function to write the variables section of the file
 
-function []=variables_write(fidi,dmeth,dvar)
+function []=variables_write(fidi,dmeth,dvar,corrin)
 
 display('Writing variables section of Dakota input file.');
@@ -198,4 +198,15 @@
     fhvar=str2func([dmeth.lcspec{i}    '.dakota_write']);
     fhvar(fidi,dvar);
+end
+
+% write correlation matrix: 
+if ~isempty(corrin),
+	fprintf(fidi,'uncertain_correlation_matrix\n');
+	for i=1:size(corrin,1),
+		for j=1:size(corrin,2),
+			fprintf(fidi,'%g ',corrin(i,j));
+		end
+		fprintf(fidi,'\n');
+	end
 end
 
Index: /issm/trunk-jpl/src/m/qmu/preqmu.m
===================================================================
--- /issm/trunk-jpl/src/m/qmu/preqmu.m	(revision 25040)
+++ /issm/trunk-jpl/src/m/qmu/preqmu.m	(revision 25041)
@@ -67,5 +67,5 @@
 
 %create in file for dakota
-dakota_in_data(md.qmu.method(imethod),variables,responses,md.qmu.params(iparams),qmufile);
+dakota_in_data(md.qmu.method(imethod),variables,responses,md.qmu.params(iparams),qmufile,md.qmu.correlation_matrix);
 
 %build a list of variables and responses descriptors. the list is not expanded.