Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 17827)
@@ -818,4 +818,13 @@
 				input=this->inputs->GetInput(output_enum);
 				break;
+			case StrainRatexxEnum:
+			case StrainRatexyEnum:
+			case StrainRatexzEnum:
+			case StrainRateyyEnum:
+			case StrainRateyzEnum:
+			case StrainRatezzEnum:
+				this->ComputeStrainRate();
+				input=this->inputs->GetInput(output_enum);
+				break;
 			case DeviatoricStressxxEnum: 
 			case DeviatoricStressxyEnum: 
Index: /issm/trunk-jpl/src/c/classes/Elements/Element.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 17827)
@@ -202,4 +202,5 @@
 		virtual void   ComputeSigmaNN(void)=0;
 		virtual void   ComputeBasalStress(Vector<IssmDouble>* sigma_b)=0;
+		virtual void   ComputeStrainRate()=0;
 		virtual void   ComputeStrainRate(Vector<IssmDouble>* eps)=0;
 		virtual void   ComputeStressTensor(void)=0;
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17827)
@@ -277,5 +277,5 @@
 /*}}}*/
 /*FUNCTION Penta::ComputeStrainRate {{{*/
-void  Penta::ComputeStrainRate(Vector<IssmDouble>* eps){
+void  Penta::ComputeStrainRate(){
 
 	IssmDouble      xyz_list[NUMVERTICES][3];
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 17827)
@@ -55,5 +55,6 @@
 		IssmDouble CharacteristicLength(void){_error_("not implemented yet");};
 		void   ComputeBasalStress(Vector<IssmDouble>* sigma_b);
-		void   ComputeStrainRate(Vector<IssmDouble>* eps);
+		void   ComputeStrainRate();
+		void   ComputeStrainRate(Vector<IssmDouble>* eps){_error_("not implemented yet");};
 		void   ComputeSigmaNN(){_error_("not implemented yet");};
 		void   ComputeStressTensor();
Index: /issm/trunk-jpl/src/c/classes/Elements/Seg.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 17827)
@@ -55,4 +55,5 @@
 		void        ComputeBasalStress(Vector<IssmDouble>* sigma_b){_error_("not implemented yet");};
 		void        ComputeSigmaNN(){_error_("not implemented yet");};
+		void        ComputeStrainRate(){_error_("not implemented yet");};
 		void        ComputeStrainRate(Vector<IssmDouble>* eps){_error_("not implemented yet");};
 		void        ComputeStressTensor(){_error_("not implemented yet");};
Index: /issm/trunk-jpl/src/c/classes/Elements/Tetra.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tetra.h	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tetra.h	(revision 17827)
@@ -55,4 +55,5 @@
 		void        ComputeBasalStress(Vector<IssmDouble>* sigma_b){_error_("not implemented yet");};
 		void        ComputeSigmaNN(){_error_("not implemented yet");};
+		void        ComputeStrainRate(){_error_("not implemented yet");};
 		void        ComputeStrainRate(Vector<IssmDouble>* eps){_error_("not implemented yet");};
 		void        ComputeStressTensor(){_error_("not implemented yet");};
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17827)
@@ -203,5 +203,5 @@
 /*}}}*/
 /*FUNCTION Tria::ComputeStrainRate {{{*/
-void  Tria::ComputeStrainRate(Vector<IssmDouble>* eps){
+void  Tria::ComputeStrainRate(){
 
 	IssmDouble      xyz_list[NUMVERTICES][3];
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 17826)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 17827)
@@ -53,5 +53,6 @@
 		void        ComputeBasalStress(Vector<IssmDouble>* sigma_b);
 		void        ComputeSigmaNN();
-		void        ComputeStrainRate(Vector<IssmDouble>* eps);
+		void		   ComputeStrainRate();
+		void		   ComputeStrainRate(Vector<IssmDouble>* eps){_error_("not implemented yet");};
 		void        ComputeStressTensor();
 		void        ComputeDeviatoricStressTensor();