Index: ../trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/FemModel.cpp	(revision 18828)
+++ ../trunk-jpl/src/c/classes/FemModel.cpp	(revision 18829)
@@ -513,8 +513,9 @@
 			break;
 
 		case TransientSolutionEnum:{
-			bool isSIA,isthermal,isenthalpy,ismasstransport,isgroundingline,isstressbalance,islevelset,ishydrology,isdamage;
+			bool isSIA,isFS,isthermal,isenthalpy,ismasstransport,isgroundingline,isstressbalance,islevelset,ishydrology,isdamage;
 			iomodel->Constant(&isSIA,FlowequationIsSIAEnum);
+			iomodel->Constant(&isFS,FlowequationIsFSEnum);
 			iomodel->Constant(&isthermal,TransientIsthermalEnum);
 			iomodel->Constant(&isenthalpy,ThermalIsenthalpyEnum);
 			iomodel->Constant(&islevelset,TransientIslevelsetEnum);
@@ -524,11 +525,19 @@
 			iomodel->Constant(&isdamage,TransientIsdamageevolutionEnum);
 			iomodel->Constant(&ishydrology,TransientIshydrologyEnum);
 			if(isstressbalance){
+				int  fe_FS;
+				iomodel->Constant(&fe_FS,FlowequationFeFSEnum);
 				analyses_temp[numanalyses++]=StressbalanceAnalysisEnum;
 				analyses_temp[numanalyses++]=StressbalanceVerticalAnalysisEnum;
 				if(isSIA){
 					analyses_temp[numanalyses++]=StressbalanceSIAAnalysisEnum;
 				}
+				analyses_temp[numanalyses++]=L2ProjectionBaseAnalysisEnum;
+				analyses_temp[numanalyses++]=ExtrudeFromBaseAnalysisEnum;
+				analyses_temp[numanalyses++]=DepthAverageAnalysisEnum;
+				if(fe_FS==LATaylorHoodEnum || fe_FS==LACrouzeixRaviartEnum){
+					analyses_temp[numanalyses++]=UzawaPressureAnalysisEnum;
+				}
 			}
 			if(isthermal && iomodel->domaintype==Domain3DEnum){
 				if(isenthalpy){
@@ -765,6 +774,7 @@
 				case ThicknessAlongGradientEnum:    ThicknessAlongGradientx(&double_result,elements,nodes,vertices,loads,materials,parameters);     break;
 				case ThicknessAcrossGradientEnum:   ThicknessAcrossGradientx(&double_result,elements,nodes,vertices,loads,materials,parameters);    break;
 				case RheologyBbarAbsGradientEnum:   RheologyBbarAbsGradientx(&double_result,elements,nodes,vertices,loads,materials,parameters);    break;
+				case RheologyBAbsGradientEnum:      RheologyBAbsGradientx(&double_result,elements,nodes,vertices,loads,materials,parameters);  break;
 				case DragCoefficientAbsGradientEnum:DragCoefficientAbsGradientx(&double_result,elements,nodes,vertices,loads,materials,parameters); break;
 				case BalancethicknessMisfitEnum:    BalancethicknessMisfitx(&double_result);                                                        break;
 				case SurfaceAbsMisfitEnum:          SurfaceAbsMisfitx(&double_result); break;