Index: /issm/trunk/src/m/classes/public/LaunchQueueJobGeneric.m
===================================================================
--- /issm/trunk/src/m/classes/public/LaunchQueueJobGeneric.m	(revision 2048)
+++ /issm/trunk/src/m/classes/public/LaunchQueueJobGeneric.m	(revision 2049)
@@ -19,5 +19,5 @@
 	packages={[md.name '.bin'],[md.name '.queue']};
 	if md.qmu_analysis,
-		packages{end+1}=[md.name '.qmu.in '];
+		packages{end+1}=[md.name '.qmu.in'];
 	end
 
Index: /issm/trunk/src/m/classes/public/importancefactors.m
===================================================================
--- /issm/trunk/src/m/classes/public/importancefactors.m	(revision 2049)
+++ /issm/trunk/src/m/classes/public/importancefactors.m	(revision 2049)
@@ -0,0 +1,57 @@
+function factors=importancefactors(md,variablename,responsename)
+%IMPORTANCEFACTORS - compute importance factors for a certain variable and response.
+%
+%   Usage:
+%      factors=importancefactors(md,variablename,responsename)
+%
+%
+%   Example: factors=importancefactors(md,'drag','max_vel');
+%
+
+
+variablenamelength=length(variablename);
+
+%go through all response functions and find the one corresponding to the correct responsename
+responsefunctions=md.dakotaresults.dresp_out;
+found=0;
+for i=1:length(responsefunctions),
+	if strcmpi(responsefunctions(i).descriptor,responsename),
+		found=i;
+		break;
+	end
+end
+if ~found,
+	error('importancefactors error message: could not find correct response function');
+end
+responsefunctions=responsefunctions(found);
+nfun=size(responsefunctions.var,1);
+
+%Now recover response to the correct design variable
+importancefactors=zeros(1,0);
+count=0;
+for i=1:nfun,
+	desvar=responsefunctions.var{i};
+	if strncmpi(desvar,variablename,variablenamelength),
+		importancefactors(end+1)=responsefunctions.impfac(i);
+		count=count+1;
+	end
+end
+
+if count==0,
+	error('importancefactors error message: could not find to response functions with corresponding design variable');
+end
+
+if count==1, %we have scalar
+	factors=importancefactors;
+	return;
+else
+	%Ok, get partitioning, unless already supplied
+	if isempty(md.part),
+		[epart npart]=MeshPartition(md,md.npart);
+	else
+		npart=md.part+1;
+	end
+
+	%distribute importance factor
+	factors=importancefactors(npart)';
+end
Index: sm/trunk/src/m/classes/public/plot/plot_qmuimportancefactors.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/plot_qmuimportancefactors.m	(revision 2048)
+++ 	(revision )
@@ -1,87 +1,0 @@
-function plot_qmuimportancefactors(md,options_structure,width,ii);
-%PLOT_IMPORTANCEFACTORS - plot importance factors
-%
-%   Usage:
-%      plot_importancefactors(md,options_structure,width,i);
-%
-%   See also: PLOTMODEL
-
-%first recover design variable descriptor
-if ~isnan(options_structure.designvariable),
-	descriptor=options_structure.designvariable;
-else
-	error('plot_importancefactors error message: Need to supply design variable descriptor');
-end
-descriptorlength=length(descriptor);
-
-%then recover responsfunction name
-if ~isnan(options_structure.responsefunction),
-	responsefunctiondescriptor=options_structure.responsefunction;
-else
-	error('plot_importancefactors error message: Need to supply response function descriptor');
-end
-
-%go through all response functions and find the one corresponding to the correct responsefunctiondescriptor
-responsefunctions=md.dakotaresults.dresp_out;
-found=0;
-for i=1:length(responsefunctions),
-	if strcmpi(responsefunctions(i).descriptor,responsefunctiondescriptor),
-		found=i;
-		break;
-	end
-end
-if ~found,
-	error('plot_importancefactors error message: could not find correct response function');
-end
-responsefunctions=responsefunctions(found);
-nfun=size(responsefunctions.var,1);
-
-%Now recover response to the correct design variable
-importancefactors=zeros(md.npart,1);
-count=0;
-for i=1:nfun,
-	var=responsefunctions.var{i};
-	if strncmpi(var,descriptor,descriptorlength),
-		count=count+1;
-		importancefactors(count)=responsefunctions.impfac(i);
-	end
-end
-if count==0,
-	error('plot_importancefactors error message: could not find to response functions with corresponding design variable');
-end
-
-%log?
-if ~isnan(options_structure.log),
-	logvalue=options_structure.log;
-	importancefactors=log(importancefactors)/log(logvalue);
-end
-
-%Ok, get partitioning.
-[epart npart]=MeshPartition(md,md.npart);
-
-%distribute importance factor
-gridimportance=importancefactors(npart);
-
-%process data and model
-[x y z elements is2d]=processmesh(md,options_structure);
-
-%edgecolor?
-if ~isnan(options_structure.edgecolor),
-	edgecolor=options_structure.edgecolor;
-else
-	edgecolor='none';
-end
-
-%standard plot:
-subplot(width,width,ii);
-
-%ok, plot gridimportance now.
-if is2d,
-	A=elements(:,1); B=elements(:,2); C=elements(:,3); 
-	patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', gridimportance,'FaceColor','interp','EdgeColor',edgecolor);
-else
-	error('plot_importancefactors error message: 3d meshes not supported yet');
-end
-
-%apply options
-applyoptions(md,[],options_structure);
Index: /issm/trunk/src/m/classes/public/plot/plotdoc.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/plotdoc.m	(revision 2048)
+++ /issm/trunk/src/m/classes/public/plot/plotdoc.m	(revision 2049)
@@ -42,5 +42,4 @@
 disp('                  - ''transient_movie_time'': time for each image (default 2 seconds)');
 disp('                  - ''thermaltransient_results'': this will display all the time steps of a thermal transient run');
-disp('                  - ''importancefactors'': qmu results');
 disp('                  - ''normplot'': norm plot');
 disp('                  - ''histnorm'': histogram normal distribution');
@@ -95,6 +94,4 @@
 disp('       ''icesheet'': ''none'' if ice sheets are to not be plotted');
 disp('       ''log'': value of log');
-disp('       ''designvariable'': name of design variable for importance factors');
-disp('       ''responsefunction'': name of response function for importance factors');
 disp('       ''latlon'': ''on'' or {latstep lonstep [resolution [color]]} where latstep,longstep and resolution are in degrees, color is a [r g b] array');
 disp('       ''latlonnumbering'': ''on'' or {latgap longap colornumber latangle lonangle} where latgap and longap are pixel gaps for the numbers,');