Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp	(revision 15625)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp	(revision 15626)
@@ -17,5 +17,5 @@
 	IssmDouble rho_ice;
 	IssmDouble FSreconditioning;
-	bool       isSSA,isL1L2,isHO,isFS,iscoupling;
+	bool       isSIA,isSSA,isL1L2,isHO,isFS,iscoupling;
 	bool       spcpresent = false;
 	int        Mx,Nx;
@@ -46,4 +46,5 @@
 	iomodel->Constant(&rho_ice,MaterialsRhoIceEnum);
 	iomodel->Constant(&FSreconditioning,DiagnosticFSreconditioningEnum);
+	iomodel->Constant(&isSIA,FlowequationIsSIAEnum);
 	iomodel->Constant(&isSSA,FlowequationIsSSAEnum);
 	iomodel->Constant(&isL1L2,FlowequationIsL1L2Enum);
@@ -61,5 +62,5 @@
 
 	/*Do we have coupling*/
-	if( (isSSA?1.:0.) + (isL1L2?1.:0.) + (isHO?1.:0.) + (isFS?1.:0.) >1.)
+	if((isSIA?1.:0.) + (isSSA?1.:0.) + (isL1L2?1.:0.) + (isHO?1.:0.) + (isFS?1.:0.) >1.)
 	 iscoupling = true;
 	else