Index: /issm/trunk-jpl/src/c/classes/Params/Parameters.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Params/Parameters.cpp	(revision 22559)
+++ /issm/trunk-jpl/src/c/classes/Params/Parameters.cpp	(revision 22560)
@@ -658,13 +658,4 @@
 		/*In case we are using PETSC, we do not rely on issmoptions. Instead, we dump issmoptions into the Petsc 
 		 * options database: */
-
-		#if _PETSC_MAJOR_ == 2 
-		PetscOptionsDestroy();
-		PetscOptionsCreate();
-		//PetscOptionsCheckInitial_Private();
-		//PetscOptionsCheckInitial_Components();
-		PetscOptionsSetFromOptions();
-		PetscOptionsInsertMultipleString(options); //our patch
-		#else
 		#if (_PETSC_MINOR_>=7)
 		PetscOptionsSetFromOptions(NULL);
@@ -674,7 +665,8 @@
 		PetscOptionsClear();
 		#endif
-		//PetscOptionsSetFromOptions();
-		PetscOptionsInsertMultipleString(options); //our patch
-		#endif
+		int ierr = PetscOptionsInsertString(NULL,options);
+		//int ierr = PetscOptionsInsertString(NULL,"-mat_type mpiaij -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package mumps -mat_mumps_icntl_14 120 -mat_mumps_icntl_28 2 -mat_mumps_icntl_29 2");
+		//int ierr = PetscOptionsInsertString(NULL,"-mat_type mpiaij -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package mumps -mat_mumps_icntl_14 120");
+		if(ierr) _error_("Could not enter PETSc options");
 
 	#endif
Index: /issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp	(revision 22559)
+++ /issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp	(revision 22560)
@@ -25,5 +25,5 @@
 
 	/*parameters: */
-	double sparsity=.001; //default
+	double sparsity=0.001; //default
 	int    m,n;
 	int    d_nz,o_nz,nnz;
@@ -114,13 +114,7 @@
 	MatGetType(outmatrix,&type);
 
-	#if _PETSC_MAJOR_ == 2 
-	if((strcmp(type,"mpiaij")==0) || (strcmp(type,"aijmumps")==0)){
-		MatMPIAIJSetPreallocation(outmatrix,d_nz,NULL,o_nz,NULL);
-	}
-	#else
 	if((strcmp(type,"mpiaij")==0) || (strcmp(type,"mpidense")==0)){
 		MatMPIAIJSetPreallocation(outmatrix,d_nz,NULL,o_nz,NULL);
 	}
-	#endif
 
 	return outmatrix;