Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 21893)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 21894)
@@ -313,4 +313,6 @@
 		ComputeBasalMeltingrate(element);
 	}
+	femmodel->parameters->SetParam(BasalforcingsGroundediceMeltingRateEnum,InputToExtrudeEnum);
+	extrudefrombase_core(femmodel);
 }/*}}}*/
 void           EnthalpyAnalysis::ComputeBasalMeltingrate(Element* element){/*{{{*/
Index: /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp	(revision 21893)
+++ /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp	(revision 21894)
@@ -246,4 +246,16 @@
 			}
 			break;
+		case P1xP3Enum:
+			if(start==-1){
+				for(int i=0+3;i<12;i+=3) this->values[i]=this->values[0];
+				for(int i=1+3;i<12;i+=3) this->values[i]=this->values[1];
+				for(int i=2+3;i<12;i+=3) this->values[i]=this->values[2];
+			}
+			else{
+				for(int i=0;i<9;i+=3) this->values[i]=this->values[9];
+				for(int i=1;i<9;i+=3) this->values[i]=this->values[10];
+				for(int i=2;i<9;i+=3) this->values[i]=this->values[11];
+			}
+			break;
 		default:
 			_error_("not supported yet for type "<<EnumToStringx(this->interpolation_type));