Index: /issm/trunk-jpl/src/m/solvers/asmoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/asmoptions.m	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/asmoptions.m	(revision 21988)
@@ -11,5 +11,5 @@
 %default asm options
 asm.toolkit='petsc';
-asm.mat_type=getfieldvalue(options,'mat_type','aij');
+asm.mat_type=getfieldvalue(options,'mat_type','mpiaij');
 asm.ksp_type=getfieldvalue(options,'ksp_type','gmres');
 asm.pc_type=getfieldvalue(options,'pc_type','asm');
Index: /issm/trunk-jpl/src/m/solvers/asmoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/asmoptions.py	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/asmoptions.py	(revision 21988)
@@ -10,5 +10,5 @@
 	arguments=pairoptions.pairoptions(*args) 
 	
-	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','lu'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-30]];
+	options=[['toolkit','petsc'],['mat_type','mpiaij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','lu'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-30]];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/iluasmoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/iluasmoptions.m	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/iluasmoptions.m	(revision 21988)
@@ -9,5 +9,5 @@
 iluasm=struct();
 iluasm.toolkit='petsc';
-iluasm.mat_type=getfieldvalue(options,'mat_type','aij');
+iluasm.mat_type=getfieldvalue(options,'mat_type','mpiaij');
 iluasm.ksp_type=getfieldvalue(options,'ksp_type','gmres');
 iluasm.pc_type=getfieldvalue(options,'pc_type','asm');
Index: /issm/trunk-jpl/src/m/solvers/iluasmoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/iluasmoptions.py	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/iluasmoptions.py	(revision 21988)
@@ -16,5 +16,5 @@
 	#default iluasm options
 	iluasm['toolkit']='petsc'
-	iluasm['mat_type']=options.getfieldvalue('mat_type','aij')
+	iluasm['mat_type']=options.getfieldvalue('mat_type','mpiaij')
 	iluasm['ksp_type']=options.getfieldvalue('ksp_type','gmres')
 	iluasm['pc_type']=options.getfieldvalue('pc_type','asm')
Index: /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.m	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.m	(revision 21988)
@@ -11,5 +11,5 @@
 %default jacobiasm options
 jacobiasm.toolkit='petsc';
-jacobiasm.mat_type=getfieldvalue(options,'mat_type','aij');
+jacobiasm.mat_type=getfieldvalue(options,'mat_type','mpiaij');
 jacobiasm.ksp_type=getfieldvalue(options,'ksp_type','gmres');
 jacobiasm.pc_type=getfieldvalue(options,'pc_type','asm');
Index: /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py	(revision 21988)
@@ -10,5 +10,5 @@
 	arguments=pairoptions.pairoptions(*args) 
 	
-	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','jacobi'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-15]];
+	options=[['toolkit','petsc'],['mat_type','mpiaij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','jacobi'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-15]];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/jacobicgoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobicgoptions.m	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/jacobicgoptions.m	(revision 21988)
@@ -11,5 +11,5 @@
 %default jacobiasm options
 jacobicg.toolkit='petsc';
-jacobicg.mat_type=getfieldvalue(options,'mat_type','aij');
+jacobicg.mat_type=getfieldvalue(options,'mat_type','mpiaij');
 jacobicg.ksp_type=getfieldvalue(options,'ksp_type','cg');
 jacobicg.ksp_max_it=getfieldvalue(options,'ksp_max_it',100);
Index: /issm/trunk-jpl/src/m/solvers/jacobicgoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobicgoptions.py	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/jacobicgoptions.py	(revision 21988)
@@ -10,5 +10,5 @@
 	arguments=pairoptions.pairoptions(*args) 
 	
-	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','cg'],['ksp_max_it',100],['ksp_rtol',1e-15]];
+	options=[['toolkit','petsc'],['mat_type','mpiaij'],['ksp_type','cg'],['ksp_max_it',100],['ksp_rtol',1e-15]];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/soroptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/soroptions.m	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/soroptions.m	(revision 21988)
@@ -11,5 +11,5 @@
 %default sor options
 sor.toolkit='petsc';
-sor.mat_type=getfieldvalue(options,'mat_type','aij');
+sor.mat_type=getfieldvalue(options,'mat_type','mpiaij');
 sor.ksp_type=getfieldvalue(options,'ksp_type','cg');
 sor.pc_type=getfieldvalue(options,'pc_type','sor');
Index: /issm/trunk-jpl/src/m/solvers/soroptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/soroptions.py	(revision 21987)
+++ /issm/trunk-jpl/src/m/solvers/soroptions.py	(revision 21988)
@@ -10,5 +10,5 @@
 	arguments=pairoptions.pairoptions(*args) 
 	
-	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','cg'],['pc_type','sor'],['pc_sor_omega',1.1],['pc_sor_its',2]];
+	options=[['toolkit','petsc'],['mat_type','mpiaij'],['ksp_type','cg'],['pc_type','sor'],['pc_sor_omega',1.1],['pc_sor_its',2]];
 
 	#now, go through our arguments, and write over default options.