Index: /issm/trunk/src/m/solutions/flaim.m
===================================================================
--- /issm/trunk/src/m/solutions/flaim.m	(revision 9420)
+++ /issm/trunk/src/m/solutions/flaim.m	(revision 9421)
@@ -32,4 +32,6 @@
 %  process file names
 
+[pathstr,name,ext,versn] = fileparts(md.fm_tracks);
+fm_tracks=fullfile('',[name ext versn]);
 [pathstr,name,ext,versn] = fileparts(options.fmfile);
 if isempty(ext)
@@ -81,6 +83,6 @@
 display('Calling flaim.');
 FLAIM_DIR=flaimdir();
-%system([FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' md.fm_tracks ' -kml ' filekml ' -gridsatequator 200000 -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints -usevalueordering -split-antimeridian >& ' filelog]);
-flaim_cmd=[FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' md.fm_tracks ' -kml ' filekml ' -gridsatequator ' int2str(md.fm_gridsatequator) ' -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints'];
+%system([FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' fm_tracks ' -kml ' filekml ' -gridsatequator 200000 -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints -usevalueordering -split-antimeridian >& ' filelog]);
+flaim_cmd=[FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' fm_tracks ' -kml ' filekml ' -gridsatequator ' int2str(md.fm_gridsatequator) ' -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints'];
 if (md.fm_usevalueordering)
 	flaim_cmd=[flaim_cmd ' -usevalueordering'];
@@ -100,6 +102,14 @@
 if exist(filelog,'file')
 	[status,fline]=system(['grep -i quality ' filelog]);
-	strings=textscan(fline,'%s','delimiter',' =,','MultipleDelimsAsOne',1);
-	md.results.FlaimSolution.quality =sscanf(strings{1}{2},'%f');
+	if ~isempty(fline)
+		strings=textscan(fline,'%s','delimiter',' =,','MultipleDelimsAsOne',1);
+		md.results.FlaimSolution.quality =sscanf(strings{1}{2},'%f');
+	else
+		warning('Could not find ''quality'' in the ''%s/%s'' flaim log file.',options.fmdir,filelog);
+		md.results.FlaimSolution.quality =NaN;
+	end
+end
+if ~exist(md.results.FlaimSolution.solution,'file')
+	md.results.FlaimSolution.solution='';
 end
 display(md.results.FlaimSolution);
Index: /issm/trunk/src/m/solutions/issm.m
===================================================================
--- /issm/trunk/src/m/solutions/issm.m	(revision 9420)
+++ /issm/trunk/src/m/solutions/issm.m	(revision 9421)
@@ -35,8 +35,31 @@
 	else
 		%launch dakota driver for diagnostic core solution
+		display('calling Qmu');
 		Qmu(femmodel,femmodel.parameters);
+		display('finished with Qmu');
 	end
+
+	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MinVx',true);
+	disp(sprintf('MinVx=%g\n',dtest));
+	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxVx',true);
+	disp(sprintf('MaxVx=%g\n',dtest));
+	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxAbsVx',true);
+	disp(sprintf('MaxAbsVx=%g\n',dtest));
+	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MinVy',true);
+	disp(sprintf('MinVy=%g\n',dtest));
+	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxVy',true);
+	disp(sprintf('MaxVy=%g\n',dtest));
+	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxAbsVy',true);
+	disp(sprintf('MaxAbsVy=%g\n',dtest));
+%	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MinVel',true);
+%	disp(sprintf('MinVel=%g\n',dtest));
+%	dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxVel',true);
+%	disp(sprintf('MaxVel=%g\n',dtest));
 
 	%stop timing
 	t2=clock;
 	disp(['solution converged in ' num2str(etime(t2,t1)) ' seconds']);
+
+	display('writing femmodel to femmodel.mat');
+	save femmodel femmodel
+