Index: /issm/trunk/test/Par/79North.par
===================================================================
--- /issm/trunk/test/Par/79North.par	(revision 4949)
+++ /issm/trunk/test/Par/79North.par	(revision 4949)
@@ -0,0 +1,45 @@
+%Start defining model parameters here
+
+%Create geometry
+load('79Northdata.mat');
+md.surface=surface;
+md.bed=bed;
+md.thickness=thickness;
+md.vx=vx;
+md.vy=vy;
+
+%Put melting on ice shelves
+md.melting_rate=zeros(md.numberofgrids,1);
+pos=find(md.gridoniceshelf);
+md.melting_rate(pos)=10;
+
+%Create drag
+md.drag_type=2;
+md.drag_coefficient=50*ones(md.numberofgrids,1);
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+%Create temperature and rheology
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md);
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+
+%Parallel options
+md.np=7;
+md.time=50;
+md.waitonlock=90;
+
+%Misc
+md.verbose=0;
+md.eps_abs=NaN;
+md.eps_rel=0.05;
Index: /issm/trunk/test/Par/Pig.par
===================================================================
--- /issm/trunk/test/Par/Pig.par	(revision 4949)
+++ /issm/trunk/test/Par/Pig.par	(revision 4949)
@@ -0,0 +1,36 @@
+%Start defining model parameters here
+
+load('Pigdata.mat');
+md.surface=surface;
+md.bed=bed;
+md.thickness=thickness;
+md.vx_obs=vx_obs;
+md.vy_obs=vy_obs;
+
+%Create drag
+md.drag_type=2;
+md.drag_coefficient=50*ones(md.numberofgrids,1);
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+%Create temperature and rheology
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md);
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=60;
+
+%Misc
+md.verbose=0;
+md.eps_abs=NaN;
+md.eps_rel=0.005;
Index: /issm/trunk/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 4949)
+++ /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 4949)
@@ -0,0 +1,36 @@
+%Start defining model parameters here
+
+%dynamics
+md.verbose=0;
+md.artificial_diffusivity=1;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+%Create drag
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+%Create temperature and rheology
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetIceSheetBC(md);
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Par/SquareShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareShelf.par	(revision 4949)
+++ /issm/trunk/test/Par/SquareShelf.par	(revision 4949)
@@ -0,0 +1,35 @@
+%Start defining model parameters here
+
+%dynamics
+md.verbose=0;
+md.artificial_diffusivity=1;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+%Create temperature and rheology
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetIceShelfBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 4949)
+++ /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 4949)
@@ -0,0 +1,35 @@
+%Start defining model parameters here
+
+%dynamics
+md.verbose=0;
+md.artificial_diffusivity=1;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+%Create temperature and rheology
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetIceSheetBC(md);
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;