Index: ../trunk-jpl/test/Archives/Archive437.nc
===================================================================
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Index: ../trunk-jpl/test/Archives/Archive432.nc
===================================================================
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Index: ../trunk-jpl/test/Archives/Archive436.nc
===================================================================
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Index: ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18621)
+++ ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18622)
@@ -561,6 +561,7 @@
 	/* implementation of the basal condition decision chart of Aschwanden 2012, Fig.5 */
 	if(!element->IsOnBase() || element->IsFloating()) return NULL;
 
+	bool isdynamicbasalspc;
 	int i, state;
 	IssmDouble  dt,Jdet,scalar;
 	IssmDouble	enthalpy, enthalpyup, pressure, pressureup, watercolumn, meltingrate;
@@ -578,6 +579,7 @@
 	/*Retrieve all inputs and parameters*/
 	element->GetVerticesCoordinatesBase(&xyz_list_base);
 	element->FindParam(&dt,TimesteppingTimeStepEnum);
+	element->FindParam(&isdynamicbasalspc,ThermalIsdynamicbasalspcEnum);
 	Input* vx_input             = element->GetInput(VxEnum);                          _assert_(vx_input);
 	Input* vy_input             = element->GetInput(VyEnum);                          _assert_(vy_input);
 	Input* vz_input             = element->GetInput(VzEnum);                          _assert_(vz_input);
@@ -600,14 +602,18 @@
 		element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
 		element->NodalFunctions(basis,gauss);
 
-		enthalpy_input->GetInputValue(&enthalpy,gauss);
-		enthalpy_input->GetInputValue(&enthalpyup,gaussup);
-		pressure_input->GetInputValue(&pressure,gauss);
-		pressure_input->GetInputValue(&pressureup,gaussup);
-		watercolumn_input->GetInputValue(&watercolumn,gauss);
-		meltingrate_input->GetInputValue(&meltingrate,gauss);
+		if(isdynamicbasalspc){
+			enthalpy_input->GetInputValue(&enthalpy,gauss);
+			enthalpy_input->GetInputValue(&enthalpyup,gaussup);
+			pressure_input->GetInputValue(&pressure,gauss);
+			pressure_input->GetInputValue(&pressureup,gaussup);
+			watercolumn_input->GetInputValue(&watercolumn,gauss);
+			meltingrate_input->GetInputValue(&meltingrate,gauss);
+			state=GetThermalBasalCondition(element, enthalpy, enthalpyup, pressure, pressureup, watercolumn, meltingrate);
+		}
+		else
+			state=0;
 
-		state=GetThermalBasalCondition(element, enthalpy, enthalpyup, pressure, pressureup, watercolumn, meltingrate);
 		switch (state) {
 			case 0:
 				// cold, dry base: apply basal surface forcing
