Index: /issm/trunk/test/Par/79North.par
===================================================================
--- /issm/trunk/test/Par/79North.par	(revision 9454)
+++ /issm/trunk/test/Par/79North.par	(revision 9455)
@@ -17,5 +17,4 @@
 %Friction
 pos=find(md.elementoniceshelf);
-md.drag_type=2;
 md.drag_coefficient=50*ones(md.numberofnodes,1);
 md.drag_coefficient(md.elements(pos,:))=0;
Index: /issm/trunk/test/Par/ISMIPA.par
===================================================================
--- /issm/trunk/test/Par/ISMIPA.par	(revision 9454)
+++ /issm/trunk/test/Par/ISMIPA.par	(revision 9455)
@@ -8,5 +8,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
 %Take care of iceshelves: no basal drag
Index: /issm/trunk/test/Par/ISMIPB.par
===================================================================
--- /issm/trunk/test/Par/ISMIPB.par	(revision 9454)
+++ /issm/trunk/test/Par/ISMIPB.par	(revision 9455)
@@ -8,5 +8,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
 %Take care of iceshelves: no basal drag
Index: /issm/trunk/test/Par/ISMIPC.par
===================================================================
--- /issm/trunk/test/Par/ISMIPC.par	(revision 9454)
+++ /issm/trunk/test/Par/ISMIPC.par	(revision 9455)
@@ -8,5 +8,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 %md.drag_coefficient=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x/2)).*sin(md.y*2*pi/max(md.x/2)))./(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)));
 md.drag_coefficient=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x)).*sin(md.y*2*pi/max(md.x))));
Index: /issm/trunk/test/Par/ISMIPD.par
===================================================================
--- /issm/trunk/test/Par/ISMIPD.par	(revision 9454)
+++ /issm/trunk/test/Par/ISMIPD.par	(revision 9455)
@@ -8,5 +8,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x))));
 %Take care of iceshelves: no basal drag
Index: /issm/trunk/test/Par/ISMIPE.par
===================================================================
--- /issm/trunk/test/Par/ISMIPE.par	(revision 9454)
+++ /issm/trunk/test/Par/ISMIPE.par	(revision 9455)
@@ -20,5 +20,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=zeros(md.numberofnodes,1);
 md.drag_p=ones(md.numberofelements,1);
Index: /issm/trunk/test/Par/ISMIPF.par
===================================================================
--- /issm/trunk/test/Par/ISMIPF.par	(revision 9454)
+++ /issm/trunk/test/Par/ISMIPF.par	(revision 9455)
@@ -10,5 +10,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=sqrt(md.yts/(2.140373*10^-7*1000))*ones(md.numberofnodes,1);
 md.drag_p=ones(md.numberofelements,1);
Index: /issm/trunk/test/Par/Pig.par
===================================================================
--- /issm/trunk/test/Par/Pig.par	(revision 9454)
+++ /issm/trunk/test/Par/Pig.par	(revision 9455)
@@ -22,5 +22,4 @@
 %Friction
 pos=find(md.elementoniceshelf);
-md.drag_type=2;
 md.drag_coefficient=50*ones(md.numberofnodes,1);
 md.drag_coefficient(md.elements(pos,:))=0;
Index: /issm/trunk/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 9454)
+++ /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 9455)
@@ -6,5 +6,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=20*ones(md.numberofnodes,1); %q=1. %no drag is specified in the analytical solution
 md.drag_p=ones(md.numberofelements,1);
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 9454)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 9455)
@@ -51,5 +51,4 @@
 %Friction
 pos=find(md.elementoniceshelf);
-md.drag_type=2;
 md.drag_coefficient=20*ones(md.numberofnodes,1);
 md.drag_coefficient(md.elements(pos,:))=0;
Index: /issm/trunk/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 9454)
+++ /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 9455)
@@ -10,5 +10,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=20*ones(md.numberofnodes,1); %q=1. %no drag is specified in the analytical solution
 %Take care of iceshelves: no basal drag
Index: /issm/trunk/test/Par/SquareEISMINT.par
===================================================================
--- /issm/trunk/test/Par/SquareEISMINT.par	(revision 9454)
+++ /issm/trunk/test/Par/SquareEISMINT.par	(revision 9455)
@@ -10,5 +10,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
 %Take care of iceshelves: no basal drag
Index: /issm/trunk/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 9454)
+++ /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 9455)
@@ -25,5 +25,4 @@
 %Friction
 pos=find(md.elementoniceshelf);
-md.drag_type=2;
 md.drag_coefficient=20*ones(md.numberofnodes,1);
 md.drag_coefficient(md.elements(pos,:))=0;
Index: /issm/trunk/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetShelf.par	(revision 9454)
+++ /issm/trunk/test/Par/SquareSheetShelf.par	(revision 9455)
@@ -32,5 +32,4 @@
 %Friction
 pos=find(md.elementoniceshelf);
-md.drag_type=2;
 md.drag_coefficient=20*ones(md.numberofnodes,1);
 md.drag_coefficient(md.elements(pos,:))=0;
Index: /issm/trunk/test/Par/SquareShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareShelf.par	(revision 9454)
+++ /issm/trunk/test/Par/SquareShelf.par	(revision 9455)
@@ -25,5 +25,4 @@
 %Friction
 pos=find(md.elementoniceshelf);
-md.drag_type=2;
 md.drag_coefficient=20*ones(md.numberofnodes,1);
 md.drag_coefficient(md.elements(pos,:))=0;
Index: /issm/trunk/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 9454)
+++ /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 9455)
@@ -29,5 +29,4 @@
 %Friction
 pos=find(md.elementoniceshelf);
-md.drag_type=2;
 md.drag_coefficient=20*ones(md.numberofnodes,1);
 md.drag_coefficient(md.elements(pos,:))=0;
Index: /issm/trunk/test/Par/SquareThermal.par
===================================================================
--- /issm/trunk/test/Par/SquareThermal.par	(revision 9454)
+++ /issm/trunk/test/Par/SquareThermal.par	(revision 9455)
@@ -16,5 +16,4 @@
 
 disp('      creating drag');
-md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
 %Take care of iceshelves: no basal drag