Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21818)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21819)
@@ -1736,5 +1736,5 @@
 		else{
 			/*last chance for the output definition, if the enum is one of Outputdefinition[1-10]Enum:*/
-			if(output_enum>=Outputdefinition1Enum && output_enum <=Outputdefinition10Enum){
+			if(output_enum>=Outputdefinition1Enum && output_enum <=Outputdefinition100Enum){
 				double_result = OutputDefinitionsResponsex(this,output_enum);
 				if(save_results){
@@ -2032,5 +2032,11 @@
 		case VelEnum:                      this->ElementResponsex(responses,VelEnum); break;
 		case FrictionCoefficientEnum:      NodalValuex(responses, FrictionCoefficientEnum,elements,nodes, vertices, loads, materials, parameters); break;
-		default: _error_("response descriptor \"" << EnumToStringx(response_descriptor_enum) << "\" not supported yet!"); break; 
+		default: 
+			if(response_descriptor_enum>=Outputdefinition1Enum && response_descriptor_enum <=Outputdefinition100Enum){
+				IssmDouble double_result = OutputDefinitionsResponsex(this,response_descriptor_enum);
+				*responses=double_result;
+			}
+			else _error_("response descriptor \"" << EnumToStringx(response_descriptor_enum) << "\" not supported yet!"); 
+			break; 
 	}
 
Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 21818)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 21819)
@@ -181,4 +181,9 @@
 			if isa(md.friction,'friction'),
 				md.friction.coefficient=project2d(md,md.friction.coefficient,1);
+				md.friction.p=project2d(md,md.friction.p,1);
+				md.friction.q=project2d(md,md.friction.q,1);
+			elseif isa(md.friction,'frictioncoulomb'),
+				md.friction.coefficient=project2d(md,md.friction.coefficient,1);
+				md.friction.coefficientcoulomb=project2d(md,md.friction.coefficientcoulomb,1);
 				md.friction.p=project2d(md,md.friction.p,1);
 				md.friction.q=project2d(md,md.friction.q,1);