Index: /issm/trunk/src/m/solutions/dakota/qmuoutput.m
===================================================================
--- /issm/trunk/src/m/solutions/dakota/qmuoutput.m	(revision 781)
+++ /issm/trunk/src/m/solutions/dakota/qmuoutput.m	(revision 781)
@@ -0,0 +1,69 @@
+function md=qmuoutput(md,qmufile,qmudir);
+
+%INPUT function md=qmuoutput(md,qmufile,qmudir)
+%Deal with dakota output results in files.
+
+global ISSM_DIR;
+
+%  check to see if dakota returned errors in the err file
+qmuerrfile=[qmufile '.err'];
+if md.qmu_analysis,
+	qmuerrfile=[qmufile '.qmu.err'];
+end
+
+if exist(qmuerrfile,'file')
+   fide=fopen(qmuerrfile,'r');
+   fline=fgetl(fide);
+   if ischar(fline)
+       while ischar(fline)
+           disp(sprintf('%s',fline));
+           fline=fgetl(fide);
+       end
+       status=fclose(fide);
+       cd ../
+       error(['Dakota returned error in ''' qmuerrfile ' file.  ''' qmudir ''' directory retained.'])
+    end
+    status=fclose(fide);
+else
+   cd ../
+   error(['Dakota did not generate ''' qmuerrfile ' file.  ''' qmudir ''' directory retained.'])
+end
+
+%parse inputs and results from dakota
+qmuinfile=[qmufile '.in'];
+qmuoutfile=[qmufile '.out'];
+if md.qmu_analysis,
+	qmuinfile=[qmufile '.qmu.in'];
+	qmuoutfile=[qmufile '.qmu.out'];
+end
+
+[method,dvar,dresp_in]=dakota_in_parse(qmuinfile);
+md.dakotaresults.method   =method;
+md.dakotaresults.dvar     =dvar;
+md.dakotaresults.dresp_in =dresp_in;
+
+[method,dresp_out,scm,pcm,srcm,prcm]=dakota_out_parse(qmuoutfile);
+md.dakotaresults.dresp_out=dresp_out;
+md.dakotaresults.scm      =scm;
+md.dakotaresults.pcm      =pcm;
+md.dakotaresults.srcm     =srcm;
+md.dakotaresults.prcm     =prcm;
+
+if exist('dakota_tabular.dat','file')
+    [method,dresp_dat                  ]=dakota_out_parse('dakota_tabular.dat');
+    md.dakotaresults.dresp_dat=dresp_dat;
+end
+
+%save input and output files into model
+%md.dakotain =readfile([qmufile '.in']);
+%md.dakotaout=readfile([qmufile '.out']);
+%if exist('dakota_tabular.dat','file')
+%	md.dakotadat=readfile('dakota_tabular.dat');
+%end
+
+%  move all the individual function evalutations into zip files
+if ~md.qmu_analysis,
+	system('zip -mq params.in.zip params.in.[1-9]*');
+	system('zip -mq results.out.zip results.out.[1-9]*');
+	system('zip -mq matlab.out.zip matlab*.out.[1-9]*');
+end