Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18621)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18622)
@@ -562,4 +562,5 @@
 	if(!element->IsOnBase() || element->IsFloating()) return NULL;
 
+	bool isdynamicbasalspc;
 	int i, state;
 	IssmDouble  dt,Jdet,scalar;
@@ -579,4 +580,5 @@
 	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);
@@ -601,12 +603,16 @@
 		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);
-
-		state=GetThermalBasalCondition(element, enthalpy, enthalpyup, pressure, pressureup, watercolumn, meltingrate);
+		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;
+
 		switch (state) {
 			case 0: