Index: /issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h
===================================================================
--- /issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h	(revision 5195)
+++ /issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h	(revision 5196)
@@ -268,4 +268,5 @@
 	PatchNodesEnum,
 	TimeEnum,
+	OutputFrequencyEnum,
 	/*}}}*/
 	/*Convergence{{{1*/
Index: /issm/trunk/src/c/EnumDefinitions/EnumToString.cpp
===================================================================
--- /issm/trunk/src/c/EnumDefinitions/EnumToString.cpp	(revision 5195)
+++ /issm/trunk/src/c/EnumDefinitions/EnumToString.cpp	(revision 5196)
@@ -238,4 +238,5 @@
 		case PatchNodesEnum : return "PatchNodes";
 		case TimeEnum : return "Time";
+		case OutputFrequencyEnum : return "OutputFrequency";
 		case RelativeEnum : return "Relative";
 		case ResidualEnum : return "Residual";
Index: /issm/trunk/src/c/EnumDefinitions/StringToEnum.cpp
===================================================================
--- /issm/trunk/src/c/EnumDefinitions/StringToEnum.cpp	(revision 5195)
+++ /issm/trunk/src/c/EnumDefinitions/StringToEnum.cpp	(revision 5196)
@@ -236,4 +236,5 @@
 	else if (strcmp(name,"PatchNodes")==0) return PatchNodesEnum;
 	else if (strcmp(name,"Time")==0) return TimeEnum;
+	else if (strcmp(name,"OutputFrequency")==0) return OutputFrequencyEnum;
 	else if (strcmp(name,"Relative")==0) return RelativeEnum;
 	else if (strcmp(name,"Residual")==0) return ResidualEnum;
Index: /issm/trunk/src/c/modules/ModelProcessorx/CreateParameters.cpp
===================================================================
--- /issm/trunk/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 5195)
+++ /issm/trunk/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 5196)
@@ -39,4 +39,5 @@
 	parameters->AddObject(new BoolParam(IsStokesEnum,iomodel->isstokes));
 	parameters->AddObject(new IntParam(VerboseEnum,iomodel->verbose));
+	parameters->AddObject(new IntParam(OutputFrequencyEnum,iomodel->output_frequency));
 	parameters->AddObject(new DoubleParam(EpsResEnum,iomodel->eps_res));
 	parameters->AddObject(new DoubleParam(EpsRelEnum,iomodel->eps_rel));
Index: /issm/trunk/src/c/objects/IoModel.cpp
===================================================================
--- /issm/trunk/src/c/objects/IoModel.cpp	(revision 5195)
+++ /issm/trunk/src/c/objects/IoModel.cpp	(revision 5196)
@@ -162,4 +162,5 @@
 	IoModelFetchData(&this->verbose,iomodel_handle,"verbose");
 	IoModelFetchData(&this->plot,iomodel_handle,"plot");
+	IoModelFetchData(&this->output_frequency,iomodel_handle,"output_frequency");
 	IoModelFetchData(&this->prognostic_DG,iomodel_handle,"prognostic_DG");
 	IoModelFetchData(&this->nsteps,iomodel_handle,"nsteps");
@@ -318,4 +319,5 @@
 	this->verbose=0;
 	this->plot=0;
+	this->output_frequency=0;
 	this->eps_res=0;
 	this->eps_rel=0;
Index: /issm/trunk/src/c/objects/IoModel.h
===================================================================
--- /issm/trunk/src/c/objects/IoModel.h	(revision 5195)
+++ /issm/trunk/src/c/objects/IoModel.h	(revision 5196)
@@ -137,4 +137,5 @@
 		int     verbose;
 		int     plot;
+		int     output_frequency;
 		double  eps_res;
 		double  eps_rel;
Index: /issm/trunk/src/c/solutions/transient2d_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/transient2d_core.cpp	(revision 5195)
+++ /issm/trunk/src/c/solutions/transient2d_core.cpp	(revision 5196)
@@ -20,4 +20,5 @@
 	double dt,yts;
 	int    solution_type;
+	int    output_frequency;
 	
 	/*intermediary: */
@@ -31,4 +32,5 @@
 	femmodel->parameters->FindParam(&yts,YtsEnum);
 	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
+	femmodel->parameters->FindParam(&output_frequency,OutputFrequencyEnum);
 
 	/*initialize: */
@@ -54,5 +56,5 @@
 
 		if(verbose)_printf_("%s\n","saving results:\n");
-		if(solution_type==Transient2DSolutionEnum){
+		if(solution_type==Transient2DSolutionEnum && step%output_frequency==0){
 			InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VxEnum,step,time); 
 			InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VyEnum,step,time);
Index: /issm/trunk/src/c/solutions/transient3d_core.cpp
===================================================================
--- /issm/trunk/src/c/solutions/transient3d_core.cpp	(revision 5195)
+++ /issm/trunk/src/c/solutions/transient3d_core.cpp	(revision 5196)
@@ -20,4 +20,5 @@
 	double dt,yts;
 	int    solution_type;
+	int    output_frequency;
 	
 	/*intermediary: */
@@ -30,5 +31,5 @@
 	femmodel->parameters->FindParam(&dt,DtEnum);
 	femmodel->parameters->FindParam(&yts,YtsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
+	femmodel->parameters->FindParam(&output_frequency,OutputFrequencyEnum);
 
 	/*initialize: */
@@ -60,5 +61,5 @@
 
 		if(verbose)_printf_("%s\n","saving results:\n");
-		if(solution_type==Transient3DSolutionEnum){
+		if(solution_type==Transient3DSolutionEnum && step%output_frequency==0){
 			InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VxEnum,step,time);
 			InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VyEnum,step,time);