Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17815)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17816)
@@ -150,5 +150,5 @@
 void  Tria::ComputeSigmaNN(){
 
-	if(~IsOnBase()){
+	if(!IsOnBase()){
 		IssmDouble sigma_nn=0;
 		this->inputs->AddInput(new TriaInput(SigmaNNEnum,&sigma_nn,P0Enum));
Index: /issm/trunk-jpl/src/c/shared/Elements/Arrhenius.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/Arrhenius.cpp	(revision 17815)
+++ /issm/trunk-jpl/src/c/shared/Elements/Arrhenius.cpp	(revision 17816)
@@ -34,4 +34,5 @@
 	/*convert temperature to absolute temperature*/
 	_assert_(depth>0);
+	printf("temperature = %g beta = %g depth = %g\n",temperature,beta,depth);
 	Tstar=temperature-beta*depth;
 	_assert_(Tstar>0);
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_newton.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_newton.cpp	(revision 17815)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_newton.cpp	(revision 17816)
@@ -68,5 +68,5 @@
 
 		/*Prepare next iteration using Newton's method*/
-		SystemMatricesx(&Kff,&Kfs,&pf,&df,NULL,femmodel);delete df;
+		SystemMatricesx(&Kff,&Kfs,&pf,&df,&kmax,femmodel);delete df;
 		CreateNodalConstraintsx(&ys,femmodel->nodes,configuration_type);
 		Reduceloadx(pf,Kfs,ys);delete Kfs;